I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

Structure of the earth's core and lowermost mantle from seismic PKP waves

This thesis addresses the fine structure, both radial and lateral, of compressional wave velocity and attenuation of the Earth's core and the lowermost mantle using waveforms, differential travel times and amplitudes of PKP waves, which penetrate the Earth's core.

The structure near the inner core boundary (ICB) is studied by analyzing waveforms of a regional sample. The waveform modeling approach is demonstrated to be an effective tool for constrainning the ICB structure. The best model features a sharp velocity jump of 0.78km/s at the ICB and a low velocity gradient at the lowermost outer core (indicating possible inhomogeneity) and high attenuation at the top of the inner core.

A spherically symmetric P-wave model of the core, is proposed from PKP differential times, waveforms and amplitudes. The ICB remains sharp with a velocity jump of 0. 78km/ s. A very low velocity gradient at the base of the fluid core is demonstrated to be a robust feature, indicating inhomogeneity is practically inevitable. The model also indicates that the attenuation in the inner core decreases with depth. The velocity at D" is smaller than PREM.

The inner core is confirmed to be very anisotropic, possessing a cylindrical symmetry around the Earth spin axis with the N-S direction 3% faster than the E-W direction. All of the N-S rays through the inner core were found to be faster than the E-W rays by 1.5 to 3.5s. Exhaustive data selection and efforts in insolating contributions from the region above ensure that this is an inner core feature.

The anisotropy at the very top of the inner core is found to be distinctly different from the deeper part. The top 60km of the inner core is not anisotropic. From 60km to 150km, there appears to be a transition from isotropy to anisotropy.

PKP differential travel times are used to study the P velocity structure in D". Systematic regional variations of up to 2s in AB-DF times were observed, attributed primarily to heterogeneities in the lower 500km of the mantle. However, direct comparisons with tomographic models are not successful.

Electron transfer in cytochrome c oxidase : the rate of electron equillibration between cytochrome a and copper A

Intramolecular electron transfer in partially reduced cytochrome c oxidase has been studied by means of perturbed equilibrium techniques. We have prepared a three electron reduced, CO inhibited form of the enzyme in which cytochrome a and copper A are partially reduced an in intramolecular redox equilibrium. When these samples were photolyzed using a nitrogen laser (0.6 µs, 1.0 mJ pulses) changes in absorbance at 598 nm and 830 nm were observed which are consistent with a fast electron from cytochrome a to copper A. The absorbance changes at 598 nm have an apparent rate of 17,200 ± 1,700 s^(-1) (1σ), at pH 7.0 and 25.5 °C. These changes were not observed in either the CO mixed valence or CO inhibited fully reduced forms of the enzyme. The rate is fastest at about pH 8.0, and falls off in either direction, and there is a small, but clear temperature dependence. The process was also observed in the cytochrome c -- cytochrome c oxidase high affinity complex.

This rate is far faster than any rate measured or inferred previously for the cytochrome a -- copper A electron equilibration, but the interpretation of these results is hampered by the fact that the relaxation could only be followed during the time before CO became rebound to the oxygen binding site. The meaning of our our measured rate is discussed, along with other reported rates for this process. In addition, a temperature-jump experiment on the same system is discussed.

We have also prepared a partially reduced, cyanide inhibited form of the enzyme in which cytochrome a, copper A and copper B are partially reduced and in redox equilibrium. Warming these samples produced absorbance changes at 605 nm which indicate that cytochrome a was becoming more oxidized, but there were no parallel changes in absorbance at 830 nm as would be expected if copper A was becoming reduced. We concluded that electrons were being redistributed from cytochrome a to copper B. The kinetics of the absorbance changes at 605 nm were investigated by temperature-jump methods. Although a rate could not be resolved, we concluded that the process must occur with an (apparent) rate larger than 10,000 s^(-1).

During the course of the temperature-jump experiments, we also found that non-redox related, temperature dependent absorbance changes in fully reduced CO inhibited cytochrome c oxidase, and in the cyanide mixed valence enzyme, took place with an (apparent) rate faster that 30,000 s^(-1).

Topics in randomized numerical linear algebra

