Mechanistic studies of reactions at the single-molecule level using microfluidics with applications in molecular diagnostics

Motivated by needs in molecular diagnostics and advances in microfabrication, researchers started to seek help from microfluidic technology, as it provides approaches to achieve high throughput, high sensitivity, and high resolution. One strategy applied in microfluidics to fulfill such requirements is to convert continuous analog signal into digitalized signal. One most commonly used example for this conversion is digital PCR, where by counting the number of reacted compartments (triggered by the presence of the target entity) out of the total number of compartments, one could use Poisson statistics to calculate the amount of input target.

However, there are still problems to be solved and assumptions to be validated before the technology is widely employed. In this dissertation, the digital quantification strategy has been examined from two angles: efficiency and robustness. The former is a critical factor for ensuring the accuracy of absolute quantification methods, and the latter is the premise for such technology to be practically implemented in diagnosis beyond the laboratory. The two angles are further framed into a “fate” and “rate” determination scheme, where the influence of different parameters is attributed to fate determination step or rate determination step. In this discussion, microfluidic platforms have been used to understand reaction mechanism at single molecule level. Although the discussion raises more challenges for digital assay development, it brings the problem to the attention of the scientific community for the first time.

This dissertation also contributes towards developing POC test in limited resource settings. On one hand, it adds ease of access to the tests by incorporating massively producible, low cost plastic material and by integrating new features that allow instant result acquisition and result feedback. On the other hand, it explores new isothermal chemistry and new strategies to address important global health concerns such as cyctatin C quantification, HIV/HCV detection and treatment monitoring as well as HCV genotyping.

Studies toward the total synthesis of ritterazine B : the investigation of alkynylation reactions for use in the synthesis of the western fragment

The ritterazine and cephalostatin natural products have biological activities and structures that are interesting to synthetic organic chemists. These products have been found to exhibit significant cytotoxicity against P388 murine leukemia cells, and therefore have the potential to be used as anticancer drugs. The ritterazines and cephalostatins are steroidal dimers joined by a central pyrazine ring. Given that the steroid halves are unsymmetrical and highly oxygenated, there are several challenges in synthesizing these compounds in an organic laboratory.

Ritterazine B is the most potent derivative in the ritterazine family. Its biological activity is comparable to drugs that are being used to treat cancer today. For this reason, and the fact that there are no reported syntheses of ritterazine B to date, our lab set out to synthesize this natural product.

Herein, efforts toward the synthesis of the western fragment of ritterazine B are described. Two different routes are explored to access a common intermediate. An alkyne conjugate addition reaction was initially investigated due to the success of this key reaction in the synthesis of the eastern fragment. However, it has been found that a propargylation reaction has greater reactivity and yields, and has the potential to reduce the step count of the synthesis of the western fragment of ritterazine B.

Advancing EDL technologies for future space missions : from ground testing facilities to ablative heatshields

Motivated by recent MSL results where the ablation rate of the PICA heatshield was over-predicted, and staying true to the objectives outlined in the NASA Space Technology Roadmaps and Priorities report, this work focuses on advancing EDL technologies for future space missions.

Due to the difficulties in performing flight tests in the hypervelocity regime, a new ground testing facility called the vertical expansion tunnel is proposed. The adverse effects from secondary diaphragm rupture in an expansion tunnel may be reduced or eliminated by orienting the tunnel vertically, matching the test gas pressure and the accelerator gas pressure, and initially separating the test gas from the accelerator gas by density stratification. If some sacrifice of the reservoir conditions can be made, the VET can be utilized in hypervelocity ground testing, without the problems associated with secondary diaphragm rupture.

The performance of different constraints for the Rate-Controlled Constrained-Equilibrium (RCCE) method is investigated in the context of modeling reacting flows characteristic to ground testing facilities, and re-entry conditions. The effectiveness of different constraints are isolated, and new constraints previously unmentioned in the literature are introduced. Three main benefits from the RCCE method were determined: 1) the reduction in number of equations that need to be solved to model a reacting flow; 2) the reduction in stiffness of the system of equations needed to be solved; and 3) the ability to tabulate chemical properties as a function of a constraint once, prior to running a simulation, along with the ability to use the same table for multiple simulations.

