I. Studies on alkaline phosphatase activity in developing sea urchin embryos. II. Studies on the activation of protein biosynthesis in sea urchin eggs at fertilization

I. Alkaline phosphatase activity in the developing sea urchin Lytechinus pictus has been investigated with respect to intensity at various stages, ionic requirements and intracellular localization. The activity per embryo remains the same in the unfertilized egg, fertilized egg and cleavage stages. At a time just prior to gastrulation (about 10 hours after fertilization) the activity per embryo begins to rise and increases after 300 times over the activity in the cleavage stages during the next 60 hours.

The optimum ionic strength for enzymatic activity shows a wide peak at 0.6 to 1.0. Calcium and magnesium show an additional optimum at a concentration in the range of 0.02 to 0.07 molar. EDTA at concentrations of 0.0001 molar and higher shows a definite inhibition of activity.

The intracellular localization of alkaline phosphatase in homogenates of 72-hour embryos has been studied employing the differential centrifugation method. The major portion of the total activity in these homogenates was found in mitochondrial and microsomal fractions with less than 5% in the nuclear fraction and less than 2% in the final supernatant. The activity could be released from all fractions by treatment with sodium deoxycholate.

II. The activation of protein biosynthesis at fertilization in eggs of the sea urchins Lytechinus pictus and Strongylocentrotus purpuratus has been studied in both intact eggs and cell-free homogenates. It is shown that homogenates from both unfertilized and fertilized eggs are dependent on potassium and magnesium ions for optimum amino acid incorporation activity and in the case of the latter the concentration range is quite narrow. Though the optimum magnesium concentrations appear to differ slightly in homogenates of unfertilized and fertilized eggs, in no case was it observed that unfertilized egg homogenates were stimulated to incorporate at a level comparable to that of the fertilized eggs.

An activation of amino acid incorporation into protein has also been shown to occur in parthenogenetically activated non-nucleate sea urchin egg fragments or homogenates thereof. This activation resembles that in the fertilized whole egg or fragment both in amount and pattern of activation. Furthermore, it is shown that polyribosomes form in these non-nucleate fragments upon artificial activation. These findings are discussed along with possible mechanisms for activation of the system at fertilization.

Imaging the earth with ambient noise and earthquakes

In this thesis, I develop the velocity and structure models for the Los Angeles Basin and Southern Peru. The ultimate goal is to better understand the geological processes involved in the basin and subduction zone dynamics. The results are obtained from seismic interferometry using ambient noise and receiver functions using earthquake- generated waves. Some unusual signals specific to the local structures are also studied. The main findings are summarized as follows:

(1) Los Angeles Basin

The shear wave velocities range from 0.5 to 3.0 km/s in the sediments, with lateral gradients at the Newport-Inglewood, Compton-Los Alamitos, and Whittier Faults. The basin is a maximum of 8 km deep along the profile, and the Moho rises to a depth of 17 km under the basin. The basin has a stretch factor of 2.6 in the center decreasing to 1.3 at the edges, and is in approximate isostatic equilibrium. This "high-density" (~1 km spacing) "short-duration" (~1.5 month) experiment may serve as a prototype experiment that will allow basins to be covered by this type of low-cost survey.

(2) Peruvian subduction zone

Two prominent mid-crust structures are revealed in the 70 km thick crust under the Central Andes: a low-velocity zone interpreted as partially molten rocks beneath the Western Cordillera – Altiplano Plateau, and the underthrusting Brazilian Shield beneath the Eastern Cordillera. The low-velocity zone is oblique to the present trench, and possibly indicates the location of the volcanic arcs formed during the steepening of the Oligocene flat slab beneath the Altiplano Plateau.

The Nazca slab changes from normal dipping (~25 degrees) subduction in the southeast to flat subduction in the northwest of the study area. In the flat subduction regime, the slab subducts to ~100 km depth and then remains flat for ~300 km distance before it resumes a normal dipping geometry. The flat part closely follows the topography of the continental Moho above, indicating a strong suction force between the slab and the overriding plate. A high-velocity mantle wedge exists above the western half of the flat slab, which indicates the lack of melting and thus explains the cessation of the volcanism above. The velocity turns to normal values before the slab steepens again, indicating possible resumption of dehydration and ecologitization.

(3) Some unusual signals

