An investigation of techniques for the measurement and interpretation of cosmic ray isotopic abundances

An instrument, the Caltech High Energy Isotope Spectrometer Telescope (HEIST), has been developed to measure isotopic abundances of cosmic ray nuclei in the charge range 3 ≤ Z ≤ 28 and the energy range between 30 and 800 MeV/nuc by employing an energy loss -- residual energy technique. Measurements of particle trajectories and energy losses are made using a multiwire proportional counter hodoscope and a stack of CsI(TI) crystal scintillators, respectively. A detailed analysis has been made of the mass resolution capabilities of this instrument.

Landau fluctuations set a fundamental limit on the attainable mass resolution, which for this instrument ranges between ~.07 AMU for z~3 and ~.2 AMU for z~2b. Contributions to the mass resolution due to uncertainties in measuring the path-length and energy losses of the detected particles are shown to degrade the overall mass resolution to between ~.1 AMU (z~3) and ~.3 AMU (z~2b).

A formalism, based on the leaky box model of cosmic ray propagation, is developed for obtaining isotopic abundance ratios at the cosmic ray sources from abundances measured in local interstellar space for elements having three or more stable isotopes, one of which is believed to be absent at the cosmic ray sources. This purely secondary isotope is used as a tracer of secondary production during propagation. This technique is illustrated for the isotopes of the elements O, Ne, S, Ar and Ca.

The uncertainties in the derived source ratios due to errors in fragmentation and total inelastic cross sections, in observed spectral shapes, and in measured abundances are evaluated. It is shown that the dominant sources of uncertainty are uncorrelated errors in the fragmentation cross sections and statistical uncertainties in measuring local interstellar abundances.

These results are applied to estimate the extent to which uncertainties must be reduced in order to distinguish between cosmic ray production in a solar-like environment and in various environments with greater neutron enrichments.

Observation and interpretation of tectonic strain release mechanisms

In four chapters various aspects of earthquake source are studied.

Chapter I

Surface displacements that followed the Parkfield, 1966, earthquakes were measured for two years with six small-scale geodetic networks straddling the fault trace. The logarithmic rate and the periodic nature of the creep displacement recorded on a strain meter made it possible to predict creep episodes on the San Andreas fault. Some individual earthquakes were related directly to surface displacement, while in general, slow creep and aftershock activity were found to occur independently. The Parkfield earthquake is interpreted as a buried dislocation.

Chapter II

The source parameters of earthquakes between magnitude 1 and 6 were studied using field observations, fault plane solutions, and surface wave and S-wave spectral analysis. The seismic moment, M_O, was found to be related to local magnitude, M_L, by log M_O = 1.7 M_L + 15.1. The source length vs magnitude relation for the San Andreas system found to be: M_L = 1.9 log L - 6.7. The surface wave envelope parameter AR gives the moment according to log M_O = log AR₃₀₀ + 30.1, and the stress drop, τ, was found to be related to the magnitude by τ = 0.54 M - 2.58. The relation between surface wave magnitude M_S and M_L is proposed to be M_S = 1.7 M_L - 4.1. It is proposed to estimate the relative stress level (and possibly the strength) of a source-region by the amplitude ratio of high-frequency to low-frequency waves. An apparent stress map for Southern California is presented.

Chapter III

Seismic triggering and seismic shaking are proposed as two closely related mechanisms of strain release which explain observations of the character of the P wave generated by the Alaskan earthquake of 1964, and distant fault slippage observed after the Borrego Mountain, California earthquake of 1968. The Alaska, 1964, earthquake is shown to be adequately described as a series of individual rupture events. The first of these events had a body wave magnitude of 6.6 and is considered to have initiated or triggered the whole sequence. The propagation velocity of the disturbance is estimated to be 3.5 km/sec. On the basis of circumstantial evidence it is proposed that the Borrego Mountain, 1968, earthquake caused release of tectonic strain along three active faults at distances of 45 to 75 km from the epicenter. It is suggested that this mechanism of strain release is best described as "seismic shaking."

Chapter IV

