Neural representation of auditory temporal structure

Neurons in the songbird forebrain nucleus HVc are highly sensitive to auditory temporal context and have some of the most complex auditory tuning properties yet discovered. HVc is crucial for learning, perceiving, and producing song, thus it is important to understand the neural circuitry and mechanisms that give rise to these remarkable auditory response properties. This thesis investigates these issues experimentally and computationally.

Extracellular studies reported here compare the auditory context sensitivity of neurons in HV c with neurons in the afferent areas of field L. These demonstrate that there is a substantial increase in the auditory temporal context sensitivity from the areas of field L to HVc. Whole-cell recordings of HVc neurons from acute brain slices are described which show that excitatory synaptic transmission between HVc neurons involve the release of glutamate and the activation of both AMPA/kainate and NMDA-type glutamate receptors. Additionally, widespread inhibitory interactions exist between HVc neurons that are mediated by postsynaptic GABA_A receptors. Intracellular recordings of HVc auditory neurons in vivo provides evidence that HV c neurons encode information about temporal structure using a variety of cellular and synaptic mechanisms including syllable-specific inhibition, excitatory post-synaptic potentials with a range of different time courses, and burst-firing, and song-specific hyperpolarization.

The final part of this thesis presents two computational approaches for representing and learning temporal structure. The first method utilizes comput ational elements that are analogous to temporal combination sensitive neurons in HVc. A network of these elements can learn using local information and lateral inhibition. The second method presents a more general framework which allows a network to discover mixtures of temporal features in a continuous stream of input.

Seismic structure above and below the core-mantle boundary

Seismic structure above and below the core-mantle boundary (CMB) has been studied through use of travel time and waveform analyses of several different seismic wave groups. Anomalous systematic trends in observables document mantle heterogeneity on both large and small scales. Analog and digital data has been utilized, and in many cases the analog data has been optically scanned and digitized prior to analysis.

Differential travel times of S - SKS are shown to be an excellent diagnostic of anomalous lower mantle shear velocity (V s) structure. Wavepath geometries beneath the central Pacific exhibit large S- SKS travel time residuals (up to 10 sec), and are consistent with a large scale 0(1000 km) slower than average V_s region (≥3%). S - SKS times for paths traversing this region exhibit smaller scale patterns and trends 0(100 km) indicating V_s perturbations on many scale lengths. These times are compared to predictions of three tomographically derived aspherical models: MDLSH of Tanimoto [1990], model SH12_WM13 of Suet al. [1992], and model SH.10c.17 of Masters et al. [1992]. Qualitative agreement between the tomographic model predictions and observations is encouraging, varying from fair to good. However, inconsistencies are present and suggest anomalies in the lower mantle of scale length smaller than the present 2000+ km scale resolution of tomographic models. 2-D wave propagation experiments show the importance of inhomogeneous raypaths when considering lateral heterogeneities in the lowermost mantle.

A dataset of waveforms and differential travel times of S, ScS, and the arrival from the D" layer, Scd, provides evidence for a laterally varying V_s velocity discontinuity at the base of the mantle. Two different localized D" regions beneath the central Pacific have been investigated. Predictions from a model having a V_s discontinuity 180 km above the CMB agree well with observations for an eastern mid-Pacific CMB region. This thickness differs from V_s discontinuity thicknesses found in other regions, such as a localized region beneath the western Pacific, which average near 280 km. The "sharpness" of the V_s jump at the top of D", i.e., the depth range over which the V_s increase occurs, is not resolved by our data, and our data can in fact may be modeled equally well by a lower mantle with the increase in V_s at the top of D" occurring over a 100 krn depth range. It is difficult at present to correlate D" thicknesses from this study to overall lower mantle heterogeneity, due to uncertainties in the 3-D models, as well as poor coverage in maps of D" discontinuity thicknesses.

P-wave velocity structure (V_p) at the base of the mantle is explored using the seismic phases SKS and SPdKS. SPdKS is formed when SKS waves at distances around 107° are incident upon the CMB with a slowness that allows for coupling with diffracted P-waves at the base of the mantle. The P-wave diffraction occurs at both the SKS entrance and exit locations of the outer core. SP_dKS arrives slightly later in time than SKS, having a wave path through the mantle and core very close to SKS. The difference time between SKS and SP_dKS strongly depends on V_p at the base of the mantle near SK Score entrance and exit points. Observations from deep focus Fiji-Tonga events recorded by North American stations, and South American events recorded by European and Eurasian stations exhibit anomalously large SP_dKS - SKS difference times. SKS and the later arriving SP_dKS phase are separated by several seconds more than predictions made by 1-D reference models, such as the global average PREM [Dziewonski and Anderson, 1981] model. Models having a pronounced low-velocity zone (5%) in V_p in the bottom 50-100 km of the mantle predict the size of the observed SP_dK S-SKS anomalies. Raypath perturbations from lower mantle V_s structure may also be contributing to the observed anomalies.

