Black holes in the early universe, in compact binaries, and as energy sources inside solar-type stars

This thesis consists of three separate studies of roles that black holes might play in our universe.

In the first part we formulate a statistical method for inferring the cosmological parameters of our universe from LIGO/VIRGO measurements of the gravitational waves produced by coalescing black-hole/neutron-star binaries. This method is based on the cosmological distance-redshift relation, with "luminosity distances" determined directly, and redshifts indirectly, from the gravitational waveforms. Using the current estimates of binary coalescence rates and projected "advanced" LIGO noise spectra, we conclude that by our method the Hubble constant should be measurable to within an error of a few percent. The errors for the mean density of the universe and the cosmological constant will depend strongly on the size of the universe, varying from about 10% for a "small" universe up to and beyond 100% for a "large" universe. We further study the effects of random gravitational lensing and find that it may strongly impair the determination of the cosmological constant.

In the second part of this thesis we disprove a conjecture that black holes cannot form in an early, inflationary era of our universe, because of a quantum-field-theory induced instability of the black-hole horizon. This instability was supposed to arise from the difference in temperatures of any black-hole horizon and the inflationary cosmological horizon; it was thought that this temperature difference would make every quantum state that is regular at the cosmological horizon be singular at the black-hole horizon. We disprove this conjecture by explicitly constructing a quantum vacuum state that is everywhere regular for a massless scalar field. We further show that this quantum state has all the nice thermal properties that one has come to expect of "good" vacuum states, both at the black-hole horizon and at the cosmological horizon.

In the third part of the thesis we study the evolution and implications of a hypothetical primordial black hole that might have found its way into the center of the Sun or any other solar-type star. As a foundation for our analysis, we generalize the mixing-length theory of convection to an optically thick, spherically symmetric accretion flow (and find in passing that the radial stretching of the inflowing fluid elements leads to a modification of the standard Schwarzschild criterion for convection). When the accretion is that of solar matter onto the primordial hole, the rotation of the Sun causes centrifugal hangup of the inflow near the hole, resulting in an "accretion torus" which produces an enhanced outflow of heat. We find, however, that the turbulent viscosity, which accompanies the convective transport of this heat, extracts angular momentum from the inflowing gas, thereby buffering the torus into a lower luminosity than one might have expected. As a result, the solar surface will not be influenced noticeably by the torus's luminosity until at most three days before the Sun is finally devoured by the black hole. As a simple consequence, accretion onto a black hole inside the Sun cannot be an answer to the solar neutrino puzzle.

On reconstruction theorems in noncommutative Riemannian geometry

We present a novel account of the theory of commutative spectral triples and their two closest noncommutative generalisations, almost-commutative spectral triples and toric noncommutative manifolds, with a focus on reconstruction theorems, viz, abstract, functional-analytic characterisations of global-analytically defined classes of spectral triples. We begin by reinterpreting Connes's reconstruction theorem for commutative spectral triples as a complete noncommutative-geometric characterisation of Dirac-type operators on compact oriented Riemannian manifolds, and in the process clarify folklore concerning stability of properties of spectral triples under suitable perturbation of the Dirac operator. Next, we apply this reinterpretation of the commutative reconstruction theorem to obtain a reconstruction theorem for almost-commutative spectral triples. In particular, we propose a revised, manifestly global-analytic definition of almost-commutative spectral triple, and, as an application of this global-analytic perspective, obtain a general result relating the spectral action on the total space of a finite normal compact oriented Riemannian cover to that on the base space. Throughout, we discuss the relevant refinements of these definitions and results to the case of real commutative and almost-commutative spectral triples. Finally, we outline progess towards a reconstruction theorem for toric noncommutative manifolds.

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

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

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

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

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

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

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

On erasure coding for distributed storage and streaming communications

The work presented in this thesis revolves around erasure correction coding, as applied to distributed data storage and real-time streaming communications.

First, we examine the problem of allocating a given storage budget over a set of nodes for maximum reliability. The objective is to find an allocation of the budget that maximizes the probability of successful recovery by a data collector accessing a random subset of the nodes. This optimization problem is challenging in general because of its combinatorial nature, despite its simple formulation. We study several variations of the problem, assuming different allocation models and access models, and determine the optimal allocation and the optimal symmetric allocation (in which all nonempty nodes store the same amount of data) for a variety of cases. Although the optimal allocation can have nonintuitive structure and can be difficult to find in general, our results suggest that, as a simple heuristic, reliable storage can be achieved by spreading the budget maximally over all nodes when the budget is large, and spreading it minimally over a few nodes when it is small. Coding would therefore be beneficial in the former case, while uncoded replication would suffice in the latter case.