Strong higher-mode Rayleigh waves due to the basin structure are observed in the periods less than 5 s. The particle motions provide a good test for distinguishing between the fundamental and higher mode. The precursor and coda waves relative to the interstation Rayleigh waves are observed, and modeled with a strong scatterer located in the active volcanic area in Southern Peru. In contrast with the usual receiver function analysis, multiples are extensively involved in this thesis. In the LA Basin, a good image is only from PpPs multiples, while in Peru, PpPp multiples contribute significantly to the final results.

I. The structure and properties of closed circular duplex DNA. II. Chemically interacting systems at equilibrium in a buoyant density gradient

I. The binding of the intercalating dye ethidium bromide to closed circular SV 40 DNA causes an unwinding of the duplex structure and a simultaneous and quantitatively equivalent unwinding of the superhelices. The buoyant densities and sedimentation velocities of both intact (I) and singly nicked (II) SV 40 DNAs were measured as a function of free dye concentration. The buoyant density data were used to determine the binding isotherms over a dye concentration range extending from 0 to 600 µg/m1 in 5.8 M CsCl. At high dye concentrations all of the binding sites in II, but not in I, are saturated. At free dye concentrations less than 5.4 µg/ml, I has a greater affinity for dye than II. At a critical amount of dye bound I and II have equal affinities, and at higher dye concentration I has a lower affinity than II. The number of superhelical turns, τ, present in I is calculated at each dye concentration using Fuller and Waring's (1964) estimate of the angle of duplex unwinding per intercalation. The results reveal that SV 40 DNA I contains about -13 superhelical turns in concentrated salt solutions.

The free energy of superhelix formation is calculated as a function of τ from a consideration of the effect of the superhelical turns upon the binding isotherm of ethidium bromide to SV 40 DNA I. The value of the free energy is about 100 kcal/mole DNA in the native molecule. The free energy estimates are used to calculate the pitch and radius of the superhelix as a function of the number of superhelical turns. The pitch and radius of the native I superhelix are 430 Å and 135 Å, respectively.

A buoyant density method for the isolation and detection of closed circular DNA is described. The method is based upon the reduced binding of the intercalating dye, ethidium bromide, by closed circular DNA. In an application of this method it is found that HeLa cells contain in addition to closed circular mitochondrial DNA of mean length 4.81 microns, a heterogeneous group of smaller DNA molecules which vary in size from 0.2 to 3.5 microns and a paucidisperse group of multiples of the mitochondrial length.

II. The general theory is presented for the sedimentation equilibrium of a macromolecule in a concentrated binary solvent in the presence of an additional reacting small molecule. Equations are derived for the calculation of the buoyant density of the complex and for the determination of the binding isotherm of the reagent to the macrospecies. The standard buoyant density, a thermodynamic function, is defined and the density gradients which characterize the four component system are derived. The theory is applied to the specific cases of the binding of ethidium bromide to SV 40 DNA and of the binding of mercury and silver to DNA.

In search of slow-moving ionizing massive particles

Many particles proposed by theories, such as GUT monopoles, nuclearites and 1/5 charge superstring particles, can be categorized as Slow-moving, Ionizing, Massive Particles (SIMPs).

Detailed calculations of the signal-to-noise ratios in vanous acoustic and mechanical methods for detecting such SIMPs are presented. It is shown that the previous belief that such methods are intrinsically prohibited by the thermal noise is incorrect, and that ways to solve the thermal noise problem are already within the reach of today's technology. In fact, many running and finished gravitational wave detection ( GWD) experiments are already sensitive to certain SIMPs. As an example, a published GWD result is used to obtain a flux limit for nuclearites.

The result of a search using a scintillator array on Earth's surface is reported. A flux limit of 4.7 x 10^(-12) cm^(-2)sr^(-1)s^(-1) (90% c.l.) is set for any SIMP with 2.7 x 10^(-4) less than β less than 5 x 10^(-3) and ionization greater than 1/3 of minimum ionizing muons. Although this limit is above the limits from underground experiments for typical supermassive particles (10^(16)GeV), it is a new limit in certain β and ionization regions for less massive ones (~10^9 GeV) not able to penetrate deep underground, and implies a stringent limit on the fraction of the dark matter that can be composed of massive electrically and/ or magnetically charged particles.

The prospect of the future SIMP search in the MACRO detector is discussed. The special problem of SIMP trigger is examined and a circuit proposed, which may solve most of the problems of the previous ones proposed or used by others and may even enable MACRO to detect certain SIMP species with β as low as the orbital velocity around the earth.

Sequence specific complexation of b DNA at sites containing g,c base pairs