Outer core structure is investigated using the family of SmKS (m=2,3,4) seismic waves. SmKS are waves that travel as S-waves in the mantle, P-waves in the core, and reflect (m-1) times on the underside of the CMB, and are well-suited for constraining outermost core V_p structure. This is due to closeness of the mantle paths and also the shallow depth range these waves travel in the outermost core. S3KS - S2KS and S4KS - S3KS differential travel times were measured using the cross-correlation method and compared to those from reflectivity synthetics created from core models of past studies. High quality recordings from a deep focus Java Sea event which sample the outer core beneath the northern Pacific, the Arctic, and northwestern North America (spanning 1/8th of the core's surface area), have SmKS wavepaths that traverse regions where lower mantle heterogeneity is pre- dieted small, and are well-modeled by the PREM core model, with possibly a small V_p decrease (1.5%) in the outermost 50 km of the core. Such a reduction implies chemical stratification in this 50 km zone, though this model feature is not uniquely resolved. Data having wave paths through areas of known D" heterogeneity (±2% and greater), such as the source-side of SmKS lower mantle paths from Fiji-Tonga to Eurasia and Africa, exhibit systematic SmKS differential time anomalies of up to several seconds. 2-D wave propagation experiments demonstrate how large scale lower mantle velocity perturbations can explain long wavelength behavior of such anomalous SmKS times. When improperly accounted for, lower mantle heterogeneity maps directly into core structure. Raypaths departing from homogeneity play an important role in producing SmKS anomalies. The existence of outermost core heterogeneity is difficult to resolve at present due to uncertainties in global lower mantle structure. Resolving a one-dimensional chemically stratified outermost core also remains difficult due to the same uncertainties. Restricting study to higher multiples of SmKS (m=2,3,4) can help reduce the affect of mantle heterogeneity due to the closeness of the mantle legs of the wavepaths. SmKS waves are ideal in providing additional information on the details of lower mantle heterogeneity.

Teleseismic array analysis of upper mantle compressional velocity structure

Large quantities of teleseismic short-period seismograms recorded at SCARLET provide travel time, apparent velocity and waveform data for study of upper mantle compressional velocity structure. Relative array analysis of arrival times from distant (30° < Δ < 95°) earthquakes at all azimuths constrains lateral velocity variations beneath southern California. We compare dT/dΔ back azimuth and averaged arrival time estimates from the entire network for 154 events to the same parameters derived from small subsets of SCARLET. Patterns of mislocation vectors for over 100 overlapping subarrays delimit the spatial extent of an east-west striking, high-velocity anomaly beneath the Transverse Ranges. Thin lens analysis of the averaged arrival time differences, called 'net delay' data, requires the mean depth of the corresponding lens to be more than 100 km. Our results are consistent with the PKP-delay times of Hadley and Kanamori (1977), who first proposed the high-velocity feature, but we place the anomalous material at substantially greater depths than their 40-100 km estimate.

Detailed analysis of travel time, ray parameter and waveform data from 29 events occurring in the distance range 9° to 40° reveals the upper mantle structure beneath an oceanic ridge to depths of over 900 km. More than 1400 digital seismograms from earthquakes in Mexico and Central America yield 1753 travel times and 58 dT/dΔ measurements as well as high-quality, stable waveforms for investigation of the deep structure of the Gulf of California. The result of a travel time inversion with the tau method (Bessonova et al., 1976) is adjusted to fit the p(Δ) data, then further refined by incorporation of relative amplitude information through synthetic seismogram modeling. The application of a modified wave field continuation method (Clayton and McMechan, 1981) to the data with the final model confirms that GCA is consistent with the entire data set and also provides an estimate of the data resolution in velocity-depth space. We discover that the upper mantle under this spreading center has anomalously slow velocities to depths of 350 km, and place new constraints on the shape of the 660 km discontinuity.

Seismograms from 22 earthquakes along the northeast Pacific rim recorded in southern California form the data set for a comparative investigation of the upper mantle beneath the Cascade Ranges-Juan de Fuca region, an ocean-continent transit ion. These data consist of 853 seismograms (6° < Δ < 42°) which produce 1068 travel times and 40 ray parameter estimates. We use the spreading center model initially in synthetic seismogram modeling, and perturb GCA until the Cascade Ranges data are matched. Wave field continuation of both data sets with a common reference model confirms that real differences exist between the two suites of seismograms, implying lateral variation in the upper mantle. The ocean-continent transition model, CJF, features velocities from 200 and 350 km that are intermediate between GCA and T7 (Burdick and Helmberger, 1978), a model for the inland western United States. Models of continental shield regions (e.g., King and Calcagnile, 1976) have higher velocities in this depth range, but all four model types are similar below 400 km. This variation in rate of velocity increase with tectonic regime suggests an inverse relationship between velocity gradient and lithospheric age above 400 km depth.

The structure and petrology of the Johnny Lyon Hills area, Cochise County, Arizona

The Johnny Lyon Hills area is located in Cochise County in southeastern Arizona. The rocks of the area include a central core of Lower pre-Cambrian igneous and metamorphic rocks surrounded by a complexly faulted and tilted section of Upper pre-Cambrian and Paleozoic strata. Limited exposures of Mesozoic and Tertiary sedimentary and volcanic rocks are present at the north end of the map area. Late Tertiary and Quaternary alluvium almost completely surrounds and overlaps upon the older rocks.

The older pre-Cambrian rocks include a section of more than 9000 feet of generally moderately metamorphosed graywackes, slates and conglomerates of the Pinal schist injected in zones by somewhat younger rnyolite sheets. The original sediments were deposited in a geosyncline whose extent probably included large parts of Arizona, New Mexico and west Texas. During the Mazatzal Revolution the Pinal schist was deformed into northeast-trending, steeply dipping and plunging structures and the entire local section was overturned steeply toward the northwest. The pre-Cambrian Johnny Lyon granodiorite was emplaced as a large epi-tectonic pluton which modified the metamorphic character of part of the Pinal schist. Larsen method determinations indicate an age of about 715 million years for this rock, which is about the minimum age compatible with the geologic relations.

The Laramide orogeny produced numerous major thrust faults in the area involving all rocks older than and including the Lower Cretaceous Bisbee group. Major compression from the southwest and subsequent superimposed thrusting from the southeast and east are indicated. Minimum thrust displacements of more than a mile are clear and the probable displacements are of much greater magnitude. The crystalline core behaved as a single structural unit and probably caused important local divergences from the regional pattern of northeast-trending compressive forces. The massif was rotated as a unit 40 degrees or more about a northwest-trending axis overturning the pre-Cambrian fold axes in the Pinal schist.

Swarms of Late Cretaceous(?) or Early Tertiary(?) lamprophyric dikes cross the Laramide structures and are probably related to the large Texas Canyon stock several miles southeast of the map area. Intermittent high angle faulting, both older and younger than the dikes, has continued since the Laramide orogeny and has been superimposed on the older structures. This steep faulting combined with the fundamental northwesterly Laramide structural grain to produce the northwesterly trends characteristic of the mountain ridges and valleys of the area.

Experimental and theoretical developments in extended x-ray absorption fine structure (EXAFS) spectroscopy

To obtain accurate information from a structural tool it is necessary to have an understanding of the physical principles which govern the interaction between the probe and the sample under investigation. In this thesis a detailed study of the physical basis for Extended X-ray Absorption Fine Structure (EXAFS) spectroscopy is presented. A single scattering formalism of EXAFS is introduced which allows a rigorous treatment of the central atom potential. A final state interaction formalism of EXAFS is also discussed. Multiple scattering processes are shown to be significant for systems of certain geometries. The standard single scattering EXAFS analysis produces erroneous results if the data contain a large multiple scattering contribution. The effect of thermal vibrations on such multiple scattering paths is also discussed. From symmetry considerations it is shown that only certain normal modes contribute to the Debye-Waller factor for a particular scattering path. Furthermore, changes in the scattering angles induced by thermal vibrations produces additional EXAFS components called modification factors. These factors are shown to be small for most systems.

A study of the physical basis for the determination of structural information from EXAFS data is also presented. An objective method of determining the background absorption and the threshold energy is discussed and involves Gaussian functions. In addition, a scheme to determine the nature of the scattering atom in EXAFS experiments is introduced. This scheme is based on the fact that the phase intercept is a measure of the type of scattering atom. A method to determine bond distances is also discussed and does not require the use of model compounds or calculated phase shifts. The physical basis for this method is the absence of a linear term in the scattering phases. Therefore, it is possible to separate these phases from the linear term containing the distance information in the total phase.

Structure and properties of an amorphous ferromagnetic alloy

The structure and the electrical and magnetic properties of an amorphous alloy containing approximately 80 at .% iron, 13 at.% phos phorus and 7 at.% carbon (Fe_(80)Fe_(13)C_7) obtained by rapid quenching from the liquid state have been studied. Transmission electron diffraction data confirm the amorphous nature of this alloy. An analysis of the radial distribution function obtained from X-ray diffraction data indicates that the number of nearest neighbors is approximately seven, at a distance of 2.6A. The structure of the alloy can be related to that of silicate glasses and is based on a random arrangement of trigonal prisms of Fe_2P and Fe_3C types in which the iron atoms have an average ligancy of seven. Electrical resistance measurements show that the alloys are metallic. A minimum in the electrical resistivity vs. temperature curve is observed between 10° K to 50° K depending on the specimen, and the temperature at which the minimum occurs is related to the degree of local ordering. The Fe-P-C amorphous alloys are ferromagnetic. The Curie temperature measured by the induction method and by Mossbauer spectroscopy is 315° C. The field dependence of the magneto-resistance at temperatures from liquid helium to room temperature is similar to that found in crystalline iron. The ordinary Hall coefficient is approximately 10^(-11) volt-cm/amp-G. The spontaneous Hall coefficient is about 0.6 x 10^(-9) volt-cm/amp-G and is practically independent of temperature from liquid helium temperature up to 300° c.

Structure-, stereochemistry-, and metal-regulated DNA binding/cleaving molecules

In the cell, the binding of proteins to specific sequences of double helical DNA is essential for controlling the processes of protein synthesis (at the level of DNA transcription) and cell proliferation (at the level of DNA replication). In the laboratory, the sequence-specific DNA binding/cleaving properties of restriction endonuclease enzymes (secreted by microorganisms to protect them from foreign DNA molecules) have helped to fuel a revolution in molecular biology. The strength and specificity of a protein:DNA interaction depend upon structural features inherent to the protein and DNA sequences, but it is now appreciated that these features (and therefore protein:DNA complexation) may be altered (regulated) by other protein:DNA complexes, or by environmental factors such as temperature or the presence of specific organic molecules or inorganic ions. It is also now appreciated that molecules much smaller than proteins (including antibiotics of molecular weight less than 2000 and oligonucleotides) can bind to double-helical DNA in sequence-specific fashion. Elucidation of structural motifs and microscopic interactions responsible for the specific molecular recognition of DNA leads to greater understanding of natural processes and provides a basis for the design of novel sequence-specific DNA binding molecules. This thesis describes the synthesis and DNA binding/cleaving characteristics of molecules designed to probe structural, stereochemical, and environmental factors that regulate sequence-specific DNA recognition.

Chapter One introduces the DNA minor groove binding antibiotics Netropsin and Distamycin A, which are di- and tri(N-methylpyrrolecarboxamide) peptides, respectively. The method of DNA affinity cleaving, which has been employed to determine DNA binding properties of designed synthetic molecules is described. The design and synthesis of a series of Netropsin dimers linked in tail-to-tail fashion (by oxalic, malonic, succinic, or fumaric acid), or in head-to-tail fashion (by glycine, β-alanine, and γ-aminobutanoic acid (Gaba)) are presented. These Bis(Netropsin)s were appended with the iron-chelating functionality EDTA in order to make use of the technique of DNA affinity cleaving. Bis(Netropsin)-EDTA compounds are analogs of penta(N-methylpyrrolecarboxamide)-EDTA (P5E), which may be considered a head-to-tail Netropsin dimer linked by Nmethylpyrrolecarboxamide. Low- and high-resolution analysis of pBR322 DNA affinity cleaving by the iron complexes of these molecules indicated that small changes in the length and nature of the linker had significant effects on DNA binding/cleaving efficiency (a measure of DNA binding affinity). DNA binding/cleaving efficiency was found to decrease with changes in the linker in the order β-alanine > succinamide > fumaramide > N-methylpyrrolecarboxamide > malonamide >glycine, γ-aminobutanamide > oxalamide. In general, the Bis(Netropsin)-EDTA:Fe compounds retained the specificity for seven contiguous A:T base pairs characteristic of P5E:Fe binding. However, Bis(Netropsin)Oxalamide- EDTA:Fe exhibited decreased specificity for A:T base pairs, and Bis(Netropsin)-Gaba-EDT A:Fe exhibited some DNA binding sites of less than seven base pairs. Bis(Netropsin)s linked with diacids have C₂-symmmetrical DNA binding subunits and exhibited little DNA binding orientation preference. Bis(Netropsin)s linked with amino acids lack C₂-symmetrical DNA binding subunits and exhibited higher orientation preferences. A model for the high DNA binding orientation preferences observed with head-to-tail DNA minor groove binding molecules is presented.

Chapter Two describes the design, synthesis, and DNA binding properties of a series of chiral molecules: Bis(Netropsin)-EDTA compounds with linkers derived from (R,R)-, (S,S)-, and (RS,SR)-tartaric acids, (R,R)-, (S,S)-, and (RS,SR)-tartaric acid acetonides, (R)- and (S)-malic acids, N ,N-dimethylaminoaspartic acid, and (R)- and (S)-alanine, as well as three constitutional isomers in which an N-methylpyrrolecarboxamide (P1) subunit and a tri(N-methylpyrrolecarboxamide)-EDTA (P3-EDTA) subunit were linked by succinic acid, (R ,R)-, and (S ,S)-tartaric acids. DNA binding/cleaving efficiencies among this series of molecules and the Bis(Netropsin)s described in Chapter One were found to decrease with changes in the linker in the order β-alanine > succinamide > P1-succinamide-P3 > fumaramide > (S)-malicamide > N-methylpyrrolecarboxamide > (R)-malicamide > malonamide > N ,N-dimethylaminoaspanamide > glycine = Gaba = (S,S)-tartaramide = P1-(S,S)-tanaramide-P3 > oxalamide > (RS,SR)-tartaramide = P1- (R,R)-tanaramide-P3 > (R,R)-tartaramide (no sequence-specific DNA binding was detected for Bis(Netropsin)s linked by (R)- or (S)-alanine or by tartaric acid acetonides). The chiral molecules retained DNA binding specificity for seven contiguous A:T base pairs. From the DNA affinity cleaving data it could be determined that: 1) Addition of one or two substituents to the linker of Bis(Netropsin)-Succinamide resulted in stepwise decreases in DNA binding affinity; 2) molecules with single hydroxyl substituents bound DNA more strongly than molecules with single dimethylamino substituents; 3) hydroxyl-substituted molecules of (S) configuration bound more strongly to DNA than molecules of (R) configuration. This stereochemical regulation of DNA binding is proposed to arise from the inherent right-handed twist of (S)-enantiomeric Bis(Netropsin)s versus the inherent lefthanded twist of (R)-enantiomeric Bis(Netropsin)s, which makes the (S)-enantiomers more complementary to the right-handed twist of B form DNA.

