^(18)O/^(16)O and D/H studies on the interactions between heated meteoric ground waters and igneous intrusions: Western San Juan Mountains, Colorado and the Isle of Skye, Scotland

Detailed oxygen, hydrogen and carbon isotope studies have been carried out on igneous and metamorphic rocks of the Stony Mountain complex, Colorado, and the Isle of Skye, Scotland, in order to better understand the problems of hydrothermal meteoric water-rock interaction.

The Tertiary Stony Mountain stock (~1.3 km in diameter), is composed of an outer diorite, a main mass of biotite gabbro, and an inner diorite. The entire complex and most of the surrounding country rocks have experienced various degrees of ¹⁸O depletion (up to 10 per mil) due to interaction with heated meteoric waters. The inner diorite apparently formed from a low-¹⁸O magma with δ¹⁸O ≃ +2.5, but most of the isotopic effects are a result of exchange between H₂O and solidified igneous rocks. The low-¹⁸O inner diorite magma was probably produced by massive assimilation and/or melting of hydrothermally altered country rocks. The δ¹⁸O values of the rocks generally increase with increasing grain size, except that quartz typically has δ¹⁸O = +6 to +8, and is more resistant to hydrothermal exchange than any other mineral studied. Based on atom % oxygen, the outer diorites, gabbros, and volcanic rocks exhibit integrated water/rock ratios of 0.3 ± 0.2, 0.15 ± 0.1, and 0.2 ± 0.1, respectively. Locally, water/rock ratios attain values greater than 1.0. Hydrogen isotopic analyses of sericites, chlorites, biotites, and amphiboles range from -117 to -150. δD in biotites varies inversely with Fe/Fe+Mg, as predicted by Suzuoki and Epstein (1974), and positively with elevation, over a range of 600 m. The calculated δD of the mid-to-late-Tertiary meteoric waters is about -100. Carbonate δ¹³C values average -5.5 (PDB), within the generally accepted range for deep-seated carbon.

Almost all the rocks within 4 km of the central Tertiary intrusive complex of Skye are depleted in ¹⁸O. Whole-rock δ¹⁸O values of basalts (-7. 1 to +8.4), Mesozoic shales (-0.6 to + 12.4), and Precambrian sandstones (-6.2 to + 10.8) systematically decrease inward towards the center of the complex. The Cuillin gabbro may have formed from a ¹⁸O-depleted magma (depleted by about 2 per mil); δ¹⁸O of plagioclase (-7.1 to + 2.5) and pyroxene (-0.5 to + 3.2) decrease outward toward the margins of the pluton. The Red Hills epigranite plutons have δ¹⁸O quartz (-2.7 to + 7.6) and feldspar (-6.7 to + 6.0) that suggest about 3/4 of the exchange took place at subsolidus temperatures; profound disequilibrium quartz-feldspar fractionations (up to 12) are characteristic. The early epigranites were intruded as low-¹⁸O melts (depletions of up to 3 per mil) with δ¹⁸O of the primary, igneous quartz decreasing progressively with time. The Southern Porphyritic Epigranite was apparently intruded as a low-¹⁸O magma with δ¹⁸O ≃ -2.6. A good correlation exists between grain size and δ¹⁸O for the unique, high-¹⁸O Beinn an Dubhaich granite which intrudes limestone having a δ¹⁸O range of +0.5 to +20.8, and δ¹³C of -4.9 to -1.0. The δD values of sericites (-104 to -107), and amphiboles, chlorites, and biotites (-105 to -128) from the igneous rocks , indicate that Eocene surface waters at Skye had δD ≃ -90. The average water/rock ratio for the Skye hydrothermal system is approximately one; at least 2000 km³ of heated meteoric waters were cycled through these rocks.

Thus these detailed isotopic studies of two widely separated areas indicate that (1) ¹⁸O-depleted magmas are commonly produced in volcanic terranes invaded by epizonal intrusions; (2) most of the ¹⁸O-depletion in such areas are a result of subsolidus exchange (particularly of feldspars); however correlation of δ¹⁸O with grain size is generally preserved only for systems that have undergone relatively minor meteoric hydrothermal exchange; (3) feldspar and calcite are the minerals mos t susceptible to oxygen isotopic exchange, whereas quartz is very resistant to oxygen isotope exchange; biotite, magnetite, and pyroxene have intermediate susceptibilities; and (4) basaltic country rocks are much more permeable to the hydrothermal convective system than shale, sandstone, or the crystalline basement complex.

Structural requirements for catalysis: studies on RTEM-1 β-lactamase

