DFT study of Co-C bond cleavage in the neutral and one-electron-reduced alkyl-cobalt(III) phthalocyanines

Density functional theory (DFT) has been applied to the analysis of the structural and electronic properties of the alkyl-cobalt(III) phthalocyanine complexes, [CoIIIPc]-R (Pc = phthalocyanine, R = Me or Et), and their pyridine adducts. The BP86/6-31G(d) level of theory shows good reliability for the optimized axial bond lengths and bond dissociation energies (BDEs). The mechanism of the reductive cleavage was probed for the [CoIIIPc]-Me complex which is known as a highly effective methyl group donor. In the present analysis, which follows a recent study on the reductive Co-C bond cleavage in methylcobalamin (J. Phys. Chem. B 2007, 111, 7638-7645), it is demonstrated that addition of an electron and formation of the pi-anion radical [CoIII(Pc*)]-Me- significantly lowers the energetic barrier required for homolytic Co-C bond dissociation. Such BDE lowering in [CoIII(Pc*)]-Me- arises from the involvement of two electronic states: upon electron addition, a quasi-degenerate pi*Pc state is initially formed, but when the cobalt-carbon bond is stretched, the unpaired electron moves to a sigma*Co-C state and the final cleavage involves the three-electron (sigma)2(sigma*)1 bond. As in corrin complexes, the pi*Pc-sigma*Co-C states crossing does not take place at the equilibrium geometry of [CoIII(Pc*)]-Me- but only when the Co-C bond is stretched to approximately 2.3 A. The DFT computed Co-C BDE of 23.3 kcal/mol in the one-electron-reduced phthalocyanine species, [CoIII(Pc*)]-Me-, is lowered by approximately 37% compared to the neutral Py-[CoIIIPc]-Me complex where BDE = 36.8 kcal/mol. A similar comparison for the corrin-containing complexes shows that a DFT computed BDE of 20.4 kcal/mol for [CoIII(corrin*)]-Me leads to approximately 45% bond strength reduction, in comparison to 37.0 kcal/mol for Im-[CoIII(corrin)]-Me+. These results suggest some preference by the alkylcorrinoids for the reductive cleavage mechanism.     

Theoretical study of M(PH3)2 complexes of C60, corannulene (C20H10), and sumanene (C21H12) (M = Pd or Pt). Unexpectedly large binding energy of M(PH3)2(C60)

DFT and MP2 to MP4(SDQ) methods were applied to M(PH3)2(C60), Pt(PH3)2(C20H10), and Pt(PH3)2(C21H12) (M = Pd or Pt, C20H10 = corannulene, and C21H12 = sumanene). The binding energy considerably fluctuates around MP2 and MP3 levels but much less upon going from MP3 to MP4(SDQ) in Pt(PH3)2(C2H4), Pt(PH3)2(C20H10), and Pt(PH3)2(C21H12). Also, the MP4(SDQ) method presents a binding energy similar to that of the CCSD(T) method in Pt(PH3)2(C2H4). Thus, it is likely that the MP4(SDQ) method is useful to evaluate binding energies of these complexes. The binding energies of Pt(PH3)2(C20H10) and Pt(PH3)2(C21H12) are evaluated to be 24.9 and 26.1 kcal/mol, respectively, by the MP4(SDQ) method and only +5.8 and -2.6 kcal/mol, respectively, by the DFT(B3LYP) method. These MP4(SDQ)-calculated binding energies of Pt(PH3)2(C20H10) and Pt(PH3)2(C21H12) are similar to that of Pt(PH3)2(C2H4), which strongly suggests that these complexes can be successfully synthesized. The binding energy of Pt(PH3)2(C60) is evaluated to be 44.8 and 45.5 kcal/mol with the ONIOM(MP4(SDQ):UFF) and ONIOM(MP4(SDQ):B3LYP) methods, respectively, and that of the Pd analogue is evaluated to be 39.9 kcal/mol with the ONIOM(MP4(SDQ):UFF) method, whereas the DFT(B3LYP), DFT(BVP86), and DFT(BPW91) methods provide much smaller binding energies. It is noted that these binding energies are much larger than those of the ethylene, corannulene, and sumanene analogues. This difference is reasonably interpreted in terms that the LUMO of C60 is at much lower energy than those of ethylene, corannulene, and sumanene. We investigated also how to separate the high level and the low level regions in the ONIOM calculation of M(PH3)2(C60) and proposed here the reasonable way to evaluate the binding energy of transition-metal complexes of C60.     