Chapter Three describes the design and synthesis of molecules for the study of metalloregulated DNA binding phenomena. Among a series of Bis(Netropsin)-EDTA compounds linked by homologous tethers bearing four, five, or six oxygen atoms, the Bis(Netropsin) linked by a pentaether tether exhibited strongly enhanced DNA binding/cleaving in the presence of strontium or barium cations. The observed metallospecificity was consistent with the known affinities of metal cations for the cyclic hexaether 18-crown-6 in water. High-resolution DNA affinity cleaving analysis indicated that DNA binding by this molecule in the presence of strontium or barium was not only stronger but of different sequence-specificity than the (weak) binding observed in the absence of metal cations. The metalloregulated binding sites were consistent with A:T binding by the Netropsin subunits and G:C binding by a strontium or barium:pentaether complex. A model for the observed positive metalloregulation and novel sequence-specificity is presented. The effects of 44 different cations on DNA affinity cleaving by P5E:Fe were examined. A series of Bis(Netropsin)-EDTA compounds linked by tethers bearing two, three, four, or five amino groups was also synthesized. These molecules exhibited strong and specific binding to A:T rich regions of DNA. It was found that the iron complexes of these molecules bound and cleaved DNA most efficiently at pH 6.0-6.5, while P5E:Fe bound and cleaved most efficiently at pH 7.5-8.0. Incubating the Bis(Netropsin) Polyamine-EDTA:Fe molecules with K₂PdCl₄ abolished their DNA binding/cleaving activity. It is proposed that the observed negative metalloregulation arises from kinetically inert Bis(Netropsin) Polyamine:Pd(II) complexes or aggregates, which are sterically unsuitable for DNA complexation. Finally, attempts to produce a synthetic metalloregulated DNA binding protein are described. For this study, five derivatives of a synthetic 52 amino acid residue DNA binding/cleaving protein were produced. The synthetic mutant proteins carried a novel pentaether ionophoric amino acid residue at different positions within the primary sequence. The proteins did not exhibit significant DNA binding/cleaving activity, but they served to illustrate the potential for introducing novel amino acid residues within DNA binding protein sequences, and for the development of the tricyclohexyl ester of EDTA as a superior reagent for the introduction of EDT A into synthetic proteins.