The changes of apparent stress with depth are studied in the South American deep seismic zone. For shallow earthquakes the apparent stress is 20 bars on the average, the same as for earthquakes in the Aleutians and on Oceanic Ridges. At depths between 50 and 150 km the apparent stresses are relatively high, approximately 380 bars, and around 600 km depth they are again near 20 bars. The seismic efficiency is estimated to be 0.1. This suggests that the true stress is obtained by multiplying the apparent stress by ten. The variation of apparent stress with depth is explained in terms of the hypothesis of ocean floor consumption.

Mechanistic and reactivity studies of bent metallocene complexes of niobium and tantalum

A variety of olefin hydride complexes of niobium and tantalum has been prepared in order to study their reactivity and to gain insight into organometallic reaction mechanisms. Examination of a series of ethylene and propylene complexes of niobocene (CP_2Nb; Cp = η^5-C_5H_5), permethylniobocene (Cp*_2Nb; Cp* = η^5-C_5(CH_3)_5), tantalocene, and permethyltantalocene has indicated that there are both large electronic and steric effects deriving from the metal (and its ancillary ligands) in the olefin insertion (β-migratory insertion) process. Furthermore, a thermodynamic and kinetic analysis has been completed for a series of substituted styrene complexes of niobocene in order to better understand the important electronic properties of the olefin. The results are in accord with a concerted four-center process with only moderate charge development.

The special case of β-migratory insertion of a hydride ligand into coordinated benzyne has also been studied for the permethyltantalocene system. The coordinatively unsaturated (sixteen electron) phenyl tautomer, which is made accessible by the facile benzyne hydride insertion reaction, readily reacts with a variety of ligands, L, to afford Cp*_2 Ta(C_6H_5)L complexes (L = CO, O_2, NC≡R, :CH_2, H_2, etc.). This family of compounds exhibits interesting reactivity (a-migratory insertion, O_2 activation, and reductive elimination) which is discussed in some detail.

Finally a series of paramagnetic seventeen electron Cp*_2 TaX_2 (X = halide, alkyl, hydride) complexes, and the corresponding cationic and anionic species, have been prepared and studied. The odd electron neutral complexes exhibit surprising thermal stability and undergo very little reactivity. While the chemistry of the anionic compounds is almost completely dominated by their potent reducing power, that of the cations is quite diverse and amenable for study. Therefore the syntheses and reactivity (1 ,2-eliminations, ligand insertions, and deprotonation reactions) of these coordinatively unsaturated sixteen electron species are presented.

Isolation and characterization of proteins that bind to yeast origins of DNA replication

Yeast chromosomes contain sequences called ARSs which function as origins of replication in vitro and in vivo. We have carried out a systematic deletion analysis of ARS1, allowing us to define three functionally distinct domains, designated A, B, and C. Domain A is a sequence of 11 to 19bp, containing the core consensus element that is required for replication. The core consensus sequence, A/TTTTATPuTTTA/T, is conserved at all ARSs sequenced to date. A fragment containing only element A and 8 flanking nucleotides enables autonomous replication of centromeric plasmids. These plasmids replicate very inefficiently, suggesting that flanking sequences must be important for ARS function. Domain B also provides important sequences needed for efficient replication. Deletion of domain B drastically increases the doubling times of transformants and reduces plasmid stability. Domain B contains a potential consensus sequence conserved at some ARSs which overlaps a region of bent DNA. Mutational analysis suggests this bent DNA may be important for ARS function. Deletion of domain C has only a slight effect on replication of plasmids carrying those deletions.

We have identified a protein called ARS binding factor I (ABF-I) that binds to the HMR-E ARS and ARS1. We have purified this protein to homogeneity using conventional and oligonucleotide affinity chromatography. The protein has an apparent molecular weight of 135kDa and is present at about 700 molecules per diploid cell, based on the yield of purified protein and in situ antibody staining. DNaseI footprinting reveals that ABF-I binds sequence-specifically to an approximately 24bp sequence that overlaps element Bat ARSl. This same protein binds to and protects a similar size region at the HMR-E ARS.

We also find evidence for another ARS binding protein, ABF-III, based on DN asei footprint analysis and gel retardation assays. The protein protects approximately 22bp adjacent to the ABF-I site. There appears to be no interaction between ABF-I and ABF-III despite the proximity of their binding sites.

To address the function of ABF-I in DNA replication, we have cloned the ABF-I gene using rabbit polyclonal anti-sera and murine monoclonal antibodies against ABF-I to screen a λgt11 expression library. Four EcoRI restriction fragments were isolated which encoded proteins that were recognized by both polyclonal and monoclonal antibodies. A gene disruption can now be constructed to determine the in vivo function of ABF-I.

