Fundamental Studies of Carbon Oxygen Bond Activation in Nickel Diphosphine Ether Complexes. And, Metallomacrocycles as Ligands: Synthesis and Characterization of Aluminum-Bridged Bisglyoximato Complexes of Iron and Cobalt

In order to develop better catalysts for the cleavage of aryl-X bonds fundamental studies of the mechanism and individual steps of the mechanism have been investigated in detail. As the described studies are difficult at best in catalytic systems, model systems are frequently used. To study aryl-oxygen bond activation, a terphenyl diphosphine scaffold containing an ether moiety in the central arene was designed. The first three chapters of this dissertation focus on the studies of the nickel complexes supported by this diphosphine backbone and the research efforts in regards to aryl-oxygen bond activation.

Chapter 2 outlines the synthesis of a variety of diphosphine terphenyl ether ligand scaffolds. The metallation of these scaffolds with nickel is described. The reactivity of these nickel(0) systems is also outlined. The systems were found to typically undergo a reductive cleavage of the aryl oxygen bond. The mechanism was found to be a subsequent oxidative addition, β-H elimination, reductive elimination and (or) decarbonylation.

Chapter 3 presents kinetic studies of the aryl oxygen bond in the systems outlined in Chapter 2. Using a series of nickel(0) diphosphine terphenyl ether complexes the kinetics of aryl oxygen bond activation was studied. The activation parameters of oxidative addition for the model systems were determined. Little variation was observed in the rate and activation parameters of oxidative addition with varying electronics in the model system. The cause of the lack of variation is due to the ground state and oxidative addition transition state being affected similarly. Attempts were made to extend this study to catalytic systems.

Chapter 4 investigates aryl oxygen bond activation in the presence of additives. It was found that the addition of certain metal alkyls to the nickel(0) model system lead to an increase in the rate of aryl oxygen bond activation. The addition of excess Grignard reagent led to an order of magnitude increase in the rate of aryl oxygen bond activation. Similarly the addition of AlMe3 led to a three order of magnitude rate increase. Addition of AlMe₃ at -80 °C led to the formation of an intermediate which was identified by NOESY correlations as a system in which the AlMe3 is coordinated to the ether moiety of the backbone. The rates and activation parameters of aryl oxygen bond activation in the presence of AlMe3 were investigated.

The last two chapters involve the study of metalla-macrocycles as ligands. Chapter 5 details the synthesis of a variety of glyoxime backbones and diphenol precursors and their metallation with aluminum. The coordination chemistry of iron on the aluminum scaffolds was investigated. Varying the electronics of the aluminum macrocycle was found to affect the observed electrochemistry of the iron center.

Chapter 6 extends the studies of chapter 5 to cobalt complexes. The synthesis of cobalt dialuminum glyoxime metal complexes is described. The electrochemistry of the cobalt complexes was investigated. The electrochemistry was compared to the observed electrochemistry of a zinc analog to identify the redox activity of the ligand. In the presence of acid the cobalt complexes were found to electrochemically reduce protons to dihydrogen. The electronics of the ancillary aluminum ligands were found to affect the potential of proton reduction in the cobalt complexes. These potentials were compared to other diglyoximate complexes.

Sequence specific recognition and cleavage of DNA by synthetic molecules

DNA recognition is an essential biological process responsible for the regulation of cellular functions including protein synthesis and cell division and is implicated in the mechanism of action of some anticancer drugs. Studies directed towards defining the elements responsible for sequence specific DNA recognition through the study of the interactions of synthetic organic ligands with DNA are described.

DNA recognition by poly-N-methylpyrrolecarboxamides was studied by the synthesis and characterization of a series of molecules where the number of contiguous N-methylpyrrolecarboxamide units was increased from 2 to 9. The effect of this incremental change in structure on DNA recognition has been investigated at base pair resolution using affinity cleaving and MPE•Fe(II) footprinting techniques. These studies led to a quantitative relationship between the number of amides in the molecule and the DNA binding site size. This relationship is called the n + 1 rule and it states that a poly-N methylpyrrolecarboxamide molecule with n amides will bind n + 1 base pairs of DNA. This rule is consistent with a model where the carboxamides of these compounds form three center bridging hydrogen bonds between adjacent base pairs on opposite strands of the helix. The poly-N methylpyrrolecarboxamide recognition element was found to preferentially bind poly dA•poly dT stretches; however, both binding site selection and orientation were found to be affected by flanking sequences. Cleavage of large DNA is also described.