Chapter Four describes the discovery and characterization of a new DNA binding/cleaving agent, [SalenMn(III)]OAc. This metal complex produces single- and double-strand cleavage of DNA, with specificity for A:T rich regions, in the presence of oxygen atom donors such as iodosyl benzene, hydrogen peroxide, or peracids. Maximal cleavage by [SalenMn(III)]OAc was produced at pH 6-7. A comparison of DNA singleand double-strand cleavage by [SalenMn(III)]⁺ and other small molecules (Methidiumpropyl-EDTA:Fe, Distamycin-EDTA:Fe, Neocarzinostatin, Bleomycin:Fe) is presented. It was found that DNA cleavage by [SalenMn(III)]⁺ did not require the presence of dioxygen, and that base treatment of DNA subsequent to cleavage by [SalenMn(III)]⁺ afforded greater cleavage and alterations in the cleavage patterns. Analysis of DNA products formed upon DNA cleavage by [SalenMn(III)] indicated that cleavage was due to oxidation of the sugar-phosphate backbone of DNA. Several mechanisms consistent with the observed products and reaction requirements are discussed.

Chapter Five describes progress on some additional studies. In one study, the DNA binding/cleaving specificities of Distamycin-EDTA derivatives bearing pyrrole N-isopropyl substituents were found to be the same as those of derivatives bearing pyrrole N-methyl substituents. In a second study, the design of and synthetic progress towards a series of nucleopeptide activators of transcription are presented. Five synthetic plasmids designed to test for activation of in vitro run-off transcription by DNA triple helix-forming oligonucleotides or nucleopeptides are described.