Synthetic and structural studies of permethylscandocene and permethyltantalocene derivatives

Cp*_2Sc-H reacts with H_2 and CO at -78°C to yield Cp*_2ScOCH_3. A stepwise reduction of CO to an alkoxide is observed when CO reacts with Cp*_2ScC_6H_4CH_3-p to give the η^2-acyl Cp*_2Sc(CO)C_6H_4CH_3-p, which then reacts with H_2 to produce Cp*_2ScOCH_2C_6H_4CH_3-p. Cp*_2ScCH_3 and Cp*_2ScH(THF) react with CO to give unchar- uncharacterizable products. Cp*_2ScH and Cp*_2ScCH_3 react with Cp_2MCO (M = Mo, W) to give scandoxycarbenes, Cp_2M=C(CH_3)OScCp*_2, while a wide variety of Cp*_2ScX (X = H, CH_3, N(CH_3)_2, CH_2CH_2C_6H_5) reacts with CpM(CO)_2 (M = Co, Rh) to yield similar carbene complexes. An x-ray crystal structure determination of Cp(CO)Co=C(CH_3)- OScCp*_2 revealed a µ^2: η^1, η^1 carbonyl interaction between the Co-CO and Sc.

CO_2 inserts cleanly into Sc-phenyl bonds at -78°C to produce a carboxylate complex, Cp*_2Sc(O_2C)C_6H_4CH_3-p. The structure of this compound was determined by x-ray crystallographic techniques.

Excess C_2H_2 reacts with Cp*_2ScR (R = H, alkyl, aryl, alkenyl, alkynyl, amide) at temperatures below -78°C to form the alkynyl species Cp*_2Sc-C≡C-H, which then reacts with the remaining acetylene to form polyacetylene. Cp*_2Sc-C≡C-H further reacts to yield Cp*_2sc-C≡C-ScCp*_2. This unusual C_2 bridged dimer was characterized by x-ray crystallography.

