Time-Domain TeraHertz Spectroscopy and Observational Probes of Prebiotic Interstellar Gas and Ice Chemistry

Understanding the origin of life on Earth has long fascinated the minds of the global community, and has been a driving factor in interdisciplinary research for centuries. Beyond the pioneering work of Darwin, perhaps the most widely known study in the last century is that of Miller and Urey, who examined the possibility of the formation of prebiotic chemical precursors on the primordial Earth [1]. More recent studies have shown that amino acids, the chemical building blocks of the biopolymers that comprise life as we know it on Earth, are present in meteoritic samples, and that the molecules extracted from the meteorites display isotopic signatures indicative of an extraterrestrial origin [2]. The most recent major discovery in this area has been the detection of glycine (NH₂CH₂COOH), the simplest amino acid, in pristine cometary samples returned by the NASA STARDUST mission [3]. Indeed, the open questions left by these discoveries, both in the public and scientific communities, hold such fascination that NASA has designated the understanding of our "Cosmic Origins" as a key mission priority.

Despite these exciting discoveries, our understanding of the chemical and physical pathways to the formation of prebiotic molecules is woefully incomplete. This is largely because we do not yet fully understand how the interplay between grain-surface and sub-surface ice reactions and the gas-phase affects astrophysical chemical evolution, and our knowledge of chemical inventories in these regions is incomplete. The research presented here aims to directly address both these issues, so that future work to understand the formation of prebiotic molecules has a solid foundation from which to work.

From an observational standpoint, a dedicated campaign to identify hydroxylamine (NH₂OH), potentially a direct precursor to glycine, in the gas-phase was undertaken. No trace of NH₂OH was found. These observations motivated a refinement of the chemical models of glycine formation, and have largely ruled out a gas-phase route to the synthesis of the simplest amino acid in the ISM. A molecular mystery in the case of the carrier of a series of transitions was resolved using observational data toward a large number of sources, confirming the identity of this important carbon-chemistry intermediate B11244 as l-C₃H⁺ and identifying it in at least two new environments. Finally, the doubly-nitrogenated molecule carbodiimide HNCNH was identified in the ISM for the first time through maser emission features in the centimeter-wavelength regime.

In the laboratory, a TeraHertz Time-Domain Spectrometer was constructed to obtain the experimental spectra necessary to search for solid-phase species in the ISM in the THz region of the spectrum. These investigations have shown a striking dependence on large-scale, long-range (i.e. lattice) structure of the ices on the spectra they present in the THz. A database of molecular spectra has been started, and both the simplest and most abundant ice species, which have already been identified, as well as a number of more complex species, have been studied. The exquisite sensitivity of the THz spectra to both the structure and thermal history of these ices may lead to better probes of complex chemical and dynamical evolution in interstellar environments.

Oxidation of Volatile Organic Compounds in Aqueous Solution and at the Air-water Interface of Aqueous Microdroplets

Isoprene (ISO),the most abundant non-methane VOC, is the major contributor to secondary organic aerosols (SOA) formation. The mechanisms involved in such transformation, however, are not fully understood. Current mechanisms, which are based on the oxidation of ISO in the gas-phase, underestimate SOA yields. The heightened awareness that ISO is only partially processed in the gas-phase has turned attention to heterogeneous processes as alternative pathways toward SOA.

During my research project, I investigated the photochemical oxidation of isoprene in bulk water. Below, I will report on the λ > 305 nm photolysis of H₂O₂ in dilute ISO solutions. This process yields C₁₀H₁₅OH species as primary products, whose formation both requires and is inhibited by O₂. Several isomers of C₁₀H₁₅OH were resolved by reverse-phase high-performance liquid chromatography and detected as MH⁺ (m/z = 153) and MH⁺-18 (m/z = 135) signals by electrospray ionization mass spectrometry. This finding is consistent with the addition of ·OH to ISO, followed by HO-ISO· reactions with ISO (in competition with O₂) leading to second generation HO(ISO)₂· radicals that terminate as C₁₀H₁₅OH via β-H abstraction by O₂.

It is not generally realized that chemistry on the surface of water cannot be deduced, extrapolated or translated to those in bulk gas and liquid phases. The water density drops a thousand-fold within a few Angstroms through the gas-liquid interfacial region and therefore hydrophobic VOCs such as ISO will likely remain in these relatively 'dry' interfacial water layers rather than proceed into bulk water. In previous experiments from our laboratory, it was found that gas-phase olefins can be protonated on the surface of pH < 4 water. This phenomenon increases the residence time of gases at the interface, an event that makes them increasingly susceptible to interaction with gaseous atmospheric oxidants such as ozone and hydroxyl radicals.

