Density Functional Study of Cobalt-59 Nuclear Magnetic Resonance Chemical Shifts and Shielding Tensor Elements in Co(III) Complexes

The Scope of the Applicability of Non-relativistic DFT Calculations of NMR Chemical Shifts in Pyridine-Metal Complexes for Applied Applications

Heavy metals are toxic, but it is impossible to stop using them. Considering the variety of molecular systems in which they can be present, the multicomponent nature and disorder of the structure of such systems, one of the most effective methods for studying them is NMR spectroscopy. This determines the need to calculate NMR chemical shifts for expected model systems. For elements beyond the third row of the periodic table, corrections for relativistic effects are necessary when calculating NMR parameters. Such corrections may be necessary even for light atoms due to the shielding effect of a neighboring heavy atom. This work examines the extent to which non-relativistic DFT calculations are able to reproduce experimental 15N and 113Cd NMR chemical shift tensors in pyridine-metal coordination complexes. It is shown that while for the calculation of 15N NMR chemical shift tensors there is no real need to consider relativistic corrections, for 113Cd, on the contrary, none of the tested calculation methods could reproduce the experimentally obtained tensor to any extent correctly.

Electron and Spin Delocalization in [Co6 Se8 (PEt3 )6 ]0/+1 Superatoms

Molecular clusters can function as nanoscale atoms/superatoms, assembling into superatomic solids, a new class of solid-state materials with designable properties through modifications on superatoms. To explore possibilities on diversifying building blocks, here we thoroughly studied one representative superatom, Co6 Se8 (PEt3 )6 . We probed its structural, electronic, and magnetic properties and revealed its detailed electronic structure as valence electrons delocalize over inorganic [Co6 Se8 ] core while ligands function as an insulated shell. 59 Co SSNMR measurements on the core and 31 P, 13 C on the ligands show that the neutral Co6 Se8 (PEt3 )6 is diamagnetic and symmetric, with all ligands magnetically equivalent. Quantum computations cross-validate NMR results and reveal degenerate delocalized HOMO orbitals, indicating aromaticity. Ligand substitution keeps the inorganic core nearly intact. After losing one electron, the unpaired electron in [Co6 Se8 (PEt3 )6 ]+1 is delocalized, causing paramagnetism and a delocalized electron spin. Notably, this feature of electron/spin delocalization over a large cluster is attractive for special single-electron devices.

A Small Lattice Change Induces Significant Dynamic Changes of CH3NH3+ Caged in Hybrid Perovskite Crystals: Toward Understanding the Interplay between Host Lattices and Guest Molecules

Two perovskite-type compounds, (MA)₂[B'Co(CN)₆] (MA = methylammonium, B' = K(I) and Na(I)), have very similar structures, but exhibit marked differences in the phase and dielectric transitions. Solid state ²H NMR studies reveal the detailed dynamic changes of the caged methylammonium (MA) cations before and after the phase transitions, which are correlated with the different dielectric states of the compounds. Using solid state ⁵⁹Co NMR, the dynamic changes of the host lattices before and after the transitions, which accompany the changes in the dynamics of the caged MA cations, are unveiled, demonstrating the intriguing interplay between the MA cations and the host lattices. On the basis of these observations, the molecular origins of the dielectric transitions are discussed in detail.

Computational, electrochemical, and spectroscopic studies of two mononuclear cobaloximes: the influence of an axial pyridine and solvent on the redox behaviour and evidence for pyridine coordination to cobalt(i) and cobalt(ii) metal centres