Second, we study how distributed storage allocations affect the recovery delay in a mobile setting. Specifically, two recovery delay optimization problems are considered for a network of mobile storage nodes: the maximization of the probability of successful recovery by a given deadline, and the minimization of the expected recovery delay. We show that the first problem is closely related to the earlier allocation problem, and solve the second problem completely for the case of symmetric allocations. It turns out that the optimal allocations for the two problems can be quite different. In a simulation study, we evaluated the performance of a simple data dissemination and storage protocol for mobile delay-tolerant networks, and observed that the choice of allocation can have a significant impact on the recovery delay under a variety of scenarios.

Third, we consider a real-time streaming system where messages created at regular time intervals at a source are encoded for transmission to a receiver over a packet erasure link; the receiver must subsequently decode each message within a given delay from its creation time. For erasure models containing a limited number of erasures per coding window, per sliding window, and containing erasure bursts whose maximum length is sufficiently short or long, we show that a time-invariant intrasession code asymptotically achieves the maximum message size among all codes that allow decoding under all admissible erasure patterns. For the bursty erasure model, we also show that diagonally interleaved codes derived from specific systematic block codes are asymptotically optimal over all codes in certain cases. We also study an i.i.d. erasure model in which each transmitted packet is erased independently with the same probability; the objective is to maximize the decoding probability for a given message size. We derive an upper bound on the decoding probability for any time-invariant code, and show that the gap between this bound and the performance of a family of time-invariant intrasession codes is small when the message size and packet erasure probability are small. In a simulation study, these codes performed well against a family of random time-invariant convolutional codes under a number of scenarios.

Finally, we consider the joint problems of routing and caching for named data networking. We propose a backpressure-based policy that employs virtual interest packets to make routing and caching decisions. In a packet-level simulation, the proposed policy outperformed a basic protocol that combines shortest-path routing with least-recently-used (LRU) cache replacement.

Algorithmic challenges in green data centers

With data centers being the supporting infrastructure for a wide range of IT services, their efficiency has become a big concern to operators, as well as to society, for both economic and environmental reasons. The goal of this thesis is to design energy-efficient algorithms that reduce energy cost while minimizing compromise to service. We focus on the algorithmic challenges at different levels of energy optimization across the data center stack. The algorithmic challenge at the device level is to improve the energy efficiency of a single computational device via techniques such as job scheduling and speed scaling. We analyze the common speed scaling algorithms in both the worst-case model and stochastic model to answer some fundamental issues in the design of speed scaling algorithms. The algorithmic challenge at the local data center level is to dynamically allocate resources (e.g., servers) and to dispatch the workload in a data center. We develop an online algorithm to make a data center more power-proportional by dynamically adapting the number of active servers. The algorithmic challenge at the global data center level is to dispatch the workload across multiple data centers, considering the geographical diversity of electricity price, availability of renewable energy, and network propagation delay. We propose algorithms to jointly optimize routing and provisioning in an online manner. Motivated by the above online decision problems, we move on to study a general class of online problem named "smoothed online convex optimization", which seeks to minimize the sum of a sequence of convex functions when "smooth" solutions are preferred. This model allows us to bridge different research communities and help us get a more fundamental understanding of general online decision problems.

Catalytic conversion of nitrogen to ammonia by an iron model complex

Threefold symmetric Fe phosphine complexes have been used to model the structural and functional aspects of biological N₂ fixation by nitrogenases. Low-valent bridging Fe-S-Fe complexes in the formal oxidation states Fe(II)Fe(II), Fe(II)/Fe(I), and Fe(I)/Fe(I) have been synthesized which display rich spectroscopic and magnetic behavior. A series of cationic tris-phosphine borane (TPB) ligated Fe complexes have been synthesized and been shown to bind a variety of nitrogenous ligands including N₂H₄, NH₃, and NH₂-. These complexes are all high spin S = 3/2 and display EPR and magnetic characteristics typical of this spin state. Furthermore, a sequential protonation and reduction sequence of a terminal amide results in loss of NH₃ and uptake of N₂. These stoichiometric transformations represent the final steps in potential N₂ fixation schemes.