Chapter Six contains the experimental documentation of the thesis work.

Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.

Full and Model-Reduced Structure-Preserving Simulation of Incompressible Fluids

This thesis outlines the construction of several types of structured integrators for incompressible fluids. We first present a vorticity integrator, which is the Hamiltonian counterpart of the existing Lagrangian-based fluid integrator. We next present a model-reduced variational Eulerian integrator for incompressible fluids, which combines the efficiency gains of dimension reduction, the qualitative robustness to coarse spatial and temporal resolutions of geometric integrators, and the simplicity of homogenized boundary conditions on regular grids to deal with arbitrarily-shaped domains with sub-grid accuracy.

Both these numerical methods involve approximating the Lie group of volume-preserving diffeomorphisms by a finite-dimensional Lie-group and then restricting the resulting variational principle by means of a non-holonomic constraint. Advantages and limitations of this discretization method will be outlined. It will be seen that these derivation techniques are unable to yield symplectic integrators, but that energy conservation is easily obtained, as is a discretized version of Kelvin's circulation theorem.

Finally, we outline the basis of a spectral discrete exterior calculus, which may be a useful element in producing structured numerical methods for fluids in the future.

Imidazolium-Based Organic Structure Directing Agents for the Synthesis of Microporous Materials

The central theme of this thesis is the use of imidazolium-based organic structure directing agents (OSDAs) in microporous materials synthesis. Imidazoliums are advantageous OSDAs as they are relatively inexpensive and simple to prepare, show robust stability under microporous material synthesis conditions, have led to a wide range of products, and have many permutations in structure that can be explored. The work I present involves the use of mono-, di-, and triquaternary imidazolium-based OSDAs in a wide variety of microporous material syntheses. Much of this work was motivated by successful computational predictions (Chapter 2) that led me to continue to explore these types of OSDAs. Some of the important discoveries with these OSDAs include the following: 1) Experimental evaluation and confirmation of a computational method that predicted a new OSDA for pure-silica STW, a desired framework containing helical pores that was previously very difficult to synthesize. 2) Discovery of a number of new imidazolium OSDAs to synthesize zeolite RTH, a zeolite desired for both the methanol-to-olefins reaction as well as NOX reduction in exhaust gases. This discovery enables the use of RTH for many additional investigations as the previous OSDA used to make this framework was difficult to synthesize, such that no large scale preparations would be practical. 3) The synthesis of pure-silica RTH by topotactic condensation from a layered precursor (denoted CIT-10), that can also be pillared to make a new framework material with an expanded pore system, denoted CIT-11, that can be calcined to form a new microporous material, denoted CIT-12. CIT-10 is also interesting since it is the first layered material to contain 8 membered rings through the layers, making it potentially useful in separations if delamination methods can be developed. 4) The synthesis of a new microporous material, denoted CIT-7 (framework code CSV) that contains a 2-dimensional system of 8 and 10 membered rings with a large cage at channel intersections. This material is especially important since it can be synthesized as a pure-silica framework under low-water, fluoride-mediated synthesis conditions, and as an aluminosilicate material under hydroxide mediated conditions. 5) The synthesis of high-silica heulandite (HEU) by topotactic condensation as well as direct synthesis, demonstrating new, more hydrothermally stable compositions of a previously known framework. 6) The synthesis of germanosilicate and aluminophosphate LTA using a triquaternary OSDA. All of these materials show the diverse range of products that can be formed from OSDAs that can be prepared by straightforward syntheses and have made many of these materials accessible for the first time under facile zeolite synthesis conditions.