Attempts were made to model the C-N bond breaking step of hydrodenitrogenation by synthesizing Cp*_2TaH(η^2-H_2C=N(C_6H_4X)) and studying its rearrangement to Cp*_2Ta(=N(C_6H_4X))(CH_3). The 1,2 addition/elimination reactions of Cp*_2Ta(η^2- H_2C=N(CH_3)H and Cp*_2Ta(=X)H (X=O, S, NH, N(C_6H_5)) were investigated. Cp*_2Ta(=NH)H was found to react with D_2 to give Cp*_2Ta(=ND)H, implying a nonsymmetric amide-dihydride intermediate for the addition/elimination process. Cp*_2Ta(=S)H and H_2O equilibrate with Cp*_2Ta(=O)H and H_2S, which allowed determination of the difference in bond strengths for Ta=O and Ta=S. Ta=O was found to be approximately 41 kcals/mole stronger than Ta=S.

Characterization and developmental regulation of a gene expressed specifically in the skeletogenic lineage of the sea urchin embryo

The sea urchin embryonic skeleton, or spicule, is deposited by mesenchymal progeny of four precursor cells, the micromeres, which are determined to the skeletogenic pathway by a process known as cytoplasmic localization. A gene encoding one of the major products of the skeletogenic mesenchyme, a prominent 50 kD protein of the spicule matrix, has been characterized in detail. cDNA clones were first isolated by antibody screening of a phage expression library, followed by isolation of homologous genomic clones. The gene, known as SM50, is single copy in the sea urchin genome, is divided into two exons of 213 and 1682 bp, and is expressed only in skeletogenic cells. Transcripts are first detectable at the 120 cell stage, shortly after the segregation of the skeletogenic precursors from the rest of the embryo. The SM50 open reading frame begins within the first exon, is 450 amino acids in length, and contains a loosely repeated 13 amino acid motif rich in acidic residues which accounts for 45% of the protein and which is possibly involved in interaction with the mineral phase of the spicule.

The important cis-acting regions of the SM50 gene necessary for proper regulation of expression were identified by gene transfer experiments. A 562 bp promoter fragment, containing 438 bp of 5' promoter sequence and 124 bp of the SM50 first exon (including the SM50 initiation codon), was both necessary and sufficient to direct high levels of expression of the bacterial chloramphenicol acetyltransferase (CAT) reporter gene specifically in the skeletogenic cells. Removal of promoter sequences between positions -2200 and -438, and of transcribed regions downstream of +124 (including the SM50 intron), had no effect on the spatial or transcriptional activity of the transgenes.

Regulatory proteins that interact with the SM50 promoter were identified by the gel retardation assay, using bulk embryo mesenchyme blastula stage nuclear proteins. Five protein binding sites were identified and mapped to various degrees of resolution. Two sites are homologous, may be enhancer elements, and at least one is required for expression. Two additional sites are also present in the promoter of the aboral ectoderm specific cytoskeletal actin gene CyIIIa; one of these is a CCAA T element, the other a putative repressor element. The fifth site overlaps the binding site of the putative repressor and may function as a positive regulator by interfering with binding of the repressor. All of the proteins are detectable in nuclear extracts prepared from 64 cell stage embryos, a stage just before expression of SM50 is initiated, as well as from blastula and gastrula stage; the putative enhancer binding protein may be maternal as well.

Electrochemical and electron transfer investigations of copper proteins

A study of the pH and temperature dependence of the redox potentials of azurins from five species of bacteria has been performed. The variations in the potentials with pH have been interpreted in terms of electrostatic interactions between the copper site and titrating histidine residues, including the effects of substitutions in the amino acid sequences of the proteins on the electrostatic interactions. A comparison of the observed pH dependences with predictions based on histidine pK_a values known for Pseudomonas aeruginosa (Pae), Alcaligenes denitrificans (Ade), and Alcaligenes faecalis (Afa) azurins indicates that the Pae and Ade redox potentials exhibit pH dependences in line with electrostatic arguments, while Afa azurin exhibits more complex behavior. Redox enthalpies and entropies for four of the azurins at low and high pH values have also been obtained. Based on these results in conjuction with the variable pH experiments, it appears that Bordetella bronchiseptica azurin may undergo a more substantial conformational change with pH than has been observed for other species of azurin.

The temperature dependence of the redox potential of bovine erythrocyte superoxide dismutase (SOD) has been determined at pH 7.0, with potassium ferricyanide as the mediator. The following thermodynamic parameters have been obtained (T = 25°C): E°' = 403±5 mV vs. NHE, ΔG°' = -9.31 kcal/mol, ΔH°' = -21.4 kcal/mol, ΔS°' = -40.7 eu, ΔS°'_(rc) = -25.1 eu. It is apparent from these results that ΔH°', rather than ΔS°', is the dominant factor in establishing the high redox potential of SOD. The large negative enthalpy of reduction may also reflect the factors which give SOD its high specificity toward reduction and oxidation by superoxide.

Preparation and characterization of an intramolecular electron transfer in a pentaammineruthenium derivative of Candida krisei cytochrome c

A semisynthetic binuclear metalloprotein has been prepared by appending the pentaammineruthenium moiety to histidine 39 of the cytochrome c from the yeast Candida krusei. The site of ruthenium binding was identified by peptide mapping. Spectroscopic and electrochemical properties of the derivative indicate the protein conformation is unperturbed by the modification. A preliminary (minimum) rate constant of 170s^(-1) has been determined for the intramolecular electron transfer from ruthenium(II) to iron(III), which occurs over a distance of at least 13Å (barring major conformational changes). Electrochemical studies indicate that this reaction should proceed with a driving force of ~170mV. The rate constant is an order of magnitude faster than that observed in horse heart cytochrome c for intramolecular electron transfer from pentaammineruthenium(II)(histidine 33) to iron(III) (over a similar distance, and with a similar driving force), suggesting a medium or orientation effect makes the Candida intramolecular electron transfer more favorable.

Structural and mechanistic motifs in membrane proteins: the three-dimensional modelling of rhodopsin, band 3, and the nicotinic acetylcholine receptor

Because so little is known about the structure of membrane proteins, an attempt has been made in this work to develop techniques by which to model them in three dimensions. The procedures devised rely heavily upon the availability of several sequences of a given protein. The modelling procedure is composed of two parts. The first identifies transmembrane regions within the protein sequence on the basis of hydrophobicity, β-turn potential, and the presence of certain amino acid types, specifically, proline and basic residues. The second part of the procedure arranges these transmembrane helices within the bilayer based upon the evolutionary conservation of their residues. Conserved residues are oriented toward other helices and variable residues are positioned to face the surrounding lipids. Available structural information concerning the protein's helical arrangement, including the lengths of interhelical loops, is also taken into account. Rhodopsin, band 3, and the nicotinic acetylcholine receptor have all been modelled using this methodology, and mechanisms of action could be proposed based upon the resulting structures.

Specific residues in the rhodopsin and iodopsin sequences were identified, which may regulate the proteins' wavelength selectivities. A hinge-like motion of helices M3, M4, and M5 with respect to the rest of the protein was proposed to result in the activation of transducin, the G-protein associated with rhodopsin. A similar mechanism is also proposed for signal transduction by the muscarinic acetylcholine and β-adrenergic receptors.

The nicotinic acetylcholine receptor was modelled with four trans-membrane helices per subunit and with the five homologous M2 helices forming the cation channel. Putative channel-lining residues were identified and a mechanism of channel-opening based upon the concerted, tangential rotation of the M2 helices was proposed.

Band 3, the anion exchange protein found in the erythrocyte membrane, was modelled with 14 transmembrane helices. In general the pathway of anion transport can be viewed as a channel composed of six helices that contains a single hydrophobic restriction. This hydrophobic region will not allow the passage of charged species, unless they are part of an ion-pair. An arginine residue located near this restriction is proposed to be responsible for anion transport. When ion-paired with a transportable anion it rotates across the barrier and releases the anion on the other side of the membrane. A similar process returns it to its original position. This proposed mechanism, based on the three-dimensional model, can account for the passive, electroneutral, anion exchange observed for band 3. Dianions can be transported through a similar mechanism with the additional participation of a histidine residue. Both residues are located on M10.

Dating and characterizing late holocene earthquakes using paleomagnetics

In this thesis I apply paleomagnetic techniques to paleoseismological problems. I investigate the use of secular-variation magnetostratigraphy to date prehistoric earthquakes; I identify liquefaction remanent magnetization (LRM), and I quantify coseismic deformation within a fault zone by measuring the rotation of paleomagnetic vectors.

In Chapter 2 I construct a secular-variation reference curve for southern California. For this curve I measure three new well-constrained paleomagnetic directions: two from the Pallett Creek paleoseismological site at A.D. 1397-1480 and A.D. 1465-1495, and one from Panum Crater at A.D. 1325-1365. To these three directions I add the best nine data points from the Sternberg secular-variation curve, five data points from Champion, and one point from the A.D. 1480 eruption of Mt. St. Helens. I derive the error due to the non-dipole field that is added to these data by the geographical correction to southern California. Combining these yields a secular variation curve for southern California covering the period A.D. 670 to 1910, with the best coverage in the range A.D. 1064 to 1505.

In Chapter 3 I apply this curve to a problem in southern California. Two paleoseismological sites in the Salton trough of southern California have sediments deposited by prehistoric Lake Cahuilla. At the Salt Creek site I sampled sediments from three different lakes, and at the Indio site I sampled sediments from four different lakes. Based upon the coinciding paleomagnetic directions I correlate the oldest lake sampled at Salt Creek with the oldest lake sampled at Indio. Furthermore, the penultimate lake at Indio does not appear to be present at Salt Creek. Using the secular variation curve I can assign the lakes at Salt Creek to broad age ranges of A.D. 800 to 1100, A.D. 1100 to 1300, and A.D. 1300 to 1500. This example demonstrates the large uncertainties in the secular variation curve and the need to construct curves from a limited geographical area.

Chapter 4 demonstrates that seismically induced liquefaction can cause resetting of detrital remanent magnetization and acquisition of a liquefaction remanent magnetization (LRM). I sampled three different liquefaction features, a sandbody formed in the Elsinore fault zone, diapirs from sediments of Mono Lake, and a sandblow in these same sediments. In every case the liquefaction features showed stable magnetization despite substantial physical disruption. In addition, in the case of the sandblow and the sandbody, the intensity of the natural remanent magnetization increased by up to an order of magnitude.

In Chapter 5 I apply paleomagnetics to measuring the tectonic rotations in a 52 meter long transect across the San Andreas fault zone at the Pallett Creek paleoseismological site. This site has presented a significant problem because the brittle long-term average slip-rate across the fault is significantly less than the slip-rate from other nearby sites. I find sections adjacent to the fault with tectonic rotations of up to 30°. If interpreted as block rotations, the non-brittle offset was 14.0+2.8, -2.1 meters in the last three earthquakes and 8.5+1.0, -0.9 meters in the last two. Combined with the brittle offset in these events, the last three events all had about 6 meters of total fault offset, even though the intervals between them were markedly different.

In Appendix 1 I present a detailed description of my standard sampling and demagnetization procedure.

In Appendix 2 I present a detailed discussion of the study at Panum Crater that yielded the well-constrained paleomagnetic direction for use in developing secular variation curve in Chapter 2. In addition, from sampling two distinctly different clast types in a block-and-ash flow deposit from Panum Crater, I find that this flow had a complex emplacement and cooling history. Angular, glassy "lithic" blocks were emplaced at temperatures above 600° C. Some of these had cooled nearly completely, whereas others had cooled only to 450° C, when settling in the flow rotated the blocks slightly. The partially cooled blocks then finished cooling without further settling. Highly vesicular, breadcrusted pumiceous clasts had not yet cooled to 600° C at the time of these rotations, because they show a stable, well clustered, unidirectional magnetic vector.

Chemical and mineralogical characterization of micro-inclusions in diamonds

Secondary-ion mass spectrometry (SIMS), electron probe analysis (EPMA), analytical scanning electron microscopy (SEM) and infrared (IR) spectroscopy were used to determine the chemical composition and the mineralogy of sub-micrometer inclusions in cubic diamonds and in overgrowths (coats) on octahedral diamonds from Zaire, Botswana, and some unknown localities.

The inclusions are sub-micrometer in size. The typical diameter encountered during transmission electron microscope (TEM) examination was 0.1-0.5 µm. The micro-inclusions are sub-rounded and their shape is crystallographically controlled by the diamond. Normally they are not associated with cracks or dislocations and appear to be well isolated within the diamond matrix. The number density of inclusions is highly variable on any scale and may reach 10^(11) inclusions/cm^3 in the most densely populated zones. The total concentration of metal oxides in the diamonds varies between 20 and 1270 ppm (by weight).

SIMS analysis yields the average composition of about 100 inclusions contained in the sputtered volume. Comparison of analyses of different volumes of an individual diamond show roughly uniform composition (typically ±10% relative). The variation among the average compositions of different diamonds is somewhat greater (typically ±30%). Nevertheless, all diamonds exhibit similar characteristics, being rich in water, carbonate, SiO_2, and K_2O, and depleted in MgO. The composition of micro-inclusions in most diamonds vary within the following ranges: SiO_2, 30-53%; K_2O, 12-30%; CaO, 8-19%; FeO, 6-11%; Al_2O_3, 3-6%; MgO, 2-6%; TiO_2, 2-4%; Na_2O, 1-5%; P_2O_5, 1-4%; and Cl, 1-3%. In addition, BaO, 1-4%; SrO, 0.7-1.5%; La_2O_3, 0.1-0.3%; Ce_2O_3, 0.3-0.5%; smaller amounts of other rare-earth elements (REE), as well as Mn, Th, and U were also detected by instrumental neutron activation analysis (INAA). Mg/(Fe+Mg), 0.40-0.62 is low compared with other mantle derived phases; K/ AI ratios of 2-7 are very high, and the chondrite-normalized Ce/Eu ratios of 10-21 are also high, indicating extremely fractionated REE patterns.

SEM analyses indicate that individual inclusions within a single diamond are roughly of similar composition. The average composition of individual inclusions as measured with the SEM is similar to that measured by SIMS. Compositional variations revealed by the SEM are larger than those detected by SIMS and indicate a small variability in the composition of individual inclusions. No compositions of individual inclusions were determined that might correspond to mono-mineralic inclusions.

IR spectra of inclusion- bearing zones exhibit characteristic absorption due to: (1) pure diamonds, (2) nitrogen and hydrogen in the diamond matrix; and (3) mineral phases in the micro-inclusions. Nitrogen concentrations of 500-1100 ppm, typical of the micro-inclusion-bearing zones, are higher than the average nitrogen content of diamonds. Only type IaA centers were detected by IR. A yellow coloration may indicate small concentration of type IB centers.

The absorption due to the micro-inclusions in all diamonds produces similar spectra and indicates the presence of hydrated sheet silicates (most likely, Fe-rich clay minerals), carbonates (most likely calcite), and apatite. Small quantities of molecular CO_2 are also present in most diamonds. Water is probably associated with the silicates but the possibility of its presence as a fluid phase cannot be excluded. Characteristic lines of olivine, pyroxene and garnet were not detected and these phases cannot be significant components of the inclusions. Preliminary quantification of the IR data suggests that water and carbonate account for, on average, 20-40 wt% of the micro-inclusions.

The composition and mineralogy of the micro-inclusions are completely different from those of the more common, larger inclusions of the peridotitic or eclogitic assemblages. Their bulk composition resembles that of potassic magmas, such as kimberlites and lamproites, but is enriched in H_2O, CO_3, K_2O, and incompatible elements, and depleted in MgO.

It is suggested that the composition of the micro-inclusions represents a volatile-rich fluid or a melt trapped by the diamond during its growth. The high content of K, Na, P, and incompatible elements suggests that the trapped material found in the micro-inclusions may represent an effective metasomatizing agent. It may also be possible that fluids of similar composition are responsible for the extreme enrichment of incompatible elements documented in garnet and pyroxene inclusions in diamonds.

The origin of the fluid trapped in the micro-inclusions is still uncertain. It may have been formed by incipient melting of a highly metasomatized mantle rocks. More likely, it is the result of fractional crystallization of a potassic parental magma at depth. In either case, the micro-inclusions document the presence of highly potassic fluids or melts at depths corresponding to the diamond stability field in the upper mantle. The phases presently identified in the inclusions are believed to be the result of closed system reactions at lower pressures.

Singularities and horizons in the collisions of gravitational waves

This thesis presents a study of the dynamical, nonlinear interaction of colliding gravitational waves, as described by classical general relativity. It is focused mainly on two fundamental questions: First, what is the general structure of the singularities and Killing-Cauchy horizons produced in the collisions of exactly plane-symmetric gravitational waves? Second, under what conditions will the collisions of almost-plane gravitational waves (waves with large but finite transverse sizes) produce singularities?

In the work on the collisions of exactly-plane waves, it is shown that Killing horizons in any plane-symmetric spacetime are unstable against small plane-symmetric perturbations. It is thus concluded that the Killing-Cauchy horizons produced by the collisions of some exactly plane gravitational waves are nongeneric, and that generic initial data for the colliding plane waves always produce "pure" spacetime singularities without such horizons. This conclusion is later proved rigorously (using the full nonlinear theory rather than perturbation theory), in connection with an analysis of the asymptotic singularity structure of a general colliding plane-wave spacetime. This analysis also proves that asymptotically the singularities created by colliding plane waves are of inhomogeneous-Kasner type; the asymptotic Kasner axes and exponents of these singularities in general depend on the spatial coordinate that runs tangentially to the singularity in the non-plane-symmetric direction.

In the work on collisions of almost-plane gravitational waves, first some general properties of single almost-plane gravitational-wave spacetimes are explored. It is shown that, by contrast with an exact plane wave, an almost-plane gravitational wave cannot have a propagation direction that is Killing; i.e., it must diffract and disperse as it propagates. It is also shown that an almost-plane wave cannot be precisely sandwiched between two null wavefronts; i.e., it must leave behind tails in the spacetime region through which it passes. Next, the occurrence of spacetime singularities in the collisions of almost-plane waves is investigated. It is proved that if two colliding, almost-plane gravitational waves are initially exactly plane-symmetric across a central region of sufficiently large but finite transverse dimensions, then their collision produces a spacetime singularity with the same local structure as in the exact-plane-wave collision. Finally, it is shown that a singularity still forms when the central regions are only approximately plane-symmetric initially. Stated more precisely, it is proved that if the colliding almost-plane waves are initially sufficiently close to being exactly plane-symmetric across a bounded central region of sufficiently large transverse dimensions, then their collision necessarily produces spacetime singularities. In this case, nothing is now known about the local and global structures of the singularities.