Treatment of an anionic FeN₂ complex with excess acid also results in the formation of some NH₃, suggesting the possibility of a catalytic cycle for the conversion of N₂ to NH₃ mediated by Fe. Indeed, use of excess acid and reductant results in the formation of seven equivalents of NH₃ per Fe center, demonstrating Fe mediated catalytic N₂ fixation with acids and protons for the first time. Numerous control experiments indicate that this catalysis is likely being mediated by a molecular species.

A number of other phosphine ligated Fe complexes have also been tested for catalysis and suggest that a hemi-labile Fe-B interaction may be critical for catalysis. Additionally, various conditions for the catalysis have been investigated. These studies further support the assignment of a molecular species and delineate some of the conditions required for catalysis.

Finally, combined spectroscopic studies have been performed on a putative intermediate for catalysis. These studies converge on an assignment of this new species as a hydrazido(2-) complex. Such species have been known on group 6 metals for some time, but this represents the first characterization of this ligand on Fe. Further spectroscopic studies suggest that this species is present in catalytic mixtures, which suggests that the first steps of a distal mechanism for N₂ fixation are feasible in this system.

Singularities and horizons in the collisions of gravitational waves

This thesis presents a study of the dynamical, nonlinear interaction of colliding gravitational waves, as described by classical general relativity. It is focused mainly on two fundamental questions: First, what is the general structure of the singularities and Killing-Cauchy horizons produced in the collisions of exactly plane-symmetric gravitational waves? Second, under what conditions will the collisions of almost-plane gravitational waves (waves with large but finite transverse sizes) produce singularities?

In the work on the collisions of exactly-plane waves, it is shown that Killing horizons in any plane-symmetric spacetime are unstable against small plane-symmetric perturbations. It is thus concluded that the Killing-Cauchy horizons produced by the collisions of some exactly plane gravitational waves are nongeneric, and that generic initial data for the colliding plane waves always produce "pure" spacetime singularities without such horizons. This conclusion is later proved rigorously (using the full nonlinear theory rather than perturbation theory), in connection with an analysis of the asymptotic singularity structure of a general colliding plane-wave spacetime. This analysis also proves that asymptotically the singularities created by colliding plane waves are of inhomogeneous-Kasner type; the asymptotic Kasner axes and exponents of these singularities in general depend on the spatial coordinate that runs tangentially to the singularity in the non-plane-symmetric direction.

In the work on collisions of almost-plane gravitational waves, first some general properties of single almost-plane gravitational-wave spacetimes are explored. It is shown that, by contrast with an exact plane wave, an almost-plane gravitational wave cannot have a propagation direction that is Killing; i.e., it must diffract and disperse as it propagates. It is also shown that an almost-plane wave cannot be precisely sandwiched between two null wavefronts; i.e., it must leave behind tails in the spacetime region through which it passes. Next, the occurrence of spacetime singularities in the collisions of almost-plane waves is investigated. It is proved that if two colliding, almost-plane gravitational waves are initially exactly plane-symmetric across a central region of sufficiently large but finite transverse dimensions, then their collision produces a spacetime singularity with the same local structure as in the exact-plane-wave collision. Finally, it is shown that a singularity still forms when the central regions are only approximately plane-symmetric initially. Stated more precisely, it is proved that if the colliding almost-plane waves are initially sufficiently close to being exactly plane-symmetric across a bounded central region of sufficiently large transverse dimensions, then their collision necessarily produces spacetime singularities. In this case, nothing is now known about the local and global structures of the singularities.

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.

^(18)O/^(16)O and D/H studies on the interactions between heated meteoric ground waters and igneous intrusions: Western San Juan Mountains, Colorado and the Isle of Skye, Scotland

Detailed oxygen, hydrogen and carbon isotope studies have been carried out on igneous and metamorphic rocks of the Stony Mountain complex, Colorado, and the Isle of Skye, Scotland, in order to better understand the problems of hydrothermal meteoric water-rock interaction.