β-lactamases are a group of enzymes that confer resistance to penam and cephem antibiotics by hydrolysis of the β-lactam ring, thereby inactivating the antibiotic. Crystallographic and computer modeling studies of RTEM-1 β-lactamase have indicated that Asp 132, a strictly conserved residue among the class A β-lactamases, appears to be involved in substrate binding, catalysis, or both. To study the contribution of residue 132 to β-lactamase function, site saturation mutagenesis was used to generate mutants coding for all 20 amino acids at position 132. Phenotypic screening of all mutants indicated that position 132 is very sensitive to amino acid changes, with only N132C, N132D, N132E, and N132Q showing any appreciable activity. Kinetic analysis of three of these mutants showed increases in K_M, along with substantial decreases in k_(cat). Efforts to trap a stable acyl-enzyme intermediate were unsuccessfuL These results indicate that residue 132 is involved in substrate binding, as well as catalysis, and supports the involvement of this residue in acylation as suggested by Strynadka et al.

Crystallographic and computer modeling studies of RTEM-1 β-lactamase have indicated that Lys 73 and Glu 166, two strictly conserved residues among the class A β-lactamases, appear to be involved in substrate binding, catalysis, or both. To study the contribution of these residues to β-lactamase function, site saturation mutagenesis was used to generate mutants coding for all 20 amino acids at positions 73 and 166. Then all 400 possible combinations of mutants were created by combinatorial mutagenesis. The colonies harboring the mutants were screened for growth in the presence of ampicillin. The competent colonys' DNA were sequenced, and kinetic parameters investigated. It was found that lysine is essential at position 73, and that position 166 only tolerated fairly conservative changes (Aspartic acid, Histidine, and Tyrosine). These functional mutants exhibited decreased kcat's, but K_M was close to wild-type levels. The results of the combinatorial mutagenesis experiments indicate that Lysis absolutely required for activity at position 73; no mutation at residue 166 can compensate for loss of the long side chain amine. The active mutants found--K73K/E166D, K73KIE166H, and K73KIE166Y were studied by kinetic analysis. These results reaffirmed the function of residue 166 as important in catalysis, specifically deacylation.

The identity of the residue responsible for enhancing the active site serine (Ser 70) in RTEM-1 β-lactamase has been disputed for some time. Recently, analysis of a crystal structure of RTEM-1 β-lactamase with covalently bound intermediate was published, and it was suggested that Lys 73, a strictly conserved residue among the class A β-lactamases, was acting as a general base, activating Ser 70. For this to be possible, the pK_a of Lys 73 would have to be depressed significantly. In an attempt to assay the pK_a of Lys 73, the mutation K73C was made. This mutant protein can be reacted with 2-bromoethylamine, and activity is restored to near wild type levels. ^(15)N-2-bromoethylamine hydrobromide and ^(13)C-2-bromoethylamine hydrobromide were synthesized. Reacting these compounds with the K73C mutant gives stable isotopic enrichment at residue 73 in the form of aminoethylcysteine, a lysine homologue. The pK_a of an amine can be determined by NMR titration, following the change in chemical shift of either the ^(15)N-amine nuclei or adjacent Be nuclei as pH is changed. Unfortunately, low protein solubility, along with probable label scrambling in the Be experiment, did not permit direct observation of either the ^(15)N or ^(13)C signals. Indirect detection experiments were used to observe the protons bonded directly to the ^(13)C atoms. Two NMR signals were seen, and their chemical shift change with pH variation was noted. The peak which was determined to correspond to the aminoethylcysteine residue shifted from 3.2 ppm down to 2.8 ppm over a pH range of 6.6 to 12.5. The pK_a of the amine at position 73 was determined to be ~10. This indicates that residue 73 does not function as a general base in the acylation step of the reaction. However the experimental measurement takes place in the absence of substrate. Since the enzyme undergoes conformational changes upon substrate binding, the measured pK_a of the free enzyme may not correspond to the pK_a of the enzyme substrate complex.

Oxygen, carbon and hydrogen isotope studies of contact metamorphism

The O¹⁸/O¹⁶, C¹³/C¹², and D/H ratios have been determined for rocks and coexisting minerals from several granitic plutons and their contact metamorphic aureoles in northern Nevada, eastern California, central Colorado, and Texas, with emphasis on oxygen isotopes. A consistent order of O¹⁸/O¹⁶, C¹³/C¹², and D/H enrichment in coexisting minerals, and a correlation between isotopic fractionations among coexisting mineral pairs are in general observed, suggesting that mineral assemblages tend to approach isotopic equilibrium during contact metamorphism. In certain cases, a correlation is observed between oxygen isotopic fractionations of a mineral pair and sample distance from intrusive contacts. Isotopic temperatures generally show good agreement with heat flow considerations. Based on the experimentally determined quartz-muscovite O¹⁸/O¹⁶ fractionation calibration curve, temperatures are estimated to be 525 to 625°C at the contacts of the granitic stocks studied.