This thesis studies three classes of randomized numerical linear algebra algorithms, namely: (i) randomized matrix sparsification algorithms, (ii) low-rank approximation algorithms that use randomized unitary transformations, and (iii) low-rank approximation algorithms for positive-semidefinite (PSD) matrices.

Randomized matrix sparsification algorithms set randomly chosen entries of the input matrix to zero. When the approximant is substituted for the original matrix in computations, its sparsity allows one to employ faster sparsity-exploiting algorithms. This thesis contributes bounds on the approximation error of nonuniform randomized sparsification schemes, measured in the spectral norm and two NP-hard norms that are of interest in computational graph theory and subset selection applications.

Low-rank approximations based on randomized unitary transformations have several desirable properties: they have low communication costs, are amenable to parallel implementation, and exploit the existence of fast transform algorithms. This thesis investigates the tradeoff between the accuracy and cost of generating such approximations. State-of-the-art spectral and Frobenius-norm error bounds are provided.

The last class of algorithms considered are SPSD "sketching" algorithms. Such sketches can be computed faster than approximations based on projecting onto mixtures of the columns of the matrix. The performance of several such sketching schemes is empirically evaluated using a suite of canonical matrices drawn from machine learning and data analysis applications, and a framework is developed for establishing theoretical error bounds.

In addition to studying these algorithms, this thesis extends the Matrix Laplace Transform framework to derive Chernoff and Bernstein inequalities that apply to all the eigenvalues of certain classes of random matrices. These inequalities are used to investigate the behavior of the singular values of a matrix under random sampling, and to derive convergence rates for each individual eigenvalue of a sample covariance matrix.

The ^(26)Al(p,γ)^(27)Si reaction : stellar origins of galactic ^(26)Al

To explain the ^(26)Mg isotopic anomaly seen in meteorites (^(26)Al daughter) as well as the observation of 1809-keV γ rays in the interstellar medium (live decay of 26Al) one must know, among other things, the destruction rate of ^(26)Al. Properties of states in ^(27)Si just above the ^(26)Al + p mass were investigated to determine the destruction rate of ^(26)Al via the ^(26)Al(p,γ)^(27)Si reaction at astrophysical temperatures.

Twenty micrograms of ^(26)Al were used to produce two types of Al_2O_3 targets by evaporation of the oxide. One was onto a thick platinum backing suitable for (p,γ) work, and the other onto a thin carbon foil for the (^3He,d) reaction.

The ^(26)Al(p,γ)^(27)Si excitation function, obtained using a germanium detector and voltage-ramped target, confirmed known resonances and revealed new ones at 770, 847, 876, 917, and 928 keV. Possible resonances below the lowest observed one at E_p = 286 keV were investigated using the ^(26)Al(^3He,d)^(27)Si proton-transfer reaction. States in 27Si corresponding to 196- and 286-keV proton resonances were observed. A possible resonance at 130 keV (postulated in prior work) was shown to have a strength of wγ less than 0.02 µeV.

By arranging four large Nal detector as a 47π calorimeter, the 196-keV proton resonance, and one at 247 keV, were observed directly, having wγ = 55± 9 and 10 ± 5 µeV, respectively.

Large uncertainties in the reaction rate have been reduced. At novae temperatures, the rate is about 100 times faster than that used in recent model calculations, casting some doubt on novae production of galactic ^(26)Al.

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

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

Methods of multiparameter inversion of seismic data using the acoustic and elastic born approximations

This thesis presents two different forms of the Born approximations for acoustic and elastic wavefields and discusses their application to the inversion of seismic data. The Born approximation is valid for small amplitude heterogeneities superimposed over a slowly varying background. The first method is related to frequency-wavenumber migration methods. It is shown to properly recover two independent acoustic parameters within the bandpass of the source time function of the experiment for contrasts of about 5 percent from data generated using an exact theory for flat interfaces. The independent determination of two parameters is shown to depend on the angle coverage of the medium. For surface data, the impedance profile is well recovered.