Finally, published physical properties of PICA are compiled, and the composition of the pyrolysis gases that form at high temperatures internal to a heatshield is investigated. A necessary link between the composition of the solid resin, and the composition of the pyrolysis gases created is provided. This link, combined with a detailed investigation into a reacting pyrolysis gas mixture, allows a much needed consistent, and thorough description of many of the physical phenomena occurring in a PICA heatshield, and their implications, to be presented.

Through the use of computational fluid mechanics and computational chemistry methods, significant contributions have been made to advancing ground testing facilities, computational methods for reacting flows, and ablation modeling.

Kinetic investigations of heterogeneously catalyzed reactions on the Ru(001) and Pt(110)-(1x2) surfaces

The initial probabilities of activated, dissociative chemisorption of methane and ethane on Pt(110)-(1 x 2) have been measured. The surface temperature was varied from 450 to 900 K with the reactant gas temperature constant at 300 K. Under these conditions, we probe the kinetics of dissociation via trapping-mediated (as opposed to 'direct') mechanism. It was found that the probabilities of dissociation of both methane and ethane were strong functions of the surface temperature with an apparent activation energies of 14.4 kcal/mol for methane and 2.8 kcal/mol for ethane, which implys that the methane and ethane molecules have fully accommodated to the surface temperature. Kinetic isotope effects were observed for both reactions, indicating that the C-H bond cleavage was involved in the rate-limiting step. A mechanistic model based on the trapping-mediated mechanism is used to explain the observed kinetic behavior. The activation energies for C-H bond dissociation of the thermally accommodated methane and ethane on the surface extracted from the model are 18.4 and 10.3 kcal/mol, respectively.

The studies of the catalytic decomposition of formic acid on the Ru(001) surface with thermal desorption mass spectrometry following the adsorption of DCOOH and HCOOH on the surface at 130 and 310 K are described. Formic acid (DCOOH) chemisorbs dissociatively on the surface via both the cleavage of its O-H bond to form a formate and a hydrogen adatom, and the cleavage of its C-O bond to form a carbon monoxide, a deuterium adatom and an hydroxyl (OH). The former is the predominant reaction. The rate of desorption of carbon dioxide is a direct measure of the kinetics of decomposition of the surface formate. It is characterized by a kinetic isotope effect, an increasingly narrow FWHM, and an upward shift in peak temperature with Ɵ_T, the coverage of the dissociatively adsorbed formic acid. The FWHM and the peak temperature change from 18 K and 326 K at Ɵ_T = 0.04 to 8 K and 395 K at Ɵ_T = 0.89. The increase in the apparent activation energy of the C-D bond cleavage is largely a result of self-poisoning by the formate, the presence of which on the surface alters the electronic properties of the surface such that the activation energy of the decomposition of formate is increased. The variation of the activation energy for carbon dioxide formation with Ɵ_T accounts for the observed sharp carbon dioxide peak. The coverage of surface formate can be adjusted over a relatively wide range so that the activation energy for C-D bond cleavage in the case of DCOOH can be adjusted to be below, approximately equal to, or well above the activation energy for the recombinative desorption of the deuterium adatoms. Accordingly, the desorption of deuterium was observed to be governed completely by the desorption kinetics of the deuterium adatoms at low Ɵ_T, jointly by the kinetics of deuterium desorption and C-D bond cleavage at intermediate Ɵ_T, and solely by the kinetics of C-D bond cleavage at high Ɵ_T. The overall branching ratio of the formate to carbon dioxide and carbon monoxide is approximately unity, regardless the initial coverage Ɵ_T, even though the activation energy for the production of carbon dioxide varies with Ɵ_T. The desorption of water, which implies C-O bond cleavage of the formate, appears at approximately the same temperature as that of carbon dioxide. These observations suggest that the cleavage of the C-D bond and that of the C-O bond of two surface formates are coupled, possibly via the formation of a short-lived surface complex that is the precursor to to the decomposition.

