Shentan Chen

Moving Protons with Pendant Amines: Proton Mobility in a Nickel Catalyst for Oxidation of Hydrogen

Proton transport is ubiquitous in chemical and biological processes, including the reduction of dioxygen to water, the reduction of CO(2) to formate, and the production/oxidation of hydrogen. In this work we describe intramolecular proton transfer between Ni and positioned pendant amines for the hydrogen oxidation electrocatalyst [Ni(P(Cy)(2)N(Bn)(2)H)(2)](2+) (P(Cy)(2)N(Bn)(2) = 1,5-dibenzyl-3,7-dicyclohexyl-1,5-diaza-3,7-diphosphacyclooctane). Rate constants are determined by variable-temperature one-dimensional NMR techniques and two-dimensional EXSY experiments. Computational studies provide insight into the details of the proton movement and energetics of these complexes. Intramolecular proton exchange processes are observed for two of the three experimentally observable isomers of the doubly protonated Ni(0) complex, [Ni(P(Cy)(2)N(Bn)(2)H)(2)](2+), which have N-H bonds but no Ni-H bonds. For these two isomers, with pendant amines positioned endo to the Ni, the rate constants for proton exchange range from 10(4) to 10(5) s(-1) at 25 °C, depending on isomer and solvent. No exchange is observed for protons on pendant amines positioned exo to the Ni. Analysis of the exchange as a function of temperature provides a barrier for proton exchange of ΔG(‡) = 11-12 kcal/mol for both isomers, with little dependence on solvent. Density functional theory calculations and molecular dynamics simulations support the experimental observations, suggesting metal-mediated intramolecular proton transfers between nitrogen atoms, with chair-to-boat isomerizations as the rate-limiting steps. Because of the fast rate of proton movement, this catalyst may be considered a metal center surrounded by a cloud of exchanging protons. The high intramolecular proton mobility provides information directly pertinent to the ability of pendant amines to accelerate proton transfers during catalysis of hydrogen oxidation. These results may also have broader implications for proton movement in homogeneous catalysts and enzymes in general, with specific implications for the proton channel in the Ni-Fe hydrogenase enzyme.

Hydrogen oxidation catalysis by a nickel diphosphine complex with pendant tert-butyl amines

A bis-diphosphine nickel complex with tert-butyl functionalized pendant amines [Ni(P(Cy)(2)N(t-Bu)(2))(2)](2+) has been synthesized. It is a highly active electrocatalyst for the oxidation of hydrogen in the presence of base. The turnover rate of 50 s(-1) under 1.0 atm H(2) at a potential of -0.77 V vs. the ferrocene couple is 5 times faster than the rate reported heretofore for any other synthetic molecular H(2) oxidation catalyst.

Synthesis, Characterization, and Reactivity of Fe Complexes Containing Cyclic Diazadiphosphine Ligands: The Role of the Pendant Base in Heterolytic Cleavage of H2

The iron complexes CpFe(P(Ph)(2)N(Bn)(2))Cl (1-Cl), CpFe(P(Ph)(2)N(Ph)(2))Cl (2-Cl), and CpFe(P(Ph)(2)C(5))Cl (3-Cl)(where P(Ph)(2)N(Bn)(2) is 1,5-dibenzyl-1,5-diaza-3,7-diphenyl-3,7-diphosphacyclooctane, P(Ph)(2)N(Ph)(2) is 1,3,5,7-tetraphenyl-1,5-diaza-3,7-diphosphacyclooctane, and P(Ph)(2)C(5) is 1,4-diphenyl-1,4-diphosphacycloheptane) have been synthesized and characterized by NMR spectroscopy, electrochemical studies, and X-ray diffraction. These chloride derivatives are readily converted to the corresponding hydride complexes [CpFe(P(Ph)(2)N(Bn)(2))H (1-H), CpFe(P(Ph)(2)N(Ph)(2))H (2-H), CpFe(P(Ph)(2)C(5))H (3-H)] and H(2) complexes [CpFe(P(Ph)(2)N(Bn)(2))(H(2))]BAr(F)(4), [1-H(2)]BAr(F)(4), (where BAr(F)(4) is B[(3,5-(CF(3))(2)C(6)H(3))(4)](-)), [CpFe(P(Ph)(2)N(Ph)(2))(H(2))]BAr(F)(4), [2-H(2)]BAr(F)(4), and [CpFe(P(Ph)(2)C(5))(H(2))]BAr(F)(4), [3-H(2)]BAr(F)(4), as well as [CpFe(P(Ph)(2)N(Bn)(2))(CO)]BAr(F)(4), [1-CO]Cl. Structural studies are reported for [1-H(2)]BAr(F)(4), 1-H, 2-H, and [1-CO]Cl. The conformations adopted by the chelate rings of the P(Ph)(2)N(Bn)(2) ligand in the different complexes are determined by attractive or repulsive interactions between the sixth ligand of these pseudo-octahedral complexes and the pendant N atom of the ring adjacent to the sixth ligand. An example of an attractive interaction is the observation that the distance between the N atom of the pendant amine and the C atom of the coordinated CO ligand for [1-CO]BAr(F)(4) is 2.848 Å, considerably shorter than the sum of the van der Waals radii of N and C atoms. Studies of H/D exchange by the complexes [1-H(2)](+), [2-H(2)](+), and [3-H(2)](+) carried out using H(2) and D(2) indicate that the relatively rapid H/D exchange observed for [1-H(2)](+) and [2-H(2)](+) compared to [3-H(2)](+) is consistent with intramolecular heterolytic cleavage of H(2) mediated by the pendant amine. Computational studies indicate a low barrier for heterolytic cleavage of H(2). These mononuclear Fe(II) dihydrogen complexes containing pendant amines in the ligands mimic crucial features of the distal Fe site of the active site of the [FeFe]-hydrogenase required for H-H bond formation and cleavage.