The second method explored is mathematically similar to iterative tomographic methods recently introduced in the geophysical literature. Its basis is an integral relation between the scattered wavefield and the medium parameters obtained after applying a far-field approximation to the first-order Born approximation. The Davidon-Fletcher-Powell algorithm is used since it converges faster than the steepest descent method. It consists essentially of successive backprojections of the recorded wavefield, with angular and propagation weighing coefficients for density and bulk modulus. After each backprojection, the forward problem is computed and the residual evaluated. Each backprojection is similar to a before-stack Kirchhoff migration and is therefore readily applicable to seismic data. Several examples of reconstruction for simple point scatterer models are performed. Recovery of the amplitudes of the anomalies are improved with successive iterations. Iterations also improve the sharpness of the images.

The elastic Born approximation, with the addition of a far-field approximation is shown to correspond physically to a sum of WKBJ-asymptotic scattered rays. Four types of scattered rays enter in the sum, corresponding to P-P, P-S, S-P and S-S pairs of incident-scattered rays. Incident rays propagate in the background medium, interacting only once with the scatterers. Scattered rays propagate as if in the background medium, with no interaction with the scatterers. An example of P-wave impedance inversion is performed on a VSP data set consisting of three offsets recorded in two wells.

Synthesis and characterization of 1,1-di-tert-butyldiazene

The synthesis and direct observation of 1,1-di-tert-butyldiazene (16) at -127°C is described. The absorption spectrum of a red solution of 1,1-diazene 16 reveals a structured absorption band with λ max at 506 run (Me_2O, -125°C). The vibrational spacing in S_1 is about 1200 cm^(-1). The excited state of 16 emits weakly with a single maximum at 715 run observed in the fluorescence spectrum (Me_2O:CD_2Cl_2, -196°C). The proton NMR spectrum of 16 occurs as a singlet at 1.41 ppm. Monitoring this NMR absorption at -94^0 ± 2°C shows that 1,1-diazene 16 decomposes with a first-order rate of 1.8 x 10^(-3) sec(-1) to form isobutane, isobutylene and hexarnethylethane. This rate is 10^8 and 10^(34) times faster than the thermal decomposition of the corresponding cis and trans 1,2-di-tert-butyldiazene isomers. The free energy of activation for decomposition of 1,1-diazene 16 is found to be 12.5 ± 0.2 kcal/mol at -94°C which is much lower than the values of 19.1 and 19.4 kcal/lmole calculated at -94°C for N-(2,2,6,6- tetramethylpiperidyl)nitrene (3) and N-(2,2,5,5- tetrarnethylpyrrolidyl)nitrene (4), respectively. This difference between 16 and the cyclic-1,1-diazenes 3 and 4 can be attributed to a large steric interaction between the tert-butyl groups in 1,1-diazene 16.

In order to investigate the nature of the singlet-triplet gap in 1,1-diazenes, 2,5-di-tert-butyl-N-pyrrolynitrene (22) was generated but was found to be too reactive towards dimerization to be persistent. In the presence of dimethylsulfoxide, however, N-pyrrolynitrene (22) can be trapped as N-(2,5-di-tert-butyl- N'-pyrrolyl)dimethylsulfoxirnine (38). N-(2,5-di-tert-butyl-N'-pyrrolyl)dimethylsulfoximine (38-d^6) exchanges with free dimethylsulfoxide at 50°C in solution, presumably by generation and retrapping of pyrrolynitrene 22.

Model studies of Ziegler-Natta olefin polymerization using group 3 and group 4 metallocenes

Evidence for the stereochemical isomerization of a variety of ansa metallocene compounds is presented. For the scandocene allyl derivatives described here, we have established that the process is promoted by a variety of salts in both ether and hydrocarbon solvents and is not accelerated by light. A plausible mechanism based on an earlier proposal by Marks, et al., is offered as an explanation of this process. It involves coordination of anions and/or donor solvents to the metal center with cation assistance to encourage metalcyclopentadienyl bond heterolysis, rotation about the Si-Cp bond of the detached cyclopentadienide and recoordination of the opposite face. Our observations in some cases of thermodynamic racemic:meso ratios under the reaction conditions commonly used for the synthesis of the metallocene chlorides suggests that the interchange is faster than metallation, such that the composition of the reaction mixture is determined by thermodynamic, not kinetic, control in these cases.

