Design, Synthesis, and Study of Novel Platforms for Iron-N2 Chemistry and Photoinduced, Copper-mediated C-N Bond Formation

Several new ligand platforms designed to support iron dinitrogen chemistry have been developed. First, we report Fe complexes of a tris(phosphino)alkyl (CP^iPr₃) ligand featuring an axial carbon donor intended to conceptually model the interstitial carbide atom of the nitrogenase iron-molybdenum cofactor (FeMoco). It is established that in this scaffold, the iron center binds dinitrogen trans to the C_alkyl anchor in three structurally characterized oxidation states. Fe-C_alkyl lengthening is observed upon reduction, reflective of significant ionic character in the Fe-C_alkyl interaction. The anionic (CP^iPr₃)FeN₂^- species can be functionalized by a silyl electrophile to generate (CP^iPr₃)Fe-N₂SiR₃. This species also functions as a modest catalyst for the reduction of N₂ to NH₃. Next, we introduce a new binucleating ligand scaffold that supports an Fe(μ-SAr)Fe diiron subunit that coordinates dinitrogen (N₂-Fe(μ-SAr)Fe-N₂) across at least three oxidation states (Fe^IIFe^II, Fe^IIFe^I, and Fe^IFe^I). Despite the sulfur-rich coordination environment of iron in FeMoco, synthetic examples of transition metal model complexes that bind N₂ and also feature sulfur donor ligands remain scarce; these complexes thus represent an unusual series of low-valent diiron complexes featuring thiolate and dinitrogen ligands. The (N₂-Fe(μ-SAr)Fe-N₂) system undergoes reduction of the bound N₂ to produce NH₃ (~50% yield) and can efficiently catalyze the disproportionation of N₂H₄ to NH₃ and N₂. The present scaffold also supports dinitrogen binding concomitant with hydride as a co-ligand. Next, inspired by the importance of secondary-sphere interactions in many metalloenzymes, we present complexes of iron in two new ligand scaffolds ([SiP^NMe₃] and [SiP^iPr₂P^NMe]) that incorporate hydrogen-bond acceptors (tertiary amines) which engage in interactions with nitrogenous substrates bound to the iron center (NH₃ and N₂H₄). Cation binding is also facilitated in anionic Fe(0)-N₂ complexes. While Fe-N₂ complexes of a related ligand ([SiP^iPr₃]) lacking hydrogen-bond acceptors produce a substantial amount of ammonia when treated with acid and reductant, the presence of the pendant amines instead facilitates the formation of metal hydride species.

Additionally, we present the development and mechanistic study of copper-mediated and copper-catalyzed photoinduced C-N bond forming reactions. Irradiation of a copper-amido complex, ((m-tol)₃P)₂Cu(carbazolide), in the presence of aryl halides furnishes N-phenylcarbazole under mild conditions. The mechanism likely proceeds via single-electron transfer from an excited state of the copper complex to the aryl halide, generating an aryl radical. An array of experimental data are consistent with a radical intermediate, including a cyclization/stereochemical investigation and a reactivity study, providing the first substantial experimental support for the viability of a radical pathway for Ullmann C-N bond formation. The copper complex can also be used as a precatalyst for Ullmann C-N couplings. We also disclose further study of catalytic C_alkyl-N couplings using a CuI precatalyst, and discuss the likely role of [Cu(carbazolide)₂]^- and [Cu(carbazolide)₃]^- species as intermediates in these reactions.