Small-scale oxygen isotope exchange effects between intrusive and country rock are observed over distances of 0.5 to 3 feet on both sides of the contacts; the isotopic gradients are typically 2 to 3 per mil per foot. The degree of oxygen isotopic exchange is essentially identical for different coexisting minerals. This presumably occurred through a diffusion-controlled recrystallization process. The size of the oxygen isotope equilibrium systems in the small-scale exchanged zones vary from about 1.5 cm to 30 cm. A xenolith and a re-entrant of country rock projecting into on intrusive hove both undergone much more extensive isotopic exchange (to hundreds of feet); they also show abnormally high isotopic temperatures. The marginal portions of most plutons have unusually high O¹⁸/O¹⁶ ratios compared to "normal" igneous rocks, presumably due to large-scale isotopic exchange with meta-sedimentary country rocks when the igneous rocks were essentially in a molten state. The isotopic data suggest that outward horizontal movement of H₂O into the contact metamorphic aureoles is almost negligible, but upward movement of H₂O may be important. Also, direct influx and absorption of water from the country rock may be significant in certain intrusive stocks.

Except in the exchanged zones, the O¹⁸/O¹⁶ ratios of pelitic rocks do not change appreciably during contact metamorphism, even in the cordierite and sillimanite grades; this is in contrast to regional metamorphic rocks which commonly decrease in O¹⁸ with increasing grade. Low O¹⁸/O¹⁶ and C¹³/C¹² ratios of the contact metamorphic marbles generally correlate well with the presence of calc-silicate minerals, indicating that the CO₂ liberated during metamorphic decarbonation reactions is enriched in both O¹⁸ and C¹³ relative to the carbonates.

The D/H ratios of biotites in the contact metamorphic rocks and their associated intrusions show a geographic correlation that is similar to that shown by the D/H ratios of meteoric surface waters, perhaps indicating that meteoric waters were present in the rocks during crystallization of the biotites.

Theory of ozone isotopic effects and various electron transfer reactions

Three separate topics, each stimulated by experiments, are treated theoretically in this dessertation: isotopic effects of ozone, electron transfer at interfaces, and intramolecular directional electron transfer in a supramolecular system.

The strange mass-independent isotope effect for the enrichment of ozone, which has been a puzzle in the literature for some 20 years, and the equally puzzling unconventional strong mass-dependent effect of individual reaction rate constants are studied as different aspects of a symmetry-driven behavior. A statistical (RRKM-based) theory with a hindered-rotor transition state is used. The individual rate constant ratios of recombination reactions at low pressures are calculated using the theory involving (1) small deviation from the statistical density of states for symmetric isotopomers, and (2) weak collisions for deactivation of the vibrationally excited ozone molecules. The weak collision and partitioning among exit channels play major roles in producing the large unconventional isotope effect in "unscrambled" systems. The enrichment studies reflect instead the non-statistical effect in "scrambled" systems. The theoretical results of low-pressure ozone enrichments and individual rate constant ratios obtained from these calculations are consistent with the corresponding experimental results. The isotopic exchange rate constant for the reaction ^(16)O + ^(18)O ^(18)O→+ ^(16)O ^(18)O + ^(18)O provides information on the nature of a variationally determined hindered-rotor transition state using experimental data at 130 K and 300 K. Pressure effects on the recombination rate constant, on the individual rate constant ratios and on the enrichments are also investigated. The theoretical results are consistent with the experimental data. The temperature dependence of the enrichment and rate constant ratios is also discussed, and experimental tests are suggested. The desirability of a more accurate potential energy surface for ozone in the transition state region is also noted.

Electron transfer reactions at semiconductor /liquid interfaces are studied using a tight-binding model for the semiconductors. The slab method and a z-transform method are employed in obtaining the tight-binding electronic structures of semiconductors having surfaces. The maximum electron transfer rate constants at Si/viologen^(2-/+) and InP /Me_(2)Fc^(+/O) interfaces are computed using the tight-binding type calculations for the solid and the extended-Huckel for the coupling to the redox agent at the interface. These electron transfer reactions are also studied using a free electron model for the semiconductor and the redox molecule, where Bardeen's method is adapted to calculate the coupling matrix element between the molecular and semiconductor electronic states. The calculated results for maximum rate constant of the electron transfer from the semiconductor bulk states are compared with the experimentally measured values of Lewis and coworkers, and are in reasonable agreement, without adjusting parameters. In the case of InP /liquid interface, the unusual current vs applied potential behavior is additionally interpreted, in part, by the presence of surface states.