[Co(dmgBF2)2(H2O)2] (where dmgBF2 = difluoroboryldimethylglyoximato) was used to synthesize [Co(dmgBF2)2(H2O)(py)]·0.5(CH3)2CO (where py = pyridine) in acetone. The formulation of complex was confirmed by elemental analysis, high resolution MS, and various spectroscopic techniques. The complex [Co(dmgBF2)2(solv)(py)] (where solv = solvent) was readily formed in situ upon the addition of pyridine to complex . A spectrophotometric titration involving complex and pyridine proved the formation of such a species, with formation constants, log K = 5.5, 5.1, 5.0, 4.4, and 3.1 in 2-butanone, dichloromethane, acetone, 1,2-difluorobenzene/acetone (4 : 1, v/v), and acetonitrile, respectively, at 20 °C. In strongly coordinating solvents, such as acetonitrile, the lower magnitude of K along with cyclic voltammetry, NMR, and UV-visible spectroscopic measurements indicated extensive dissociation of the axial pyridine. In strongly coordinating solvents, [Co(dmgBF2)2(solv)(py)] can only be distinguished from [Co(dmgBF2)2(solv)2] upon addition of an excess of pyridine, however, in weakly coordinating solvents the distinctions were apparent without the need for excess pyridine. The coordination of pyridine to the cobalt(ii) centre diminished the peak current at the Epc value of the Co(I/0) redox couple, which was indicative of the relative position of the reaction equilibrium. Herein we report the first experimental and theoretical (59)Co NMR spectroscopic data for the formation of Co(i) species of reduced cobaloximes in the presence and absence of py (and its derivatives) in CD3CN. From spectroelectrochemical studies, it was found that pyridine coordination to a cobalt(i) metal centre is more favourable than coordination to a cobalt(ii) metal centre as evident by the larger formation constant, log K = 4.6 versus 3.1, respectively, in acetonitrile at 20 °C. The electrosynthesis of hydrogen by complexes and in various solvents demonstrated the dramatic effects of the axial ligand and the solvent on the turnover number of the respective catalyst.

Computational rationalization of the selective C–H and C–F activations of fluoroaromatic imines and ketones by cobalt complexes

The selective C–H and C–F activations of fluoroaromatic imines and ketones by cobalt complexes have been rationalized well by performing DFT calculations.

NMR, IR, Mössbauer and quantum chemical investigations of metalloporphyrins and metalloproteins

We review contributions made towards the elucidation of CO and O2binding geometries in respiratory proteins. Nuclear magnetic resonance, infrared spectroscopy, Mössbauer spectroscopy, X-ray crystallography and quantum chemistry have all been used to investigate the Fe –ligand interactions. Early experimental results showed linear correlations between17O chemical shifts and the infrared stretching frequency (νCO) of the CO ligand in carbonmonoxyheme proteins and between the17O chemical shift and the13CO shift. These correlations led to early theoretical investigations of the vibrational frequency of carbon monoxide and of the13C and17O NMR chemical shifts in the presence of uniform and non-uniform electric fields. Early success in modeling these spectroscopic observables then led to the use of computational methods, in conjunction with experiment, to evaluate ligand-binding geometries in heme proteins. Density functional theory results are described which predict57Fe chemical shifts and Mössbauer electric field gradient tensors,17O NMR isotropic chemical shifts, chemical shift tensors and nuclear quadrupole coupling constants (e2qQ/h) as well as13C isotropic chemical shifts and chemical shift tensors in organometallic clusters, heme model metalloporphyrins and in metalloproteins. A principal result is that CO in most heme proteins has an essentially linear and untilted geometry (τ = 4 °, β = 7 °) which is in extremely good agreement with a recently published X-ray synchrotron structure. CO / O2discrimination is thus attributable to polar interactions with the distal histidine residue, rather than major Fe–C–O geometric distortions.

Solid-state 13C and 59Co NMR spectroscopy of 13C-methylcobalt(iii) complexes with amine ligands

