Testing inflationary cosmology with the BICEP1 and BICEP2 experiments

Recent observations of the temperature anisotropies of the cosmic microwave background (CMB) favor an inflationary paradigm in which the scale factor of the universe inflated by many orders of magnitude at some very early time. Such a scenario would produce the observed large-scale isotropy and homogeneity of the universe, as well as the scale-invariant perturbations responsible for the observed (10 parts per million) anisotropies in the CMB. An inflationary epoch is also theorized to produce a background of gravitational waves (or tensor perturbations), the effects of which can be observed in the polarization of the CMB. The E-mode (or parity even) polarization of the CMB, which is produced by scalar perturbations, has now been measured with high significance. Con- trastingly, today the B-mode (or parity odd) polarization, which is sourced by tensor perturbations, has yet to be observed. A detection of the B-mode polarization of the CMB would provide strong evidence for an inflationary epoch early in the universe’s history.

In this work, we explore experimental techniques and analysis methods used to probe the B- mode polarization of the CMB. These experimental techniques have been used to build the Bicep2 telescope, which was deployed to the South Pole in 2009. After three years of observations, Bicep2 has acquired one of the deepest observations of the degree-scale polarization of the CMB to date. Similarly, this work describes analysis methods developed for the Bicep1 three-year data analysis, which includes the full data set acquired by Bicep1. This analysis has produced the tightest constraint on the B-mode polarization of the CMB to date, corresponding to a tensor-to-scalar ratio estimate of r = 0.04±0.32, or a Bayesian 95% credible interval of r < 0.70. These analysis methods, in addition to producing this new constraint, are directly applicable to future analyses of Bicep2 data. Taken together, the experimental techniques and analysis methods described herein promise to open a new observational window into the inflationary epoch and the initial conditions of our universe.

Tools for the study of dynamical spacetimes

This thesis covers a range of topics in numerical and analytical relativity, centered around introducing tools and methodologies for the study of dynamical spacetimes. The scope of the studies is limited to classical (as opposed to quantum) vacuum spacetimes described by Einstein's general theory of relativity. The numerical works presented here are carried out within the Spectral Einstein Code (SpEC) infrastructure, while analytical calculations extensively utilize Wolfram's Mathematica program.

We begin by examining highly dynamical spacetimes such as binary black hole mergers, which can be investigated using numerical simulations. However, there are difficulties in interpreting the output of such simulations. One difficulty stems from the lack of a canonical coordinate system (henceforth referred to as gauge freedom) and tetrad, against which quantities such as Newman-Penrose Psi_4 (usually interpreted as the gravitational wave part of curvature) should be measured. We tackle this problem in Chapter 2 by introducing a set of geometrically motivated coordinates that are independent of the simulation gauge choice, as well as a quasi-Kinnersley tetrad, also invariant under gauge changes in addition to being optimally suited to the task of gravitational wave extraction.

Another difficulty arises from the need to condense the overwhelming amount of data generated by the numerical simulations. In order to extract physical information in a succinct and transparent manner, one may define a version of gravitational field lines and field strength using spatial projections of the Weyl curvature tensor. Introduction, investigation and utilization of these quantities will constitute the main content in Chapters 3 through 6.

For the last two chapters, we turn to the analytical study of a simpler dynamical spacetime, namely a perturbed Kerr black hole. We will introduce in Chapter 7 a new analytical approximation to the quasi-normal mode (QNM) frequencies, and relate various properties of these modes to wave packets traveling on unstable photon orbits around the black hole. In Chapter 8, we study a bifurcation in the QNM spectrum as the spin of the black hole a approaches extremality.

An experimental investigation of the heat transfer from a buoyant gas plume to a horizontal unobstructed ceiling

This thesis presents an experimental investigation of the axisymmetric heat transfer from a small scale fire and resulting buoyant plume to a horizontal, unobstructed ceiling during the initial stages of development. A propane-air burner yielding a heat source strength between 1.0 kW and 1.6 kW was used to simulate the fire, and measurements proved that this heat source did satisfactorily represent a source of buoyancy only. The ceiling consisted of a 1/16" steel plate of 0.91 m. diameter, insulated on the upper side. The ceiling height was adjustable between 0.5 m and 0.91 m. Temperature measurements were carried out in the plume, ceiling jet, and on the ceiling.