In order to test this hypothesis, I carried out experiments in which ISO(g) collides with the surface of aqueous microdroplets of various compositions. Herein I report that ISO(g) is oxidized into soluble species via Fenton chemistry on the surface of aqueous Fe(II)Cl₂ solutions simultaneously exposed to H₂O₂(g). Monomer and oligomeric species (ISO)1-8H⁺ were detected via online electrospray ionization mass spectrometry (ESI-MS) on the surface of pH ~ 2 water, and were then oxidized into a suite of products whose combined yields exceed ~ 5% of (ISO)1-8H⁺. MS/MS analysis revealed that products mainly consisted of alcohols, ketones, epoxides and acids. Our experiments demonstrated that olefins in ambient air may be oxidized upon impact on the surface of Fe-containing aqueous acidic media, such as those of typical to tropospheric aerosols.

Related experiments involving the reaction of ISO(g) with ·OH radicals from the photolysis of dissolved H₂O₂ were also carried out to test the surface oxidation of ISO(g) by photolyzing H₂O₂(aq) at 266 nm at various pH. The products were analyzed via online electrospray ionization mass spectrometry. Similar to our Fenton experiments, we detected (ISO)1-7H⁺ at pH < 4, and new m/z⁺ = 271 and m/z^- = 76 products at pH > 5.

Field and Laboratory Studies of Atmospheric Organic Aerosol

This thesis is the culmination of field and laboratory studies aimed at assessing processes that affect the composition and distribution of atmospheric organic aerosol. An emphasis is placed on measurements conducted using compact and high-resolution Aerodyne Aerosol Mass Spectrometers (AMS). The first three chapters summarize results from aircraft campaigns designed to evaluate anthropogenic and biogenic impacts on marine aerosol and clouds off the coast of California. Subsequent chapters describe laboratory studies intended to evaluate gas and particle-phase mechanisms of organic aerosol oxidation.

The 2013 Nucleation in California Experiment (NiCE) was a campaign designed to study environments impacted by nucleated and/or freshly formed aerosol particles. Terrestrial biogenic aerosol with > 85% organic mass was observed to reside in the free troposphere above marine stratocumulus. This biogenic organic aerosol (BOA) originated from the Northwestern United States and was transported to the marine atmosphere during periodic cloud-clearing events. Spectra recorded by a cloud condensation nuclei counter demonstrated that BOA is CCN active. BOA enhancements at latitudes north of San Francisco, CA coincided with enhanced cloud water concentrations of organic species such as acetate and formate.

Airborne measurements conducted during the 2011 Eastern Pacific Emitted Aerosol Cloud Experiment (E-PEACE) were aimed at evaluating the contribution of ship emissions to the properties of marine aerosol and clouds off the coast of central California. In one study, analysis of organic aerosol mass spectra during periods of enhanced shipping activity yielded unique tracers indicative of cloud-processed ship emissions (m/z 42 and 99). The variation of their organic fraction (f₄₂ and f₉₉) was found to coincide with periods of heavy (f₄₂ > 0.15; f₉₉ > 0.04), moderate (0.05 < f₄₂ < 0.15; 0.01 < f₉₉ < 0.04), and negligible (f₄₂ < 0.05; f₉₉ < 0.01) ship influence. Application of these conditions to all measurements conducted during E-PEACE demonstrated that a large fraction of cloud droplet (72%) and dry aerosol mass (12%) sampled in the California coastal study region was heavily or moderately influenced by ship emissions. Another study investigated the chemical and physical evolution of a controlled organic plume emitted from the R/V Point Sur. Under sunny conditions, nucleated particles composed of oxidized organic compounds contributed nearly an order of magnitude more cloud condensation nuclei (CCN) than less oxidized particles formed under cloudy conditions. The processing time necessary for particles to become CCN active was short ( < 1 hr) compared to the time needed for particles to become hygroscopic at sub-saturated humidity ( > 4 hr).

Laboratory chamber experiments were also conducted to evaluate particle-phase processes influencing aerosol phase and composition. In one study, ammonium sulfate seed was coated with a layer of secondary organic aerosol (SOA) from toluene oxidation followed by a layer of SOA from α-pinene oxidation. The system exhibited different evaporative properties than ammonium sulfate seed initially coated with α-pinene SOA followed by a layer of toluene SOA. This behavior is consistent with a shell-and-core model and suggests limited mixing among different SOA types. Another study investigated the reactive uptake of isoprene epoxy diols (IEPOX) onto non-acidified aerosol. It was demonstrated that particle acidity has limited influence on organic aerosol formation onto ammonium sulfate seed, and that the chemical system is limited by the availability of nucleophiles such as sulfate.

Flow tube experiments were conducted to examine the role of iron in the reactive uptake and chemical oxidation of glycolaldehyde. Aerosol particles doped with iron and hydrogen peroxide were mixed with gas-phase glycolaldehyde and photochemically aged in a custom-built flow reactor. Compared to particles free of iron, iron-doped aerosols significantly enhanced the oxygen to carbon (O/C) ratio of accumulated organic mass. The primary oxidation mechanism is suggested to be a combination of Fenton and photo-Fenton reactions which enhance particle-phase OH radical concentrations.

I. The attenuation and dispersion of sound in a condensing medium. II. The flow of a gas-particle mixture downstream of a normal shock in a nozzle