The measurement of steady-state rate is demonstrated here to be valuable in determining kinetics associated with short-lived, molecularly adsorbed precursor to further reactions on the surface, by determining the kinetic parameters of the molecular precursor of formaldehyde to its dissociation on the Pt(110)-(1 x 2) surface.

Overlayers of nitrogen adatoms on Ru(001) have been characterized both by thermal desorption mass spectrometry and low-energy electron diffraction, as well as chemically via the postadsorption and desorption of ammonia and carbon monoxide.

The nitrogen-adatom overlayer was prepared by decomposing ammonia thermally on the surface at a pressure of 2.8 x 10^(-6) Torr and a temperature of 480 K. The saturated overlayer prepared under these conditions has associated with it a (√247/10 x √247/10)R22.7° LEED pattern, has two peaks in its thermal desorption spectrum, and has a fractional surface coverage of 0.40. Annealing the overlayer to approximately 535 K results in a rather sharp (√3 x √3)R30° LEED pattern with an associated fractional surface coverage of one-third. Annealing the overlayer further to 620 K results in the disappearance of the low-temperature thermal desorption peak and the appearance of a rather fuzzy p(2x2) LEED pattern with an associated fractional surface coverage of approximately one-fourth. In the low coverage limit, the presence of the (√3 x √3)R30° N overlayer alters the surface in such a way that the binding energy of ammonia is increased by 20% relative to the clean surface, whereas that of carbon monoxide is reduced by 15%.

A general methodology for the indirect relative determination of the absolute fractional surface coverages has been developed and was utilized to determine the saturation fractional coverage of hydrogen on Ru(001). Formaldehyde was employed as a bridge to lead us from the known reference point of the saturation fractional coverage of carbon monoxide to unknown reference point of the fractional coverage of hydrogen on Ru(001), which is then used to determine accurately the saturation fractional coverage of hydrogen. We find that ƟSAT/H = 1.02 (±0.05), i.e., the surface stoichiometry is Ru : H = 1 : 1. The relative nature of the method, which cancels systematic errors, together with the utilization of a glass envelope around the mass spectrometer, which reduces spurious contributions in the thermal desorption spectra, results in high accuracy in the determination of absolute fractional coverages.

Development of Asymmetric Protonation Reactions for the Synthesis of Indoline Alkaloids

Nitrogen-containing heterocycles, such as indolines and pyrroloindolines, are prevalent in a variety of diverse natural products, many of which exhibit remarkable biological activities. These frameworks have inspired innovative research aimed at discovering novel methods for their stereoselective preparation.

We have developed an enantioselective synthesis of pyrroloindolines based on a formal (3 + 2) cycloaddition of indoles and 2-amidoacrylates. This reaction is promoted by (R)-BINOL•SnCl₄; this complex is a Lewis acid-assisted Brønsted acid that effects a highly face-selective catalyst-controlled protonation of an enolate. Mechanistic studies also determined that the initial product of this reaction is an indolinium ion, which upon aqueous workup undergoes cyclization to the pyrroloindoline.

Based on this result, we investigated alternative nucleophiles to trap the indolinium ion. First, addition of sodium borohydride to the optimized reaction conditions yields indoline-containing amino acid derivatives.

Next, carbon nucleophiles were explored. Indole substrates incorporating a tethered alkene were exposed to the conditions for the formal (3 + 2) cycloaddition, resulting in a conjugate addition/asymmetric protonation/Prins cyclization cascade. In this transformation, the indolinium ion is attacked by the olefin, and the resulting carbocation is quenched by a chloride ion. Zirconium tetrachloride was found to be the optimal Lewis acid. Stoichiometric proton and chloride sources were also found to be crucial for reactivity.

Study of Solid-Phase Reactions of metals on GaAs