Heat transfer data were obtained by using the transient method and applying corrections for the radial conduction along the ceiling and losses through the insulation material. The ceiling heat transfer coefficient was based on the adiabatic ceiling jet temperature (recovery temperature) reached after a long time. A parameter involving the source strength Q and ceiling height H was found to correlate measurements of this temperature and its radial variation. A similar parameter for estimating the ceiling heat transfer coefficient was confirmed by the experimental results.

This investigation therefore provides reasonable estimates for the heat transfer from a buoyant gas plume to a ceiling in the axisymmetric case, for the stagnation region where such heat transfer is a maximum and for the ceiling jet region (r/H ≤ 0.7). A comparison with data from experiments which involved larger heat sources indicates that the predicted scaling of temperatures and heat transfer rates for larger scale fires is adequate.

Wnt and FGF signaling in C. elegans vulval cell lineage polarity

The interpretation of extracellular cues leading to the polarization of intracellular components and asymmetric cell divisions is a fundamental part of metazoan organogenesis. The C. elegans vulva, with its invariant cell lineage and interaction of multiple cell signaling pathways, provides an excellent model for the study of cell polarity within an organized epithelial tissue. Herein I discuss the interaction of Wnt and FGF signaling in controlling vulval cell lineage polarity with emphasis on the posterior-most cell that forms the vulva, P7.p.

The mirror symmetry of the C. elegans vulva is achieved by the opposite division orientation of the vulval precursor cells (VPCs) flanking the axis of symmetry. Opposing Wnt signals control the division patterns of the VPCs by controlling the localization of SYS-1/ β-catenin toward the direction of the Wnt gradient. Multiple Wnt signals, expressed at the axis of symmetry, promote the wild-type, anterior-facing, P7.p orientation, whereas Wnts EGL-20 and CWN-1 from the tail and posterior body wall muscle, respectively, promote the daughter cells of P7.p to face the posterior. EGL-20 acts through a member of the LDL receptor superfamily, LRP-2, along with Ror/CAM-1 and Van Gogh/VANG-1. All three transmembrane proteins control orientation through the localization of the SYS-1.

The Fibroblast Growth Factor (FGF) pathway acts in concert with LIN-17/Frizzled to regulate the localization of SYS-1. The source of the FGF ligand is the 1° VPC, P6.p, which controls the polarity of the neighboring 2° VPC, P7.p, by signaling through the sex myoblasts (SMs), activating the FGF pathway. The Wnt, cwn-1, is expressed in the posterior body wall muscle of the worm as well as the SMs, making it the only Wnt expressed on the posterior and anterior sides of P7.p at the time of the polarity decision. Both sources of cwn-1 act instructively to influence P7.p polarity in the direction of the Wnt gradient. The FGF pathway leads to the regulation of cwn-1 transcripts in the SMs. These results illustrate the first evidence of the interaction between FGF and Wnt in C. elegans development and vulval cell lineage polarity as well as highlight the promiscuous nature of Wnt signaling within C. elegans.

Topics in LIGO-related physics: interferometric speed meters and tidal work

In the quest to develop viable designs for third-generation optical interferometric gravitational-wave detectors, one strategy is to monitor the relative momentum or speed of the test-mass mirrors, rather than monitoring their relative position. The most straightforward design for a speed-meter interferometer that accomplishes this is described and analyzed in Chapter 2. This design (due to Braginsky, Gorodetsky, Khalili, and Thorne) is analogous to a microwave-cavity speed meter conceived by Braginsky and Khalili. A mathematical mapping between the microwave speed meter and the optical interferometric speed meter is developed and used to show (in accord with the speed being a quantum nondemolition observable) that in principle the interferometric speed meter can beat the gravitational-wave standard quantum limit (SQL) by an arbitrarily large amount, over an arbitrarily wide range of frequencies . However, in practice, to reach or beat the SQL, this specific speed meter requires exorbitantly high input light power. The physical reason for this is explored, along with other issues such as constraints on performance due to optical dissipation.