Five octahedral Co(iii) cations, [trans-Co(en)(2)(X)((13)CH(3))](n+) where en = ethylenediamine, X = CN(-), N(3)(-), NH(3), NO(2)(-) or H(2)O and n = 1 or 2, as well as [Co(NH(3))(5)(13)CH(3)](2+), have been investigated by solid-state (13)C and (59)Co NMR spectroscopy. We show that the determination of the (59)Co nuclear quadrupolar parameters both directly via(59)Co NMR and indirectly via(13)C NMR provide complementary information that is unavailable if one investigates only one nucleus. Specifically, (1)J((59)Co,(13)C) and the orientation of the largest component of the EFG were determined via(13)C NMR spectroscopy, which also established the negative sign of C(Q)((59)Co). Cobalt-59 NMR spectroscopy was used to characterize the cobalt magnetic shielding tensor, to verify the magnitudes of C(Q)((59)Co) and to establish the value of eta(Q), which is difficult to determine indirectly. The measurements show that the EFG tensors are either axially symmetric or close to being so, but there is a wide range of C(Q) values, from -40 MHz for the complex with X = H(2)O to -105 MHz with X = CN(-). The Co chemical shift tensors are approximately axially symmetric with the spans, delta(11)-delta(33), ranging from 3700 to 5600 ppm for X = H(2)O and CN(-), respectively. The latter measurements also established the relative orientations of the Co EFG and chemical shift tensors. Density functional theory calculations of the (59)Co EFG and magnetic shielding tensors as well as of (1)J((59)Co,(13)C) for the NO(2)(-) and N(3)(-) complexes were undertaken. These calculations confirm the experimental observation that the sign of C(Q) is negative and that the largest component of the EFG is along the Co-methyl-carbon bond.

Dinuclear fluoro-peroxovanadium(v) complexes with symmetric and asymmetric peroxo bridges: syntheses, structures and DFT studies

Two new dinuclear fluoro peroxovanadium(v) complexes, Cs3[V2O2(O2)4F] x H2O (1) and Cs3[V2O2(O2)3F3] x 2HF x H2O (2), were prepared and characterized by elemental analysis, IR spectroscopy, thermal analysis and X-ray crystallography. While the anion in possesses an asymmetric structure with a micro-eta1:eta2 bridging peroxo group, the [V2O2(O2)3F3]3- ion in exhibits a symmetrical structure with a unique mu-fluoro and micro-eta2:eta2 peroxo double bridge. The X-ray structure data were compared with equilibrium and vibrationally-averaged (effective) DFT calculated geometries. The decomposition reactions of and in aqueous solution were studied by 51V NMR spectroscopy. The calculations of vibrationally averaged NMR chemical shifts (DFT-GIAO) were used to support the empirical assignment of NMR signals and afforded excellent agreement with experimental values for the studied peroxovanadium species. The ESI mass spectra of the prepared compounds are in accordance with the assignment of NMR spectra and with DFT study.

Computational study of structures and properties of metallaboranes: cobalt bis(dicarbollide)

A density functional study at the BP86/AE1 level is presented for the cobalt bis(dicarbollide) ion [3-Co-(1,2-C2B9H11)2]- (1) and selected isomers and rotamers thereof. Rotation of the two dicarbollide moieties with respect to each other is facile, as judged by the small energetic separation of the three rotamers located (within 11 kJ mol(-1)) and by the low barriers for their interconversion (at most 41 kJ mol(-1)). Among the isomers differing in carbon atom positions that contain two equivalent dicarbollide ligands, the 1,7 ("carbon apart") form [2-Co-(1,7-C2B9H11)2]- is the most stable, 121 kJ mol(-1) below 1. The electronic structure of 1 is characterized in terms of molecular orbitals, population analysis, and excitation energies from time-dependent density functional theory, relevant to UV/Vis spectroscopy. Experimental 11B NMR chemical shifts of 1 are reproduced to better than 5 ppm at the GIAO-B3LYP/II' level, and the computed delta(11B) values are only little affected by rotational averaging or the presence of a polarizable continuum. Larger such effects are found for the as-yet unknown 59Co chemical shift, for which a value in the range between -1800 and -2400 ppm is predicted. Even though the accuracy achieved for the theoretical delta(11B) values is somewhat lower than that for heteroboranes at conventional ab initio levels, the level of density functional employed can afford qualitatively reliable chemical shifts, which can be useful in assignments and structural refinements of heteroboranes containing transition metal.