A fully-nonlocal energy-based formulation and high-performance realization of the quasicontinuum method

The quasicontinuum (QC) method was introduced to coarse-grain crystalline atomic ensembles in order to bridge the scales from individual atoms to the micro- and mesoscales. Though many QC formulations have been proposed with varying characteristics and capabilities, a crucial cornerstone of all QC techniques is the concept of summation rules, which attempt to efficiently approximate the total Hamiltonian of a crystalline atomic ensemble by a weighted sum over a small subset of atoms. In this work we propose a novel, fully-nonlocal, energy-based formulation of the QC method with support for legacy and new summation rules through a general energy-sampling scheme. Our formulation does not conceptually differentiate between atomistic and coarse-grained regions and thus allows for seamless bridging without domain-coupling interfaces. Within this structure, we introduce a new class of summation rules which leverage the affine kinematics of this QC formulation to most accurately integrate thermodynamic quantities of interest. By comparing this new class of summation rules to commonly-employed rules through analysis of energy and spurious force errors, we find that the new rules produce no residual or spurious force artifacts in the large-element limit under arbitrary affine deformation, while allowing us to seamlessly bridge to full atomistics. We verify that the new summation rules exhibit significantly smaller force artifacts and energy approximation errors than all comparable previous summation rules through a comprehensive suite of examples with spatially non-uniform QC discretizations in two and three dimensions. Due to the unique structure of these summation rules, we also use the new formulation to study scenarios with large regions of free surface, a class of problems previously out of reach of the QC method. Lastly, we present the key components of a high-performance, distributed-memory realization of the new method, including a novel algorithm for supporting unparalleled levels of deformation. Overall, this new formulation and implementation allows us to efficiently perform simulations containing an unprecedented number of degrees of freedom with low approximation error.