Chapter 3 proposes a more sophisticated version of a speed meter. This new design requires only a modest input power and appears to be a fully practical candidate for third-generation LIGO. It can beat the SQL (the approximate sensitivity of second-generation LIGO interferometers) over a broad range of frequencies (~ 10 to 100 Hz in practice) by a factor h/h_SQL ~ √W^(SQL)_(circ)/W_circ. Here W_circ is the light power circulating in the interferometer arms and W_SQL ≃ 800 kW is the circulating power required to beat the SQL at 100 Hz (the LIGO-II power). If squeezed vacuum (with a power-squeeze factor e^-2R) is injected into the interferometer's output port, the SQL can be beat with a much reduced laser power: h/h_SQL ~ √W^(SQL)_(circ)/W_circe^-2R. For realistic parameters (e^-2R ≃ 10 and W_circ ≃ 800 to 2000 kW), the SQL can be beat by a factor ~ 3 to 4 from 10 to 100 Hz. [However, as the power increases in these expressions, the speed meter becomes more narrow band; additional power and re-optimization of some parameters are required to maintain the wide band.] By performing frequency-dependent homodyne detection on the output (with the aid of two kilometer-scale filter cavities), one can markedly improve the interferometer's sensitivity at frequencies above 100 Hz.

Chapters 2 and 3 are part of an ongoing effort to develop a practical variant of an interferometric speed meter and to combine the speed meter concept with other ideas to yield a promising third- generation interferometric gravitational-wave detector that entails low laser power.

Chapter 4 is a contribution to the foundations for analyzing sources of gravitational waves for LIGO. Specifically, it presents an analysis of the tidal work done on a self-gravitating body (e.g., a neutron star or black hole) in an external tidal field (e.g., that of a binary companion). The change in the mass-energy of the body as a result of the tidal work, or "tidal heating," is analyzed using the Landau-Lifshitz pseudotensor and the local asymptotic rest frame of the body. It is shown that the work done on the body is gauge invariant, while the body-tidal-field interaction energy contained within the body's local asymptotic rest frame is gauge dependent. This is analogous to Newtonian theory, where the interaction energy is shown to depend on how one localizes gravitational energy, but the work done on the body is independent of that localization. These conclusions play a role in analyses, by others, of the dynamics and stability of the inspiraling neutron-star binaries whose gravitational waves are likely to be seen and studied by LIGO.

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

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

Measurement of the neutron (^3He) spin structure function at low Q^2 : a connection between the Bjorken and Drell-Hearn-Gerasimov sum rules

The spin dependent cross sections, σT_1/2 and σT_3/2 , and asymmetries, A_∥ and A_⊥ for ³He have been measured at the Jefferson Lab's Hall A facility. The inclusive scattering process ³He(e,e)X was performed for initial beam energies ranging from 0.86 to 5.1 GeV, at a scattering angle of 15.5°. Data includes measurements from the quasielastic peak, resonance region, and the deep inelastic regime. An approximation for the extended Gerasimov-Drell-Hearn integral is presented at a 4-momentum transfer Q² of 0.2-1.0 GeV².

Also presented are results on the performance of the polarized ³He target. Polarization of ³He was achieved by the process of spin-exchange collisions with optically pumped rubidium vapor. The ³He polarization was monitored using the NMR technique of adiabatic fast passage (AFP). The average target polarization was approximately 35% and was determined to have a systematic uncertainty of roughly ±4% relative.

A study of a new electromechanical energy conversion process using superconducting frequency modulated resonators

Experimental demonstrations and theoretical analyses of a new electromechanical energy conversion process which is made feasible only by the unique properties of superconductors are presented in this dissertation. This energy conversion process is characterized by a highly efficient direct energy transformation from microwave energy into mechanical energy or vice versa and can be achieved at high power level. It is an application of a well established physical principle known as the adiabatic theorem (Boltzmann-Ehrenfest theorem) and in this case time dependent superconducting boundaries provide the necessary interface between the microwave energy on one hand and the mechanical work on the other. The mechanism which brings about the conversion is another known phenomenon - the Doppler effect. The resonant frequency of a superconducting resonator undergoes continuous infinitesimal shifts when the resonator boundaries are adiabatically changed in time by an external mechanical mechanism. These small frequency shifts can accumulate coherently over an extended period of time to produce a macroscopic shift when the resonator remains resonantly excited throughout this process. In addition, the electromagnetic energy in s ide the resonator which is proportional to the oscillation frequency is al so accordingly changed so that a direct conversion between electromagnetic and mechanical energies takes place. The intrinsically high efficiency of this process is due to the electromechanical interactions involved in the conversion rather than a process of thermodynamic nature and therefore is not limited by the thermodynamic value.