Homogeneous Ni catalysts for H2 oxidation and production: An assessment of theoretical methods, from density functional theory to post hartree-fock correlated wave-function theory

A systematic assessment of theoretical methods applicable to the accurate characterization of catalytic cycles of homogeneous catalysts for H(2) oxidation and evolution is reported. The key elementary steps involve heterolytic cleavage of the H-H bond and formation/cleavage of Ni-H and N-H bonds. In the context of density functional theory (DFT), we investigated the use of functionals in the generalized gradient approximation (GGA) as well as hybrid functionals. We compared the results with wave-function theories based on perturbation theory (MP2 and MP4) and on coupled-cluster expansions [CCD, CCSD, and CCSD(T)]. Our findings indicate that DFT results based on Perdew correlation functionals are in semiquantitative agreement with the CCSD(T) results, with deviations of only a few kilocalories/mole. On the other hand, the B3LYP functional is not even in qualitative agreement with CCSD(T). Surprisingly, the MP2 results are found to be extremely poor, in particular for the diproton Ni(0) and dihydride Ni(IV) species on the reaction potential energy surface. The Hartree-Fock reference wave function in MP2 theory gives a poor reference state description for these states that are electron rich on Ni, giving rise to erroneous MP2 energies. We present a detailed potential-energy diagram for the oxidation of H(2) by these catalysts after accounting for the effects of solvation, as modeled by a polarizable continuum, and of free energy estimated at the harmonic level of theory.

Dinitrogen Reduction by a Chromium(0) Complex Supported by a 16-Membered Phosphorus Macrocycle

We report a rare example of a Cr-N2 complex supported by a 16-membered phosphorus macrocycle containing pendant amine bases. Reactivity with acid afforded hydrazinium and ammonium, representing the first example of N2 reduction by a Cr-N2 complex. Computational analysis examined the thermodynamically favored protonation steps of N2 reduction with Cr leading to the formation of hydrazine.

Toward molecular catalysts by computer.

CONSPECTUS: Rational design of molecular catalysts requires a systematic approach to designing ligands with specific functionality and precisely tailored electronic and steric properties. It then becomes possible to devise computer protocols to design catalysts by computer. In this Account, we first review how thermodynamic properties such as redox potentials (E°), acidity constants (pKa), and hydride donor abilities (ΔGH(-)) form the basis for a framework for the systematic design of molecular catalysts for reactions that are critical for a secure energy future. We illustrate this for hydrogen evolution and oxidation, oxygen reduction, and CO conversion, and we give references to other instances where it has been successfully applied. The framework is amenable to quantum-chemical calculations and conducive to predictions by computer. We review how density functional theory allows the determination and prediction of these thermodynamic properties within an accuracy relevant to experimentalists (∼0.06 eV for redox potentials, ∼1 pKa unit for pKa values, and 1-2 kcal/mol for hydricities). Computation yielded correlations among thermodynamic properties as they reflect the electron population in the d shell of the metal center, thus substantiating empirical correlations used by experimentalists. These correlations point to the key role of redox potentials and other properties (pKa of the parent aminium for the proton-relay-based catalysts designed in our laboratory) that are easily accessible experimentally or computationally in reducing the parameter space for design. These properties suffice to fully determine free energies maps and profiles associated with catalytic cycles, i.e., the relative energies of intermediates. Their prediction puts us in a position to distinguish a priori between desirable and undesirable pathways and mechanisms. Efficient catalysts have flat free energy profiles that avoid high activation barriers due to low- and high-energy intermediates. The criterion of a flat energy profile can be mathematically resolved in a functional in the reduced parameter space that can be efficaciously calculated by means of the correlation expressions. Optimization of the functional permits the prediction by computer of design points for optimum catalysts. Specifically, the optimization yields the values of the thermodynamic properties for efficient (high rate and low overpotential) catalysts. We are on the verge of design of molecular electrocatalysts by computer. Future efforts must focus on identifying actual ligands that possess these properties. We believe that this can also be achieved through computation, using Taft-like relationships linking molecular composition and structure with electron-donating ability and steric effects. We note also that the approach adopted here of using free energy maps to decipher catalytic pathways and mechanisms does not account for kinetic barriers associated with elementary steps along the catalytic pathway, which may make thermodynamically accessible intermediates kinetically inaccessible. Such an extension of the approach will require further computations that, however, can take advantage of Polanyi-like linear free energy relationships linking activation barriers and reaction free energies.