Electroreduction of a series of alkylcobalamins: mechanism of stepwise reductive cleavage of the Co-C bond

The electrochemical (EC) reduction mechanism of methylcobalamin (Me-Cbl) in a mixed DMF/MeOH solvent in 0.2 M tetrabutylammonium fluoroborate electrolyte was studied as a function of temperature and solvent ratio vs a nonaqueous Ag/AgCl/Cl(-) reference electrode. Double-potential-step chronoamperometry allowed the rate constant of the subsequent homogeneous reaction to be measured over the temperature range from 0 to -80 degrees C in 40:60 and 50:50 DMF:MeOH ratios. Activation enthalpies are 5.8 +/- 0.5 and 7.6 +/- 0.3 kcal/mol in the 40:60 and 50:50 mixtures of DMF/MeOH, respectively. Digital simulation and curve-fitting for an EC mechanism using a predetermined homogeneous rate constant of 5.5 x 10(3) s(-1) give E degrees' = -1.466 V, k degrees = 0.016 cm/s, and alpha = 0.77 at 20 degrees C for a quasi-reversible electrode process. Digital simulation of the results of Lexa and Savéant (J. Am. Chem. Soc. 1978, 100, 3220-3222) shows that the mechanism is a series of stepwise homogeneous equilibrium processes with an irreversible step following the initial electron transfer (ET) and allows estimation of the equilibrium and rate constants of these reactions. An electron coupling matrix element of H(kA) = (4.7 +/- 1.1) x 10(-4) eV ( approximately 46 J/mol) is calculated for the nonadiabatic ET step for reduction to the radical anion. A reversible bond dissociation enthalpy for homolytic cleavage of Me-Cbl is calculated as 31 +/- 2 kcal/mol. The voltammetry of the ethyl-, n-propyl-, n-butyl-, isobutyl-, and adenosyl-substituted cobalamin was studied, and estimated reversible redox potentials were correlated with Co-C bond distances as determined by DFT (B3LYP/ LANL2DZ) calculations.     

Comparative study of metal-porphyrins, -porphyrazines,and -phthalocyanines

A theoretical comparative study of complexes of porphyrin (P), porphyrazine (Pz), and phthalocyanine (Pc) with metal (M) = Fe, Co, Ni, Cu, and Zn has been carried out using a DFT method. The calculations provide a clear elucidation of the ground states for the MP/Pz/Pc molecules and for a series of [MP/Pz/Pc](x-) and [MP/Pz/Pc](y+) ions (x = 1, 2, 3, 4; y = 1, 2). There are significant differences among MP, MPz, and MPc in the electronic structure and other calculated properties. For FeP/Pz and CoP/Pz, the first oxidation occurs at the central metal, while it is the macroring of FePc and CoPc that is the site of oxidation. The smaller coordination cavity results in a stronger ligand field in Pz than in P. However, the benzo annulation produces a surprisingly strong destabilizing effect on the metal-macrocycle bonding. The effects of Cl axial bonding upon the electronic structures of the iron(III) complexes of P, Pz, and Pc were examined, as was the bonding of pyridine (py) to NiP, NiPz, and NiPc. The porphinato core size plays a crucial role in controlling the spin state of Fe(III) in these complexes. FePc(Cl) is predicted to be a pure intermediate-spin system, whereas NiPz(py)(2) and NiPc(py)(2) are metastable in high-spin (S = 1) states. The NiPz/Pc-(py)(2) binding energy curve has only a shallow well that facilitates decomposition of the complex. The NiP-(py)(2) bond energy is small, but the relatively deep well in the binding energy curve ought to make this system stable to decomposition.     

Quantum-chemical modeling of the molecular structures of (555)macrotricyclic chelates in M(II) ion–thiooxamide–glyoxal ternary systems (M = Mn,Fe, Co,Ni, Cu,Zn)