The Tertiary Stony Mountain stock (~1.3 km in diameter), is composed of an outer diorite, a main mass of biotite gabbro, and an inner diorite. The entire complex and most of the surrounding country rocks have experienced various degrees of ¹⁸O depletion (up to 10 per mil) due to interaction with heated meteoric waters. The inner diorite apparently formed from a low-¹⁸O magma with δ¹⁸O ≃ +2.5, but most of the isotopic effects are a result of exchange between H₂O and solidified igneous rocks. The low-¹⁸O inner diorite magma was probably produced by massive assimilation and/or melting of hydrothermally altered country rocks. The δ¹⁸O values of the rocks generally increase with increasing grain size, except that quartz typically has δ¹⁸O = +6 to +8, and is more resistant to hydrothermal exchange than any other mineral studied. Based on atom % oxygen, the outer diorites, gabbros, and volcanic rocks exhibit integrated water/rock ratios of 0.3 ± 0.2, 0.15 ± 0.1, and 0.2 ± 0.1, respectively. Locally, water/rock ratios attain values greater than 1.0. Hydrogen isotopic analyses of sericites, chlorites, biotites, and amphiboles range from -117 to -150. δD in biotites varies inversely with Fe/Fe+Mg, as predicted by Suzuoki and Epstein (1974), and positively with elevation, over a range of 600 m. The calculated δD of the mid-to-late-Tertiary meteoric waters is about -100. Carbonate δ¹³C values average -5.5 (PDB), within the generally accepted range for deep-seated carbon.

Almost all the rocks within 4 km of the central Tertiary intrusive complex of Skye are depleted in ¹⁸O. Whole-rock δ¹⁸O values of basalts (-7. 1 to +8.4), Mesozoic shales (-0.6 to + 12.4), and Precambrian sandstones (-6.2 to + 10.8) systematically decrease inward towards the center of the complex. The Cuillin gabbro may have formed from a ¹⁸O-depleted magma (depleted by about 2 per mil); δ¹⁸O of plagioclase (-7.1 to + 2.5) and pyroxene (-0.5 to + 3.2) decrease outward toward the margins of the pluton. The Red Hills epigranite plutons have δ¹⁸O quartz (-2.7 to + 7.6) and feldspar (-6.7 to + 6.0) that suggest about 3/4 of the exchange took place at subsolidus temperatures; profound disequilibrium quartz-feldspar fractionations (up to 12) are characteristic. The early epigranites were intruded as low-¹⁸O melts (depletions of up to 3 per mil) with δ¹⁸O of the primary, igneous quartz decreasing progressively with time. The Southern Porphyritic Epigranite was apparently intruded as a low-¹⁸O magma with δ¹⁸O ≃ -2.6. A good correlation exists between grain size and δ¹⁸O for the unique, high-¹⁸O Beinn an Dubhaich granite which intrudes limestone having a δ¹⁸O range of +0.5 to +20.8, and δ¹³C of -4.9 to -1.0. The δD values of sericites (-104 to -107), and amphiboles, chlorites, and biotites (-105 to -128) from the igneous rocks , indicate that Eocene surface waters at Skye had δD ≃ -90. The average water/rock ratio for the Skye hydrothermal system is approximately one; at least 2000 km³ of heated meteoric waters were cycled through these rocks.

Thus these detailed isotopic studies of two widely separated areas indicate that (1) ¹⁸O-depleted magmas are commonly produced in volcanic terranes invaded by epizonal intrusions; (2) most of the ¹⁸O-depletion in such areas are a result of subsolidus exchange (particularly of feldspars); however correlation of δ¹⁸O with grain size is generally preserved only for systems that have undergone relatively minor meteoric hydrothermal exchange; (3) feldspar and calcite are the minerals mos t susceptible to oxygen isotopic exchange, whereas quartz is very resistant to oxygen isotope exchange; biotite, magnetite, and pyroxene have intermediate susceptibilities; and (4) basaltic country rocks are much more permeable to the hydrothermal convective system than shale, sandstone, or the crystalline basement complex.

Theory of ozone isotopic effects and various electron transfer reactions

Three separate topics, each stimulated by experiments, are treated theoretically in this dessertation: isotopic effects of ozone, electron transfer at interfaces, and intramolecular directional electron transfer in a supramolecular system.