I. The attenuation of sound due to particles suspended in a gas was first calculated by Sewell and later by Epstein in their classical works on the propagation of sound in a two-phase medium. In their work, and in more recent works which include calculations of sound dispersion, the calculations were made for systems in which there was no mass transfer between the two phases. In the present work, mass transfer between phases is included in the calculations.

The attenuation and dispersion of sound in a two-phase condensing medium are calculated as functions of frequency. The medium in which the sound propagates consists of a gaseous phase, a mixture of inert gas and condensable vapor, which contains condensable liquid droplets. The droplets, which interact with the gaseous phase through the interchange of momentum, energy, and mass (through evaporation and condensation), are treated from the continuum viewpoint. Limiting cases, for flow either frozen or in equilibrium with respect to the various exchange processes, help demonstrate the effects of mass transfer between phases. Included in the calculation is the effect of thermal relaxation within droplets. Pressure relaxation between the two phases is examined, but is not included as a contributing factor because it is of interest only at much higher frequencies than the other relaxation processes. The results for a system typical of sodium droplets in sodium vapor are compared to calculations in which there is no mass exchange between phases. It is found that the maximum attenuation is about 25 per cent greater and occurs at about one-half the frequency for the case which includes mass transfer, and that the dispersion at low frequencies is about 35 per cent greater. Results for different values of latent heat are compared.

II. In the flow of a gas-particle mixture through a nozzle, a normal shock may exist in the diverging section of the nozzle. In Marble’s calculation for a shock in a constant area duct, the shock was described as a usual gas-dynamic shock followed by a relaxation zone in which the gas and particles return to equilibrium. The thickness of this zone, which is the total shock thickness in the gas-particle mixture, is of the order of the relaxation distance for a particle in the gas. In a nozzle, the area may change significantly over this relaxation zone so that the solution for a constant area duct is no longer adequate to describe the flow. In the present work, an asymptotic solution, which accounts for the area change, is obtained for the flow of a gas-particle mixture downstream of the shock in a nozzle, under the assumption of small slip between the particles and gas. This amounts to the assumption that the shock thickness is small compared with the length of the nozzle. The shock solution, valid in the region near the shock, is matched to the well known small-slip solution, which is valid in the flow downstream of the shock, to obtain a composite solution valid for the entire flow region. The solution is applied to a conical nozzle. A discussion of methods of finding the location of a shock in a nozzle is included.

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.

Proton transfers at the air-water interface

Proton transfer reactions at the interface of water with hydrophobic media, such as air or lipids, are ubiquitous on our planet. These reactions orchestrate a host of vital phenomena in the environment including, for example, acidification of clouds, enzymatic catalysis, chemistries of aerosol and atmospheric gases, and bioenergetic transduction. Despite their importance, however, quantitative details underlying these interactions have remained unclear. Deeper insight into these interfacial reactions is also required in addressing challenges in green chemistry, improved water quality, self-assembly of materials, the next generation of micro-nanofluidics, adhesives, coatings, catalysts, and electrodes. This thesis describes experimental and theoretical investigation of proton transfer reactions at the air-water interface as a function of hydration gradients, electrochemical potential, and electrostatics. Since emerging insights hold at the lipid-water interface as well, this work is also expected to aid understanding of complex biological phenomena associated with proton migration across membranes.

Based on our current understanding, it is known that the physicochemical properties of the gas-phase water are drastically different from those of bulk water. For example, the gas-phase hydronium ion, H₃O⁺(g), can protonate most (non-alkane) organic species, whereas H₃O⁺(aq) can neutralize only relatively strong bases. Thus, to be able to understand and engineer water-hydrophobe interfaces, it is imperative to investigate this fluctuating region of molecular thickness wherein the ‘function’ of chemical species transitions from one phase to another via steep gradients in hydration, dielectric constant, and density. Aqueous interfaces are difficult to approach by current experimental techniques because designing experiments to specifically sample interfacial layers (< 1 nm thick) is an arduous task. While recent advances in surface-specific spectroscopies have provided valuable information regarding the structure of aqueous interfaces, but structure alone is inadequate to decipher the function. By similar analogy, theoretical predictions based on classical molecular dynamics have remained limited in their scope.

Recently, we have adapted an analytical electrospray ionization mass spectrometer (ESIMS) for probing reactions at the gas-liquid interface in real time. This technique is direct, surface-specific,and provides unambiguous mass-to-charge ratios of interfacial species. With this innovation, we have been able to investigate the following:

1. How do anions mediate proton transfers at the air-water interface?

2. What is the basis for the negative surface potential at the air-water interface?

3. What is the mechanism for catalysis ‘on-water’?

In addition to our experiments with the ESIMS, we applied quantum mechanics and molecular dynamics to simulate our experiments toward gaining insight at the molecular scale. Our results unambiguously demonstrated the role of electrostatic-reorganization of interfacial water during proton transfer events. With our experimental and theoretical results on the ‘superacidity’ of the surface of mildly acidic water, we also explored implications on atmospheric chemistry and green chemistry. Our most recent results explained the basis for the negative charge of the air-water interface and showed that the water-hydrophobe interface could serve as a site for enhanced autodissociation of water compared to the condensed phase.