A series of eight related analogs of distamycin A has been synthesized. Footprinting and affinity cleaving reveal that only two of the analogs, pyridine-2- car box amide-netropsin (2-Py N) and 1-methylimidazole-2-carboxamide-netrops in (2-ImN), bind to DNA with a specificity different from that of the parent compound. A new class of sites, represented by a TGACT sequence, is a strong site for 2-PyN binding, and the major recognition site for 2-ImN on DNA. Both compounds recognize the G•C bp specifically, although A's and T's in the site may be interchanged without penalty. Additional A•T bp outside the binding site increase the binding affinity. The compounds bind in the minor groove of the DNA sequence, but protect both grooves from dimethylsulfate. The binding evidence suggests that 2-PyN or 2-ImN binding induces a DNA conformational change.

In order to understand this sequence specific complexation better, the Ackers quantitative footprinting method for measuring individual site affinity constants has been extended to small molecules. MPE•Fe(II) cleavage reactions over a 10^5 range of free ligand concentrations are analyzed by gel electrophoresis. The decrease in cleavage is calculated by densitometry of a gel autoradiogram. The apparent fraction of DNA bound is then calculated from the amount of cleavage protection. The data is fitted to a theoretical curve using non-linear least squares techniques. Affinity constants at four individual sites are determined simultaneously. The distamycin A analog binds solely at A•T rich sites. Affinities range from 10^(6)- 10^(7)M^(-1) The data for parent compound D fit closely to a monomeric binding curve. 2-PyN binds both A•T sites and the TGTCA site with an apparent affinity constant of 10^(5) M^(-1). 2-ImN binds A•T sites with affinities less than 5 x 10^(4) M^(-1). The affinity of 2-ImN for the TGTCA site does not change significantly from the 2-PyN value. At the TGTCA site, the experimental data fit a dimeric binding curve better than a monomeric curve. Both 2-PyN and 2-ImN have substantially lower DNA affinities than closely related compounds.

In order to probe the requirements of this new binding site, fourteen other derivatives have been synthesized and tested. All compounds that recognize the TGTCA site have a heterocyclic aromatic nitrogen ortho to the N or C-terminal amide of the netropsin subunit. Specificity is strongly affected by the overall length of the small molecule. Only compounds that consist of at least three aromatic rings linked by amides exhibit TGTCA site binding. Specificity is only weakly altered by substitution on the pyridine ring, which correlates best with steric factors. A model is proposed for TGTCA site binding that has as its key feature hydrogen bonding to both G's by the small molecule. The specificity is determined by the sequence dependence of the distance between G's.

One derivative of 2-PyN exhibits pH dependent sequence specificity. At low pH, 4-dimethylaminopyridine-2-carboxamide-netropsin binds tightly to A•T sites. At high pH, 4-Me_(2)NPyN binds most tightly to the TGTCA site. In aqueous solution, this compound protonates at the pyridine nitrogen at pH 6. Thus presence of the protonated form correlates with A•T specificity.

The binding site of a class of eukaryotic transcriptional activators typified by yeast protein GCN4 and the mammalian oncogene Jun contains a strong 2-ImN binding site. Specificity requirements for the protein and small molecule are similar. GCN4 and 2-lmN bind simultaneously to the same binding site. GCN4 alters the cleavage pattern of 2-ImN-EDTA derivative at only one of its binding sites. The details of the interaction suggest that GCN4 alters the conformation of an AAAAAAA sequence adjacent to its binding site. The presence of a yeast counterpart to Jun partially blocks 2-lmN binding. The differences do not appear to be caused by direct interactions between 2-lmN and the proteins, but by induced conformational changes in the DNA protein complex. It is likely that the observed differences in complexation are involved in the varying sequence specificity of these proteins.

Scanning tunneling spectroscopic studies on high-temperature superconductors and Dirac materials

This thesis details the investigations of the unconventional low-energy quasiparticle excitations in electron-type cuprate superconductors and electron-type ferrous superconductors as well as the electronic properties of Dirac fermions in graphene and three-dimensional strong topological insulators through experimental studies using spatially resolved scanning tunneling spectroscopy (STS) experiments.