The thermodynamic and geometric parameters of isomeric macrotricyclic Mn(II), Fe(II), Co(II), Ni(II), Cu(II), and Zn(II) complexes that can form upon the complexation of the corresponding hexacyanoferrates( II) with thiooxamide H₂N–C(=S)–C(=O)–NH₂ and glyoxal HC(=O)–CH(=O) in gelatin-immobilized matrices have been calculated by the OPBE/TZVP DFT method with the use of the Gaussian09 program package. It has been found that a complex with the MN₄ chelate core is most stable for M = Mn, Fe, Co, Ni, and Zn, and the MN₂S₂ core is most stable for M = Cu. Bond lengths and bond angles have been reported, and it has been noted that in all complexes, except the Zn(II) one, the chelate core and three fivemembered chelate rings are almost planar.     

Nickel dithiolenes revisited: structures and electron distribution from density functional theory for the three-member electron-transfer series [Ni(S2C2Me2)2]0,1-,2-

The complexes [Ni(S2C2Me2)2](z) (z = 0, 1-, 2-) have been isolated for the purpose of investigating their electronic structures in a reversible three-member electron-transfer series. Members are interrelated by reversible redox reactions with E(1/2)(0/1-) = -0.15 V and E(1/2)(1-/2-) = -1.05 V versus SCE in acetonitrile. The three complexes have nearly planar structures of idealized D(2)(h) symmetry. As the series is traversed in the reducing direction, Ni-S and C-S bond lengths increase; the chelate ring C-C bond length decreases from the neutral complex to the monoanion and does not change significantly in the dianion. Structural trends are compared with previous results for [Ni(S2C2R2)2)](1-,2-). Following the geometrical changes, values of nu(Ni)(-)(S) and nu(C)(-)(S) decrease, while the value of nu(C)(-)(C) increases with increased reduction. Geometry optimizations at the density functional theory (DFT) level were performed for all members of the series. Geometrical parameters obtained from the calculations are in good agreement with the experimental findings. The 5b(2g) orbital was identified as the LUMO in [Ni(S2C2Me2)2], the SOMO in [Ni(S2C2Me2)2](1-), and the HOMO in [Ni(S2C2Me2)2]2-. Unlike in the situation in the [M(CO)2-(S2C2Me2)2]z series (M = Mo, W; z = 0, 1-, 2-), the apparent contribution from the metal d orbital in the electroactive orbital is not constant. In the present series, the d(xz) contribution increases from 13 to 20 to 39% upon passing from the neutral to the monoanionic to the dianionic complex. Accurate calculation of EPR g-values of [Ni(S2C2Me2)2]1- by DFT serves as a test for the reliability of the electronic structure calculations.     

60]Fullerene metal complexes with large effective two-photon absorption cross-section

Fourteen novel complexes of fullerene C(60) with metal dialkyldithiophosphate, {[(RO)(2)PS(2)](2)M}·2C(60) (R = Me and Et; M = Mn, Fe, Co, Ni, Cu, Zn and Pd), can be obtained in high yield by the reaction of metal(II) dialkyldithiophosphate complexes with C(60) fullerene. Their structures are determined by (1)H, (13)C, and (31)P NMR spectroscopy, elemental analysis, FT-IR and UV-vis, and are supplied by single crystal X-ray data. The compound {[(EtO)(2)PS(2)](2)Pd}·2C(60) was also studied by the HRTEM image and SAED patterns. Studies of the optical properties of the complexes of fullerene C(60) with metal dialkyldithiophosphate show that these compounds all have very strong nonlinear optical absorption effects. The nonlinear absorption of fourteen complexes of fullerene C(60) with metal dialkyldithiophosphate (in o-dichlorobenzene and in the solid state) was measured by the open-aperture Z-scan technique at a wavelength of 532 nm. DFT calculations are also used to discuss the stability of these complexes and to confirm the structural assignments.     

N4 ring as a square planar ligand in novel MN4 species