Studies of ambient organic and inorganic aerosol in Southern California

The negative impacts of ambient aerosol particles, or particulate matter (PM), on human health and climate are well recognized. However, owing to the complexity of aerosol particle formation and chemical evolution, emissions control strategies remain difficult to develop in a cost effective manner. In this work, three studies are presented to address several key issues currently stymieing California's efforts to continue improving its air quality.

Gas-phase organic mass (GPOM) and CO emission factors are used in conjunction with measured enhancements in oxygenated organic aerosol (OOA) relative to CO to quantify the significant lack of closure between expected and observed organic aerosol concentrations attributable to fossil-fuel emissions. Two possible conclusions emerge from the analysis to yield consistency with the ambient organic data: (1) vehicular emissions are not a dominant source of anthropogenic fossil SOA in the Los Angeles Basin, or (2) the ambient SOA mass yields used to determine the SOA formation potential of vehicular emissions are substantially higher than those derived from laboratory chamber studies. Additional laboratory chamber studies confirm that, owing to vapor-phase wall loss, the SOA mass yields currently used in virtually all 3D chemical transport models are biased low by as much as a factor of 4. Furthermore, predictions from the Statistical Oxidation Model suggest that this bias could be as high as a factor of 8 if the influence of the chamber walls could be removed entirely.

Once vapor-phase wall loss has been accounted for in a new suite of laboratory chamber experiments, the SOA parameterizations within atmospheric chemical transport models should also be updated. To address the numerical challenges of implementing the next generation of SOA models in atmospheric chemical transport models, a novel mathematical framework, termed the Moment Method, is designed and presented. Assessment of the Moment Method strengths and weaknesses provide valuable insight that can guide future development of SOA modules for atmospheric CTMs.

Finally, regional inorganic aerosol formation and evolution is investigated via detailed comparison of predictions from the Community Multiscale Air Quality (CMAQ version 4.7.1) model against a suite of airborne and ground-based meteorological measurements, gas- and aerosol-phase inorganic measurements, and black carbon (BC) measurements over Southern California during the CalNex field campaign in May/June 2010. Results suggests that continuing to target sulfur emissions with the hopes of reducing ambient PM concentrations may not the most effective strategy for Southern California. Instead, targeting dairy emissions is likely to be an effective strategy for substantially reducing ammonium nitrate concentrations in the eastern part of the Los Angeles Basin.

Geochemistry and resonance ionization of platinum-group elements

Experimental studies were conducted with the goals of 1) determining the origin of Pt- group element (PGE) alloys and associated mineral assemblages in refractory inclusions from meteorites and 2) developing a new ultrasensitive method for the in situ chemical and isotopic analysis of PGE. A general review of the geochemistry and cosmochemistry of the PGE is given, and specific research contributions are presented within the context of this broad framework.

An important step toward understanding the cosmochemistry of the PGE is the determination of the origin of POE-rich metallic phases (most commonly εRu-Fe) that are found in Ca, AJ-rich refractory inclusions (CAI) in C3V meteorites. These metals occur along with γNi-Fe metals, Ni-Fe sulfides and Fe oxides in multiphase opaque assemblages. Laboratory experiments were used to show that the mineral assemblages and textures observed in opaque assemblages could be produced by sulfidation and oxidation of once homogeneous Ni-Fe-PGE metals. Phase equilibria, partitioning and diffusion kinetics were studied in the Ni-Fe-Ru system in order to quantify the conditions of opaque assemblage formation. Phase boundaries and tie lines in the Ni-Fe-Ru system were determined at 1273, 1073 and 873K using an experimental technique that allowed the investigation of a large portion of the Ni-Fe-Ru system with a single experiment at each temperature by establishing a concentration gradient within which local equilibrium between coexisting phases was maintained. A wide miscibility gap was found to be present at each temperature, separating a hexagonal close-packed εRu-Fe phase from a face-centered cubic γNi-Fe phase. Phase equilibria determined here for the Ni-Fe-Ru system, and phase equilibria from the literature for the Ni-Fe-S and Ni-Fe-O systems, were compared with analyses of minerals from opaque assemblages to estimate the temperature and chemical conditions of opaque assemblage formation. It was determined that opaque assemblages equilibrated at a temperature of ~770K, a sulfur fugacity 10 times higher than an equilibrium solar gas, and an oxygen fugacity 10⁶ times higher than an equilibrium solar gas.

Diffusion rates between -γNi-Fe and εRu-Fe metal play a critical role in determining the time (with respect to CAI petrogenesis) and duration of the opaque assemblage equilibration process. The diffusion coefficient for Ru in Ni (D^Ru_Ni) was determined as an analog for the Ni-Fe-Ru system by the thin-film diffusion method in the temperature range of 1073 to 1673K and is given by the expression:

D^Ru_Ni (cm² sec^-1) = 5.0(±0.7) x 10^-3 exp(-2.3(±0.1) x 10¹² erg mole^-1/RT) where R is the gas constant and T is the temperature in K. Based on the rates of dissolution and exsolution of metallic phases in the Ni-Fe-Ru system it is suggested that opaque assemblages equilibrated after the melting and crystallization of host CAI during a metamorphic event of ≥ 10³ years duration. It is inferred that opaque assemblages originated as immiscible metallic liquid droplets in the CAI silicate liquid. The bulk compositions of PGE in these precursor alloys reflects an early stage of condensation from the solar nebula and the partitioning of V between the precursor alloys and CAI silicate liquid reflects the reducing nebular conditions under which CAI were melted. The individual mineral phases now observed in opaque assemblages do not preserve an independent history prior to CAI melting and crystallization, but instead provide important information on the post-accretionary history of C3V meteorites and allow the quantification of the temperature, sulfur fugacity and oxygen fugacity of cooling planetary environments. This contrasts with previous models that called upon the formation of opaque assemblages by aggregation of phases that formed independently under highly variable conditions in the solar nebula prior to the crystallization of CAI.

Analytical studies were carried out on PGE-rich phases from meteorites and the products of synthetic experiments using traditional electron microprobe x-ray analytical techniques. The concentrations of PGE in common minerals from meteorites and terrestrial rocks are far below the ~100 ppm detection limit of the electron microprobe. This has limited the scope of analytical studies to the very few cases where PGE are unusually enriched. To study the distribution of PGE in common minerals will require an in situ analytical technique with much lower detection limits than any methods currently in use. To overcome this limitation, resonance ionization of sputtered atoms was investigated for use as an ultrasensitive in situ analytical technique for the analysis of PGE. The mass spectrometric analysis of Os and Re was investigated using a pulsed primary Ar⁺ ion beam to provide sputtered atoms for resonance ionization mass spectrometry. An ionization scheme for Os that utilizes three resonant energy levels (including an autoionizing energy level) was investigated and found to have superior sensitivity and selectivity compared to nonresonant and one and two energy level resonant ionization schemes. An elemental selectivity for Os over Re of ≥ 10³ was demonstrated. It was found that detuning the ionizing laser from the autoionizing energy level to an arbitrary region in the ionization continuum resulted in a five-fold decrease in signal intensity and a ten-fold decrease in elemental selectivity. Osmium concentrations in synthetic metals and iron meteorites were measured to demonstrate the analytical capabilities of the technique. A linear correlation between Os⁺ signal intensity and the known Os concentration was observed over a range of nearly 10⁴ in Os concentration with an accuracy of ~ ±10%, a millimum detection limit of 7 parts per billion atomic, and a useful yield of 1%. Resonance ionization of sputtered atoms samples the dominant neutral-fraction of sputtered atoms and utilizes multiphoton resonance ionization to achieve high sensitivity and to eliminate atomic and molecular interferences. Matrix effects should be small compared to secondary ion mass spectrometry because ionization occurs in the gas phase and is largely independent of the physical properties of the matrix material. Resonance ionization of sputtered atoms can be applied to in situ chemical analysis of most high ionization potential elements (including all of the PGE) in a wide range of natural and synthetic materials. The high useful yield and elemental selectivity of this method should eventually allow the in situ measurement of Os isotope ratios in some natural samples and in sample extracts enriched in PGE by fire assay fusion.

Phase equilibria and diffusion experiments have provided the basis for a reinterpretation of the origin of opaque assemblages in CAI and have yielded quantitative information on conditions in the primitive solar nebula and cooling planetary environments. Development of the method of resonance ionization of sputtered atoms for the analysis of Os has shown that this technique has wide applications in geochemistry and will for the first time allow in situ studies of the distribution of PGE at the low concentration levels at which they occur in common minerals.

Force chains, friction, and flow: Behavior of granular media across length scales

We study the behavior of granular materials at three length scales. At the smallest length scale, the grain-scale, we study inter-particle forces and "force chains". Inter-particle forces are the natural building blocks of constitutive laws for granular materials. Force chains are a key signature of the heterogeneity of granular systems. Despite their fundamental importance for calibrating grain-scale numerical models and elucidating constitutive laws, inter-particle forces have not been fully quantified in natural granular materials. We present a numerical force inference technique for determining inter-particle forces from experimental data and apply the technique to two-dimensional and three-dimensional systems under quasi-static and dynamic load. These experiments validate the technique and provide insight into the quasi-static and dynamic behavior of granular materials.

At a larger length scale, the mesoscale, we study the emergent frictional behavior of a collection of grains. Properties of granular materials at this intermediate scale are crucial inputs for macro-scale continuum models. We derive friction laws for granular materials at the mesoscale by applying averaging techniques to grain-scale quantities. These laws portray the nature of steady-state frictional strength as a competition between steady-state dilation and grain-scale dissipation rates. The laws also directly link the rate of dilation to the non-steady-state frictional strength.