The strange mass-independent isotope effect for the enrichment of ozone, which has been a puzzle in the literature for some 20 years, and the equally puzzling unconventional strong mass-dependent effect of individual reaction rate constants are studied as different aspects of a symmetry-driven behavior. A statistical (RRKM-based) theory with a hindered-rotor transition state is used. The individual rate constant ratios of recombination reactions at low pressures are calculated using the theory involving (1) small deviation from the statistical density of states for symmetric isotopomers, and (2) weak collisions for deactivation of the vibrationally excited ozone molecules. The weak collision and partitioning among exit channels play major roles in producing the large unconventional isotope effect in "unscrambled" systems. The enrichment studies reflect instead the non-statistical effect in "scrambled" systems. The theoretical results of low-pressure ozone enrichments and individual rate constant ratios obtained from these calculations are consistent with the corresponding experimental results. The isotopic exchange rate constant for the reaction ^(16)O + ^(18)O ^(18)O→+ ^(16)O ^(18)O + ^(18)O provides information on the nature of a variationally determined hindered-rotor transition state using experimental data at 130 K and 300 K. Pressure effects on the recombination rate constant, on the individual rate constant ratios and on the enrichments are also investigated. The theoretical results are consistent with the experimental data. The temperature dependence of the enrichment and rate constant ratios is also discussed, and experimental tests are suggested. The desirability of a more accurate potential energy surface for ozone in the transition state region is also noted.

Electron transfer reactions at semiconductor /liquid interfaces are studied using a tight-binding model for the semiconductors. The slab method and a z-transform method are employed in obtaining the tight-binding electronic structures of semiconductors having surfaces. The maximum electron transfer rate constants at Si/viologen^(2-/+) and InP /Me_(2)Fc^(+/O) interfaces are computed using the tight-binding type calculations for the solid and the extended-Huckel for the coupling to the redox agent at the interface. These electron transfer reactions are also studied using a free electron model for the semiconductor and the redox molecule, where Bardeen's method is adapted to calculate the coupling matrix element between the molecular and semiconductor electronic states. The calculated results for maximum rate constant of the electron transfer from the semiconductor bulk states are compared with the experimentally measured values of Lewis and coworkers, and are in reasonable agreement, without adjusting parameters. In the case of InP /liquid interface, the unusual current vs applied potential behavior is additionally interpreted, in part, by the presence of surface states.

Photoinduced electron transfer reactions in small supramolecular systems, such as 4-aminonaphthalimide compounds, are interesting in that there are, in principle, two alternative pathways (directions) for the electron transfer. The electron transfer, however, is unidirectional, as deduced from pH-dependent fluorescence quenching studies on different compounds. The role of electronic coupling matrix element and the charges in protonation are considered to explain the directionality of the electron transfer and other various results. A related mechanism is proposed to interpret the fluorescence behavior of similar molecules as fluorescent sensors of metal ions.

Studies directed toward the synthesis of palau'amine and axinellamines A-D

The use of spiro [2.4]hepta-4,6-diene-1-methanol 7 as a general precursor for the synthesis of highly functionalized cyclopentyl rings is described. Diene 7 was converted to its silyl protected 4-nitrile derivative 24 in 46% overall yield. The cyclopropyl ring of 24 reacted with soft carbanionic nucleophiles to give ring opened homo-conjugate addition products 25a-h in 76-97% yield without loss of optical purity. The addition products could be further manipulated by selective mono-hydrogenation to give 1,2 substituted cyclopentenes 26a-e in 85-96% yield.

Diene 7 was used as a starting material for studies directed toward the synthesis of the stereochemically dense chloro-cyclopentyl core of palau'amine 1. Two advanced intermediates 50 and 72 were synthesized. Attempts to effect intramolecular chlorine transfer with 50 were unsuccessful. Attempted intramolecular chlorine transfer with 72 led, instead, to an oxygenated species resulting from oxygen radical trapping.

The enantioselective synthesis of the stereochemically dense chloro-cyclopenty l core of axinellamines A-D 2-5 starting from 7 is also described. The core is synthesized in 4.6% yield over 24 steps. Nakamura's radical dehalogenative hydroxylation is applied for the first time to a cyclopropyl carbonyl iodide to give the ring-opened product in 86% yield. Bolm's meso-anhydride desymmetrization is used to introduce asymmetry in a norbornene intermediate. The final step is a diastereoselective intermolecular chlorination using Barton's methodology to achieve chlorine transfer in 76% yield.

Investigations of the origin of stereocontrol in syndiospecific Ziegler-Natta polymerizations