To make stable and reproducible contacts to GaAs, metals which react with GaAs in the solid-phase should be favored. In this study, contacts formed employing Pd/TiN/Pd/Ag, Pd:Mg/TiN/Pd:Mg/Ag and Ru/TiN/Ru/Ag are studied. The TiN layer is included to investigate its application as diffusion barrier in these metallizations. Contacts to n-GaAs are rectifying and the value of barrier height is modified upon annealing. Contacts to p-GaAs are initially rectifying but exhibit ohmic behaviour after annealing. The modifications in the electrical properties are attributed to the solid-phase reaction of metal and GaAs. The integrity of the contacts relies critically on the success of TiN to prevent the intermixing of Ag overlayer and the underlying layers. At elevated annealing temperatures (450°C), TiN fails to function as a diffusion barrier. As a result, the properties of the contact deteriorates.

Engineering and Characterization of Cytochrome P450 Enzymes for Nitrogen-Atom Transfer Reactions

The creation of novel enzyme activity is a great challenge to protein engineers, but nature has done so repeatedly throughout the process of natural selection. I begin by outlining the multitude of distinct reactions catalyzed by a single enzyme class, cytochrome P450 monooxygenases. I discuss the ability of cytochrome P450 to generate reactive intermediates capable of diverse reactivity, suggesting this enzyme can also be used to generate novel reactive intermediates in the form of metal-carbenoid and nitrenoid species. I then show that cytochrome P450 from Bacillus megaterium (P450_BM3) and its isolated cofactor can catalyze metal-nitrenoid transfer in the form of intramolecular C–H bond amination. Mutations to the protein sequence can enhance the reactivity and selectivity of this transformation significantly beyond that of the free cofactor. Next, I demonstrate an intermolecular nitrene transfer reaction catalyzed by P450_BM3 in the form of sulfide imidation. Understanding that sulfur heteroatoms are strong nucleophiles, I show that increasing the sulfide nucleophilicity through substituents on the aryl sulfide ring can dramatically increase reaction productivity. To explore engineering nitrenoid transfer in P450BM3, active site mutagenesis is employed to tune the regioselectivity intramolecular C–H amination catalysts. The solution of the crystal structure of a highly selective variant demonstrates that hydrophobic residues in the active site strongly modulate reactivity and regioselectivity. Finally, I use a similar strategy to develop P450-based catalysts for intermolecular olefin aziridination, demonstrating that active site mutagenesis can greatly enhance this nitrene transfer reaction. The resulting variant can catalyze intermolecular aziridination with more than 1000 total turnovers and enantioselectivity of up to 99% ee.

Investigations of group IVA transition metal mediated carbon-carbon bond forming reactions

Zirconocene aldehyde and ketone complexes were synthesized in high yield by treatment of zirconocene acyl complexes with trimethylaluminum or diisobutylaluminum hydride. These complexes, which are activated by dialkylaluminum chloride ligands, inserted unsaturated substrates such as alkynes, allenes, ethylene, nitriles, ketenes, aldehydes, ketones, lactones, and acid chlorides with moderate to high conversion. Insertion of aldehyde substrates yielded zirconocene diolate complexes with up to 20:1 (anti:syn) diastereoselectivity. The zirconocene diolates were hydrolyzed to afford unsymmetrical 1,2-diols in 40-80% isolated yield. Unsymmetrical ketones gave similar insertion yields with little or no diastereoselectivity. A high yielding one-pot method was developed that coupled carbonyl substrates with zirconocene aldehyde complexes that were derived from olefins by hydrozirconation and carbonylation. The zirconocene aldehyde complexes also inserted carbon monoxide and gave acyloins in 50% yield after hydrolysis.

The insertion reaction of aryl epoxides with the trimethylphoshine adduct of titanocene methylidene was examined. The resulting oxytitanacyclopentanes were carbonylated and oxidatively cleaved with dioxygen to afford y-lactones in moderate yields. Due to the instability and difficult isolation of titanocene methylidene trimethylphoshine adducts, a one-pot method involving the addition of catalytic amounts of trimethylphosphine to β,β-dimethyltitanacyclobutane was developed. A series of disubstituted aryl epoxides were examined which gave mixtures of diastereomeric insertion products. Based on these results, as well as earlier Hammett studies and labeling experiments, a biradical transition state intermediate is proposed. The method is limited to aryl substituted epoxide substrates with aliphatic examples showing no insertion reactivity.