One approach towards the design of molecules that bind large sequences of double helical DNA sequence specifically is to couple DNA binding subunits of similar or diverse base pair specificity. Bis-EDTA-distamycin-fumaramide (BEDF) is an octaamide dimer of two tri-N methylpyrrolecarboxamide subunits linked by fumaramide. DNA recognition by BEDF was compared to P7E, an octaamide molecule containing seven consecutive pyrroles. These two compounds were found to recognize the same sites on pBR322 with approximately the same affinities demonstrating that fumaramide is an effective linking element for Nmethylpyrrolecarboxamide recognition subunits. Further studies involved the synthesis and characterization of a trimer of tetra-N-methylpyrrolecarboxamide subunits linked by β-alanine ((P4)_(3)E). This trimerization produced a molecule which is capable of recognizing 16 base pairs of A•T DNA, more than a turn and a half of the DNA helix.

DNA footprinting is a powerful direct method for determining the binding sites of proteins and small molecules on heterogeneous DNA. It was found that attachment of EDTA•Fe(II) to spermine creates a molecule, SE•Fe(II), which binds and cleaves DNA sequence neutrally. This lack of specificity provides evidence that at the nucleotide level polyamines recognize heterogeneous DNA independent of sequence and allows SE•Fe(II) to be used as a footprinting reagent. SE•Fe(II) was compared with two other small molecule footprinting reagents, EDTA•Fe(II) and MPE•Fe(II).

Part I. A study of the velocity structure of the earth by the use of core phases. Part II. The 1971 San Fernando earthquake series focal mechanisms and tectonics

The initial objective of Part I was to determine the nature of upper mantle discontinuities, the average velocities through the mantle, and differences between mantle structure under continents and oceans by the use of P'dP', the seismic core phase P'P' (PKPPKP) that reflects at depth d in the mantle. In order to accomplish this, it was found necessary to also investigate core phases themselves and their inferences on core structure. P'dP' at both single stations and at the LASA array in Montana indicates that the following zones are candidates for discontinuities with varying degrees of confidence: 800-950 km, weak; 630-670 km, strongest; 500-600 km, strong but interpretation in doubt; 350-415 km, fair; 280-300 km, strong, varying in depth; 100-200 km, strong, varying in depth, may be the bottom of the low-velocity zone. It is estimated that a single station cannot easily discriminate between asymmetric P'P' and P'dP' for lead times of about 30 sec from the main P'P' phase, but the LASA array reduces this uncertainty range to less than 10 sec. The problems of scatter of P'P' main-phase times, mainly due to asymmetric P'P', incorrect identification of the branch, and lack of the proper velocity structure at the velocity point, are avoided and the analysis shows that one-way travel of P waves through oceanic mantle is delayed by 0.65 to 0.95 sec relative to United States mid-continental mantle.

A new P-wave velocity core model is constructed from observed times, dt/dΔ's, and relative amplitudes of P'; the observed times of SKS, SKKS, and PKiKP; and a new mantle-velocity determination by Jordan and Anderson. The new core model is smooth except for a discontinuity at the inner-core boundary determined to be at a radius of 1215 km. Short-period amplitude data do not require the inner core Q to be significantly lower than that of the outer core. Several lines of evidence show that most, if not all, of the arrivals preceding the DF branch of P' at distances shorter than 143° are due to scattering as proposed by Haddon and not due to spherically symmetric discontinuities just above the inner core as previously believed. Calculation of the travel-time distribution of scattered phases and comparison with published data show that the strongest scattering takes place at or near the core-mantle boundary close to the seismic station.