Photoinduced electron transfer reactions in small supramolecular systems, such as 4-aminonaphthalimide compounds, are interesting in that there are, in principle, two alternative pathways (directions) for the electron transfer. The electron transfer, however, is unidirectional, as deduced from pH-dependent fluorescence quenching studies on different compounds. The role of electronic coupling matrix element and the charges in protonation are considered to explain the directionality of the electron transfer and other various results. A related mechanism is proposed to interpret the fluorescence behavior of similar molecules as fluorescent sensors of metal ions.

I. Electron spin relaxation studies of manganese(II) complexes in acetonitrile. II. Nuclear spin relaxation studies of bromine in tetrabromomanganese(II) in acetonitrile

Magnetic resonance techniques have given us a powerful means for investigating dynamical processes in gases, liquids and solids. Dynamical effects manifest themselves in both resonance line shifts and linewidths, and, accordingly, require detailed analyses to extract desired information. The success of a magnetic resonance experiment depends critically on relaxation mechanisms to maintain thermal equilibrium between spin states. Consequently, there must be an interaction between the excited spin states and their immediate molecular environment which promote changes in spin orientation while excess magnetic energy is coupled into other degrees of freedom by non-radiative processes. This is well known as spin-lattice relaxation. Certain dynamical processes cause fluctuations in the spin state energy levels leading to spin-spin relaxation and, here again, the environment at the molecular level plays a significant role in the magnitude of interaction. Relatively few electron spin relaxation studies of solutions have been conducted and the present work is addressed toward the extension of our knowledge in this area and the retrieval of dynamical information from line shape analyses on a time scale comparable to diffusion controlled phenomena.

Specifically, the electron spin relaxation of three Mn⁺²3d⁵ complexes, Mn(CH₃CN)₆⁺², MnCl₄^-2 in acetonitrile has been studied in considerable detail. The effective spin Hamiltonian constants were carefully evaluated under a wide range of experimental conditions. Resonance widths of these Mn⁺² complexes were studied in the presence of various excess ligand ions and as a function of concentration, viscosity, temperature and frequency (X-band, ~9.5 Ԍ Hz and K-band, ~35 Ԍ Hz).

A number of interesting conclusions were drawn from these studies. For the Et₄NCl-₄^-2 system several relaxation mechanisms leading to resonance broadening were observed. One source appears to arise through spin-orbit interactions caused by modulation of the ligand field resulting from transient distortions of the complex imparted by solvent fluctuations in the immediate surroundings of the paramagnetic ion. An additional spin relaxation was assigned to the formation of ion pairs [Et₄N⁺…MnCl₄^-2] and it was possible to estimate the dissociation constant for this specie in acetonitrile.

The Bu₄NBr-MnBr₄^-2 study was considerably more interesting. As in the former case, solvent fluctuations and ion-pairing of the paramagnetic complex [Bu₄N⁺…MnBr₄^-2] provide significant relaxation for the electronic spin system. Most interesting, without doubt, is the onset of a new relaxation mechanism leading to resonance broadening which is best interpreted as chemical exchange. Thus, assuming that resonance widths were simply governed by electron spin state lifetimes, we were able to extract dynamical information from an interaction in which the initial and final states are the same

MnBr₄^-2 + Br^- = MnBr₄^-2 + Br^-.

The bimolecular rate constants were obtained at six different temperatures and their magnitudes suggested that the exchange is probably diffusion controlled with essentially a zero energy of activation. The most important source of spin relaxation in this system stems directly from dipolar interactions between the manganese 3d⁵ electrons. Moreover, the dipolar broadening is strongly frequency dependent indicating a deviation between the transverse and longitudinal relaxation times. We are led to the conclusion that the 3d⁵ spin states of ion-paired MnBr₄^-2 are significantly correlated so that dynamical processes are also entering the picture. It was possible to estimate the correlation time, ^Td, characterizing this dynamical process.