A highly reentrant superconducting resonator resonating in the range of 90 to 160 MHz was used for demonstrating this new conversion technique. The resonant frequency was mechanically modulated at a rate of two kilohertz. Experimental results showed that the time evolution of the electromagnetic energy inside this frequency modulated (FM) superconducting resonator indeed behaved as predicted and thus demonstrated the unique features of this process. A proposed usage of FM superconducting resonators as electromechanical energy conversion devices is given along with some practical design considerations. This device seems to be very promising in producing high power (~10W/cm^3) microwave energy at 10 - 30 GHz.

Weakly coupled FM resonator system is also analytically studied for its potential applications. This system shows an interesting switching characteristic with which the spatial distribution of microwave energies can be manipulated by external means. It was found that if the modulation was properly applied, a high degree (>95%) of unidirectional energy transfer from one resonator to the other could be accomplished. Applications of this characteristic to fabricate high efficiency energy switching devices and high power microwave pulse generators are also found feasible with present superconducting technology.

Magnetospheric access of solar particles and the configuration of the distant geomagnetic field

The access of 1.2-40 MeV protons and 0.4-1.0 MeV electrons from interplanetary space to the polar cap regions has been investigated with an experiment on board a low altitude, polar orbiting satellite (OG0-4).

A total of 333 quiet time observations of the electron polar cap boundary give a mapping of the boundary between open and closed geomagnetic field lines which is an order of magnitude more comprehensive than previously available.

Persistent features (north/south asymmetries) in the polar cap proton flux, which are established as normal during solar proton events, are shown to be associated with different flux levels on open geomagnetic field lines than on closed field lines. The pole in which these persistent features are observed is strongly correlated to the sector structure of the interplanetary magnetic field and uncorrelated to the north/south component of this field. The features were observed in the north (south) pole during a negative (positive) sector 91% of the time, while the solar field had a southward component only 54% of the time. In addition, changes in the north/south component have no observable effect on the persistent features.

Observations of events associated with co-rotating regions of enhanced proton flux in interplanetary space are used to establish the characteristics of the 1.2 - 40 MeV proton access windows: the access window for low polar latitudes is near the earth, that for one high polar latitude region is ~250 R_⊕ behind the earth, while that for the other high polar latitude region is ~1750 R_⊕ behind the earth. All of the access windows are of approximately the same extent (~120 R_⊕). The following phenomena contribute to persistent polar cap features: limited interplanetary regions of enhanced flux propagating past the earth, radial gradients in the interplanetary flux, and anisotropies in the interplanetary flux.

These results are compared to the particle access predictions of the distant geomagnetic tail configurations proposed by Michel and Dessler, Dungey, and Frank. The data are consistent with neither the model of Michel and Dessler nor that of Dungey. The model of Frank can yield a consistent access window configuration provided the following constraints are satisfied: the merging rate for open field lines at one polar neutral point must be ~5 times that at the other polar neutral point, related to the solar magnetic field configuration in a consistent fashion, the migration time for open field lines to move across the polar cap region must be the same in both poles, and the open field line merging rate at one of the polar neutral points must be at least as large as that required for almost all the open field lines to have merged in 0 (one hour). The possibility of satisfying these constraints is investigated in some detail.

The role played by interplanetary anisotropies in the observation of persistent polar cap features is discussed. Special emphasis is given to the problem of non-adiabatic particle entry through regions where the magnetic field is changing direction. The degree to which such particle entry can be assumed to be nearly adiabatic is related to the particle rigidity, the angle through which the field turns, and the rate at which the field changes direction; this relationship is established for the case of polar cap observations.

A quantitative investigation of the solar modulation of cosmic-ray protons and helium nuclei

The differential energy spectra of cosmic-ray protons and He nuclei have been measured at energies up to 315 MeV/nucleon using balloon- and satellite-borne instruments. These spectra are presented for solar quiet times for the years 1966 through 1970. The data analysis is verified by extensive accelerator calibrations of the detector systems and by calculations and measurements of the production of secondary protons in the atmosphere.