In Part II, the largest events in the San Fernando earthquake series, initiated by the main shock at 14 00 41.8 GMT on February 9, 1971, were chosen for analysis from the first three months of activity, 87 events in all. The initial rupture location coincides with the lower, northernmost edge of the main north-dipping thrust fault and the aftershock distribution. The best focal mechanism fit to the main shock P-wave first motions constrains the fault plane parameters to: strike, N 67° (± 6°) W; dip, 52° (± 3°) NE; rake, 72° (67°-95°) left lateral. Focal mechanisms of the aftershocks clearly outline a downstep of the western edge of the main thrust fault surface along a northeast-trending flexure. Faulting on this downstep is left-lateral strike-slip and dominates the strain release of the aftershock series, which indicates that the downstep limited the main event rupture on the west. The main thrust fault surface dips at about 35° to the northeast at shallow depths and probably steepens to 50° below a depth of 8 km. This steep dip at depth is a characteristic of other thrust faults in the Transverse Ranges and indicates the presence at depth of laterally-varying vertical forces that are probably due to buckling or overriding that causes some upward redirection of a dominant north-south horizontal compression. Two sets of events exhibit normal dip-slip motion with shallow hypocenters and correlate with areas of ground subsidence deduced from gravity data. Several lines of evidence indicate that a horizontal compressional stress in a north or north-northwest direction was added to the stresses in the aftershock area 12 days after the main shock. After this change, events were contained in bursts along the downstep and sequencing within the bursts provides evidence for an earthquake-triggering phenomenon that propagates with speeds of 5 to 15 km/day. Seismicity before the San Fernando series and the mapped structure of the area suggest that the downstep of the main fault surface is not a localized discontinuity but is part of a zone of weakness extending from Point Dume, near Malibu, to Palmdale on the San Andreas fault. This zone is interpreted as a decoupling boundary between crustal blocks that permits them to deform separately in the prevalent crustal-shortening mode of the Transverse Ranges region.

I. A seismotectonic study of the Middle America Subduction Zone. II. Lithosphere and upper mantle structure of the Canadian Shield and Eastern North America

This thesis consists of two separate parts. Part I (Chapter 1) is concerned with seismotectonics of the Middle America subduction zone. In this chapter, stress distribution and Benioff zone geometry are investigated along almost 2000 km of this subduction zone, from the Rivera Fracture Zone in the north to Guatemala in the south. Particular emphasis is placed on the effects on stress distribution of two aseismic ridges, the Tehuantepec Ridge and the Orozco Fracture Zone, which subduct at seismic gaps. Stress distribution is determined by studying seismicity distribution, and by analysis of 190 focal mechanisms, both new and previously published, which are collected here. In addition, two recent large earthquakes that have occurred near the Tehuantepec Ridge and the Orozco Fracture Zone are discussed in more detail. A consistent stress release pattern is found along most of the Middle America subduction zone: thrust events at shallow depths, followed down-dip by an area of low seismic activity, followed by a zone of normal events at over 175 km from the trench and 60 km depth. The zone of low activity is interpreted as showing decoupling of the plates, and the zone of normal activity as showing the breakup of the descending plate. The portion of subducted lithosphere containing the Orozco Fracture Zone does not differ significantly, in Benioff zone geometry or in stress distribution, from adjoining segments. The Playa Azul earthquake of October 25, 1981, M_s=7.3, occurred in this area. Body and surface wave analysis of this event shows a simple source with a shallow thrust mechanism and gives M_o=1.3x10²⁷ dyne-cm. A stress drop of about 45 bars is calculated; this is slightly higher than that of other thrust events in this subduction zone. In the Tehuantepec Ridge area, only minor differences in stress distribution are seen relative to adjoining segments. For both ridges, the only major difference from adjoining areas is the infrequency or lack of occurrence of large interplate thrust events.