Ab initio (MP2) and density functional theory (DFT) methods are used to examine a series of MN4 compounds, where M is an alkaline-earth cation (Mg2+, Ca2+, Sr2+, Ba2+), and N4(2-) is a six--electron ring. All pyramidal structures except MgN4 are the most energetically favored for all singlet MN4 systems considered here. For MgN4, the CS structure with dicoordinated Mg out of the N4 ring plane is the most stable of all. Among these systems, the pyramidal CaN4, SrN4, BaN4 and the planar C(S) structure containing dicoordinated Ba are stable as singlet molecules due to their significant isomerization or dissociation barriers (21.3-94.1 kcal/mol). Structural, natural bond orbital (NBO), and molecular orbital (MO) analyses indicate that the bonding in the BaN4 system has a larger covalent character as compared with other MN4 systems. In addition, substantial d character is found in the bonding of the MN4 (M = Ca2+, Sr2+, Ba2+) species.     

Experimental and theoretical studies on the thermal decomposition of heterocyclic nitrosimines

A series of substituted 2-nitrosiminobenzothiazolines (2) were synthesized by the nitrosation of the corresponding 2-iminobenzothiazolines (6). Thermal decomposition of 2a--f and of the seleno analogue 7 in methanol and of 3-methyl-2-nitrosobenzothiazoline (2a) in acetonitrile, 1,4-dioxane, and cyclohexane followed first-order kinetics. The activation parameters for thermal deazetization of 2a were measured in cyclohexane (Delta H(++) = 25.3 +/- 0.5 kcal/mol, Delta S(++) = 1.3 +/- 1.5 eu) and in methanol (Delta H(++) = 22.5 +/- 0.7 kcal/mol, Delta S(++) = -12.9 +/- 2.1 eu). These results indicate a unimolecular decomposition and are consistent with a proposed stepwise mechanism involving cyclization of the nitrosimine followed by loss of N(2). The ground-state conformations of the parent nitrosiminothiazoline (9a) and transition states for rotation around the exocyclic C==N bond, electrocyclic ring closure, and loss of N(2) were calculated using ab initio molecular orbital theory at the MP2/6-31G* level. The calculated gas-phase barrier height for the loss of N(2) from 9a (25.2 kcal/mol, MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE) compares favorably with the experimental barrier for 2a of 25.3 kcal/mol in cyclohexane. The potential energy surface is unusual; the rotational transition state 9a-rot-ts connects directly to the orthogonal transition state for ring-closure 9aTS. The decoupling of rotational and pseudopericyclic bond-forming transition states is contrasted with the single pericyclic transition state (15TS) for the electrocyclic ring-opening of oxetene (15) to acrolein (16). For comparison, the calculated homolytic strength of the N--NO bond is 40.0 kcal/mol (MP4(SDQ, FC)/6-31G*//MP2/6-31G* + ZPE).     

Interactions of metal ions with water: ab initio molecular orbital studies of structure, vibrational frequencies, charge distributions, bonding enthalpies, and deprotonation enthalpies. 2. Monohydroxides

The formation and properties of a wide range of metal ion monohydroxides, M(n)(+)[OH(-)], where n = 1 and 2, have been studied by ab initio molecular orbital calculations at the MP2(FULL)/6-311++G**//MP2(FULL)/6-311++G** and CCSD(T)(FULL)/6-311++G**//MP2(FULL)/6-311++G** computational levels. The ions M(n)()(+) are from groups 1A, 2A, 3A, and 4A in the second, third, and fourth periods of the Periodic Table and from the first transition series. Geometrical parameters, vibrational frequencies, atomic charge distributions, orbital occupancies, and bonding enthalpies are reported. The M(n)(+)-O distances are shorter in the hydroxides than in the corresponding hydrates (published previously as Part 1, Inorg. Chem. 1998, 37, 4421-4431) due to a greater electrostatic interaction in the hydroxides. The natural bond orbitals for most of the first-row transition metal ion hydroxides do not contain a formal metal-oxygen bonding orbital; nevertheless the atomic charge distributions show that for both n = 1 and 2 a significant amount of electron density is consistently transferred from the hydroxide ion to the bound metal ion. Deprotonation enthalpies for the hydrates have been evaluated according to the simple dissociation process, M(n)(+)[OH(2)] --> M(n)(+)[OH(-)] + H(+), and also via proton transfer to another water molecule, M(n)(+)[OH(2)] + H(2)O --> M(n)(+)[OH(-)] + H(3)O(+). The drastic reduction in these deprotonation enthalpies as H(2)O molecules are sequentially bonded in the first coordination shell of the metal ion (amounting to 71, 64, 85, and 91 kcal/mol for the bonding of six water molecules to Mg(2+), Ca(2+), Mn(2+), and Zn(2+), respectively) is found to be due to the greater decrease in the bonding enthalpies for the hydroxides relative to the hydrates. Proton transfer to bases other than water, for example side chain groups of certain amino acids, could more than offset the decrease in deprotonation energy due to the filling of the first coordination shell. Linear relationships have been found between the pK(a) values for ionization of the Mg(2+), Ca(2+), Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), and Zn(2+) aquo ions, and Delta for the bonding of the first water molecule, for the bonding of the hydroxide ion, and for proton dissociation from the monohydrate. Similar relationships have also been found between the pK(a) values and the reciprocal of the M-O bond lengths in both the monohydrates and hydroxides. Thus the ionization of metal hydrates in water echoes the properties of the monomeric species M(n)(+)[OH(2)].     

A computational characterization of boron-oxygen multiple bonding in HN=CH-CH=CH-NH-BO

Structures, relative energies, and bonding characteristics for various conformers of 3-imino-N-(oxoboryl)prop-1-en-1-amine, HN=CH-CH=CH-NH-BO, and the corresponding borocycle (-HN=CH-CH=CH-NH-B-)O are discussed using results from second-order M?ller-Plesset (MP2) perturbation theory with the Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. These MP2 results are compared to those from computationally efficient density functional theory (DFT) calculations using the LDA, PBE, TPSS, BLYP, B3LYP, BVP86, OLYP, O3LYP, and PBE1PBE functionals in conjunction with the economical Pople-type 6-311++G(d,p) basis set to evaluate the suitability of these DFT/6-311++G(d,p) levels for use with larger boron-containing systems. The effects of an aqueous environment were incorporated into the calculations using COSMO methodology. The calculated boron-oxygen bond lengths, orbital compositions, and bond orders in all the (acyclic) HN=CH-CH=CH-NH-BO conformers were consistent with the presence of a boron-oxygen triple bond, similar to that found in H-BO and H2N-BO. The (-HN=CH-CH=CH-NH-B-)O borocycle is predicted to be planar (C2v symmetry), and it is approximately 30 kcal/mol lower in energy than any of the (acyclic) HN=CH-CH=CH-NH-BO conformers; the boron-oxygen bond in this borocycle has significant double bond character, a bonding scheme for which there has been only one experimental structure reported in the literature (Vidovic, D. ; et al. J. Am. Chem. Soc. 2005, 127, 4566- 4569).     

Mapping the d-d excited-state manifolds of transition metal beta-diiminato-imido complexes. Comparison of density functional theory and CASPT2 energetics

Trigonal-planar, middle transition metal diiminato-imido complexes do not exhibit high-spin states, as might be naively expected on the basis of their low coordination numbers. Instead, the known Fe(III), Co(III), and Ni(III) complexes exhibit S = 3/2, S = 0, and S = 1/2 ground states, respectively. Kohn-Sham DFT calculations have provided a basic molecular orbital picture of these compounds as well as a qualitative rationale for the observed spin states. Reported herein are ab initio multiconfiguration second-order perturbation theory (CASPT2) calculations, which provide a relatively detailed picture of the d-d excited-state manifolds of these complexes. Thus, for a C(2v) Fe(III)(diiminato)(NPh) model complex, two near-degenerate states ((4)B(2) and (4)B(1)) compete as contenders for the ground state. Moreover, the high-spin sextet, two additional quartets and even a low-spin doublet all occur at <0.5 eV, relative to the ground state. For the Co(III) system, although CASPT2 reproduces an S = 0 ground state, as observed experimentally for a related complex, the calculations also predict two exceedingly low-energy triplet states; there are, however, no other particularly low-energy d-d excited states. In contrast to the Fe(III) and Co(III) cases, the Ni(III) complex has a clearly nondegenerate (2)B(2) ground state. The CASPT2 energetics provide benchmarks against which we can evaluate the performance of several common DFT methods. Although none of the functionals examined perform entirely satisfactorily, the B3LYP hybrid functional provides the best overall spin-state energetics.     

Theoretical study of C-H and N-H sigma-bond activation reactions by titinium(IV)-imido complex. Good understanding based on orbital interaction and theoretical proposal for N-H sigma-bond activation of ammonia

The C-H sigma-bond activation of methane and the N-H sigma-bond activation of ammonia by (Me3SiO)2Ti(=NSiMe3) 1 were theoretically investigated with DFT, MP2 to MP4(SDQ), and CCSD(T) methods. The C-H sigma-bond activation of methane takes place with an activation barrier (Ea) of 14.6 (21.5) kcal/mol and a reaction energy (DeltaE) of -22.7 (-16.5) kcal/mol to afford (Me3SiO)2Ti(Me)[NH(SiMe3)], where DFT- and MP4(SDQ)-calculated values are given without and in parentheses, respectively, hereafter. The electron population of the CH3 group increases, but the H atomic population decreases upon going to the transition state from the precursor complex, which indicates that the C-H sigma-bond activation occurs in heterolytic manner unlike the oxidative addition. The Ti atomic population considerably increases upon going to the transition state from the precursor complex, which indicates that the charge transfer (CT) occurs from methane to Ti. These population changes are induced by the orbital interactions among the d(pi)-p(pi) bonding orbital of the Ti=NSiMe3 moiety, the Ti d(z2) orbital and the C-H sigma-bonding and sigma*-antibonding orbitals of methane. The reverse regioselective C-H sigma-bond activation which leads to formation of (Me3SiO)2Ti(H)[NMe(SiMe3)] takes place with a larger Ea value and smaller exothermicity. The reasons are discussed in terms of Ti-H, Ti-CH3, Ti-NH3, N-H, and N-CH3 bond energies and orbital interactions in the transition state. The N-H sigma-bond activation of ammonia takes place in a heterolytic manner with a larger Ea value of 19.0 (27.9) kcal/mol and considerably larger exothermicity of -45.0 (-39.4) kcal/mol than those of the C-H sigma-bond activation. The N-H sigma-bond activation of ammonia by a Ti-alkylidyne complex, [(PNP)Ti(CSiMe3)] 3 (PNP = N-[2-(PH2)2-phenyl]2-]) ,was also investigated. This reaction takes place with a smaller E(a) value of 7.5 (15.3) kcal/mol and larger exothermicity of -60.2 (-56.1) kcal/mol. These results lead us to predict that the N-H sigma-bond activation of ammonia can be achieved by these complexes.     