Two new ansa-scandocene alkenyl compounds react with olefins resulting in the formation of η³-allyl complexes. Kinetics and labeling experiments indicate a tuck-in intermediate on the reaction pathway; in this intermediate the metal is bound to the carbon adjacent to the silyllinker in the rear of the metallocene wedge. In contrast, reaction of permethylscandocene alkenyl compounds with olefins results, almost exclusively, in vinylic C-H bond activation. It is proposed that relieving transition state steric interactions between the cyclopentadienyl rings and the olefin by either linking the rings together or using a larger lanthanide metal may allow for olefin coordination, stabilizing the transition state for allylic σ-bond metathesis.

A selectively isotopically labeled propylene, CH₂CD(¹³CH₃), was synthesized and its polymerization was carried out at low concentration in toluene solution using isospecific metallocene catalysts. Analysis of the NMR spectra (¹³C, ¹H, and ²H) of the resultant polymers revealed that the production of stereoerrors through chain epimerization proceeds exclusively by the tertiaryalkyl mechanism. Additionally, enantiofacial inversion of the terminally unsaturated polymer chain occurs by a non-dissociative process. The implications of these results on the mechanism of olefin polymerization with these catalysts is discussed.

Random propagation in complex systems : nonlinear matrix recursions and epidemic spread

This dissertation studies long-term behavior of random Riccati recursions and mathematical epidemic model. Riccati recursions are derived from Kalman filtering. The error covariance matrix of Kalman filtering satisfies Riccati recursions. Convergence condition of time-invariant Riccati recursions are well-studied by researchers. We focus on time-varying case, and assume that regressor matrix is random and identical and independently distributed according to given distribution whose probability distribution function is continuous, supported on whole space, and decaying faster than any polynomial. We study the geometric convergence of the probability distribution. We also study the global dynamics of the epidemic spread over complex networks for various models. For instance, in the discrete-time Markov chain model, each node is either healthy or infected at any given time. In this setting, the number of the state increases exponentially as the size of the network increases. The Markov chain has a unique stationary distribution where all the nodes are healthy with probability 1. Since the probability distribution of Markov chain defined on finite state converges to the stationary distribution, this Markov chain model concludes that epidemic disease dies out after long enough time. To analyze the Markov chain model, we study nonlinear epidemic model whose state at any given time is the vector obtained from the marginal probability of infection of each node in the network at that time. Convergence to the origin in the epidemic map implies the extinction of epidemics. The nonlinear model is upper-bounded by linearizing the model at the origin. As a result, the origin is the globally stable unique fixed point of the nonlinear model if the linear upper bound is stable. The nonlinear model has a second fixed point when the linear upper bound is unstable. We work on stability analysis of the second fixed point for both discrete-time and continuous-time models. Returning back to the Markov chain model, we claim that the stability of linear upper bound for nonlinear model is strongly related with the extinction time of the Markov chain. We show that stable linear upper bound is sufficient condition of fast extinction and the probability of survival is bounded by nonlinear epidemic map.

Investigation of bulk indium antimonide as heterodyne detector for the submillimeter wavelength region

Bulk n-lnSb is investigated at a heterodyne detector for the submillimeter wavelength region. Two modes or operation are investigated: (1) the Rollin or hot electron bolometer mode (zero magnetic field), and (2) the Putley mode (quantizing magnetic field). The highlight of the thesis work is the pioneering demonstration or the Putley mode mixer at several frequencies. For example, a double-sideband system noise temperature of about 510K was obtained using a 812 GHz methanol laser for the local oscillator. This performance is at least a factor or 10 more sensitive than any other performance reported to date at the same frequency. In addition, the Putley mode mixer achieved system noise temperatures of 250K at 492 GHz and 350K at 625 GHz. The 492 GHz performance is about 50% better and the 625 GHz is about 100% better than previous best performances established by the Rollin-mode mixer. To achieve these results, it was necessary to design a totally new ultra-low noise, room-temperature preamp to handle the higher source impedance imposed by the Putley mode operation. This preamp has considerably less input capacitance than comparably noisy, ambient designs.