Part II involves upper mantle P wave structure studies, for the Canadian shield and eastern North America. In Chapter 2, the P wave structure of the Canadian shield is determined through forward waveform modeling of the phases P_nl, P, and PP. Effects of lateral heterogeneity are kept to a minimum by using earthquakes just outside the shield as sources, with propagation paths largely within the shield. Previous mantle structure studies have used recordings of P waves in the upper mantle triplication range of 15-30°; however, the lack of large earthquakes in the shield region makes compilation of a complete P wave dataset difficult. By using the phase PP, which undergoes triplications at 30-60°, much more information becomes available. The WKBJ technique is used to calculate synthetic seismograms for PP, and these records are modeled almost as well as the P. A new velocity model, designated S25, is proposed for the Canadian shield. This model contains a thick, high-Q, high-velocity lid to 165 km and a deep low-velocity zone. These features combine to produce seismograms that are markedly different from those generated by other shield structure models. The upper mantle discontinuities in S25 are placed at 405 and 660 km, with a simple linear gradient in velocity between them. Details of the shape of the discontinuities are not well constrained. Below 405 km, this model is not very different from many proposed P wave models for both shield and tectonic regions.

Chapter 3 looks in more detail at recordings of P_nl in eastern North America. First, seismograms from four eastern North American earthquakes are analyzed, and seismic moments for the events are calculated. These earthquakes are important in that they are among the largest to have occurred in eastern North America in the last thirty years, yet in some cases were not large enough to produce many good long-period teleseismic records. A simple layer-over-a-halfspace model is used for the initial modeling, and is found to provide an excellent fit for many features of the observed waveforms. The effects on P_nl of varying lid structure are then investigated. A thick lid with a positive gradient in velocity, such as that proposed for the Canadian shield in Chapter 2, will have a pronounced effect on the waveforms, beginning at distances of 800 or 900 km. P_nl records from the same eastern North American events are recalculated for several lid structure models, to survey what kinds of variations might be seen. For several records it is possible to see likely effects of lid structure in the data. However, the dataset is too sparse to make any general observations about variations in lid structure. This type of modeling is expected to be important in the future, as the analysis is extended to more recent eastern North American events, and as broadband instruments make more high-quality regional recordings available.

Electrocatalytic activity of transition metal substituted heteropolytungstates

The electrochemical and electrocatalytic behavior of a series of heteropolytungstate anions in which a tungsten atom in the well known Keggin structure has been replaced by an iron atom is described. All of the iron substituted ions exhibit a one electron reversible couple associated with the Fe³⁺ center and a pair of two electron waves attributed to electron addition and removal from the tungsten oxo framework. The pH and ionic strength effects upon the various electrochemical processes are discussed and interpreted in terms of a competition between protonation and ion pairing of the highly negatively charged ions.

The anions are efficient catalysts for the electroreduction of hydrogen peroxide. A catalytic mechanism involving a formally Fe(IV) intermediate is proposed. Pulse radiolysis experiments were employed to detect the intermediate and evaluate the rate constants for the reactions in which it is formed and decomposed. A chain mechanism for the catalytic decomposition of hydrogen peroxide in which the Fe center shuttles between the +2, +3, and +4 oxidation states is proposed to explain the non-integral stoichiometry observed for the iron substituted polytungstate catalyzed electroreduction of hydrogen peroxide.

The anions are also efficient electrocatalyst for the electrochemical conversion of nitric oxide to ammonia. The catalyzed reduction does not produce hydroxylamine as an intermediate and appears to depend upon the ability of the multiply reduced heteropolytungstates to deliver several electrons to the bound NO group in a concerted step. A valuable feature of the heteropolytungstates is the ease at which the formal potentials of the several redox couples they exhibit may be shifted by changing the identity of the central heteroatom. Exploitation of this feature provided diagnostic information that was decisive in establishing the mechanism of electrocatalytic reduction.

The iron substituted heteropolytungstates are not degraded by repeated cycling between their oxidized and reduced states. They also show superior activity compared to their unsubstituted analogues, indicating that the Fe center acts as a binding site that facilitates inner-sphere electron transfer processes. The basic electrochemistry of several other transition metal substituted Keggin ions is also described.

Electrocatalytic reduction of dioxygen at chemically modified electrodes