Accurate calculation of the heats of formation for large main group compounds with spin-component scaled MP2 methods

Three MP2-type electron correlation treatments and standard density functional theory (DFT) approaches are used to predict the heats of formation for a wide variety of different molecules. The SCF and MP2 calculations are performed efficiently using the resolution-of-the-identity (RI) approximation such that large basis set (i.e., polarized valence quadruple-zeta quality) treatments become routinely possible for systems with 50-100 atoms. An atom equivalent scheme that corrects the calculated atomic energies is applied to extract the "real" accuracy of the methods for chemically relevant problems. It is found that the spin-component-scaled MP2 method (SCS-MP2, J. Chem. Phys, 2003, 118, 9095) performs best and provides chemical accuracy (MAD of 1.18 kcal/mol) for a G2/97 test set of molecules. The computationally more economical SOS-MP2 variant, which retains only the opposite-spin part of the correlation energy, is slightly less accurate (MAD of 1.36 kcal/mol) than SCS-MP2. Both spin-component-scaled MP2 treatments perform significantly better than standard MP2 (MAD of 1.77 kcal/mol) and DFT-B3LYP (MAD of 2.12 kcal/mol). These conclusions are supported by results obtained for a second test set of complex systems containing 70 molecules, including charged, strained, polyhalogenated, hypervalent, and large unsaturated species (e.g. C60). For this set, DFT-B3LYP performs badly (MAD of 8.6 kcal/mol) with many errors >10-20 kcal/mol while the spin-component-scaled MP2 methods are still very accurate (MAD of 2.8 and 3.7 kcal/mol, respectively). DFT-B3LYP shows an obvious tendency to underestimate molecular stability as the system size increases. Out of six density functionals tested, the hybrid functional PBE0 performs best. All in all, the SCS-MP2 method, together with large AO basis sets, clearly outperforms current DFT approaches and seems to be the most accurate quantum chemical model that routinely can predict the thermodynamic properties of large main group compounds.     