Magnetic-field- and temperature-dependent evolution of the spatially resolved quasiparticle spectra in the electron-type cuprate La_0.1Sr_0.9CuO₂ (La-112) T_C = 43 K, are investigated experimentally. For temperature (T) less than the superconducting transition temperature (TC), and in zero field, the quasiparticle spectra of La-112 exhibits gapped behavior with two coherence peaks and no satellite features. For magnetic field measurements at T < TC, first ever observation of vortices in La-112 are reported. Moreover, pseudogap-like spectra are revealed inside the core of vortices, where superconductivity is suppressed. The intra-vortex pseudogap-like spectra are characterized by an energy gap of V_PG = 8.5 ± 0.6 meV, while the inter-vortex quasiparticle spectra shows larger peak-to-peak gap values characterized by Δ_pk-pk(H) >V_PG, and Δ_pk-pk (0)=12.2 ± 0.8 meV > Δ_pk-pk (H > 0). The quasiparticle spectra are found to be gapped at all locations up to the highest magnetic field examined (H = 6T) and reveal an apparent low-energy cutoff at the V_PG energy scale.

Magnetic-field- and temperature-dependent evolution of the spatially resolved quasiparticle spectra in the electron-type "122" iron-based Ba(Fe_1-xCo_x)₂A_s2 are investigated for multiple doping levels (x = 0.06, 0.08, 0.12 with T_C= 14 K, 24 K, and 20 K). For all doping levels and the T < T_C, two-gap superconductivity is observed. Both superconducting gaps decrease monotonically in size with increasing temperature and disappear for temperatures above the superconducting transition temperature, T_C. Magnetic resonant modes that follow the temperature dependence of the superconducting gaps have been identified in the tunneling quasiparticle spectra. Together with quasiparticle interference (QPI) analysis and magnetic field studies, this provides strong evidence for two-gap sign-changing s-wave superconductivity.

Additionally spatial scanning tunneling spectroscopic studies are performed on mechanically exfoliated graphene and chemical vapor deposition grown graphene. In all cases lattice strain exerts a strong influence on the electronic properties of the sample. In particular topological defects give rise to pseudomagnetic fields (B ~ 50 Tesla) and charging effects resulting in quantized conductance peaks associated with the integer and fractional Quantum Hall States.

Finally, spectroscopic studies on the 3D-STI, Bi₂Se₃ found evidence of impurity resonance in the surface state. The impurities are in the unitary limit and the spectral resonances are localized spatially to within ~ 0.2 nm of the impurity. The spectral weight of the impurity resonance diverges as the Fermi energy approaches the Dirac point and the rapid recovery of the surface state suggests robust topological protection against perturbations that preserve time reversal symmetry.

The electric field of a weakly electric fish

Freshwater fish of the genus Apteronotus (family Gymnotidae) generate a weak, high frequency electric field (< 100 mV/cm, 0.5-10 kHz) which permeates their local environment. These nocturnal fish are acutely sensitive to perturbations in their electric field caused by other electric fish, and nearby objects whose impedance is different from the surrounding water. This thesis presents high temporal and spatial resolution maps of the electric potential and field on and near Apteronotus. The fish's electric field is a complicated and highly stable function of space and time. Its characteristics, such as spectral composition, timing, and rate of attenuation, are examined in terms of physical constraints, and their possible functional roles in electroreception.

Temporal jitter of the periodic field is less than 1 µsec. However, electrocyte activity is not globally synchronous along the fish 's electric organ. The propagation of electrocyte activation down the fish's body produces a rotation of the electric field vector in the caudal part of the fish. This may assist the fish in identifying nonsymmetrical objects, and could also confuse electrosensory predators that try to locate Apteronotus by following its fieldlines. The propagation also results in a complex spatiotemporal pattern of the EOD potential near the fish. Visualizing the potential on the same and different fish over timescales of several months suggests that it is stable and could serve as a unique signature for individual fish.

Measurements of the electric field were used to calculate the effects of simple objects on the fish's electric field. The shape of the perturbation or "electric image" on the fish's skin is relatively independent of a simple object's size, conductivity, and rostrocaudal location, and therefore could unambiguously determine object distance. The range of electrolocation may depend on both the size of objects and their rostrocaudal location. Only objects with very large dielectric constants cause appreciable phase shifts, and these are strongly dependent on the water conductivity.

Studies of Zγ production and constraints on anomalous triple gauge couplings in pp collisions at √s = 7 TeV

In this thesis, we test the electroweak sector of the Standard Model of particle physics through the measurements of the cross section of the simultaneous production of the neutral weak boson Z and photon γ, and the limits on the anomalous Zγγ and ZZγ triple gauge couplings h3 and h4 with the Z decaying to leptons (electrons and muons). We analyze events collected in proton-proton collisions at center of mass energy of sqrt(s) = 7 TeV corresponding to an integrated luminosity of 5.0 inverse femtobarn. The analyzed events were recorded by the Compact Muon Solenoid detector at the Large Hadron Collider in 2011.