The kinetics of the reduction of O₂ by Ru(NH₃)₆⁺² as catalyzed by cobalt(II) tetrakis(4-N-methylpyridyl)porphyrin are described both in homogeneous solution and when the reactants are confined to Nafion coatings on graphite electrodes. The catalytic mechanism is determined and the factors that can control the total reduction currents at Nafion-coated electrodes are specified. A kinetic zone diagram for analyzing the behavior of catalyst-mediator-substrate systems at polymer coated electrodes is presented and utilized in identifying the current-limiting processes. Good agreement is demonstrated between calculated and measured reduction currents at rotating disk electrodes. The experimental conditions that will yield the optimum performance of coated electrodes are discussed, and a relationship is derived for the optimal coating thickness.

The relation between the reduction potentials of adsorbed and unadsorbed cobalt(III) tetrakis(4-N-methylpyridyl)porphyrin and those where it catalyzes the electroreduction of dioxygen is described. There is an unusually large change in the formal potential of the Co(III) couple upon the adsorption of the porphyrin on the graphite electrode surface. The mechanism in which the (inevitably) adsorbed porphyrin catalyzes the reduction of O₂ is in accord with a general mechanistic scheme proposed for most monomeric cobalt porphyrins.

Four new dimeric metalloporphyrins (prepared in the laboratory of Professor C. K. Chang) have the two porphyrin rings linked by an anthracene bridge attached to meso positions. The electrocatalytic behavior of the diporphyrins towards the reduction of O₂ at graphite electrodes has been examined for the following combination of metal centers: Co-Cu, Co-Fe, Fe-Fe, Fe-H₂. The Co-Cu diporphyrin catalyzes the reduction of O₂ to H₂O₂ but no further. The other three catalysts all exhibit mixed reduction pathways leading to both H₂O₂ and H₂O. However, the pathways that lead to H₂O do not involve H₂O₂ as an intermediate. A possible mechanistic scheme is offered to account for the observed behavior.

Advancements in jet turbulence and noise modeling: accurate one-way solutions and empirical evaluation of the nonlinear forcing of wavepackets

Jet noise reduction is an important goal within both commercial and military aviation. Although large-scale numerical simulations are now able to simultaneously compute turbulent jets and their radiated sound, lost-cost, physically-motivated models are needed to guide noise-reduction efforts. A particularly promising modeling approach centers around certain large-scale coherent structures, called wavepackets, that are observed in jets and their radiated sound. The typical approach to modeling wavepackets is to approximate them as linear modal solutions of the Euler or Navier-Stokes equations linearized about the long-time mean of the turbulent flow field. The near-field wavepackets obtained from these models show compelling agreement with those educed from experimental and simulation data for both subsonic and supersonic jets, but the acoustic radiation is severely under-predicted in the subsonic case. This thesis contributes to two aspects of these models. First, two new solution methods are developed that can be used to efficiently compute wavepackets and their acoustic radiation, reducing the computational cost of the model by more than an order of magnitude. The new techniques are spatial integration methods and constitute a well-posed, convergent alternative to the frequently used parabolized stability equations. Using concepts related to well-posed boundary conditions, the methods are formulated for general hyperbolic equations and thus have potential applications in many fields of physics and engineering. Second, the nonlinear and stochastic forcing of wavepackets is investigated with the goal of identifying and characterizing the missing dynamics responsible for the under-prediction of acoustic radiation by linear wavepacket models for subsonic jets. Specifically, we use ensembles of large-eddy-simulation flow and force data along with two data decomposition techniques to educe the actual nonlinear forcing experienced by wavepackets in a Mach 0.9 turbulent jet. Modes with high energy are extracted using proper orthogonal decomposition, while high gain modes are identified using a novel technique called empirical resolvent-mode decomposition. In contrast to the flow and acoustic fields, the forcing field is characterized by a lack of energetic coherent structures. Furthermore, the structures that do exist are largely uncorrelated with the acoustic field. Instead, the forces that most efficiently excite an acoustic response appear to take the form of random turbulent fluctuations, implying that direct feedback from nonlinear interactions amongst wavepackets is not an essential noise source mechanism. This suggests that the essential ingredients of sound generation in high Reynolds number jets are contained within the linearized Navier-Stokes operator rather than in the nonlinear forcing terms, a conclusion that has important implications for jet noise modeling.

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.