At the macro-scale, we investigate continuum modeling techniques capable of simulating the distinct solid-like, liquid-like, and gas-like behaviors exhibited by granular materials in a single computational domain. We propose a Smoothed Particle Hydrodynamics (SPH) approach for granular materials with a viscoplastic constitutive law. The constitutive law uses a rate-dependent and dilation-dependent friction law. We provide a theoretical basis for a dilation-dependent friction law using similar analysis to that performed at the mesoscale. We provide several qualitative and quantitative validations of the technique and discuss ongoing work aiming to couple the granular flow with gas and fluid flows.

Secondary Organic Aerosol Composition Studies Using Mass Spectrometry

Trace volatile organic compounds emitted by biogenic and anthropogenic sources into the atmosphere can undergo extensive photooxidation to form species with lower volatility. By equilibrium partitioning or reactive uptake, these compounds can nucleate into new aerosol particles or deposit onto already-existing particles to form secondary organic aerosol (SOA). SOA and other atmospheric particulate matter have measurable effects on global climate and public health, making understanding SOA formation a needed field of scientific inquiry. SOA formation can be done in a laboratory setting, using an environmental chamber; under these controlled conditions it is possible to generate SOA from a single parent compound and study the chemical composition of the gas and particle phases. By studying the SOA composition, it is possible to gain understanding of the chemical reactions that occur in the gas phase and particle phase, and identify potential heterogeneous processes that occur at the surface of SOA particles. In this thesis, mass spectrometric methods are used to identify qualitatively and qualitatively the chemical components of SOA derived from the photooxidation of important anthropogenic volatile organic compounds that are associated with gasoline and diesel fuels and industrial activity (C12 alkanes, toluene, and o-, m-, and p-cresols). The conditions under which SOA was generated in each system were varied to explore the effect of NO_x and inorganic seed composition on SOA chemical composition. The structure of the parent alkane was varied to investigate the effect on the functionalization and fragmentation of the resulting oxidation products. Relative humidity was varied in the alkane system as well to measure the effect of increased particle-phase water on condensed-phase reactions. In all systems, oligomeric species, resulting potentially from particle-phase and heterogeneous processes, were identified. Imines produced by reactions between (NH₄)₂SO₄ seed and carbonyl compounds were identified in all systems. Multigenerational photochemistry producing low- and extremely low-volatility organic compounds (LVOC and ELVOC) was reflected strongly in the particle-phase composition as well.

Theory of the positive column of a gas discharge

An attempt is made to provide a theoretical explanation of the effect of the positive column on the voltage-current characteristic of a glow or an arc discharge. Such theories have been developed before, and all are based on balancing the production and loss of charged particles and accounting for the energy supplied to the plasma by the applied electric field. Differences among the theories arise from the approximations and omissions made in selecting processes that affect the particle and energy balances. This work is primarily concerned with the deviation from the ambipolar description of the positive column caused by space charge, electron-ion volume recombination, and temperature inhomogeneities.

The presentation is divided into three parts, the first of which involved the derivation of the final macroscopic equations from kinetic theory. The final equations are obtained by taking the first three moments of the Boltzmann equation for each of the three species in the plasma. Although the method used and the equations obtained are not novel, the derivation is carried out in detail in order to appraise the validity of numerous approximations and to justify the use of data from other sources. The equations are applied to a molecular hydrogen discharge contained between parallel walls. The applied electric field is parallel to the walls, and the dependent variables—electron and ion flux to the walls, electron and ion densities, transverse electric field, and gas temperature—vary only in the direction perpendicular to the walls. The mathematical description is given by a sixth-order nonlinear two-point boundary value problem which contains the applied field as a parameter. The amount of neutral gas and its temperature at the walls are held fixed, and the relation between the applied field and the electron density at the center of the discharge is obtained in the process of solving the problem. This relation corresponds to that between current and voltage and is used to interpret the effect of space charge, recombination, and temperature inhomogeneities on the voltage-current characteristic of the discharge.

The complete solution of the equations is impractical both numerically and analytically, and in Part II the gas temperature is assumed uniform so as to focus on the combined effects of space charge and recombination. The terms representing these effects are treated as perturbations to equations that would otherwise describe the ambipolar situation. However, the term representing space charge is not negligible in a thin boundary layer or sheath near the walls, and consequently the perturbation problem is singular. Separate solutions must be obtained in the sheath and in the main region of the discharge, and the relation between the electron density and the applied field is not determined until these solutions are matched.

In Part III the electron and ion densities are assumed equal, and the complicated space-charge calculation is thereby replaced by the ambipolar description. Recombination and temperature inhomogeneities are both important at high values of the electron density. However, the formulation of the problem permits a comparison of the relative effects, and temperature inhomogeneities are shown to be important at lower values of the electron density than recombination. The equations are solved by a direct numerical integration and by treating the term representing temperature inhomogeneities as a perturbation.