The third study involved the use of magnesium chloride supported titanium catalysts for the Lewis acid catalyzed silyl group transfer condensation of enol silanes with aldehydes. The reaction resulted in silylated aldol products with as many as 140 catalytic turnovers before catalyst inactivation. Low diastereoselectivities favoring the anti-isomer were consistent with an open transition state involving a titanium atom bound to the catalyst surface. The catalysts were also used for the aldol group transfer polymerization of t-butyldimethylsilyloxy-1-ethene resulting in polymers with molecular weights of 5000-31,000 and molar mass dispersities of 1.5-2.8. Attempts to polymerize methylmethacrylate using GTP proved unsuccessful with these catalysts.

Neutron-neutron scattering parameters from the ^3He(d,t)2p and ^3H(d,^3He)2n reactions

The energy spectra of tritons and Helium-3 nuclei from the reactions ³He(d,t)2p, ³H(d,³He)2n, ³He(d,³He)pn, and ³H(d,t)pn were measured between 6° and 20° at a bombarding energy of 10.9 MeV. An upper limit of 5 μb/sr. was obtained for producing a bound di-neutron at 6° and 7.5°. The ³He(d,t)2p and ³H(d,³He)2n data, together with previous measurements at higher energies, have been used to investigate whether one can unambiguously extract information on the two-nucleon system from these three-body final state reactions. As an aid to these theoretical investigations, Born approximation calculations were made employing realistic nucleon-nucleon potentials and an antisymmetrized final state wave function for the five-particle system. These calculations reproduce many of the features observed in the experimental data and indicate that the role of exchange processes cannot be ignored. The results show that previous attempts to obtain information on the neutron-neutron scattering length from the ³H(d,³He)2n reaction may have seriously overestimated the precision that could be attained.

Excitation energies and decay properties of T = 3/2 states in ^(17)O, ^(17)F and ^(21)Na

The two lowest T = 3/2 levels in ²¹Na have been studied in the ¹⁹F(³He, n), ²⁰Ne (p,p) and ²⁰Ne (p,p’) reactions, and their excitation energies, spins, parities and widths have been determined. In a separate investigation, branching ratios were measured for the isospin-nonconserving particle decays of the lowest T = 3/2 levels in ¹⁷O and ¹⁷F to the ground state and first two excited states of ¹⁶O, by studying the ¹⁵N(³He,n) ¹⁷F*(p) ¹⁶O and ¹⁸O(³He, α)¹⁷O*(n) ¹⁶O reactions.

The ¹⁹F(³He,n) ²¹Na reaction was studied at incident energies between 4.2 and 5.9 MeV using a pulsed-beam neutron-time-of-flight spectrometer. Two T = 3/2 levels were identified at excitation energies of 8.99 ± 0.05 MeV (J > ½) and 9.22 ± 0.015 MeV (J ^π = ½⁺, Γ ˂ 40 keV). The spins and parities were determined by a comparison of the measured angular distributions with the results of DWBA calculations.

These two levels were also obsesrved as isospin-forbidden resonances in the ²⁰Ne(p,p) and ²⁰Ne(p,p’) reactions. Excitation energies were measured and spins, parities, and widths were determined from a single level dispersion theory analysis. The following results were obtained:

E_x = 8.973 ± 0.007 MeV, J ^π = 5/2 ⁺ or 3/2⁺, Γ ≤ 1.2 keV,

^Γp_o = 0.1 ± 0.05 keV; E_x = 9.217 ± 0.007 MeV, J^π = ½ ⁺,

Γ = 2.3 ± 0.5 keV, ^Γp_o = 1.1 ± 0.3 keV.

Isospin assignments were made on the basis of excitation energies, spins, parities, and widths.