In addition to advancing receiver technology, this thesis also presents several novel results regarding the physics of n-lnSb at low temperatures. A Fourier transform spectrometer was constructed and used to measure the submillimeter wave absorption coefficient of relatively pure material at liquid helium temperatures and in zero magnetic field. Below 4.2K, the absorption coefficient was found to decrease with frequency much faster than predicted by Drudian theory. Much better agreement with experiment was obtained using a quantum theory based on inverse-Bremmstrahlung in a solid. Also the noise of the Rollin-mode detector at 4.2K was accurately measured and compared with theory. The power spectrum is found to be well fit by a recent theory of non- equilibrium noise due to Mather. Surprisingly, when biased for optimum detector performance, high purity lnSb cooled to liquid helium temperatures generates less noise than that predicted by simple non-equilibrium Johnson noise theory alone. This explains in part the excellent performance of the Rollin-mode detector in the millimeter wavelength region.

Again using the Fourier transform spectrometer, spectra are obtained of the responsivity and direct detection NEP as a function of magnetic field in the range 20-110 cm^-1. The results show a discernable peak in the detector response at the conduction electron cyclotron resonance frequency tor magnetic fields as low as 3 KG at bath temperatures of 2.0K. The spectra also display the well-known peak due to the cyclotron resonance of electrons bound to impurity states. The magnitude of responsivity at both peaks is roughly constant with magnet1c field and is comparable to the low frequency Rollin-mode response. The NEP at the peaks is found to be much better than previous values at the same frequency and comparable to the best long wavelength results previously reported. For example, a value NEP=4.5x10^-13/Hz^1/2 is measured at 4.2K, 6 KG and 40 cm^-1. Study of the responsivity under conditions of impact ionization showed a dramatic disappearance of the impurity electron resonance while the conduction electron resonance remained constant. This observation offers the first concrete evidence that the mobility of an electron in the N=0 and N=1 Landau levels is different. Finally, these direct detection experiments indicate that the excellent heterodyne performance achieved at 812 GHz should be attainable up to frequencies of at least 1200 GHz.

Investigations of the conductor-semiconductor interface

I. Schottky barriers produced by polymeric sulfur nitride, (SN)_x, on nine common III-V and II-VI compound semiconductors are compared to barriers formed by Au. The conductor (SN)_x produces significantly higher barriers to n-type semiconductors and lower barriers to p-type semiconductors than Au, the most electronegative elemental metal. The barrier height improvement, defined as ɸ(SN)_x - ɸ(Au), is smaller on covalent semiconductors than on ionic semiconductors; (SN)x barriers follow the ionic-covalent transition. Details of (SN)_x film deposition, samples preparation, and barrier height measurements are described.

II. The rate of dissolution of amorphous Si into solid Al is measured. The rate of movement of the amorphous Si/Al interface is found to be much faster than predicted by a simple model of the transport of Si through Al. This result is related to defects in the growth of epitaxial Si using the solid phase epitaxy process.

Chemical studies of viral entry mechanisms: I. Hydrophobic protein-lipid interactions during sendai virus membrane fusion. II. Kinetics of bacteriophage λ DNA injection.

Viruses possess very specific methods of targeting and entering cells. These methods would be extremely useful if they could also be applied to drug delivery, but little is known about the molecular mechanisms of the viral entry process. In order to gain further insight into mechanisms of viral entry, chemical and spectroscopic studies in two systems were conducted, examining hydrophobic protein-lipid interactions during Sendai virus membrane fusion, and the kinetics of bacteriophage λ DNA injection.