Electronic structure study of seven-coordinate first-row transition metal complexes derived from 1,10-diaza-15-crown-5: a successful marriage of theory with experiment

A detailed study of the electronic structure of seven-coordinate Mn(II), Co(II), and Ni(II) complexes with the lariat ether N,N'-bis(2-aminobenzyl)-1,10-diaza-15-crown-5 (L(1)) is presented. These complexes represent new examples of structurally characterized seven-coordinate (pentagonal bipyramidal) complexes for the Mn(II), Co(II), and Ni(II) ions. The X-ray crystal structures of the Mn(II) and Co(II) complexes show C(2) symmetries for the [M(L(1))](2+) cations, whereas the structures of the Ni(II) complexes show a more distorted coordination environment. The magnetic properties of the Mn(II) complex display a characteristic Curie law, whereas those of the Co(II) and Ni(II) ions show the occurrence of zero-field splitting of the S = 3/2 and 1 ground states, respectively. Geometry optimizations of the [M(L(1))](2+) systems (M = Mn, Co, or Ni) at the DFT (B3LYP) level of theory provide theoretical structures in good agreement with the experimental data. Electronic structure calculations predict a similar ordering of the metal-based beta spin frontier MO for the Mn(II) and Co(II) complexes. This particular ordering of the frontier MO leads to a pseudodegenerate ground state for the d(8) Ni(II) ion. The distortion of the C(2) symmetry in [Ni(L(1))](2+) is consistent with a Jahn-Teller effect that removes this pseudodegeneracy. Our electronic structure calculations predict that the binding strength of L(1) should follow the trend Co(II) approximately Mn(II) > Ni(II), in agreement with experimental data obtained from spectrophotometric titrations.     