Branching ratios for the isospin-nonconserving proton decays of the 11.20 MeV, T = 3/2 level in ¹⁷F were measured by the ¹⁵N(³He,n) ¹⁷ F*(p) ¹⁶O reaction to be 0.088 ± 0.016 to the ground state of ¹⁶O and 0.22 ± 0.04 to the unresolved 6.05 and 6.13 MeV levels of ¹⁶O. Branching ratios for the neutron decays of the analogous T = 3/2 level, at 11.08 MeV in ¹⁷O, were measured by the ¹⁶O(³He, α)¹⁷O*(n)¹⁶O reaction to be 0.91 ± 0.15 to the ground state of ¹⁶O and 0.05 ± 0.02 to the unresolved 6.05 and 6.13 MeV states. By comparing the ratios of reduced widths for the mirror decays, the form of the isospin impurity in the T = 3/2 levels is shown to depend on T_z.

Studies of response to earthquake ground motion

A study is made of the accuracy of electronic digital computer calculations of ground displacement and response spectra from strong-motion earthquake accelerograms. This involves an investigation of methods of the preparatory reduction of accelerograms into a form useful for the digital computation and of the accuracy of subsequent digital calculations. Various checks are made for both the ground displacement and response spectra results, and it is concluded that the main errors are those involved in digitizing the original record. Differences resulting from various investigators digitizing the same experimental record may become as large as 100% of the maximum computed ground displacements. The spread of the results of ground displacement calculations is greater than that of the response spectra calculations. Standardized methods of adjustment and calculation are recommended, to minimize such errors.

Studies are made of the spread of response spectral values about their mean. The distribution is investigated experimentally by Monte Carlo techniques using an electric analog system with white noise excitation, and histograms are presented indicating the dependence of the distribution on the damping and period of the structure. Approximate distributions are obtained analytically by confirming and extending existing results with accurate digital computer calculations. A comparison of the experimental and analytical approaches indicates good agreement for low damping values where the approximations are valid. A family of distribution curves to be used in conjunction with existing average spectra is presented. The combination of analog and digital computations used with Monte Carlo techniques is a promising approach to the statistical problems of earthquake engineering.

Methods of analysis of very small earthquake ground motion records obtained simultaneously at different sites are discussed. The advantages of Fourier spectrum analysis for certain types of studies and methods of calculation of Fourier spectra are presented. The digitizing and analysis of several earthquake records is described and checks are made of the dependence of results on digitizing procedure, earthquake duration and integration step length. Possible dangers of a direct ratio comparison of Fourier spectra curves are pointed out and the necessity for some type of smoothing procedure before comparison is established. A standard method of analysis for the study of comparative ground motion at different sites is recommended.

A study of T=1, T=3/2, and T=2 states in some light nuclei using (He^3,n) reactions

The (He³, n) reactions on B¹¹, N¹⁵, O¹⁶, and O¹⁸ targets have been studied using a pulsed-beam time-of-flight spectrometer. Special emphasis was placed upon the determination of the excitation energies and properties of states with T = 1 (in Ne¹⁸), T = 3/2 (in N¹³ and F¹⁷) and T = 2 (in Ne²⁰). The identification of the T = 3/2 and T = 2 levels is based on the structure of these states as revealed by intensities and shapes of angular distributions. The reactions are interpreted in terms of double stripping theory. Angular distributions have been compared with plane and distorted wave stripping theories. Results for the four reactions are summarized below:

1) O¹⁶ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV and two previously unreported levels in Ne¹⁸ were observed at E_x = 4.55 ± .015 MeV (Γ = 70 ± 30 keV) and E_x = 5.14 ± .018 MeV (Γ = 100 ± 40 keV).

2) B¹¹ (He³, n). The reaction has been studied at incident energies up to 13.5 MeV. Three T = 3/2 levels in N¹³ have been identified at E_x = 15.068 ± .008 MeV (Γ ˂ 15 keV), E_x = 18.44 ± .04, and E_x 18.98 ± .02 MeV (Γ = 40 ± 20 keV).