I. The crystal structure of a new dimer of triphenylfluorocyclobutadiene. II. A low temperature refinement of the cyanuric triazide structure

I. THE CRYSTAL STRUCTURE OF A NEW DIMER OF TRIPHENYLFLUOROCYCLOBUTADIENE

The crystal structure of thermal isomer of the “head-to-head” dimer of triphenylfluorocyclobutadiene was determined by the direct method. The Σ₂ relationship involving the low angle reflections with the largest E’s were found and solved for the signs by the symbolic method of Zachariasen. The structure was seen in the electron density map and the E-map, and was refined antisotropically by the method of least squares. The residual R was 0.065.

The structure is a gem-difluorohexaphenyldihydropentalene. All of the phenyl groups are planar as it is the cyclopentadiene ring of the dihydropentalene skeleton. Overcrowding at the position of the flourines causes some deviations from the normal bond angles in the cyclopentene ring.

The list of observed and calculated structure factors on pages 32-34 will not be legible on the microfilm. Photographic copies may be obtained from the California Institute of Technology.

II. A LOW TEMPERATURE REFINEMENT OF THE CYANURIC TRIAZIDE STRUCTURE

The structure of cyanuric triazide was refined anisotropically by the method of least squares. Three-dimensional intensity data, which has been collected photographically with MoK_α radiation at -110˚C, were used in the refinement. The residual R was reduced to 0.081.

The structure is completely planar, and there is no significant bond alternation in the cyanuric ring. The packing of the molecules causes the azide groups to deviate from linearity by 8 degrees.

Regional variations in upper mantle structure beneath North America

Several types of seismological data, including surface wave group and phase velocities, travel times from large explosions, and teleseismic travel time anomalies, have indicated that there are significant regional variations in the upper few hundred kilometers of the mantle beneath continental areas. Body wave travel times and amplitudes from large chemical and nuclear explosions are used in this study to delineate the details of these variations beneath North America.

As a preliminary step in this study, theoretical P wave travel times, apparent velocities, and amplitudes have been calculated for a number of proposed upper mantle models, those of Gutenberg, Jeffreys, Lehman, and Lukk and Nersesov. These quantities have been calculated for both P and S waves for model CIT11GB, which is derived from surface wave dispersion data. First arrival times for all the models except that of Lukk and Nersesov are in close agreement, but the travel time curves for later arrivals are both qualitatively and quantitatively very different. For model CIT11GB, there are two large, overlapping regions of triplication of the travel time curve, produced by regions of rapid velocity increase near depths of 400 and 600 km. Throughout the distance range from 10 to 40 degrees, the later arrivals produced by these discontinuities have larger amplitudes than the first arrivals. The amplitudes of body waves, in fact, are extremely sensitive to small variations in the velocity structure, and provide a powerful tool for studying structural details.

Most of eastern North America, including the Canadian Shield has a Pn velocity of about 8.1 km/sec, with a nearly abrupt increase in compressional velocity by ~ 0.3 km/sec near at a depth varying regionally between 60 and 90 km. Variations in the structure of this part of the mantle are significant even within the Canadian Shield. The low-velocity zone is a minor feature in eastern North America and is subject to pronounced regional variations. It is 30 to 50 km thick, and occurs somewhere in the depth range from 80 to 160 km. The velocity decrease is less than 0.2 km/sec.

Consideration of the absolute amplitudes indicates that the attenuation due to anelasticity is negligible for 2 hz waves in the upper 200 km along the southeastern and southwestern margins of the Canadian Shield. For compressional waves the average Q for this region is > 3000. The amplitudes also indicate that the velocity gradient is at least 2 x 10^-3 both above and below the low-velocity zone, implying that the temperature gradient is < 4.8°C/km if the regions are chemically homogeneous.

