Some modified bifurcation problems with application to imperfection sensitivity in buckling

The branching theory of solutions of certain nonlinear elliptic partial differential equations is developed, when the nonlinear term is perturbed from unforced to forced. We find families of branching points and the associated nonisolated solutions which emanate from a bifurcation point of the unforced problem. Nontrivial solution branches are constructed which contain the nonisolated solutions, and the branching is exhibited. An iteration procedure is used to establish the existence of these solutions, and a formal perturbation theory is shown to give asymptotically valid results. The stability of the solutions is examined and certain solution branches are shown to consist of minimal positive solutions. Other solution branches which do not contain branching points are also found in a neighborhood of the bifurcation point.

The qualitative features of branching points and their associated nonisolated solutions are used to obtain useful information about buckling of columns and arches. Global stability characteristics for the buckled equilibrium states of imperfect columns and arches are discussed. Asymptotic expansions for the imperfection sensitive buckling load of a column on a nonlinearly elastic foundation are found and rigorously justified.

Preparation, characterization and intramolecular electron-transfer studies of ruthenium-modified cytochromes c

Redox-active ruthenium complexes have been covalently attached to the surface of a series of natural, semisynthetic and recombinant cytochromes c. The protein derivatives were characterized by a variety of spectroscopic techniques. Distant Fe^(2+) - Ru^(3+) electronic couplings were extracted from intramolecular electron-transfer rates in Ru(bpy)_2(im)HisX (where X= 33, 39, 62, and 72) derivatives of cyt c. The couplings increase according to 62 (0.0060) < 72 (0.057) < 33 (0.097) < 39 (0.11 cm^(-1)); however, this order is incongruent with histidine to heme edge-edge distances [62 (14.8) > 39 (12.3) > 33 (11.1) > =72 (8.4 Å)]. These results suggest the chemical nature of the intervening medium needs to be considered for a more precise evaluation of couplings. The rates (and couplings) correlate with the lengths of a-tunneling pathways comprised of covalent bonds, hydrogen bonds and through-space jumps from the histidines to the heme group. Space jumps greatly decrease couplings: one from Pro71 to Met80 extends the σ-tunneling length of the His72 pathway by roughly 10 covalent bond units. Experimental couplings also correlate well with those calculated using extended Hiickel theory to evaluate the contribution of the intervening protein medium.