Finally, we report a series of four-coordinate, pseudotetrahedral P₃Fe^II-X complexes supported by tris(phosphine)borate ([PhBP₃Fe^R]^-) and phosphiniminato X-type ligands (-N=PR'₃) that in combination tune the spin-crossover behavior of the system. Low-coordinate transition metal complexes such as these that undergo reversible spin-crossover remain rare, and the spin equilibria of these systems have been studied in detail by a suite of spectroscopic techniques.

Z-selective olefin metathesis using chelating ruthenium alkylidene catalysts

Publications about olefin metathesis will generally discuss how the discovery and development of well-defined catalysts to carry out this unique transformation have revolutionized many fields, from natural product and materials chemistry, to green chemistry and biology. However, until recently, an entire manifestation of this methodology had been inaccessible. Except for a few select examples, metathesis catalysts favor the thermodynamic trans- or E-olefin products in cross metathesis (CM), macrocyclic ring closing metathesis (mRCM), ring opening metathesis polymerization (ROMP), and many other types of reactions. Judicious choice of substrates had allowed for the direct synthesis of cis- or Z-olefins or species that could be converted upon further reaction, however the catalyst controlled synthesis of Z-olefins was not possible until very recently.

Research into the structure and stability of metallacyclobutane intermediates has led to the proposal of models to impart Z-selectivity in metathesis reactions. Having the ability to influence the orientation of metallacyclobutane substituents to cause productive formation of Z- double bonds using steric and electronic effects was highly desired. The first successful realization of this concept was by Schrock and Hoveyda et al. who synthesized monoaryloxide pyrolidine (MAP) complexes of tungsten and molybdenum that promoted Z-selective CM. The Z-selectivity of these catalysts was attributed to the difference in the size of the two axial ligands. This size difference influences the orientation of the substituents on the forming/incipient metallacyclobutane intermediate to a cis-geometry and leads to productive formation of Z-olefins. These catalysts have shown great utility in the synthesis of complicated natural product precursors and stereoregular polymers. More recently, ruthenium catalysts capable of promoting Z-selective metathesis have been reported by our group and others. This thesis will discuss the development of ruthenium-based NHC chelated Z-selective catalysts, studies probing their unique metathesis mechanism, and synthetic applications that have been investigated thus far.

Chapter 1 will focus on studies into the stability of NHC chelated complexes and the synthesis of new and improved stable chelating architectures. Chapter 2 will discuss applications of the highly active and Z-selective developed in Chapter 1, including the formation of lepidopteran female sex pheromones using olefin cross metathesis and highly Z- and highly E-macrocycles using macrocyclic ring closing metathesis and Z-selective ethenolysis. Chapter 3 will explore studies into the unique mechanism of olefin metathesis reactions catalyzed by these NHC chelated, highly Z-selective catalysts, explaining observed trends by investigating the stability of relevant, substituted metallacyclobutane intermediates.

Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

Metallogenic Provinces of the Southwestern United States and Northern Mexico

Some of the metallogenic provinces of the southwestern United States and northern Mexico are defined by the geographic distribution of trace elements in the primary sulfide minerals chalcopyrite and sphalerite. The elements investigated include antimony, arsenic, bismuth, cadmium, cobalt, gallium, germanium, indium, manganese, molybdenum, nickel, silver, tellurium, thallium, and tin. Of these elements, cobalt, gallium, germanium, indium, nickel, silver, and tin exhibit the best defined geographic distribution.

The data indicate that chalcopyrite is the preferred host for tin and perhaps molybdenum; sphalerite is the preferred host for cadmium, gallium, germanium, indium, and manganese; galena is the preferred host for antimony, bismuth, silver, tellurium, and thallium; and pyrite is the preferred host for cobalt, nickel, and perhaps arsenic. With respect to the two minerals chalcopyrite and sphalerite, antimony, arsenic, molybdenum, nickel, silver, and tin prefer chalcopyrite; and bismuth, cadmium, cobalt, gallium, germanium, indium, manganese, and thallium prefer sphalerite. This distribution probably is the result of the interaction of several factors, among which are these: the various radii of the elements, the association due to chemical similarities of the major and trace elements, and the degree of ionic versus covalent and metallic character of the metal-sulfur bonds in chalcopyrite and sphalerite. The type of deposit, according to a temperature classification, appears to be of minor importance in determining the trace element content of chalcopyrite and sphalerite.

A preliminary investigation of large single crystals of sphalerite and chalcopyrite indicates that the distribution within a single crystal of some elements such as cadmium in sphalerite and indium and silver in chalcopyrite is relatively uniform, whereas the distribution of some other elements such as cobalt and manganese in sphalerite is somewhat less uniform and the distribution of tin in sphalerite is extremely erratic. The variations in trace element content probably are due largely to variations in the composition of the fluids during the growth of the crystals, but the erratic behavior of tin in sphalerite perhaps is related to the presence of numerous cavities and inclusions in the crystal studied.

Maps of the geographic distribution of trace elements in chalcopyrite and sphalerite exhibit three main belts of greater than average trace element content, which are called the Eastern, Central, and Western belts. These belts are consistent in trend and position with a beltlike distribution of copper, gold, lead, zinc, silver, and tungsten deposits and with most of the major tectonic features. However, there appear to be no definite time relationships, for as many as four metallogenic epochs, from Precambrian to late Tertiary, are represented by ore deposits within the Central belt.

The evidence suggests that the beltlike features have a deep seated origin, perhaps in the sub-crust or outer parts of the mantle, and that the deposits within each belt might be genetically related through a beltlike compositional heterogeneity in the source regions of the ores. Hence, the belts are regarded as metallogenic provinces.

Late Transition Metals Supported by Aryl Ethers and Phenoxides Bearing Pendant Phosphines: Mechanistic Insights Relevant to Ether C-O Bond Cleavage

Terphenyl diphosphines bearing pendant ethers were prepared to provide mechanistic insight into the mechanism of activation of aryl C–O bonds with Group 9 and Group 10 transition metals. Chapters 2 and 3 of this dissertation describe the reactivity of compounds supported by the model phosphine and extension of this chemistry to heterogenous C–O bond activation.

Chapter 2 describes the synthesis and reactivity of aryl-methyl and aryl-aryl model systems. The metallation of these compounds with Ni, Pd, Pt, Co, Rh, and Ir is described. Intramolecular bond activation pathways are described. In the case of the aryl-methyl ether, aryl C–O bond activation was observed only for Ni, Rh, and Ir.

Chapter 3 outlines the reactivity of heterogenous Rh and Ir catalysts for aryl ether C–O bond cleavage. Using Rh/C and an organometallic Ir precursor, aryl ethers were treated with H2 and heat to afford products of hydrogenolysis and hydrogenation. Conditions were modified to optimize the yield of hydrogenolysis product. Hydrogenation could not be fully suppressed in these systems.

Appendix A describes initial investigations of bisphenoxyiminoquinoline dichromium compounds for selective C2H4 oligomerization to afford α-olefins. The synthesis of monometallic and bimetallic Cr complexes is described. These compounds are compared to literature examples and found to be less active and non-selective for production of α-olefins.

Appendix B describes the coordination chemistry of terphenyl diphosphines, terphenyl bisphosphinophenols, and biphenyl phosphinophenols proligands with molybdenum, cobalt, and nickel. Since their synthesis, terphenyl diphosphine molybdenum compounds have been reported to be good catalysts for the dehydrogenation of ammonia borane. Biphenyl phosphinophenols are demonstrated provide both phosphine and arene donors to transition metals while maintaining a sterically accessible coordination sphere. Such ligands may be promising in the context of the activation of other small molecules.

Appendix C contains relevant NMR spectra for the compounds presented in the preceding sections.