Theoretical 49Ti NMR chemical shifts

49Ti chemical shifts for a total of 20 titanium complexes are reported, and several levels of theory are evaluated in order to identify a reliable approach for the calculation of titanium NMR data. The popular B3LYP/6-31G(d)//B3LYP/6-31G(d) proves to give very good agreement with experimental data over a range from 1,400 to -1,300 ppm. The MP2/6-31G(d)//MP2/6-31G(d) level computes even smaller average deviations but fails for TiI(4). This behavior together with its huge demand for computational resources requires careful handling of this theoretical level. In addition, NMR data for five titanium fulvene (or related) complexes are given.

Simulation of59Co NMR Chemical Shifts in Aqueous Solution

59Co chemical shifts were computed at the GIAO-B3LYP level for [Co(CN)6]3-, [Co(H2O)6]3+, [Co(NH3)6]3+, and [Co(CO)4]- in water. The aqueous solutions were modeled by Car-Parrinello molecular dynamics (CPMD) simulations, or by propagation on a hybrid quantum-mechanical/molecular-mechanical Born-Oppenheimer surface (QM/MM-BOMD). Mean absolute deviations from experiment obtained with these methods are on the order of 400 and 600 ppm, respectively, over a total delta(59Co) range of about 18,000 ppm. The effect of the solvent on delta(59Co) is mostly indirect, resulting primarily from substantial metal-ligand bond contractions on going from the gas phase to the bulk. The simulated solvent effects on geometries and delta(59Co) values are well reproduced by using a polarizable continuum model (PCM), based on optimization and perturbational evaluation of quantum-mechanical zero-point corrections.

Quantum chemical studies of protein structure

Quantum chemical methods now permit the prediction of many spectroscopic observables in proteins and related model systems, in addition to electrostatic properties, which are found to be in excellent accord with those determined from experiment. I discuss the developments over the past decade in these areas, including predictions of nuclear magnetic resonance chemical shifts, chemical shielding tensors, scalar couplings and hyperfine (contact) shifts, the isomer shifts and quadrupole splittings in Mössbauer spectroscopy, molecular energies and conformations, as well as a range of electrostatic properties, such as charge densities, the curvatures, Laplacians and Hessians of the charge density, electrostatic potentials, electric field gradients and electrostatic field effects. The availability of structure/spectroscopic correlations from quantum chemistry provides a basis for using numerous spectroscopic observables in determining aspects of protein structure, in determining electrostatic properties which are not readily accessible from experiment, as well as giving additional confidence in the use of these techniques to investigate questions about chemical bonding and chemical reactions.

Thermal Effects and Vibrational Corrections to Transition Metal NMR Chemical Shifts

Both zero-point and classical thermal effects on the chemical shift of transition metals have been calculated at appropriate levels of density functional theory for a number of complexes of titanium, vanadium, manganese and iron. The zero-point effects were computed by applying a perturbational approach, whereas classical thermal effects were probed by Car-Parrinello molecular dynamics simulations. The systematic investigation shows that both procedures lead to a deshielding of the magnetic shielding constants evaluated at the GIAO-B3 LYP level, which in general also leads to a downfield shift in the relative chemical shifts, delta. The effect is small for the titanium and vanadium complexes, where it is typically on the order of a few dozen ppm, and is larger for the manganese and iron complexes, where it can amount to several hundred ppm. Zero-point corrections are usually smaller than the classical thermal effect. The pronounced downfield shift is due to the sensitivity of the shielding of the metal centre with regard to the metal-ligand bond length, which increase upon vibrational averaging. Both applied methods improve the accuracy of the chemical shifts in some cases, but not in general.