Two horse heart cyt c variants incorporating the unnatural amino acids (S)-2- amino-3-(2,2'-bipyrid-6-yl)-propanoic acid (6Bpa) and (S)-2-amino-3-(2,2'-bipyrid-4-yl)propanoic acid ( 4Bpa) at position 72 have been prepared using semisynthetic protocols. Negligible perturbation of the protein structure results from this introduction of unnatural amino acids. Redox-active Ru(2,2'-bipyridine)_2^(2+) binds to 4Bpa72 cyt c but not to the 6Bpa protein. Enhanced ET rates were observed in the Ru(bpy)_2^(2+)-modified 4Bpa72 cyt c relative to the analogous His72 derivative. The rapid (< 60 nanosecond) photogeneration of ferrous Ru-modified 4Bpa72 cyt c in the conformationally altered alkaline state demonstrates that laser-induced ET can be employed to study submicrosecond protein-folding events.

Fracture of materials undergoing solid-solid phase transformation

A large number of technologically important materials undergo solid-solid phase transformations. Examples range from ferroelectrics (transducers and memory devices), zirconia (Thermal Barrier Coatings) to nickel superalloys and (lithium) iron phosphate (Li-ion batteries). These transformations involve a change in the crystal structure either through diffusion of species or local rearrangement of atoms. This change of crystal structure leads to a macroscopic change of shape or volume or both and results in internal stresses during the transformation. In certain situations this stress field gives rise to cracks (tin, iron phosphate etc.) which continue to propagate as the transformation front traverses the material. In other materials the transformation modifies the stress field around cracks and effects crack growth behavior (zirconia, ferroelectrics). These observations serve as our motivation to study cracks in solids undergoing phase transformations. Understanding these effects will help in improving the mechanical reliability of the devices employing these materials.

In this thesis we present work on two problems concerning the interplay between cracks and phase transformations. First, we consider the directional growth of a set of parallel edge cracks due to a solid-solid transformation. We conclude from our analysis that phase transformations can lead to formation of parallel edge cracks when the transformation strain satisfies certain conditions and the resulting cracks grow all the way till their tips cross over the phase boundary. Moreover the cracks continue to grow as the phase boundary traverses into the interior of the body at a uniform spacing without any instabilities. There exists an optimal value for the spacing between the cracks. We ascertain these conclusion by performing numerical simulations using finite elements.

Second, we model the effect of the semiconducting nature and dopants on cracks in ferroelectric perovskite materials, particularly barium titanate. Traditional approaches to model fracture in these materials have treated them as insulators. In reality, they are wide bandgap semiconductors with oxygen vacancies and trace impurities acting as dopants. We incorporate the space charge arising due the semiconducting effect and dopant ionization in a phase field model for the ferroelectric. We derive the governing equations by invoking the dissipation inequality over a ferroelectric domain containing a crack. This approach also yields the driving force acting on the crack. Our phase field simulations of polarization domain evolution around a crack show the accumulation of electronic charge on the crack surface making it more permeable than was previously believed so, as seen in recent experiments. We also discuss the effect the space charge has on domain formation and the crack driving force.

Sterochemically modified polyamides for recognition in the minor groove of DNA

The design of synthetic molecules that recognize specific sequences of DNA is an ongoing challenge in molecular medicine. Cell-permeable small molecules targeting predetermined DNA sequences offer a potential approach for offsetting the abnormal effects of misregulated gene-expression. Over the past twenty years, Professor Peter B. Dervan has developed a set of pairing rules for the rational design of minor groove binding polyamides containing pyrrole (Py), imidazole (Im), and hydroxypyrrole (Hp). Polyamides have illustrated the capability to permeate cells and inhibit transcription of specific genes in vivo. This provides impetus to identify structural elements that expand the repetoire of polyamide motifs with recognition properties comparable to naturally occurring DNA binding proteins. Through the introduction of chiral amino acids, we have developed chiral polyamides with stereochemically regulated binding characteristics. In addition, chiral substituents have facilitated the development of new polyamide motifs that broaden binding site sizes targetable by this class of ligands.

Synthesis and spectroscopy of ruthenium-modified nucleic acids

Redox-active probes are designed and prepared for use in DNA-mediated electron transfer studies. These probes consist of ruthenium(II) complexes bound to nucleosides that possess metal-binding ligands. Low- and high-potential oxidants are synthesized from these modified nucleosides and display reversible one-electron electrochemical behavior. The ruthenium-modified nucleosides exhibit distinct charge-transfer transitions in the visible region that resemble those of appropriate model complexes. Resonance Raman and time-resolved emission spectroscopy are used to characterize the nature of these transitions.

The site-specific incorporation of these redox-active probes into oligonucleotides is explored using post-synthetic modification and solid-phase synthetic methods. The preparation of the metal-binding nucleosides, their incorporation into oligonucleotides, and characterization of the resulting oligonucleotides is described. Because the insertion of these probes into modified oligonucleotides using post-synthetic modification is unsuccessful, solid-phase synthetic methods are explored. These efforts lead to the first report of 3'-metallated oligonucleotides prepared completely by automated solid-phase synthesis. Preliminary efforts to prepare a bis-metallated oligonucleotide by automated synthesis are described.

The electrochemical, absorption, and emissive features of the ruthenium-modified oligonucleotides are unchanged from those of the precursor metallonucleoside. The absence of any change in these properties upon incorporation into oligonucleotides and subsequent hybridization suggests that the incorporated ruthenium(II) complex is a valuable probe for DNA-mediated electron transfer studies.

Interaction of aqueous inorganic orthophosphate and phosphate rock

An experimental study was made of the interaction of phosphate rock and aqueous inorganic orthophosphate, calcium, and hydroxyl ions. A model of the reaction was developed by observing electron diffraction patterns in conjunction with concentration changes of chemical components. The model was applied in explaining the performance of batch suspensions of powdered phosphate rock and packed columns of granular phosphate rock. In both cases the reaction consisted initially of a rapid nucleation phase that occurred in a time period of minutes. In the batch system the calcium phosphate nuclei then ripened into larger micro-crystals of hydroxyapatite, which eventually became indistinguishable from the original phosphate rock surface. During column operation the high supersaturation ratio that existed after the rapid nucleation phase resulted in a layer of small nuclei that covered a slowly growing hydroxyapatite crystal.

The column steady-state rate constant was found to increase with increasing temperature, pH, and fluoride concentration, and to decrease with increasing concentrations of magnesium sulfate, ammonium chloride, and bicarbonate ion.

An engineering feasibility study indicated that, based on economic considerations, nucleation of apatite on phosphate rock ore has a potential use as a wastewater phosphate removal treatment process.

New insights into the formation and modification of carbonate-bearing minerals and methane gas in geological systems using multiply substituted isotopologues

This thesis describes the use of multiply-substituted stable isotopologues of carbonate minerals and methane gas to better understand how these environmentally significant minerals and gases form and are modified throughout their geological histories. Stable isotopes have a long tradition in earth science as a tool for providing quantitative constraints on how molecules, in or on the earth, formed in both the present and past. Nearly all studies, until recently, have only measured the bulk concentrations of stable isotopes in a phase or species. However, the abundance of various isotopologues within a phase, for example the concentration of isotopologues with multiple rare isotopes (multiply substituted or 'clumped' isotopologues) also carries potentially useful information. Specifically, the abundances of clumped isotopologues in an equilibrated system are a function of temperature and thus knowledge of their abundances can be used to calculate a sample’s formation temperature. In this thesis, measurements of clumped isotopologues are made on both carbonate-bearing minerals and methane gas in order to better constrain the environmental and geological histories of various samples.

Clumped-isotope-based measurements of ancient carbonate-bearing minerals, including apatites, have opened up paleotemperature reconstructions to a variety of systems and time periods. However, a critical issue when using clumped-isotope based measurements to reconstruct ancient mineral formation temperatures is whether the samples being measured have faithfully recorded their original internal isotopic distributions. These original distributions can be altered, for example, by diffusion of atoms in the mineral lattice or through diagenetic reactions. Understanding these processes quantitatively is critical for the use of clumped isotopes to reconstruct past temperatures, quantify diagenesis, and calculate time-temperature burial histories of carbonate minerals. In order to help orient this part of the thesis, Chapter 2 provides a broad overview and history of clumped-isotope based measurements in carbonate minerals.

In Chapter 3, the effects of elevated temperatures on a sample’s clumped-isotope composition are probed in both natural and experimental apatites (which contain structural carbonate groups) and calcites. A quantitative model is created that is calibrated by the experiments and consistent with the natural samples. The model allows for calculations of the change in a sample’s clumped isotope abundances as a function of any time-temperature history.

In Chapter 4, the effects of diagenesis on the stable isotopic compositions of apatites are explored on samples from a variety of sedimentary phosphorite deposits. Clumped isotope temperatures and bulk isotopic measurements from carbonate and phosphate groups are compared for all samples. These results demonstrate that samples have experienced isotopic exchange of oxygen atoms in both the carbonate and phosphate groups. A kinetic model is developed that allows for the calculation of the amount of diagenesis each sample has experienced and yields insight into the physical and chemical processes of diagenesis.

The thesis then switches gear and turns its attention to clumped isotope measurements of methane. Methane is critical greenhouse gas, energy resource, and microbial metabolic product and substrate. Despite its importance both environmentally and economically, much about methane’s formational mechanisms and the relative sources of methane to various environments remains poorly constrained. In order to add new constraints to our understanding of the formation of methane in nature, I describe the development and application of methane clumped isotope measurements to environmental deposits of methane. To help orient the reader, a brief overview of the formation of methane in both high and low temperature settings is given in Chapter 5.

In Chapter 6, a method for the measurement of methane clumped isotopologues via mass spectrometry is described. This chapter demonstrates that the measurement is precise and accurate. Additionally, the measurement is calibrated experimentally such that measurements of methane clumped isotope abundances can be converted into equivalent formational temperatures. This study represents the first time that methane clumped isotope abundances have been measured at useful precisions.

In Chapter 7, the methane clumped isotope method is applied to natural samples from a variety of settings. These settings include thermogenic gases formed and reservoired in shales, migrated thermogenic gases, biogenic gases, mixed biogenic and thermogenic gas deposits, and experimentally generated gases. In all cases, calculated clumped isotope temperatures make geological sense as formation temperatures or mixtures of high and low temperature gases. Based on these observations, we propose that the clumped isotope temperature of an unmixed gas represents its formation temperature — this was neither an obvious nor expected result and has important implications for how methane forms in nature. Additionally, these results demonstrate that methane-clumped isotope compositions provided valuable additional constraints to studying natural methane deposits.

Electron transfer in chemically and genetically modified myoglobins

The temperature dependences of the reduction potentials (E^o') of wildtype human myoglobin (Mb) and three site-directed mutants have been measured by using thin-layer spectroelectrochemistry. Residue Val68, which is in van der Waals contact with the heme in Mb, has been replaced by Glu, Asp, and Asn. At pH 7.0, reduction of the heme iron (III) in the former two proteins is accompanied by uptake of a proton by the protein. The changes in E^o', and the standard entropy (ΔS^o') and enthalpy (ΔH^o') of reduction in the mutant proteins were determined relative to values for wild-type; the change in E^o' at 25°C was about -200 millivolts for the Glu and Asp mutants, and about -80 millivolts for the Asn mutant. Reduction of Fe(III) to Fe(II) in the Glu and Asp mutants is accompanied by uptake of a proton. These studies demonstrate that Mb can tolerate substitution of a buried hydrophobic group by potentially charged and polar residues, and that such amino acid replacements can lead to substantial changes in the redox thermodynamics of the protein.

Through analysis of the temperature dependence and shapes of NMR dispersion signals, it is determined that a water molecule is bound to the sixth coordination site of the ferric heme in the Val68Asp and in the Val68Asn recombinant proteins while the carboxyl group of the sidechain of Glu68 occupies this position in Val68Glu. The relative rhombic distortions in the ESR spectra of these mutant proteins combined with H₂¹⁷O and spin interconversion experiments performed on them confirm the conclusions of the NMRD study.

The rates of intramolecular electron transfer (ET) of (NH₃)₅Ru-His48 (Val68Asp, His81GIn, Cys110AIa)Mb and (NH₃)₅Ru-His48 (Val68GIu,His81GIn,Cys110Ala)Mb were measured to be .85(3)s^-1 and .30(2)s^-1, respectively. This data supports the hypothesis that entropy of 111 reduction and reorganization energy of ET are inversely related. The rates of forward and reverse ET for (NH₃)₅ Ru-His48 (Val68GIu, His81 GIn, Cys110AIa)ZnMb -7.2(5)•10⁴s^-1and 1.4(2)•10⁵s^-1, respectively- demonstrate that the placement of a highly polar residue nearby does not significantly change the reorganization energy of the photoactive Zn porphyrin.

The distal histidine imidazoles of (NH₃)₄isnRu-His48 SWMb and (NH₃)₅Ru-His48 SWMb were cyanated with BrCN. The intramolecular ET rates of these BrCN-modified Mb derivatives are 5.5(6)s^-1 and 3.2(5)s^-1, respectively. These respective rates are 20 and 10 times faster than those of their noncyanated counterparts after the differences in ET rate from driving force are scaled according to the Marcus equation. This increase in ET rate of the cyanated Mb derivatives is attributed to lower reorganization energy since the cyanated Mb heme is pentacoordinate in both oxidation states; whereas, the native Mb heme loses a water molecule upon reduction so that it changes from six to five coordinate. The reorganization energy from Fe-OH₂ dissociation is estimated to be .2eV. This conclusion is used to reconcile data from previous experiments in our lab. ET in photoactive porphyrin-substituted myoglobins proceed faster than predicted by Marcus Theory when it is assumed that the only difference in ET parameters between photoactive porphyrins and native heme systems is driving force. However, the data can be consistently fit to Marcus Theory if one corrects for the smaller reorganization in the photoactive porphyrin systems since they do not undergo a coordination change upon ET.

Finally, the intramolecular ET rate of (NH₃)₄isnRu-His48 SWMb was measured to be 3.0(4)s^-1. This rate is within experimental error of that for (NH₃)₄pyrRu-His48 SWMb even though the former has 80mV more driving force. One likely possibility for this observation is that the tetraamminepyridineruthenium group undergoes less reorganization upon ET than the tetraammineisonicotinamideruthenium group. Moreover, analysis of the (NH₃)₄isnRu-His48 SWMb experimental system gives a likely explanation of why ET was not observed previously in (NH₃)₄isnRu-Cytochrome C.

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.