3) N¹⁵ (He³, n). The reaction has been studied at incident energies up to 11.88 MeV. T = 3/2 levels in F¹⁷ have been identified at E_x = 11.195 ± .007 MeV (Γ ˂ 20 keV), E_x = 12.540 ± .010 MeV (Γ ˂ 25 keV), and E_x = 13.095 ± .009 MeV (Γ ˂ 25 keV).

4) O¹⁸ (He³, n). The reaction has been studied at incident energies up to 9.0 MeV. The excitation energy of the lowest T = 2 level in Ne²⁰ has been found to be 16.730 ± .006 MeV (Γ ˂ 20 keV).

Angular distributions of the transitions leading to the above higher isospin states are well described by double stripping theory. Analog correspondences are established by comparing the present results with recent studies (t, p) and (He³, p) reactions on the same targets.

Reggeization of some unequal mass processes involving higher spins : the reactions pion-nucleon going to (V,Nucleon), pion-nucleon going to (V,delta-hyperon), and photon-proton going to (positive pion,neutron)

Recent theoretical developments in the reggeization of inelastic processes involving particles with high spin are incorporated into a model of vector meson production. A number of features of experimental differential cross sections and density matrices are interpreted in terms of this model.

The method chosen for reggeization of helicity amplitudes first separates kinematic zeros and singularities from the parity-conserving amplitudes and then applies results of Freedman and Wang on daughter trajectories to the remaining factors. Kinematic constraints on helicity amplitudes at t = 0 and t = (M – M_Δ)² are also considered.

It is found that data for reactions of types πN→VN and πN→VΔ are consistent with a model of this type in which all kinematic constraints at t = 0 are satisfied by evasion (vanishing of residue functions). As a quantitative test of the parametrization, experimental differential cross sections of vector meson production reactions dominated by pion trajectory exchange are compared with the theory. It is found that reduced residue functions are approximately constant, once the kinematic behavior near t = (M – M_Δ)² has been removed.

The alternative possibility of conspiracy between amplitudes is also discussed; and it is shown that unless conspiracy is present, some amplitudes allowed by angular momentum conservation will not contribute with full strength in the forward direction. An example, γp→π⁺n in which the data for dσ/dt indicate conspiracy, is studied in detail.

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.

Investigation of strong earthquake ground motion

The pattern of energy release during the Imperial Valley, California, earthquake of 1940 is studied by analysing the El Centro strong motion seismograph record and records from the Tinemaha seismograph station, 546 km from the epicenter. The earthquake was a multiple event sequence with at least 4 events recorded at El Centro in the first 25 seconds, followed by 9 events recorded in the next 5 minutes. Clear P, S and surface waves were observed on the strong motion record. Although the main part of the earthquake energy was released during the first 15 seconds, some of the later events were as large as M = 5.8 and thus are important for earthquake engineering studies. The moment calculated using Fourier analysis of surface waves agrees with the moment estimated from field measurements of fault offset after the earthquake. The earthquake engineering significance of the complex pattern of energy release is discussed. It is concluded that a cumulative increase in amplitudes of building vibration resulting from the present sequence of shocks would be significant only for structures with relatively long natural period of vibration. However, progressive weakening effects may also lead to greater damage for multiple event earthquakes.

The model with surface Love waves propagating through a single layer as a surface wave guide is studied. It is expected that the derived properties for this simple model illustrate well several phenomena associated with strong earthquake ground motion. First, it is shown that a surface layer, or several layers, will cause the main part of the high frequency energy, radiated from the nearby earthquake, to be confined to the layer as a wave guide. The existence of the surface layer will thus increase the rate of the energy transfer into the man-made structures on or near the surface of the layer. Secondly, the surface amplitude of the guided SH waves will decrease if the energy of the wave is essentially confined to the layer and if the wave propagates towards an increasing layer thickness. It is also shown that the constructive interference of SH waves will cause the zeroes and the peaks in the Fourier amplitude spectrum of the surface ground motion to be continuously displaced towards the longer periods as the distance from the source of the energy release increases.