Sendai virus glycoprotein interactions with target membranes during the early stages of fusion were examined using time-resolved hydrophobic photoaffinity labeling with the lipid-soluble carbene generator3-(trifluoromethyl)-3-(m-^(125 )I] iodophenyl)diazirine (TID). The probe was incorporated in target membranes prior to virus addition and photolysis. During Sendai virus fusion with liposomes composed of cardiolipin (CL) or phosphatidylserine (PS), the viral fusion (F) protein is preferentially labeled at early time points, supporting the hypothesis that hydrophobic interaction of the fusion peptide at the N-terminus of the F_1 subunit with the target membrane is an initiating event in fusion. Correlation of the hydrophobic interactions with independently monitored fusion kinetics further supports this conclusion. Separation of proteins after labeling shows that the F_1 subunit, containing the putative hydrophobic fusion sequence, is exclusively labeled, and that the F_2 subunit does not participate in fusion. Labeling shows temperature and pH dependence consistent with a need for protein conformational mobility and fusion at neutral pH. Higher amounts of labeling during fusion with CL vesicles than during virus-PS vesicle fusion reflects membrane packing regulation of peptide insertion into target membranes. Labeling of the viral hemagglutinin/neuraminidase (HN) at low pH indicates that HN-mediated fusion is triggered by hydrophobic interactions, after titration of acidic amino acids. HN labeling under nonfusogenic conditions reveals that viral binding may involve hydrophobic as well as electrostatic interactions. Controls for diffusional labeling exclude a major contribution from this source. Labeling during reconstituted Sendai virus envelope-liposome fusion shows that functional reconstitution involves protein retention of the ability to undergo hydrophobic interactions.

Examination of Sendai virus fusion with erythrocyte membranes indicates that hydrophobic interactions also trigger fusion between biological membranes, and that HN binding may involve hydrophobic interactions as well. Labeling of the erythrocyte membranes revealed close membrane association of spectrin, which may play a role in regulating membrane fusion. The data show that hydrophobic fusion protein interaction with both artificial and biological membranes is a triggering event in fusion. Correlation of these results with earlier studies of membrane hydration and fusion kinetics provides a more detailed view of the mechanism of fusion.

The kinetics of DNA injection by bacteriophage λ. into liposomes bearing reconstituted receptors were measured using fluorescence spectroscopy. LamB, the bacteriophage receptor, was extracted from bacteria and reconstituted into liposomes by detergent removal dialysis. The DNA binding fluorophore ethidium bromide was encapsulated in the liposomes during dialysis. Enhanced fluorescence of ethidium bromide upon binding to injected DNA was monitored, and showed that injection is a rapid, one-step process. The bimolecular rate law, determined by the method of initial rates, revealed that injection occurs several times faster than indicated by earlier studies employing indirect assays.

It is hoped that these studies will increase the understanding of the mechanisms of virus entry into cells, and to facilitate the development of virus-mimetic drug delivery strategies.

Electron transfer in chemically and genetically modified myoglobins