In order to expand our understanding of the mechanism of stereocontrol in syndiospecific α-olefin polymerization, a family of Cs-symmetric, ansa-group 3 metallocenes was targeted as polymerization catalysts. The syntheses of new ansa-yttrocene and scandocene derivatives that employ the doubly [SiMe₂]- bridged ligand array (1,2-SiMe₂)₂{C₅H-3,5-(CHMe₂)₂} (where R = t- butyl, tBuThp; where R = i-propyl, iPrThp) are described. The structures of tBuThpY(µ-Cl)₂K(THF)₂, tBuThpSc(µ-Cl)₂K(Et₂O)₂, tBuThpYCH(SiMe₃)₂, Y₂{µ₂-(tBuThp)₂}(µ₂-H)₂, and tBuThpSc(µ-CH₃)₂ have been examined by single crystal X-ray diffraction methods. Ansa-yttrocenes and scandocenes that incorporate the singly [CPh₂]-bridged ligand array (CPh₂)(C₅H₄)(C₁₃H₈)(where C₅H₄ = Cp, cyclopentadienyl; where C₁₃H₈ = Flu, fluourenyl) have also been prepared. Select meallocene alkyl complexes are active single component catalysts for homopolymerization of propylene and 1-pentene. The scandocene tetramethylaluminate complexes generate polymers with the highes molecular weights of the series. Under all conditions examined atactic polymer microstructures are observed, suggesting a chain-end mechanism for stereocontrol.

A series of ansa-tantalocenes have been prepared as models for Ziegler-Natta polymerization catalysts. A singly bridged ansa-tantalocene trimethyl complex, Me₂Si(η⁵-C₅H₄)₂TaMe₃, has been prepared and used for the synthesis of a tantalocene ethylene-methyl complex. Addition of propylene to this ethylene-methyl adduct results in olefin exchange to give a mixture of endo and exo propylene isomers. Doubly-silylene bridged ansa-tantalocene complexes have been prepared with the tBuThp ligand; a tantalocene trimethyl complex and a tantalocene methylidene-methyl complex have been synthesized and characterized by X-ray diffraction. Thermolysis of the methylidene-methyl complex affords the corresponding ethylene-hydride complex. Addition of either propylene or styrene to this ethylene-hydride compound results in olefin exchange. In both cases, only one product isomer is observed. Studies of olefin exchange with ansa-tantalocene olefin-hydride and olefin-methyl complexes have provided information about the important steric influences for olefin coordination in Ziegler-Natta polymerization.

Ancillary ligand effects : from zirconium(IV)-catalyzed homogeneous propylene polymerization to platinum(II)-mediated C-H bond activation