The production cross section has been measured for hard photons with transverse momentum greater than 15 GeV that are separated from the the final state leptons in the eta-phi plane by Delta R greater than 0.7, whose sum of the transverse energy of hadrons over the transverse energy of the photon in a cone around the photon with Delta R less than 0.3 is less than 0.5, and with the invariant mass of the dilepton system greater than 50 GeV. The measured cross section value is 5.33 +/- 0.08 (stat.) +/- 0.25 (syst.) +/- 0.12 (lumi.) picobarn. This is compatible with the Standard Model prediction that includes next-to-leading-order QCD contributions: 5.45 +/- 0.27 picobarn.

The measured 95 % confidence-level upper limits on the absolute values of the anomalous couplings h3 and h4 are 0.01 and 8.8E-5 for the Zγγ interactions, and, 8.6E-3 and 8.0E-5 for the ZZγ interactions. These values are also compatible with the Standard Model where they vanish in the tree-level approximation. They extend the sensitivity of the 2012 results from the ATLAS collaboration based on 1.02 inverse femtobarn of data by a factor of 2.4 to 3.1.

Reactivity of titanocene methylidene with metal halides, alkene sulfides, and alkene oxides

Titanocene metallacyclobutanes show a wide variety of reactivites with organic and inorganic reagents. Their reactions include methylene transfer to organic carbonyls, formation of enolates, electron transfer from activated alkyl chlorides, olefin metathesis, ring opening polymerization. Recently, preparations of heterobinuclear µ-methylene complexes were reported. In this thesis, mechanistic, synthetic, and structural studies of the heterobinuclear µ-methylene complexes will be described. Also, the reaction of titanocene methylidene trimethylphosphine complex with alkene sulfide and styrene sulfide will be presented.

Heterobinuclear µ-methylene-µ-methyl complexes C_(p2)Ti(µ-CH_2)( µ-CH_3)M(1,5-COD) have been prepared (M = Rh, Ir). X-ray crystallography showed that the methyl group of the complex was bonded to the rhodium and bridges to the titanium through an agostic bond. The ^(1)H,^(13)CNMR, IR spectra along with partial deuteration studies supported the structure in both solution and solid state. Activation of the agostic bond is demonstrated by the equilibration of the µ-CH_3 and µ-CH_2 groups. A nonlinear Arrhenius plot, an unusually large kinetic isotope effect (24(5)), and a large negative activation entropy (-64(3)eu) can be explained by the quantum-mechanical tunneling. Calculated rate constants with Bell-type barrier fitted well with the observed one. This equilibration was best explained by a 4e-4c mechanism (or σ bond metathesis) with the character of quantum-mechanical tunneling.

Heterobinuclear µ-methylene-µ-phenyl complexes were synthesized. Structural study of C_(p2)Ti(µ-CH_(2))(µ-p-Me_(2)NC_(6)H_(4))Rh(l,5-COD) showed that the two metal atoms are bridged by the methylene carbon and the ipso carbon of the p-N,N-dimethylarninophenyl group. The analogous structure of C_(p2))Ti(µ-CH_(2))(µ-o-MeOC_(6)H_(4))Rh(1,5-COD) has been verified by the differential NOE. The aromaticity of the phenyl group observed by ^(1)H NMR, was confirmed by the comparison of the C-C bond lengths in the crystallographic structure. The unusual downfield shifts of the ipso carbon in the ^(13)C NMR are assumed to be an indication of the interaction between the ipso carbon and electron-deficient titanium.

Titanium-platinum heterobinuclear µ-methylene complexes C_(p2)Ti(µ-CH_(2))(µ -X)Pt(Me)(PM_(2)Ph) have been prepared (X= Cl, Me). Structural studies indicate the following:(1) the Ti-CH2 bond possesses residual double bond character, (2) there is a dative Pt→Ti interaction which may be regarded as a π back donation from the platinum atom to the 'Ti=CH_(2)'' group, and (3) the µ-CH_3 group is bound to the titanium atom through a three-center, two-electron agostic bond.