The temperature dependences of the reduction potentials (E^o') of wildtype human myoglobin (Mb) and three site-directed mutants have been measured by using thin-layer spectroelectrochemistry. Residue Val68, which is in van der Waals contact with the heme in Mb, has been replaced by Glu, Asp, and Asn. At pH 7.0, reduction of the heme iron (III) in the former two proteins is accompanied by uptake of a proton by the protein. The changes in E^o', and the standard entropy (ΔS^o') and enthalpy (ΔH^o') of reduction in the mutant proteins were determined relative to values for wild-type; the change in E^o' at 25°C was about -200 millivolts for the Glu and Asp mutants, and about -80 millivolts for the Asn mutant. Reduction of Fe(III) to Fe(II) in the Glu and Asp mutants is accompanied by uptake of a proton. These studies demonstrate that Mb can tolerate substitution of a buried hydrophobic group by potentially charged and polar residues, and that such amino acid replacements can lead to substantial changes in the redox thermodynamics of the protein.

Through analysis of the temperature dependence and shapes of NMR dispersion signals, it is determined that a water molecule is bound to the sixth coordination site of the ferric heme in the Val68Asp and in the Val68Asn recombinant proteins while the carboxyl group of the sidechain of Glu68 occupies this position in Val68Glu. The relative rhombic distortions in the ESR spectra of these mutant proteins combined with H₂¹⁷O and spin interconversion experiments performed on them confirm the conclusions of the NMRD study.

The rates of intramolecular electron transfer (ET) of (NH₃)₅Ru-His48 (Val68Asp, His81GIn, Cys110AIa)Mb and (NH₃)₅Ru-His48 (Val68GIu,His81GIn,Cys110Ala)Mb were measured to be .85(3)s^-1 and .30(2)s^-1, respectively. This data supports the hypothesis that entropy of 111 reduction and reorganization energy of ET are inversely related. The rates of forward and reverse ET for (NH₃)₅ Ru-His48 (Val68GIu, His81 GIn, Cys110AIa)ZnMb -7.2(5)•10⁴s^-1and 1.4(2)•10⁵s^-1, respectively- demonstrate that the placement of a highly polar residue nearby does not significantly change the reorganization energy of the photoactive Zn porphyrin.

The distal histidine imidazoles of (NH₃)₄isnRu-His48 SWMb and (NH₃)₅Ru-His48 SWMb were cyanated with BrCN. The intramolecular ET rates of these BrCN-modified Mb derivatives are 5.5(6)s^-1 and 3.2(5)s^-1, respectively. These respective rates are 20 and 10 times faster than those of their noncyanated counterparts after the differences in ET rate from driving force are scaled according to the Marcus equation. This increase in ET rate of the cyanated Mb derivatives is attributed to lower reorganization energy since the cyanated Mb heme is pentacoordinate in both oxidation states; whereas, the native Mb heme loses a water molecule upon reduction so that it changes from six to five coordinate. The reorganization energy from Fe-OH₂ dissociation is estimated to be .2eV. This conclusion is used to reconcile data from previous experiments in our lab. ET in photoactive porphyrin-substituted myoglobins proceed faster than predicted by Marcus Theory when it is assumed that the only difference in ET parameters between photoactive porphyrins and native heme systems is driving force. However, the data can be consistently fit to Marcus Theory if one corrects for the smaller reorganization in the photoactive porphyrin systems since they do not undergo a coordination change upon ET.

Finally, the intramolecular ET rate of (NH₃)₄isnRu-His48 SWMb was measured to be 3.0(4)s^-1. This rate is within experimental error of that for (NH₃)₄pyrRu-His48 SWMb even though the former has 80mV more driving force. One likely possibility for this observation is that the tetraamminepyridineruthenium group undergoes less reorganization upon ET than the tetraammineisonicotinamideruthenium group. Moreover, analysis of the (NH₃)₄isnRu-His48 SWMb experimental system gives a likely explanation of why ET was not observed previously in (NH₃)₄isnRu-Cytochrome C.

Toward Understanding Astrophysical Phenomena

Fast radio bursts (FRBs), a novel type of radio pulse, whose physics is not yet understood at all. Only a handful of FRBs had been detected when we started this project. Taking account of the scant observations, we put physical constraints on FRBs. We excluded proposals of a galactic origin for their extraordinarily high dispersion measures (DM), in particular stellar coronas and HII regions. Therefore our work supports an extragalactic origin for FRBs. We show that the resolved scattering tail of FRB 110220 is unlikely to be due to propagation through the intergalactic plasma. Instead the scattering is probably caused by the interstellar medium in the FRB's host galaxy, and indicates that this burst sits in the central region of that galaxy. Pulse durations of order $\ms$ constrain source sizes of FRBs implying enormous brightness temperatures and thus coherent emission. Electric fields near FRBs at cosmological distances would be so strong that they could accelerate free electrons from rest to relativistic energies in a single wave period. When we worked on FRBs, it was unclear whether they were genuine astronomical signals as distinct from `perytons', clearly terrestrial radio bursts, sharing some common properties with FRBs. Recently, in April 2015, astronomers discovered that perytons were emitted by microwave ovens. Radio chirps similar to FRBs were emitted when their doors opened while they were still heating. Evidence for the astronomical nature of FRBs has strengthened since our paper was published. Some bursts have been found to show linear and circular polarizations and Faraday rotation of the linear polarization has also been detected. I hope to resume working on FRBs in the near future. But after we completed our FRB paper, I decided to pause this project because of the lack of observational constraints.

The pulsar triple system, J0733+1715, has its orbital parameters fitted to high accuracy owing to the precise timing of the central $\ms$ pulsar. The two orbits are highly hierarchical, namely $P_{\mathrm{orb,1}}\ll P_{\mathrm{orb,2}}$, where 1 and 2 label the inner and outer white dwarf (WD) companions respectively. Moreover, their orbital planes almost coincide, providing a unique opportunity to study secular interaction associated purely with eccentricity beyond the solar system. Secular interaction only involves effect averaged over many orbits. Thus each companion can be represented by an elliptical wire with its mass distributed inversely proportional to its local orbital speed. Generally there exists a mutual torque, which vanishes only when their apsidal lines are parallel or anti-parallel. To maintain either mode, the eccentricity ratio, $e_1/e_2$, must be of the proper value, so that both apsidal lines precess together. For J0733+1715, $e_1\ll e_2$ for the parallel mode, while $e_1\gg e_2$ for the anti-parallel one. We show that the former precesses $\sim 10$ times slower than the latter. Currently the system is dominated by the parallel mode. Although only a little anti-parallel mode survives, both eccentricities especially $e_1$ oscillate on $\sim 10^3\yr$ timescale. Detectable changes would occur within $\sim 1\yr$. We demonstrate that the anti-parallel mode gets damped $\sim 10^4$ times faster than its parallel brother by any dissipative process diminishing $e_1$. If it is the tidal damping in the inner WD, we proceed to estimate its tidal quantity parameter ($Q$) to be $\sim 10^6$, which was poorly constrained by observations. However, tidal damping may also happen during the preceding low-mass X-ray binary (LMXB) phase or hydrogen thermal nuclear flashes. But, in both cases, the inner companion fills its Roche lobe and probably suffers mass/angular momentum loss, which might cause $e_1$ to grow rather than decay.

Several pairs of solar system satellites occupy mean motion resonances (MMRs). We divide these into two groups according to their proximity to exact resonance. Proximity is measured by the existence of a separatrix in phase space. MMRs between Io-Europa, Europa-Ganymede and Enceladus-Dione are too distant from exact resonance for a separatrix to appear. A separatrix is present only in the phase spaces of the Mimas-Tethys and Titan-Hyperion MMRs and their resonant arguments are the only ones to exhibit substantial librations. When a separatrix is present, tidal damping of eccentricity or inclination excites overstable librations that can lead to passage through resonance on the damping timescale. However, after investigation, we conclude that the librations in the Mimas-Tethys and Titan-Hyperion MMRs are fossils and do not result from overstability.

Rubble piles are common in the solar system. Monolithic elements touch their neighbors in small localized areas. Voids occupy a significant fraction of the volume. In a fluid-free environment, heat cannot conduct through voids; only radiation can transfer energy across them. We model the effective thermal conductivity of a rubble pile and show that it is proportional the square root of the pressure, $P$, for $P\leq \epsy^3\mu$ where $\epsy$ is the material's yield strain and $\mu$ its shear modulus. Our model provides an excellent fit to the depth dependence of the thermal conductivity in the top $140\,\mathrm{cm}$ of the lunar regolith. It also offers an explanation for the low thermal inertias of rocky asteroids and icy satellites. Lastly, we discuss how rubble piles slow down the cooling of small bodies such as asteroids.

Electromagnetic (EM) follow-up observations of gravitational wave (GW) events will help shed light on the nature of the sources, and more can be learned if the EM follow-ups can start as soon as the GW event becomes observable. In this paper, we propose a computationally efficient time-domain algorithm capable of detecting gravitational waves (GWs) from coalescing binaries of compact objects with nearly zero time delay. In case when the signal is strong enough, our algorithm also has the flexibility to trigger EM observation {\it before} the merger. The key to the efficiency of our algorithm arises from the use of chains of so-called Infinite Impulse Response (IIR) filters, which filter time-series data recursively. Computational cost is further reduced by a template interpolation technique that requires filtering to be done only for a much coarser template bank than otherwise required to sufficiently recover optimal signal-to-noise ratio. Towards future detectors with sensitivity extending to lower frequencies, our algorithm's computational cost is shown to increase rather insignificantly compared to the conventional time-domain correlation method. Moreover, at latencies of less than hundreds to thousands of seconds, this method is expected to be computationally more efficient than the straightforward frequency-domain method.