A rare terminal dinitrogen complex of chromium.

Cis and trans-Cr-N(2) complexes supported by the diphosphine ligand P(Ph)(2)N(Bn)(2) have been prepared. Positioned pendant amines in the second coordination sphere influence the thermodynamically preferred geometric isomer. Electronic structure calculations indicate negligible Cr-N(2) back-bonding; rather, electronic polarization of N(2) ligand is thought to stabilize Cr-N(2) binding.

Protonation of Ferrous Dinitrogen Complexes Containing a Diphosphine Ligand with a Pendent Amine

The addition of acids to ferrous dinitrogen complexes [FeX(N2)(P(Et)N(Me)P(Et))(dmpm)](+) (X = H, Cl, or Br; P(Et)N(Me)P(Et) = Et2PCH2N(Me)CH2PEt2; and dmpm = Me2PCH2PMe2) gives protonation at the pendent amine of the diphosphine ligand rather than at the dinitrogen ligand. This protonation increased the νN2 band of the complex by 25 cm(-1) and shifted the Fe(II/I) couple by 0.33 V to a more positive potential. A similar IR shift and a slightly smaller shift of the Fe(II/I) couple (0.23 V) was observed for the related carbonyl complex [FeH(CO)(P(Et)N(Me)P(Et))(dmpm)](+). [FeH(P(Et)N(Me)P(Et))(dmpm)](+) was found to bind N2 about three times more strongly than NH3. Computational analysis showed that coordination of N2 to Fe(II) centers increases the basicity of N2 (vs free N2) by 13 and 20 pKa units for the trans halides and hydrides, respectively. Although the iron center increases the basicity of the bound N2 ligand, the coordinated N2 is not sufficiently basic to be protonated. In the case of ferrous dinitrogen complexes containing a pendent methylamine, the amine site was determined to be the most basic site by 30 pKa units compared to the N2 ligand. The chemical reduction of these ferrous dinitrogen complexes was performed in an attempt to increase the basicity of the N2 ligand enough to promote proton transfer from the pendent amine to the N2 ligand. Instead of isolating a reduced Fe(0)-N2 complex, the reduction resulted in isolation and characterization of HFe(Et2PC(H)N(Me)CH2PEt2)(P(Et)N(Me)P(Et)), the product of oxidative addition of the methylene C-H bond of the P(Et)N(Me)P(Et) ligand to Fe.

Protonation Studies of a Mono-Dinitrogen Complex of Chromium Supported by a 12-Membered Phosphorus Macrocycle Containing Pendant Amines

The reduction of fac-[CrCl3(P(Ph)3N(Bn)3)], (1(Cl3)), (P(Ph)3N(Bn)3 = 1,5,9-tribenzyl-3,7,11-triphenyl-1,5,9-triaza-3,7,11-triphosphacyclododecane) with Mg in the presence of dmpe (dmpe = 1,2-bis(dimethylphosphino)ethane) affords the first example of a monodinitrogen Cr(0) complex, Cr(N2)(dmpe)(P(Ph)3N(Bn)3), (2(N2)), containing a pentaphosphine coordination environment. 2(N2) is supported by a unique facially coordinating 12-membered phosphorus macrocycle containing pendant amine groups in the second coordination sphere. Treatment of 2(N2) at -78 °C with 1 equiv of [H(OEt2)2][B(C6F5)4] results in protonation of the metal center, generating the seven-coordinate Cr(II)-N2 hydride complex, [Cr(H)(N2)(dmpe)(P(Ph)3N(Bn)3)][B(C6F5)4], [2(H)(N2)](+). Treatment of 2((15)N2) with excess triflic acid at -50 °C afforded a trace amount of (15)NH4(+) from the reduction of the coordinated (15)N2 ligand (electrons originate from Cr). Electronic structure calculations were employed to evaluate the pKa values of three protonated sites of 2(N2) (metal center, pendant amine, and N2 ligand) and were used to predict the thermodynamically preferred Cr-NxHy intermediates in the N2 reduction pathway for 2(N2) and the recently published complex trans-[Cr(N2)2(P(Ph)4N(Bn)4)] upon the addition of protons and electrons.