Titanocene (η^(2)-thioformaldehyde)•PMe_3 was prepared from C_(p2)Ti=CH_(2)•PMe_3 and sulfur-containing organic compounds (e.g. alkene sulfide, triphenylphosphine sulfide) including elemental sulfur. Mechanistic studies utilizing trans-styrene sulfide-d_1 suggested the stepwise reaction to explain equimolar mixture of trans- and cis-styrene-d_1 as by-products. The product reacted with methyl iodide to produce cationic titanocene (η_(2)-thiomethoxymethyl) complex. Complexes having less coordinating anion like BF_4 or BPh_4 could be obtained through metathesis. Together with structural analyses, the further reactivities of the complexes have been explored.

The complex C_(p2)TiOCH_(2)CH(Ph)CH_2 was prepared from the compound C_(p2)Ti=CH_(2)-PMe_3 and styrene oxide. The product was characterized with ^(1)H-^(1)H correlated 2-dimensional NMR, selective decoupling of ^(1)H NMR, and differential NOE. Stereospecificity of deuterium in the product was lost when trans-styrene oxide-d_1 was allowed to react. Relative rates of the reaction were measured with varying substituents on the phenyl ring. Better linearity (r = -0.98, p^(+) = -0.79) was observed with σ_(p)^(+)than σ(r = -0.87, p = -1.26). The small magnitude of p^+ value and stereospecificity loss during the formation of product were best explained by the generation of biradicals, but partial generation of charge cannot be excluded. Carbonylation of the product followed by exposure to iodine yields the corresponding β-phenyl γ-lactone.

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.

Dating and characterizing late holocene earthquakes using paleomagnetics

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

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

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

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

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

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

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

Part I. A study of the velocity structure of the earth by the use of core phases. Part II. The 1971 San Fernando earthquake series focal mechanisms and tectonics

The initial objective of Part I was to determine the nature of upper mantle discontinuities, the average velocities through the mantle, and differences between mantle structure under continents and oceans by the use of P'dP', the seismic core phase P'P' (PKPPKP) that reflects at depth d in the mantle. In order to accomplish this, it was found necessary to also investigate core phases themselves and their inferences on core structure. P'dP' at both single stations and at the LASA array in Montana indicates that the following zones are candidates for discontinuities with varying degrees of confidence: 800-950 km, weak; 630-670 km, strongest; 500-600 km, strong but interpretation in doubt; 350-415 km, fair; 280-300 km, strong, varying in depth; 100-200 km, strong, varying in depth, may be the bottom of the low-velocity zone. It is estimated that a single station cannot easily discriminate between asymmetric P'P' and P'dP' for lead times of about 30 sec from the main P'P' phase, but the LASA array reduces this uncertainty range to less than 10 sec. The problems of scatter of P'P' main-phase times, mainly due to asymmetric P'P', incorrect identification of the branch, and lack of the proper velocity structure at the velocity point, are avoided and the analysis shows that one-way travel of P waves through oceanic mantle is delayed by 0.65 to 0.95 sec relative to United States mid-continental mantle.

A new P-wave velocity core model is constructed from observed times, dt/dΔ's, and relative amplitudes of P'; the observed times of SKS, SKKS, and PKiKP; and a new mantle-velocity determination by Jordan and Anderson. The new core model is smooth except for a discontinuity at the inner-core boundary determined to be at a radius of 1215 km. Short-period amplitude data do not require the inner core Q to be significantly lower than that of the outer core. Several lines of evidence show that most, if not all, of the arrivals preceding the DF branch of P' at distances shorter than 143° are due to scattering as proposed by Haddon and not due to spherically symmetric discontinuities just above the inner core as previously believed. Calculation of the travel-time distribution of scattered phases and comparison with published data show that the strongest scattering takes place at or near the core-mantle boundary close to the seismic station.

In Part II, the largest events in the San Fernando earthquake series, initiated by the main shock at 14 00 41.8 GMT on February 9, 1971, were chosen for analysis from the first three months of activity, 87 events in all. The initial rupture location coincides with the lower, northernmost edge of the main north-dipping thrust fault and the aftershock distribution. The best focal mechanism fit to the main shock P-wave first motions constrains the fault plane parameters to: strike, N 67° (± 6°) W; dip, 52° (± 3°) NE; rake, 72° (67°-95°) left lateral. Focal mechanisms of the aftershocks clearly outline a downstep of the western edge of the main thrust fault surface along a northeast-trending flexure. Faulting on this downstep is left-lateral strike-slip and dominates the strain release of the aftershock series, which indicates that the downstep limited the main event rupture on the west. The main thrust fault surface dips at about 35° to the northeast at shallow depths and probably steepens to 50° below a depth of 8 km. This steep dip at depth is a characteristic of other thrust faults in the Transverse Ranges and indicates the presence at depth of laterally-varying vertical forces that are probably due to buckling or overriding that causes some upward redirection of a dominant north-south horizontal compression. Two sets of events exhibit normal dip-slip motion with shallow hypocenters and correlate with areas of ground subsidence deduced from gravity data. Several lines of evidence indicate that a horizontal compressional stress in a north or north-northwest direction was added to the stresses in the aftershock area 12 days after the main shock. After this change, events were contained in bursts along the downstep and sequencing within the bursts provides evidence for an earthquake-triggering phenomenon that propagates with speeds of 5 to 15 km/day. Seismicity before the San Fernando series and the mapped structure of the area suggest that the downstep of the main fault surface is not a localized discontinuity but is part of a zone of weakness extending from Point Dume, near Malibu, to Palmdale on the San Andreas fault. This zone is interpreted as a decoupling boundary between crustal blocks that permits them to deform separately in the prevalent crustal-shortening mode of the Transverse Ranges region.