Structural studies of nonclassical cyclobutylmethyl cations by the ab initio method

Ab initio calculations at the MP2/cc-pVTZ level show that the cyclobutylmethyl cation is a nonclassical sigma-delocalized species, which is distinct from the global minimum C2-symmetric cyclopentyl cation (Schleyer, P. v. R.; Carneiro, J. W. de M.; Koch, W.; Raghavachari, K. J. Am. Chem. Soc. 1989, 111, 5475). Relatively lower level DFT calculations, on the other hand, show that the primary cyclobutylmethyl cation spontaneously collapses into the cyclopentyl cation (Prakash, G. K. S.; Reddy, V. P.; Rasul, G.; Casanova, J.; Olah, G. A. J. Am. Chem. Soc. 1998, 120, 13362). The secondary 1-cyclobutylethyl cation is also a nonclassical carbocation, as shown by calculations at the MP2/cc-pVTZ level. Two structures having energy minima are identified for the latter cation on the potential energy surface. The conformer in which the methyl group is in the exo orientation is a global minimum and is favored over the corresponding endo conformer by 1.2 kcal/mol at the MP2/cc-pVTZ//MP2/cc-pVTZ +ZPE level of calculations. The tertiary 1-cyclobutyl-1-methylethyl cation, at this level of calculations, also involves substantial nonclassical sigma-delocalization, showing that the nonclassical stabilization is more important for cyclobutylmethyl cations relative to the cyclopropylmethyl cations. The 13C NMR chemical shifts obtained from GIAO-CCSD(T)/tzp/dz calculations further substantiate the nonclassical structures for these carbocations.     

H2S键合金属(II)咪唑卟啉配合物的结构、光谱和反应活性