The spectra of protons and He nuclei in this energy range are dominated by the solar modulation of the local interstellar spectra. The transport equation governing this process includes as parameters the solar-wind velocity, V, and a diffusion coefficient, K(r,R), which is assumed to be a scalar function of heliocentric radius, r, and magnetic rigidity, R. The interstellar spectra, j_D, enter as boundary conditions on the solutions to the transport equation. Solutions to the transport equation have been calculated for a broad range of assumed values for K(r,R) and j_D and have been compared with the measured spectra.

It is found that the solutions may be characterized in terms of a dimensionless parameter, ψ(r,R) = ^∞∫_r V dr'/(K(r',R). The amount of modulation is roughly proportional to ψ. At high energies or far from the Sun, where the modulation is weak, the solution is determined primarily by the value of ψ (and the interstellar spectrum) and is not sensitive to the radial dependence of the diffusion coefficient. At low energies and for small r, where the effects of adiabatic deceleration are found to be large, the spectra are largely determined by the radial dependence of the diffusion coefficient and are not very sensitive to the magnitude of ψ or to the interstellar spectra. This lack of sensitivity to j_D implies that the shape of the spectra at Earth cannot be used to determine the interstellar intensities at low energies.

Values of ψ determined from electron data were used to calculate the spectra of protons and He nuclei near Earth. Interstellar spectra of the form j_D α (W - 0.25m)^-2.65 for both protons and He nuclei were found to yield the best fits to the measured spectra for these values of ψ, where W is the total energy and m is the rest energy. A simple model for the diffusion coefficient was used in which the radial and rigidity dependence are separable and K is independent of radius inside a modulation region which has a boundary at a distance D. Good agreement was found between the measured and calculated spectra for the years 1965 through 1968, using typical boundary distances of 2.7 and 6.1 A.U. The proton spectra observed in 1969 and 1970 were flatter than in previous years. This flattening could be explained in part by an increase in D, but also seemed to require that a noticeable fraction of the observed protons at energies as high at 50 to 100 MeV be attributed to quiet-time solar emission. The turnup in the spectra at low energies observed in all years was also attributed to solar emission. The diffusion coefficient used to fit the 1965 spectra is in reasonable agreement with that determined from the power spectra of the interplanetary magnetic field (Jokipii and Coleman, 1968). We find a factor of roughly 3 increase in ψ from 1965 to 1970, corresponding to the roughly order of magnitude decrease in the proton intensity at 250 MeV. The change in ψ might be attributed to a decrease in the diffusion coefficient, or, if the diffusion coefficient is essentially unchanged over that period (Mathews et al., 1971), might be attributed to an increase in the boundary distance, D.

Strength, Deformation and Fracture in Metallic Nanostructures

An understanding of the mechanics of nanoscale metals and semiconductors is necessary for the safe and prolonged operation of nanostructured devices from transistors to nanowire- based solar cells to miniaturized electrodes. This is a fascinating but challenging pursuit because mechanical properties that are size-invariant in conventional materials, such as strength, ductility and fracture behavior, can depend critically on sample size when materials are reduced to sub- micron dimensions. In this thesis, the effect of nanoscale sample size, microstructure and structural geometry on mechanical strength, deformation and fracture are explored for several classes of solid materials. Nanocrystalline platinum nano-cylinders with diameters of 60 nm to 1 μm and 12 nm sized grains are fabricated and tested in compression. We find that nano-sized metals containing few grains weaken as sample diameter is reduced relative to grain size due to a change from deformation governed by internal grains to surface grain governed deformation. Fracture at the nanoscale is explored by performing in-situ SEM tension tests on nanocrystalline platinum and amorphous, metallic glass nano-cylinders containing purposely introduced structural flaws. It is found that failure location, mechanism and strength are determined by the stress concentration with the highest local stress whether this is at the structural flaw or a microstructural feature. Principles of nano-mechanics are used to design and test mechanically robust hierarchical nanostructures with structural and electrochemical applications. 2-photon lithography and electroplating are used to fabricate 3D solid Cu octet meso-lattices with micron- scale features that exhibit strength higher than that of bulk Cu. An in-situ SEM lithiation stage is developed and used to simultaneously examine morphological and electrochemical changes in Si-coated Cu meso-lattices that are of interest as high energy capacity electrodes for Li-ion batteries.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.

Effect of Microstructural Interfaces on the Mechanical Response of Crystalline Metallic Materials