Metasomatism and magmatic assimilation at a gabbro-limestone contact, Christmas Mountains, Big Bend region, Texas

A composite stock of alkaline gabbro and syenite is intrusive into limestone of the Del Carmen, Sue Peake and Santa Elena Formations at the northwest end of the Christmas Mountains. There is abundant evidence of solution of wallrock by magma but nowhere are gabbro and limestone in direct contact. The sequence of lithologies developed across the intrusive contact and across xenoliths is gabbro, pyroxenite, calc-silicate skarn, marble. Pyroxenite is made up of euhedral crystals of titanaugite and sphene in a leucocratic matrix of nepheline, Wollastonite and alkali feldspar. The uneven modal distribution of phases in pyroxenite and the occurrence' of nepheline syenite dikes, intrusive into pyroxenite and skarn, suggest that pyroxenite represents an accumulation of clinopyroxene "cemented" together by late-solidifying residual magma of nepheline syenite composition. Assimilation of limestone by gabbroic magma involves reactions between calcite and magma and/or crystals in equilibrium with magma and crystallization of phases in which the magma is saturated, to supply energy for the solution reaction. Gabbroic magma was saturated with plagioclase and clinopyroxene at the time of emplacement. The textural and mineralogic features of pyroxenite can be produced by the reaction 2( 1-X) CALCITE + AN_XAB_l-X = (1-X) NEPHELINE+ 2(1-X) WOLLASTONITE+ X ANORTHITE+ 2(1-X) CO₂. Plagioclase in pyroxenite has corroded margins and is rimmed by nepheline, suggestive of resorption by magma. Anorthite and wollastonite enter solid solution in titanaugite. For each mole of calcite dissolved, approximately one mole of clinopyroxene was crystallized. Thus the amount of limestone that may be assimilated is limited by the concentration of potential clinopyroxene in the magma. Wollastonite appears as a phase when magma has been depleted in iron and magnesium by crystallization of titanaugite. The predominance of mafic and ultramafic compositions among contaminated rocks and their restriction to a narrow zone along the intrusive contact provides little evidence for the generation of a significant volume of desilicated magma as a result of limestone assimilation.

Within 60 m of the intrusive contact with the gabbro, nodular chert in the Santa Elena Limestone reacted with the enveloping marble to form spherical nodules of high-temperature calc-silicate minerals. The phases wollastonite, rankinite, spurrite, tilleyite and calcite, form a series of sharply-bounded, concentric monomineralic and two-phase shells which record a step-wise decrease in silica content from the core of a nodule to its rim. Mineral zones in the nodules vary 'with distance from the gabbro as follows:

0-5 m CALCITE + SPURRITE + RANKINITE + WOLLASTONITE
5-16 m CALCITE + TILLEYITE ± SPURRITE + RANKINITE + WOLLASTONITE
16-31 m CALCITE + TILLEYITE + WOLLASTONITE
31-60 m CALCITE + WOLLASTONITE
60-plus CALCITE + QUARTZ

The mineral of a one-phase zone is compatible with the phases bounding it on either side but these phases are incompatible in the same volume of P-T-X_CO2.