In western North America, the low-velocity zone is a pronounced feature, extending to the base of the crust and having minimum velocities of 7.7 to 7.8 km/sec. Beneath the Colorado Plateau and Southern Rocky Mountains provinces, there is a rapid velocity increase of about 0.3 km/sec, similar to that observed in eastern North America, but near a depth of 100 km.

Complicated travel time curves observed on profiles with stations in both eastern and western North America can be explained in detail by a model taking into account the lateral variations in the structure of the low-velocity zone. These variations involve primarily the velocity within the zone and the depth to the top of the zone; the depth to the bottom is, for both regions, between 140 and 160 km.

The depth to the transition zone near 400 km also varies regionally, by about 30-40 km. These differences imply variations of 250 °C in the temperature or 6 % in the iron content of the mantle, if the phase transformation of olivine to the spinel structure is assumed responsible. The structural variations at this depth are not correlated with those at shallower depths, and follow no obvious simple pattern.

The computer programs used in this study are described in the Appendices. The program TTINV (Appendix IV) fits spherically symmetric earth models to observed travel time data. The method, described in Appendix III, resembles conventional least-square fitting, using partial derivatives of the travel time with respect to the model parameters to perturb an initial model. The usual ill-conditioned nature of least-squares techniques is avoided by a technique which minimizes both the travel time residuals and the model perturbations.

Spherically symmetric earth models, however, have been found inadequate to explain most of the observed travel times in this study. TVT4, a computer program that performs ray theory calculations for a laterally inhomogeneous earth model, is described in Appendix II. Appendix I gives a derivation of seismic ray theory for an arbitrarily inhomogeneous earth model.

I. The crystal structure of potassium bis(tricyanovinyl)amine. II. Nuclear magnetic resonance studies of alkali halide solutions. III. Elucidation of the intramolecular hydrogen bond in acetylacetone by the deuterium isotope effect on the chemical shift of the bridge hydrogen

Part I

Potassium bis-(tricyanovinyl) amine, K⁺N[C(CN)=C(CN)₂]₂^-, crystallizes in the monoclinic system with the space group Cc and lattice constants, a = 13.346 ± 0.003 Å, c = 8.992 ± 0.003 Å, B = 114.42 ± 0.02°, and Z = 4. Three dimensional intensity data were collected by layers perpendicular to b* and c* axes. The crystal structure was refined by the least squares method with anisotropic temperature factor to an R value of 0.064.

The average carbon-carbon and carbon-nitrogen bond distances in –C-CΞN are 1.441 ± 0.016 Å and 1.146 ± 0.014 Å respectively. The bis-(tricyanovinyl) amine anion is approximately planar. The coordination number of the potassium ion is eight with bond distances from 2.890 Å to 3.408 Å. The bond angle C-N-C of the amine nitrogen is 132.4 ± 1.9°. Among six cyano groups in the molecule, two of them are bent by what appear to be significant amounts (5.0° and 7.2°). The remaining four are linear within the experimental error. The bending can probably be explained by molecular packing forces in the crystals.

Part II

The nuclear magnetic resonance of ⁸¹Br and ¹²⁷I in aqueous solutions were studied. The cation-halide ion interactions were studied by studying the effect of the Li⁺, Na⁺, K⁺, Mg⁺⁺, Cs⁺ upon the line width of the halide ions. The solvent-halide ion interactions were studied by studying the effects of methanol, acetonitrile, and acetone upon the line width of ⁸¹Br and ¹²⁷I in the aqueous solutions. It was found that the viscosity plays a very important role upon the halide ions line width. There is no specific cation-halide ion interaction for those ions such as Mg⁺⁺, Di⁺, Na⁺, and K⁺, whereas the Cs⁺ - halide ion interaction is strong. The effect of organic solvents upon the halide ion line width in aqueous solutions is in the order acetone ˃ acetonitrile ˃ methanol. It is suggested that halide ions do form some stable complex with the solvent molecules and the reason Cs⁺ can replace one of the ligands in the solvent-halide ion complex.

Part III

An unusually large isotope effect on the bridge hydrogen chemical shift of the enol form of pentanedione-2, 4(acetylacetone) and 3-methylpentanedione-2, 4 has been observed. An attempt has been made to interpret this effect. It is suggested from the deuterium isotope effect studies, temperature dependence of the bridge hydrogen chemical shift studies, IR studies in the OH, OD, and C=O stretch regions, and the HMO calculations, that there may probably be two structures for the enol form of acetylacetone. The difference between these two structures arises mainly from the electronic structure of the π-system. The relative population of these two structures at various temperatures for normal acetylacetone and at room temperature for the deuterated acetylacetone were calculated.

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.