NMR shielding calculations across the periodic table: diamagnetic uranium compounds. 2. Ligand and metal NMR

In this and a previous article (J. Phys. Chem. A 2000, 104, 8244), the range of application for relativistic density functional theory (DFT) is extended to the calculation of nuclear magnetic resonance (NMR) shieldings and chemical shifts in diamagnetic actinide compounds. Two relativistic DFT methods are used, ZORA ("zeroth-order regular approximation") and the quasirelativistic (QR) method. In the given second paper, NMR shieldings and chemical shifts are calculated and discussed for a wide range of compounds. The molecules studied comprise uranyl complexes, [UO(2)L(n)](+/-)(q); UF(6); inorganic UF(6) derivatives, UF(6-n)Cl(n), n = 0-6; and organometallic UF(6) derivatives, UF(6-n)(OCH(3))(n), n = 0-5. Uranyl complexes include [UO(2)F(4)](2-), [UO(2)Cl(4)](2-), [UO(2)(OH)(4)](2-), [UO(2)(CO(3))(3)](4-), and [UO(2)(H(2)O)(5)](2+). For the ligand NMR, moderate (e.g., (19)F NMR chemical shifts in UF(6-n)Cl(n)) to excellent agreement [e.g., (19)F chemical shift tensor in UF(6) or (1)H NMR in UF(6-n)(OCH(3))(n)] has been found between theory and experiment. The methods have been used to calculate the experimentally unknown (235)U NMR chemical shifts. A large chemical shift range of at least 21,000 ppm has been predicted for the (235)U nucleus. ZORA spin-orbit appears to be the most accurate method for predicting actinide metal chemical shifts. Trends in the (235)U NMR chemical shifts of UF(6-n)L(n) molecules are analyzed and explained in terms of the calculated electronic structure. It is argued that the energy separation and interaction between occupied and virtual orbitals with f-character are the determining factors.

(57)Fe Mössbauer isomer shifts of heme protein model systems: electronic structure calculations

We report the results of density functional theory (DFT) calculations of the (57)Fe Mössbauer isomer shifts (delta(Fe)) for a series of 24 inorganic, organometallic, and metalloprotein/metalloporphyrin model systems in S = 0, (1)/(2), 1, (3)/(2), 2, and (5)/(2) spin states. We find an excellent correlation between calculation and experiment over the entire 2.34 mm s(-1) range of isomer shifts: a 0.07-0.08 mm s(-1) rms deviation between calculation and experiment (corresponding to 3-4% of the total delta(Fe) range, depending on the functionals used) with R(2) values of 0.973 and 0.981 (p < 0.0001). The best results are obtained by using the hybrid exchange-correlation functional B3LYP, used previously for (57)Fe Mössbauer quadrupole splittings and (57)Fe NMR chemical shifts and chemical shielding anisotropies. The relativistically corrected value of alpha, alpha(rel), converges with the large basis set used in this work, but the exact values vary somewhat with the methods used: -0.253 a(0)(3) mm s(-1) (Hartree-Fock; HF); -0.316 a(0)(3) mm s(-1) (hybrid HF-DFT; B3LYP), or -0.367 a(0)(3) mm s(-1) (pure DFT; BPW91). Both normal and intermediate spin state isomer shifts are well reproduced by the calculations, as is the broad range of delta(Fe) values: from [Fe(VI)O(4)](2-) (-0.90 mm s(-1) expt; -1.01 mm s(-1) calc) to KFe(II)F(3) (1.44 mm s(-1) expt; 1.46 mm s(-1) calc). Molecular orbital analyses of all inorganic solids as well as all organometallic and metalloporphyrin systems studied reveal that there are three major core MO contributions to rho(tot)(0), the total charge density at the iron nucleus (and hence delta(Fe)), that do not vary with changes in chemistry, while the valence MO contributions are highly correlated with delta(Fe) (R(2) = 0.915-0.938, depending on the functionals used), and the correlation between the valence MO contributions and the total MO contribution is even better (R(2) = 0.965-0.976, depending on the functionals used). These results are of general interest since they demonstrate that DFT methods now enable the accurate prediction of delta(Fe) values in inorganic, organometallic, and metalloporphyrin systems in all spin states and over a very wide range of delta(Fe) values with a very small rms error.