Growth of a monomineralio zone is initiated by reaction between minerals of adjacent one-phase zones which become unstable with rising temperature to form a thin layer of a new single phase that separates the reactants and is compatible with both of them. Because the mineral of the new zone is in equilibrium with the phases at both of its contacts, gradients in the chemical potentials of the exchangeable components are established across it. Although zone boundaries mark discontinuities in the gradients of bulk composition, two-phase equilibria at the contacts demonstrate that the chemical potentials are continuous. Hence, Ca, Si and CO₂ were redistributed in the growing nodule by diffusion. A monomineralic zone grows at the expense of an adjacent zone by reaction between diffusing components and the mineral of the adjacent zone. Equilibria between two phases at zone boundaries buffers the chemical potentials of the diffusing species. Thus, within a monomineralic zone, the chemical potentials of the diffusing components are controlled external to the local assemblage by the two-phase equilibria at the zone boundaries.

Mineralogically zoned calc-silicate skarn occurs as a narrow band that separates pyroxenite and marble along the intrusive contact and forms a rim on marble xenoliths in gabbro. Skarn consists of melilite or idocrase pseudomorphs of melili te, one or two . stoichiometric calcsilicate phases and accessory Ti-Zr garnet, perovskite and magnetite. The sequence of mineral zones from pyroxenite to marble, defined by a characteristic calc-silicate, is wollastonite, rankinite, spurrite, calcite. Mineral assemblages of adjacent skarn zones are compatible and the set of zones in a skarn band defines a facies type, indicating that the different mineral assemblages represent different bulk compositions recrystallized under identical conditions. The number of phases in each zone is less than the number that might be expected to result from metamorphism of a general bulk composition under conditions of equilibrium, trivariant in P, T and u_CO2. The "special" bulk composition of each zone is controlled by reaction between phases of the zones bounding it on either side. The continuity of the gradients of composition of melilite and garnet solid solutions across the skarn is consistent with the local equilibrium hypothesis and verifies that diffusion was the mechanism of mass transport. The formula proportions of Ti and Zr in garnet from skarn vary antithetically with that of Si Which systematically decreases from pyroxenite to marble. The chemical potential of Si in each skarn zone was controlled by the coexisting stoichiometric calc-silicate phases in the assemblage. Thus the formula proportion of Si in garnet is a direct measure of the chemical potential of Si from point to point in skarn. Reaction between gabbroic magma saturated with plagioclase and clinopyroxene produced nepheline pyroxenite and melilite-wollastonite skarn. The calcsilicate zones result from reaction between calcite and wollastonite to form spurrite and rankinite.

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Origin of Th/U variations in chondritic meteorites

Isotope dilution thorium and uranium analyses of the Harleton chondrite show a larger scatter than previously observed in equilibrated ordinary chondrites (EOC). The linear correlation of Th/U with 1/U in Harleton (and all EOC data) is produced by variation in the chlorapatite to merrillite mixing ratio. Apatite variations control the U concentrations. Phosphorus variations are compensated by inverse variations in U to preserve the Th/U vs. 1/U correlation. Because the Th/U variations reflect phosphate ampling, a weighted Th/U average should converge to an improved solar system Th/U. We obtain Th/U=3.53 (1_-mean=0.10), significantly lower and more precise than previous estimates.

To test whether apatite also produces Th/U variation in CI and CM chondrites, we performed P analyses on the solutions from leaching experiments of Orgueil and Murchison meteorites.

A linear Th/U vs. 1/U correlation in CI can be explained by redistribution of hexavalent U by aqueous fluids into carbonates and sulfates.

Unlike CI and EOC, whole rock Th/U variations in CMs are mostly due to Th variations. A Th/U vs. 1/U linear correlation suggested by previous data for CMs is not real. We distinguish 4 components responsible for the whole rock Th/U variations: (1) P and actinide-depleted matrix containing small amounts of U-rich carbonate/sulfate phases (similar to CIs); (2) CAIs and (3) chondrules are major reservoirs for actinides, (4) an easily leachable phase of high Th/U. likely carbonate produced by CAI alteration. Phosphates play a minor role as actinide and P carrier phases in CM chondrites.

Using our Th/U and minimum galactic ages from halo globular clusters, we calculate relative supernovae production rates for ²³²Th/²³⁸U and ²³⁵U/²³⁸U for different models of r-process nucleosynthesis. For uniform galactic production, the beginning of the r-process nucleosynthesis must be less than 13 Gyr. Exponentially decreasing production is also consistent with a 13 Gyr age, but very slow decay times are required (less than 35 Gyr), approaching the uniform production. The 15 Gyr Galaxy requires either a fast initial production growth (infall time constant less than 0.5 Gyr) followed by very low decrease (decay time constant greater than 100 Gyr), or the fastest possible decrease (≈8 Gyr) preceded by slow in fall (≈7.5 Gyr).