In Part II we study nuclear magnetic relaxation of bromine ions in the MnBr4-2-Bu4NBr-acetonitrile system. Essentially we monitor the 79Br and 81Br linewidths in response to the [MnBr4-2]/[Br-] ratio with the express purpose of supporting our contention that exchange is occurring between "free" bromine ions in the solvent and bromine in the first coordination sphere of the paramagnetic anion. The complexity of the system elicited a two-part study: (1) the linewidth behavior of Bu4NBr in anhydrous CH3CN in the absence of MnBr4-2 and (2) in the presence of MnBr4-2. It was concluded in study (1) that dynamical association, Bu4NBr k1= Bu4N+ + Br-, was modulating field-gradient interactions at frequencies high enough to provide an estimation of the unimolecular rate constant, k1. A comparison of the two isotopic bromine linewidth-mole fraction results led to the conclusion that quadrupole interactions provided the dominant relaxation mechanism. In study (2) the "residual" bromine linewidths for both 79Br and 81Br are clearly controlled by quadrupole interactions which appear to be modulated by very rapid dynamical processes other than molecular reorientation. We conclude that the "residual" linewidth has its origin in chemical exchange and that bromine nuclei exchange rapidly between a "free" solvated ion and the paramagnetic complex, MnBr4-2.

Studies of noble gases in meteorites and in the earth

The isotopic and elemental abundances of noble gases in the solar system are investigated, using simple mixing models and mass-spectrometric measurements of the noble gases in meteorites and terrestrial rocks and minerals.

Primordial neon is modeled by two isotopically distinct components from the interstellar gas and dust. Neon from the gas dominates solar neon, which contains about ten times more ²⁰Ne than ²²Ne. Neon from the dust is represented in meteorites by neon-E, with ²⁰Ne/²²Ne less than 0.6. Isotopic variations in meteorites require neon from both dust and gas to be present. Mixing dust and gas without neon loss generates linear correlation lines on three-isotope and composition-concentration diagrams. A model for solar wind implantation predicts small deviations from linear mixing, due to preferential sputtering of the lighter neon isotopes.

Neon in meteorites consists of galactic cosmic ray spallation neon and at least two primordial components, neon-E and neon-S. Neon was measured in several meteorites to investigate these end- members. Cosmogenic neon produced from sodium is found to be strongly enriched in ²²Ne. Neon measurements on sodium-rich samples must be interpreted with care so not to confuse this source of ²²Ne with neon-E, which is also rich in ²²Ne.

Neon data for the carbonaceous chondrite Mokoia show that the end member composition of neon-Si in meteorites is ²⁰Ne/²²Ne = 13.7, the same as the present solar wind. The solar wind composition evidently has remained constant since before the compaction of Mokoia.

Ca, Al-rich inclusions from the Allende meteorite were examined for correlation between neon-E and oxygen or magnesium isotopic anomalies. ²²Ne and ³⁶Ar enrichments found in some inclusions are attributed to cosmic- ray-induced reactions on Na and Cl, not to a primordial component. Neon-E is not detectably enriched in Allende.

Measurements were made to determine the noble gas contents of various terrestrial rocks and minerals, and to investigate the cycling of noble gases between different terrestrial reservoirs. Beryl crystals contain a characteristic suite of magmatic gases including nucleogenic ²¹Ne and ²²Ne from (α,n) reactions, radiogenic ⁴⁰Ar, and fissiogenic ^131-136Xe from the decay of K and U in the continental crust. Significant concentrations of atmospheric noble gases are also present in beryl.

Both juvenile and atmospheric noble gases are found in rocks from the Skaergaard intrusion. The ratio ⁴⁰Ar/³⁶Ar (corrected for in situ decay of ⁴⁰K) correlates with δ¹⁸O in plagioclase. Atmospheric argon has been introduced into samples that have experienced oxygen-isotope exchange with circulating meteoric hydrothermal fluids. Unexchanged samples contain juvenile argon with ⁴⁰Ar/³⁶Ar greater than 6000 that was trapped from the Skaergaard magma.

Juvenile and atmospheric gases have been measured in the glassy rims of mid-ocean ridge (MOR) pillow basalts. Evidence is presented that three samples contain excess radiogenic ¹²⁹Xe and fission xenon, in addition to the excess radiogenic ⁴⁰Ar found in all samples. These juvenile gases are being outgassed from the upper-mantle source region of the MOR magma. No isotopic evidence has been found here for juvenile primordial noble gases accompanying the juvenile radiogenic gases in the MOR glasses. Large argon isotopic variations in a single specimen provide a clear indication of the late-stage addition of atmospheric argon, probably from seawater.

The Skaergaard data demonstrate that atmospheric noble gases dissolved in ground water can be transferred into crustal rocks. Subduction of oceanic crust altered by seawater can transport atmospheric noble gases into the upper mantle. A substantial portion of the noble gases in mantle derived rocks may represent subducted gases, not a primordial component as is often assumed.