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.

Part I. Energy straggling of ^4He below 2.0 MeV in Al, Ni, Au and Pt. Part II. Studies of the Ti-W metallization system on Si

Part I.

In recent years, backscattering spectrometry has become an important tool for the analysis of thin films. An inherent limitation, though, is the loss of depth resolution due to energy straggling of the beam. To investigate this, energy straggling of ⁴He has been measured in thin films of Ni, Al, Au and Pt. Straggling is roughly proportional to square root of thickness, appears to have a slight energy dependence and generally decreases with decreasing atomic number of the adsorber. The results are compared with predictions of theory and with previous measurements. While Ni measurements are in fair agreement with Bohr's theory, Al measurements are 30% above and Au measurements are 40% below predicted values. The Au and Pt measurements give straggling values which are close to one another.

Part II.

MeV backscattering spectrometry and X-ray diffraction are used to investigate the behavior of sputter-deposited Ti-W mixed films on Si substrates. During vacuum anneals at temperatures near 700°C for several hours, the metallization layer reacts with the substrate. Backscattering analysis shows that the resulting compound layer is uniform in composition and contains Ti, Wand Si. The Ti:W ratio in the compound corresponds to that of the deposited metal film. X-ray analyses with Reed and Guinier cameras reveal the presence of the ternary Ti_xW_(1-x)Si₂ compound. Its composition is unaffected by oxygen contamination during annealing, but the reaction rate is affected. The rate measured on samples with about 15% oxygen contamination after annealing is linear, of the order of 0.5 Å per second at 725°C, and depends on the crystallographic orientation of the substrate and the dc bias during sputter-deposition of the Ti-W film.

Au layers of about 1000 Å thickness were deposited onto unreacted Ti-W films on Si. When annealed at 400°C these samples underwent a color change,and SEM micrographs of the samples showed that an intricate pattern of fissures which were typically 3µm wide had evolved. Analysis by electron microprobe revealed that Au had segregated preferentially into the fissures. This result suggests that Ti-W is not a barrier to Au-Si intermixing at 400°C.

Studies of the equation of state and elasticity of mantle minerals

Because the Earth’s upper mantle is inaccessible to us, in order to understand the chemical and physical processes that occur in the Earth’s interior we must rely on both experimental work and computational modeling. This thesis addresses both of these geochemical methods. In the first chapter, I develop an internally consistent comprehensive molar volume model for spinels in the oxide system FeO-MgO-Fe₂O₃-Cr₂O₃-Al₂O₃-TiO₂. The model is compared to the current MELTS spinel model with a demonstration of the impact of the model difference on the estimated spinel-garnet lherzolite transition pressure. In the second chapter, I calibrate a molar volume model for cubic garnets in the system SiO₂-Al₂O₃-TiO₂-Fe₂O₃-Cr₂O₃-FeO-MnO-MgO-CaO-Na₂O. I use the method of singular value analysis to calibrate excess volume of mixing parameters for the garnet model. The implications the model has for the density of the lithospheric mantle are explored. In the third chapter, I discuss the nuclear inelastic X-ray scattering (NRIXS) method, and present analysis of three orthopyroxene samples with different Fe contents. Longitudinal and shear wave velocities, elastic parameters, and other thermodynamic information are extracted from the raw NRIXS data.

Structure, chemistry, and energetics of organic and inorganic adsorbates on Ga-rich GaAs and GaP(001) surfaces

The work described in this dissertation includes fundamental investigations into three surface processes, namely inorganic film growth, water-induced oxidation, and organic functionalization/passivation, on the GaP and GaAs(001) surfaces. The techniques used to carry out this work include scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. Atomic structure, electronic structure, reaction mechanisms, and energetics related to these surface processes are discussed at atomic or molecular levels.

First, we investigate epitaxial Zn₃P₂ films grown on the Ga-rich GaAs(001)(6×6) surface. The film growth mechanism, electronic properties, and atomic structure of the Zn₃P₂/GaAs(001) system are discussed based on experimental and theoretical observations. We discover that a P-rich amorphous layer covers the crystalline Zn3P2 film during and after growth. We also propose more accurate picture of the GaP interfacial layer between Zn₃P₂ and GaAs, based on the atomic structure, chemical bonding, band diagram, and P-replacement energetics, than was previously anticipated.

Second, DFT calculations are carried out in order to understand water-induced oxidation mechanisms on the Ga-rich GaP(001)(2×4) surface. Structural and energetic information of every step in the gaseous water-induced GaP oxidation reactions are elucidated at the atomic level in great detail. We explore all reasonable ground states involved in most of the possible adsorption and decomposition pathways. We also investigate structures and energies of the transition states in the first hydrogen dissociation of a water molecule on the (2×4) surface.

Finally, adsorption structures and thermal decomposition reactions of 1-propanethiol on the Ga-rich GaP(001)(2×4) surface are investigated using high resolution STM, XPS, and DFT simulations. We elucidate adsorption locations and their associated atomic structures of a single 1-propanethiol molecule on the (2×4) surface as a function of annealing temperature. DFT calculations are carried out to optimize ground state structures and search transition states. XPS is used to investigate variations of the chemical bonding nature and coverage of the adsorbate species.

The needle in the hundred square degree haystack: The hunt for binary neutron star mergers with LIGO and Palomar Transient Factory

The Advanced LIGO and Virgo experiments are poised to detect gravitational waves (GWs) directly for the first time this decade. The ultimate prize will be joint observation of a compact binary merger in both gravitational and electromagnetic channels. However, GW sky locations that are uncertain by hundreds of square degrees will pose a challenge. I describe a real-time detection pipeline and a rapid Bayesian parameter estimation code that will make it possible to search promptly for optical counterparts in Advanced LIGO. Having analyzed a comprehensive population of simulated GW sources, we describe the sky localization accuracy that the GW detector network will achieve as each detector comes online and progresses toward design sensitivity. Next, in preparation for the optical search with the intermediate Palomar Transient Factory (iPTF), we have developed a unique capability to detect optical afterglows of gamma-ray bursts (GRBs) detected by the Fermi Gamma-ray Burst Monitor (GBM). Its comparable error regions offer a close parallel to the Advanced LIGO problem, but Fermi's unique access to MeV-GeV photons and its near all-sky coverage may allow us to look at optical afterglows in a relatively unexplored part of the GRB parameter space. We present the discovery and broadband follow-up observations (X-ray, UV, optical, millimeter, and radio) of eight GBM-IPTF afterglows. Two of the bursts (GRB 130702A / iPTF13bxl and GRB 140606B / iPTF14bfu) are at low redshift (z=0.145 and z = 0.384, respectively), are sub-luminous with respect to "standard" cosmological bursts, and have spectroscopically confirmed broad-line type Ic supernovae. These two bursts are possibly consistent with mildly relativistic shocks breaking out from the progenitor envelopes rather than the standard mechanism of internal shocks within an ultra-relativistic jet. On a technical level, the GBM--IPTF effort is a prototype for locating and observing optical counterparts of GW events in Advanced LIGO with the Zwicky Transient Facility.