选取8个典型的二价金属咪唑卟啉MP(M=Ca, Mg, Zn, Cu, Ni, Fe, Co, Mn; P代表咪唑卟啉)与H₂S(L)形成轴向金属配合物(L-MP; L-MP^*-L, P^*代表卟啉), 应用轨道和自旋概念密度泛函工具, 在优化构型的基础上, 通过自然键轨道(NBO)方法和前线轨道能级研究了它们的分子结构、光谱性质和反应活性. 模拟结果揭示L-MP和L-MP^*-L结构、光谱及其反应活性不同于其前体MP. MP排斥钙而选择镁; L对MP的结构影响较少, 与咪唑铁卟啉(FeP)能形成最稳定的单轴配合物(L-FeP), 其电子吸收光谱较前体FeP有显著的变化; 铁的亲核Fukui轨道指数值(f_Fe⁺)大于其他原子的Fukui指数, 且发生符号改变. 铁体系的自旋极化Fukui密度图也支持以上结论. 在这些典型的赤道键合配合物中, 金属M与N(S)原子之间的二级微扰相互作用能、自然电荷以及概念密度泛函指数等存在一系列线性关系. 以上结果可为理解内源性H₂S与血管性物质的相互作用机理提供启示.     

Characterizing the later 3d cyanides: the submillimeter spectrum of CoCN(X 3Phi i)

The pure rotational spectrum of the CoCN radical has been recorded in the frequency range 350-500 GHz using direct absorption techniques. This study is the first spectroscopic observation of this molecule by any experimental technique. Spectra of Co (13)CN have been measured as well. These data indicate that this species is linear in its ground electronic state and has the cyanide, as opposed to the isocyanide, geometry. The ground state term has been assigned as (3)Phi(i), based on the measurement of three spin components (Omega=4, 3, and 2) and in analogy to other isovalent cobalt-bearing species. Hyperfine splittings resulting from the (59)Co nuclear spin of I=7/2 were observed in every transition, each of which exhibited an octet pattern. For the lowest energy spin component, Omega=4, vibrational satellite features were also identified arising from the first quantum of the Co-C (v(1)=1) stretch and the v(2)=1 and v(2)=2 quanta of the bending mode, which were split by Renner-Teller interactions. The ground state measurements of CoCN were analyzed with a case a(beta) Hamiltonian, establishing rotational, fine structure, and hyperfine parameters. The vibrational and Co (13)CN spectra for the Omega=4 component were fit as well. An r(0) structure was also calculated, providing estimates of the Co-C and C-N bond distances, based on the Omega=4 transitions. CoCN is the fourth molecule in the 3d transition metal series to exhibit the linear cyanide structure, along with the Zn, Cu, and Ni analogs. The preference for this geometry, as opposed to the isocyanide form, may indicate a greater degree of covalent bonding in these species.     

Manganese,Iron, Cobalt,and Nickel Oxo‐, Peroxo‐, and Superoxoclusters: A Density Functional Theory Study

The 3d‐transition‐metal dioxo‐, peroxo‐, and superoxoclusters with the general composition MO₂, M(O₂), and MOO (M=Mn, Fe, Co, and Ni) were studied by DFT by the B1LYP functional. The dioxides in their ground states represent the global minima for the M+O₂ system. Both ground‐state dioxides and the lowest‐energy peroxides are in their (d‐only) highest spin states. The ⁶A₁ state of Co(O₂) exceeds the d‐only spin‐multiplicity value (quartet), being nearly isoenergetic with the ⁴A₁ state of Co(O₂). The energy gain on transforming the peroxides to the corresponding dioxides decreases in the order Mn(O₂)>Fe(O₂)>Co(O₂)>Ni(O₂) and varies in the range 0.27–1.8 eV. The dissociation energy to M+O₂ for all studied peroxides is less than 1 eV being the lowest (0.47 eV) for Mn(O₂). The Mn and Fe peroxides need less than 0.3 eV to rupture one of the MO bonds to form the corresponding superoxide. Mn and Fe superoxides are less stable than the corresponding peroxides; the superoxide of Co is more stable than its peroxide, while Ni superoxide is unstable—its energy is above the limit of dissociation to Ni+O₂. According to the electrostatic potential maps, the oxygen atoms in the peroxides are more nucleophilic than those in the dioxides and superoxides, in which the terminal oxygen atom is more nucleophilic than the M‐bonded oxygen atom. This result differs from the expectations based on charge‐distribution analysis.