Advances in nano-scale mechanical testing have brought about progress in the understanding of physical phenomena in materials and a measure of control in the fabrication of novel materials. In contrast to bulk materials that display size-invariant mechanical properties, sub-micron metallic samples show a critical dependence on sample size. The strength of nano-scale single crystalline metals is well-described by a power-law function, σαD^-n, where D is a critical sample size and n is a experimentally-fit positive exponent. This relationship is attributed to source-driven plasticity and demonstrates a strengthening as the decreasing sample size begins to limit the size and number of dislocation sources. A full understanding of this size-dependence is complicated by the presence of microstructural features such as interfaces that can compete with the dominant dislocation-based deformation mechanisms. In this thesis, the effects of microstructural features such as grain boundaries and anisotropic crystallinity on nano-scale metals are investigated through uniaxial compression testing. We find that nano-sized Cu covered by a hard coating displays a Bauschinger effect and the emergence of this behavior can be explained through a simple dislocation-based analytic model. Al nano-pillars containing a single vertically-oriented coincident site lattice grain boundary are found to show similar deformation to single-crystalline nano-pillars with slip traces passing through the grain boundary. With increasing tilt angle of the grain boundary from the pillar axis, we observe a transition from dislocation-dominated deformation to grain boundary sliding. Crystallites are observed to shear along the grain boundary and molecular dynamics simulations reveal a mechanism of atomic migration that accommodates boundary sliding. We conclude with an analysis of the effects of inherent crystal anisotropy and alloying on the mechanical behavior of the Mg alloy, AZ31. Through comparison to pure Mg, we show that the size effect dominates the strength of samples below 10 μm, that differences in the size effect between hexagonal slip systems is due to the inherent crystal anisotropy, suggesting that the fundamental mechanism of the size effect in these slip systems is the same.

Norbornenyl- and nortricyclyl free radicals

The number, symmetry, and product-forming capabilities of the intermediates in the photoinitiated reductions of endo- and exo-5- bromonorbornene and 2-bromonortricyclene with tri-n-butyltin hydride at temperatures between -10° and 22° were investigated.

Three mechanisms were evaluated:

1. The 5-norbornenyl- and 2-nortricyclyl radicals isomerize reversibly with the former producing nortricyclene by abstraction of hydrogen from tri-n-butyltin hydride.

2. The 5-norbornenyl- and 2-nortricyclyl radicals isomerize reversibly, but some norbornene can be formed from the 2-nortricyclyl radical or some nortricyclene can be formed from the 5-norbornenyl radical by abstraction of hydrogen.

3. There is intervention of a "bridged" radical which may be for med reversibly or irreversibly from the 5-norbornenyl- and 2-nortricyclyl radicals.

Within small error limits, the ratios of norbornene to nortricyclene as a function of the concentration of tri-n-butyltin hydride are consistent with the first mechanism.

In the reductions with tri-n-butyltin deuteride, primary deuterium isotope effects of 2. 3 and 2. 1 for the abstraction of deuterium by the 2-nortricyclyl- and 5-norbornenyl radicals, respectively, were found. The primary deuterium isotope effects were invariant with the concentration of tri-n-butyltin deuteride, although the ratios of norbornene to nortricyclene changed appreciably over this range. This is consistent with the first mechanism, and can accommodate the formation of either product from more than one intermediate only if the primary kinetic deuterium isotope effects are nearly equal for all reactions leading to the single product.

The reduction of endo-5-bromonorbornene-5, 6, 6-d₃ with tri-n-butyltin hydride or tri-n-butyltin deuteride leads to both unrearranged and rearranged norbornenes. The ratios of unrearranged to rearranged norbornene require that the 5-norbornenyl-5, 6, 6-d₃ radical isomerize to an intermediate with the symmetry expected of a nortricyclyl free radical. The results are consistent with mechanism 1, but imply a surprising normal secondary kinetic deuterium isotope effect of about 1.25 for the abstraction of hydrogen by the 5-norbornenyl- 5, 6, 6-d₃ radical.

Approximate calculations show that there does not appear to be any substantial difference in the stabilities of the 5-norbornenyl and 2-nortricyclyl radicals.

Although the results can not exclude a small contribution by a mechanism other than mechanism 1, no such contribution is required to adequately explain the results.