The conclusions reached in the study are primarily concerned with the association of the relation between electron density and axial field with the voltage-current characteristic. It is known that the effect of space charge can account for the subnormal glow discharge and that the normal glow corresponds to a close approach to an ambipolar situation. The effect of temperature inhomogeneities helps explain the decreasing characteristic of the arc, and the effect of recombination is not expected to appear except at very high electron densities.

A study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas

This thesis is a study of nonlinear phenomena in the propagation of electromagnetic waves in a weakly ionized gas externally biased with a magnetostatic field. The present study is restricted to the nonlinear phenomena rising from the interaction of electromagnetic waves in the ionized gas. The important effects of nonlinearity are wave-form distortion leads to cross modulation of one wave by a second amplitude-modulated wave.

The nonlinear effects are assumed to be small so that a perturbation method can be used. Boltzmann’s kinetic equation with an appropriate expression for the collision term is solved by expanding the electron distribution function into spherical harmonics in velocity space. In turn, the electron convection current density and the conductivity tensors of the nonlinear ionized gas are found from the distribution function. Finally, the expression for the current density and Maxwell’s equations are employed to investigate the effects of nonlinearity on the propagation of electromagnetic waves in the ionized gas, and also on the reflection of waves from an ionized gas of semi-infinite extent.

Chemistry of secondary organic aerosol

The photooxidation of volatile organic compounds (VOCs) in the atmosphere can lead to the formation of secondary organic aerosol (SOA), a major component of fine particulate matter. Improvements to air quality require insight into the many reactive intermediates that lead to SOA formation, of which only a small fraction have been measured at the molecular level. This thesis describes the chemistry of secondary organic aerosol (SOA) formation from several atmospherically relevant hydrocarbon precursors. Photooxidation experiments of methoxyphenol and phenolic compounds and C₁₂ alkanes were conducted in the Caltech Environmental Chamber. These experiments include the first photooxidation studies of these precursors run under sufficiently low NO_x levels, such that RO₂ + HO₂ chemistry dominates, an important chemical regime in the atmosphere. Using online Chemical Ionization Mass Spectrometery (CIMS), key gas-phase intermediates that lead to SOA formation in these systems were identified. With complementary particle-phase analyses, chemical mechanisms elucidating the SOA formation from these compounds are proposed.

Three methoxyphenol species (phenol, guaiacol, and syringol) were studied to model potential photooxidation schemes of biomass burning intermediates. SOA yields (ratio of mass of SOA formed to mass of primary organic reacted) exceeding 25% are observed. Aerosol growth is rapid and linear with the organic conversion, consistent with the formation of essentially non-volatile products. Gas and aerosol-phase oxidation products from the guaiacol system show that the chemical mechanism consists of highly oxidized aromatic species in the particle phase. Syringol SOA yields are lower than that of phenol and guaiacol, likely due to unique chemistry dependent on methoxy group position.

The photooxidation of several C₁₂ alkanes of varying structure n-dodecane, 2-methylundecane, cyclododecane, and hexylcyclohexane) were run under extended OH exposure to investigate the effect of molecular structure on SOA yields and photochemical aging. Peroxyhemiacetal formation from the reactions of several multifunctional hydroperoxides and aldehyde intermediates was found to be central to organic growth in all systems, and SOA yields increased with cyclic character of the starting hydrocarbon. All of these studies provide direction for future experiments and modeling in order to lessen outstanding discrepancies between predicted and measured SOA.

Detailed properties of high redshift galaxies

Galaxies evolve throughout the history of the universe from the first star-forming sources, through gas-rich asymmetric structures with rapid star formation rates, to the massive symmetrical stellar systems observed at the present day. Determining the physical processes which drive galaxy formation and evolution is one of the most important questions in observational astrophysics. This thesis presents four projects aimed at improving our understanding of galaxy evolution from detailed measurements of star forming galaxies at high redshift.

We use resolved spectroscopy of gravitationally lensed z ≃ 2 - 3 star forming galaxies to measure their kinematic and star formation properties. The combination of lensing with adaptive optics yields physical resolution of ≃ 100 pc, sufficient to resolve giant Hii regions. We find that ~ 70 % of galaxies in our sample display ordered rotation with high local velocity dispersion indicating turbulent thick disks. The rotating galaxies are gravitationally unstable and are expected to fragment into giant clumps. The size and dynamical mass of giant Hii regions are in agreement with predictions for such clumps indicating that gravitational instability drives the rapid star formation. The remainder of our sample is comprised of ongoing major mergers. Merging galaxies display similar star formation rate, morphology, and local velocity dispersion as isolated sources, but their velocity fields are more chaotic with no coherent rotation.

We measure resolved metallicity in four lensed galaxies at z = 2.0 − 2.4 from optical emission line diagnostics. Three rotating galaxies display radial gradients with higher metallicity at smaller radii, while the fourth is undergoing a merger and has an inverted gradient with lower metallicity at the center. Strong gradients in the rotating galaxies indicate that they are growing inside-out with star formation fueled by accretion of metal-poor gas at large radii. By comparing measured gradients with an appropriate comparison sample at z = 0, we demonstrate that metallicity gradients in isolated galaxies must flatten at later times. The amount of size growth inferred by the gradients is in rough agreement with direct measurements of massive galaxies. We develop a chemical evolution model to interpret these data and conclude that metallicity gradients are established by a gradient in the outflow mass loading factor, combined with radial inflow of metal-enriched gas.