Calculation of electronic g-tensors for transition metal complexes using hybrid density functionals and atomic meanfield spin-orbit operators

We report the first implementation of the calculation of electronic g-tensors by density functional methods with hybrid functionals. Spin-orbit coupling is treated by the atomic meanfield approximation. g-Tensors for a set of small main group radicals and for a series of ten 3d and two 4d transition metal complexes have been compared using the local density approximation (VWN functional), the generalized gradient approximation (BP86 functional), as well as B3-type (B3PW91) and BH-type (BHPW91) hybrid functionals. For main group radicals, the effect of exact-exchange mixing is small. In contrast, significant differences between the various functionals arise for transition metal complexes. As has been shown previously, local and in particular gradient-corrected functionals tend to underestimate the "paramagnetic" contributions to the g-tensors in these cases and thereby recover only about 40-50% of the range of experimental g-tensor components. This is improved to ca. 60% by the B3PW91 functional, which also gives slightly reduced standard deviations. The range increases to almost 100% using the half-and-half functional BHPW91. However, the quality of the correlation with experimental data worsens due to a significant overestimate of some intermediate g-tensor values. The worse performance of the BHPW91 functional in these cases is accompanied by spin contamination. Although none of the functionals tested thus appears to be ideal for the treatment of electronic g-tensors in transition metal complexes, the B3PW91 hybrid functional exhibited the overall most satisfactory performance. Apart from the validation of hybrid functionals, some aspects in the treatment of spin-orbit contributions to the g-tensor are discussed.

Medium effects on 51V NMR chemical shifts: a density functional study

Car-Parrinello molecular dynamics simulations were performed for [H2VO4], [VO2(OH2)4]+, and [VO(O2)2(OH2)]- in periodic boxes with 30, 28, and 29 water molecules, respectively, employing the BLYP density functional. On the timescale of the simulations, up to 2 ps, well-structured first solvation spheres are discernible for [H2VO4]- and [VO(O2)2(OH2)]- containing, on average, eight and ten water molecules, respectively. One of the four water molecules directly attached to the metal in [VO2(OH2)4]+ is only loosely bound, and the average coordination number of vanadium in aqueous VO2+ is between five and six. 51V chemical shifts were evaluated at the B3LYP level for representative snapshots along the trajectories, including the water molecules of the solvent by means of point charges. The resulting averaged delta(51V) values are proposed to model the combined effects of temperature (dynamic averaging) and solvent (charge polarization). Both effects are shown to be rather small, of the order of a few dozen ppm. The observed shielding of 51V in the bis(peroxo) complex with respect to the vanadate species is not reproduced computationally.

Cobalt-59 chemical shift and quadrupolar tensors of simple octahedral cobalt(III) complexes

Analysis of the solid-state powder 59Co NMR spectra of ten simple inorganic cobalt(III) complexes at 11.75, and in most cases, 4.7 T have permitted the assignment of specific ligand planes to ranges of values of the observed chemical shift principal components. The relevant chemical shift components were determined from the simulations of the powder line shapes. These simulations also provided the relative orientations of the chemical shift (CS) and electric field gradient (efg) tensors, as well as magnitude and asymmetry of the 59Co quadrupolar coupling. Using symmetry arguments and ab initio calculations, as appropriate or necessary, the orientations of the efg tensors in the molecular frame were deduced. This allowed the determination of the CS tensors in the molecular frame and thus assignment of the ligand planes responsible for the observed values of chemical shifts.Key words: cobalt, chemical shift, quadrupolar coupling, solid state NMR.