A series of C_s- and C₁-symmetric doubly-linked ansa-metallocenes of the general formula {1,1'-SiMe₂-2,2'-E-('ƞ⁵-C₅H₂-4-R¹)-(ƞ⁵-C₅H-3',5'-(CHMe₂)₂)}ZrC₂ (E = SiMe₂ (1), SiPh₂ (2), SiMe₂ -SiMe₂ (3); R¹ = H, CHMe₂, C₅H₉, C₆H₁₁, C₆H₅) has been prepared. When activated by methylaluminoxane, these are active propylene polymerization catalysts. 1 and 2 produce syndiotactic polypropylenes, and 3 produces isotactic polypropylenes. Site epimerization is the major pathway for stereoerror formation for 1 and 2. In addition, the polymer chain has slightly stronger steric interaction with the diphenylsilylene linker than with the dimethylsilylene linker. This results in more frequent site epimerization and reduced syndiospecificity for 2 compared to 1.

C₁-Symmetric ansa-zirconocenes [1,1 '-SiMe₂-(C₅H₄)-(3-R-C₅H₃)]ZrCl₂ (4), [1,1 '-SiMe₂-(C₅H₄)-(2,4-R₂-C₅H₂)]ZrCl₂ (5) and [1,1 '-SiMe₂-2,2 '-(SiMe₂-SiMe₂)-(C₅H₃)-( 4-R-C₅H₂)]ZrCl₂ (6) have been prepared to probe the origin of isospecificity in 3. While 4 and 3 produce polymers with similar isospecificity, 5 and 6 give mostly hemi-isotactic-like polymers. It is proposed that the facile site epimerization via an associative pathway allows rapid equilibration of the polymer chain between the isospecific and aspecific insertion sites. This results in more frequent insertion from the isospecific site, which has a lower kinetic barrier for chain propagation. On the other hand, site epimerization for 5 and 6 is slow. This leads to mostly alternating insertion from the isospecific and aspecific sites, and consequently, a hemi-isotactic-like polymers. In comparison, site epimerization is even slower for 3, but enchainment from the aspecific site has an extremely high kinetic barrier for monomer coordination. Therefore, enchainment occurs preferentially from the isospecific site to produce isotactic polymers.

A series of cationic complexes [(ArN=CR-CR=NAr)PtMe(L)]⁺[BF₄]⁺ (Ar = aryl; R = H, CH₃; L = water, trifluoroethanol) has been prepared. They react smoothly with benzene at approximately room temperature in trifluoroethanol solvent to yield methane and the corresponding phenyl Pt(II) cations, via Pt(IV)-methyl-phenyl-hydride intermediates. The reaction products of methyl-substituted benzenes suggest an inherent reactivity preference for aromatic over benzylic C-H bond activation, which can however be overridden by steric effects. For the reaction of benzene with cationic Pt(II) complexes, in which the diimine ligands bear 3,5-disubstituted aryl groups at the nitrogen atoms, the rate-determining step is C-H bond activation. For the more sterically crowded analogs with 2,6-dimethyl-substituted aryl groups, benzene coordination becomes rate-determining. The more electron-rich the ligand, as reflected by the CO stretching frequency in the IR spectrum of the corresponding cationic carbonyl complex, the faster the rate of C-H bond activation. This finding, however, does not reflect the actual C-H bond activation process, but rather reflects only the relative ease of solvent molecules displacing water molecules to initiate the reaction. That is, the change in rates is mostly due to a ground state effect. Several lines of evidence suggest that associative substitution pathways operate to get the hydrocarbon substrate into, and out of, the coordination sphere; i.e., that benzene substitution proceeds by a solvent- (TFE-) assisted associative pathway.

Baryon number violation beyond the standard model

This thesis describes simple extensions of the standard model with new sources of baryon number violation but no proton decay. The motivation for constructing such theories comes from the shortcomings of the standard model to explain the generation of baryon asymmetry in the universe, and from the absence of experimental evidence for proton decay. However, lack of any direct evidence for baryon number violation in general puts strong bounds on the naturalness of some of those models and favors theories with suppressed baryon number violation below the TeV scale. The initial part of the thesis concentrates on investigating models containing new scalars responsible for baryon number breaking. A model with new color sextet scalars is analyzed in more detail. Apart from generating cosmological baryon number, it gives nontrivial predictions for the neutron-antineutron oscillations, the electric dipole moment of the neutron, and neutral meson mixing. The second model discussed in the thesis contains a new scalar leptoquark. Although this model predicts mainly lepton flavor violation and a nonzero electric dipole moment of the electron, it includes, in its original form, baryon number violating nonrenormalizable dimension-five operators triggering proton decay. Imposing an appropriate discrete symmetry forbids such operators. Finally, a supersymmetric model with gauged baryon and lepton numbers is proposed. It provides a natural explanation for proton stability and predicts lepton number violating processes below the supersymmetry breaking scale, which can be tested at the Large Hadron Collider. The dark matter candidate in this model carries baryon number and can be searched for in direct detection experiments as well. The thesis is completed by constructing and briefly discussing a minimal extension of the standard model with gauged baryon, lepton, and flavor symmetries.

Geometric discretization through primal-dual meshes

This thesis introduces new tools for geometric discretization in computer graphics and computational physics. Our work builds upon the duality between weighted triangulations and power diagrams to provide concise, yet expressive discretization of manifolds and differential operators. Our exposition begins with a review of the construction of power diagrams, followed by novel optimization procedures to fully control the local volume and spatial distribution of power cells. Based on this power diagram framework, we develop a new family of discrete differential operators, an effective stippling algorithm, as well as a new fluid solver for Lagrangian particles. We then turn our attention to applications in geometry processing. We show that orthogonal primal-dual meshes augment the notion of local metric in non-flat discrete surfaces. In particular, we introduce a reduced set of coordinates for the construction of orthogonal primal-dual structures of arbitrary topology, and provide alternative metric characterizations through convex optimizations. We finally leverage these novel theoretical contributions to generate well-centered primal-dual meshes, sphere packing on surfaces, and self-supporting triangulations.