We present the first rest-frame optical spectroscopic survey of a large sample of low-luminosity galaxies at high redshift (L < L*, 1.5 < z < 3.5). This population dominates the star formation density of the universe at high redshifts, yet such galaxies are normally too faint to be studied spectroscopically. We take advantage of strong gravitational lensing magnification to compile observations for a sample of 29 galaxies using modest integration times with the Keck and Palomar telescopes. Balmer emission lines confirm that the sample has a median SFR ∼ 10 M_sun yr^−1 and extends to lower SFR than has been probed by other surveys at similar redshift. We derive the metallicity, dust extinction, SFR, ionization parameter, and dynamical mass from the spectroscopic data, providing the first accurate characterization of the star-forming environment in low-luminosity galaxies at high redshift. For the first time, we directly test the proposal that the relation between galaxy stellar mass, star formation rate, and gas phase metallicity does not evolve. We find lower gas phase metallicity in the high redshift galaxies than in local sources with equivalent stellar mass and star formation rate, arguing against a time-invariant relation. While our result is preliminary and may be biased by measurement errors, this represents an important first measurement that will be further constrained by ongoing analysis of the full data set and by future observations.

We present a study of composite rest-frame ultraviolet spectra of Lyman break galaxies at z = 4 and discuss implications for the distribution of neutral outflowing gas in the circumgalactic medium. In general we find similar spectroscopic trends to those found at z = 3 by earlier surveys. In particular, absorption lines which trace neutral gas are weaker in less evolved galaxies with lower stellar masses, smaller radii, lower luminosity, less dust, and stronger Lyα emission. Typical galaxies are thus expected to have stronger Lyα emission and weaker low-ionization absorption at earlier times, and we indeed find somewhat weaker low-ionization absorption at higher redshifts. In conjunction with earlier results, we argue that the reduced low-ionization absorption is likely caused by lower covering fraction and/or velocity range of outflowing neutral gas at earlier epochs. This result has important implications for the hypothesis that early galaxies were responsible for cosmic reionization. We additionally show that fine structure emission lines are sensitive to the spatial extent of neutral gas, and demonstrate that neutral gas is concentrated at smaller galactocentric radii in higher redshift galaxies.

The results of this thesis present a coherent picture of galaxy evolution at high redshifts 2 ≲ z ≲ 4. Roughly 1/3 of massive star forming galaxies at this period are undergoing major mergers, while the rest are growing inside-out with star formation occurring in gravitationally unstable thick disks. Star formation, stellar mass, and metallicity are limited by outflows which create a circumgalactic medium of metal-enriched material. We conclude by describing some remaining open questions and prospects for improving our understanding of galaxy evolution with future observations of gravitationally lensed galaxies.

Grain growth in protoplanetary disks

The majority of young, low-mass stars are surrounded by optically thick accretion disks. These circumstellar disks provide large reservoirs of gas and dust that will eventually be transformed into planetary systems. Theory and observations suggest that the earliest stage toward planet formation in a protoplanetary disk is the growth of particles, from sub-micron-sized grains to centimeter- sized pebbles. Theory indicates that small interstellar grains are well coupled into the gas and are incorporated to the disk during the proto-stellar collapse. These dust particles settle toward the disk mid-plane and simultaneously grow through collisional coagulation in a very short timescale. Observationally, grain growth can be inferred by measuring the spectral energy distribution at long wavelengths, which traces the continuum dust emission spectrum and hence the dust opacity. Several observational studies have indicated that the dust component in protoplanetary disks has evolved as compared to interstellar medium dust particles, suggesting at least 4 orders of magnitude in particle- size growth. However, the limited angular resolution and poor sensitivity of previous observations has not allowed for further exploration of this astrophysical process.

As part of my thesis, I embarked in an observational program to search for evidence of radial variations in the dust properties across a protoplanetary disk, which may be indicative of grain growth. By making use of high angular resolution observations obtained with CARMA, VLA, and SMA, I searched for radial variations in the dust opacity inside protoplanetary disks. These observations span more than an order of magnitude in wavelength (from sub-millimeter to centimeter wavelengths) and attain spatial resolutions down to 20 AU. I characterized the radial distribution of the circumstellar material and constrained radial variations of the dust opacity spectral index, which may originate from particle growth in these circumstellar disks. Furthermore, I compared these observational constraints with simple physical models of grain evolution that include collisional coagulation, fragmentation, and the interaction of these grains with the gaseous disk (the radial drift problem). For the parameters explored, these observational constraints are in agreement with a population of grains limited in size by radial drift. Finally, I also discuss future endeavors with forthcoming ALMA observations.