Seismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan

The fine-scale seismic structure of the central Mexico, southern Peru, and southwest Japan subduction zones is studied using intraslab earthquakes recorded by temporary and permanent regional seismic arrays. The morphology of the transition from flat to normal subduction is explored in central Mexico and southern Peru, while in southwest Japan the spatial coincidence of a thin ultra-slow velocity layer (USL) atop the flat slab with locations of slow slip events (SSEs) is explored. This USL is also observed in central Mexico and southern Peru, where its lateral extent is used as one constraint on the nature of the flat-to-normal transitions.

In western central Mexico, I find an edge to this USL which is coincident with the western boundary of the projected Orozco Fracture Zone (OFZ) region. Forward modeling of the 2D structure of the subducted Cocos plate using a finite-difference algorithm provides constraints on the velocity and geometry of the slab’s seismic structure in this region and confirms the location of the USL edge. I propose that the Cocos slab is currently fragmenting into a North Cocos plate and a South Cocos plate along the projection of the OFZ, by a process analogous to that which occurred when the Rivera plate separated from the proto-Cocos plate 10 Ma.

In eastern central Mexico, observations of a sharp transition in slab dip near the abrupt end of the Trans Mexican Volcanic Belt (TMVB) suggest a possible slab tear located within the subducted South Cocos plate. The eastern lateral extent of the USL is found to be coincident with these features and with the western boundary of a zone of decreased seismicity, indicating a change in structure which I interpret as evidence of a possible tear. Analysis of intraslab seismicity patterns and focal mechanism orientations and faulting types provides further support for a possible tear in the South Cocos slab. This potential tear, together with the tear along the projection of the OFZ to the northwest, indicates a slab rollback mechanism in which separate slab segments move independently, allowing for mantle flow between the segments.

In southern Peru, observations of a gradual increase in slab dip coupled with a lack of any gaps or vertical offsets in the intraslab seismicity suggest a smooth contortion of the slab. Concentrations of focal mechanisms at orientations which are indicative of slab bending are also observed along the change in slab geometry. The lateral extent of the USL atop the horizontal Nazca slab is found to be coincident with the margin of the projected linear continuation of the subducting Nazca Ridge, implying a causal relationship, but not a slab tear. Waveform modeling of the 2D structure in southern Peru provides constraints on the velocity and geometry of the slab’s seismic structure and confirms the absence of any tears in the slab.

In southwest Japan, I estimate the location of a possible USL along the Philippine Sea slab surface and find this region of low velocity to be coincident with locations of SSEs that have occurred in this region. I interpret the source of the possible USL in this region as fluids dehydrated from the subducting plate, forming a high pore-fluid pressure layer, which would be expected to decrease the coupling on the plate interface and promote SSEs.

Variable-angle electron spectroscopic studies of various polyatomic molecular systems

The complementary techniques of low-energy, variable-angle electron-impact spectroscopy and ultraviolet variable-angle photoelectron spectroscopy have been used to study the electronic spectroscopy and structure of several series of molecules. Electron-impact studies were performed at incident beam energies between 25 eV and 100 eV and at scattering angles ranging from 0° to 90°. The energy-loss regions from 0 eV to greater than 15 eV were studied. Photoelectron spectroscopic studies were conducted using a HeI radiation source and spectra were measured at scattering angles from 45° to 90°. The molecules studied were chosen because of their spectroscopic, chemical, and structural interest. The operation of a new electron-impact spectrometer with multiple-mode target source capability is described. This spectrometer has been used to investigate the spin-forbidden transitions in a number of molecular systems.

The electron-impact spectroscopy of the six chloro-substituted ethylenes has been studied over the energy-loss region from 0-15 eV. Spin-forbidden excitations corresponding to the π → π*, N → T transition have been observed at excitation energies ranging from 4.13 eV in vinyl chloride to 3.54 eV in tetrachloroethylene. Symmetry-forbidden transitions of the type π → np have been oberved in trans-dichloroethyene and tetrachlor oethylene. In addition, transitions to many states lying above the first ionization potential were observed for the first time. Many of these bands have been assigned to Rydberg series converging to higher ionization potentials. The trends observed in the measured transition energies for the π → π*, N → T, and N → V as well as the π → 3s excitation are discussed and compared to those observed in the methyl- and fluoro- substituted ethylenes.

The electron energy-loss spectra of the group VIb transition metal hexacarbonyls have been studied in the 0 eV to 15 eV region. The differential cross sections were obtained for several features in the 3-7 eV energy-loss region. The symmetry-forbidden nature of the ¹A_1g → ¹A_1g, 2t_2g(π) → 3t_2g(π*) transition in these compounds was confirmed by the high-energy, low-angle behavior of their relative intensities. Several low lying transitions have been assigned to ligand field transitions on the basis of the energy and angular behavior of the differential cross sections for these transitions. No transitions which could clearly be assigned to singlet → triplet excitations involving metal orbitals were located. A number of states lying above the first ionization potential have been observed for the first time. A number of features in the 6-14 eV energy-loss region of the spectra of these compounds correspond quite well to those observed in free CO.

A number of exploratory studies have been performed. The π → π*, N → T, singlet → triplet excitation has been located in vinyl bromide at 4.05 eV. We have also observed this transition at approximately 3.8 eV in a cis-/trans- mixture of the 1,2-dibromoethylenes. The low-angle spectrum of iron pentacarbonyl was measured over the energy-loss region extending from 2-12 eV. A number of transitions of 8 eV or greater excitation energy were observed for the first time. Cyclopropane was also studied at both high and low angles but no clear evidence for any spin- forbidden transitions was found. The electron-impact spectrum of the methyl radical resulting from the pyrolysis of tetramethyl tin was obtained at 100 eV incident energy and at 0° scattering angle. Transitions observed at 5.70 eV and 8.30 eV agree well with the previous optical results. In addition, a number of bands were observed in the 8-14 eV region which are most likely due to Rydberg transitions converging to the higher ionization potentials of this molecule. This is the first reported electron-impact spectrum of a polyatomic free radical.

Variable-angle photoelectron spectroscopic studies were performed on a series of three-membered-ring heterocyclic compounds. These compounds are of great interest due to their highly unusual structure. Photoelectron angular distributions using HeI radiation have been measured for the first time for ethylene oxide and ethyleneimine. The measured anisotropy parameters, β, along with those measured for cyclopropane were used to confirm the orbital correlations and photoelectron band assignments. No high values of β similar to those expected for alkene π orbitals were observed for the Walsh or Forster-Coulson-Moffit type orbitals.

Spin dynamics of pulsed nuclear magnetic double resonance in solids

In the first part of this thesis, experiments utilizing an NMR phase interferometric concept are presented. The spinor character of two-level systems is explicitly demonstrated by using this concept. Following this is the presentation of an experiment which uses this same idea to measure relaxation times of off-diagonal density matrix elements corresponding to magnetic-dipole-forbidden transitions in a ^(13)C-^1H, AX spin system. The theoretical background for these experiments and the spin dynamics of the interferometry are discussed also.

The second part of this thesis deals with NMR dipolar modulated chemical shift spectroscopy, with which internuclear bond lengths and bond angles with respect to the chemical shift principal axis frame are determined from polycrystalline samples. Experiments using benzene and calcium formate verify the validity of the technique in heteronuclear (^(13)C-^1H) systems. Similar experiments on powdered trichloroacetic acid confirm the validity in homonuclear (^1H- ^1H) systems. The theory and spin dynamics are explored in detail, and the effects of a number of multiple pulse sequences are discussed.

The last part deals with an experiment measuring the ^(13)C chemical shift tensor in K_2Pt(CN)_4Br_(0.3) • 3H_2O, a one-dimensional conductor. The ^(13)C spectra are strongly affected by ^(14)N quadrupolar interactions via the ^(13)C - ^(14)N dipolar interaction. Single crystal rotation spectra are shown.

An appendix discussing the design, construction, and performance of a single-coil double resonance NMR sample probe is included.

Determination of electronic energy levels of molecules by 90° low energy electron scattering

A review of the theory of electron scattering indicates that low incident beam energies and large scattering angles are the favorable conditions for the observation of optically forbidden transitions in atoms and molecules.

An apparatus capable of yielding electron impact spectra at 90° with incident electron beam energies between 30 and 50 electron volts is described. The resolution of the instrument is about 1 electron volt.

Impact spectra of thirteen molecules have been obtained. Known forbidden transitions to the helium 2³S, the hydrogen b³Ʃ⁺_u, the nitrogen A³Ʃ⁺_u, B³π_g, a’π_g, and C³π_u, the carbon monoxide a³π, the ethylene ᾶ³B_1u, and the benzene ᾶ³B_1u states from the corresponding ground states have been observed.

In addition, singlet-triplet vertical transitions in acetylene, propyne, propadiene, norbornadiene and quadricyclene, peaking at 5.9, 5.9, 4.5, 3.8, and 4.0 ev (±0.2 ev), respectively, have been observed and assigned for the first time.

Three-dimensional models of mantle convection : influence of phase transitions and temperature-dependent viscosity

Two of the most important questions in mantle dynamics are investigated in three separate studies: the influence of phase transitions (studies 1 and 2), and the influence of temperature-dependent viscosity (study 3).

(1) Numerical modeling of mantle convection in a three-dimensional spherical shell incorporating the two major mantle phase transitions reveals an inherently three-dimensional flow pattern characterized by accumulation of cold downwellings above the 670 km discontinuity, and cylindrical 'avalanches' of upper mantle material into the lower mantle. The exothermic phase transition at 400 km depth reduces the degree of layering. A region of strongly-depressed temperature occurs at the base of the mantle. The temperature field is strongly modulated by this partial layering, both locally and in globally-averaged diagnostics. Flow penetration is strongly wavelength-dependent, with easy penetration at long wavelengths but strong inhibition at short wavelengths. The amplitude of the geoid is not significantly affected.

(2) Using a simple criterion for the deflection of an upwelling or downwelling by an endothermic phase transition, the scaling of the critical phase buoyancy parameter with the important lengthscales is obtained. The derived trends match those observed in numerical simulations, i.e., deflection is enhanced by (a) shorter wavelengths, (b) narrower up/downwellings (c) internal heating and (d) narrower phase loops.

(3) A systematic investigation into the effects of temperature-dependent viscosity on mantle convection has been performed in three-dimensional Cartesian geometry, with a factor of 1000-2500 viscosity variation, and Rayleigh numbers of 10^5-10^7. Enormous differences in model behavior are found, depending on the details of rheology, heating mode, compressibility and boundary conditions. Stress-free boundaries, compressibility, and temperature-dependent viscosity all favor long-wavelength flows, even in internally heated cases. However, small cells are obtained with some parameter combinations. Downwelling plumes and upwelling sheets are possible when viscosity is dependent solely on temperature. Viscous dissipation becomes important with temperature-dependent viscosity.

The sensitivity of mantle flow and structure to these various complexities illustrates the importance of performing mantle convection calculations with rheological and thermodynamic properties matching as closely as possible those of the Earth.

Fermionic quantum systems. Part I: Phase transitions in quantum dots. Part II: Nuclear matter on a lattice

In the first part I perform Hartree-Fock calculations to show that quantum dots (i.e., two-dimensional systems of up to twenty interacting electrons in an external parabolic potential) undergo a gradual transition to a spin-polarized Wigner crystal with increasing magnetic field strength. The phase diagram and ground state energies have been determined. I tried to improve the ground state of the Wigner crystal by introducing a Jastrow ansatz for the wave function and performing a variational Monte Carlo calculation. The existence of so called magic numbers was also investigated. Finally, I also calculated the heat capacity associated with the rotational degree of freedom of deformed many-body states and suggest an experimental method to detect Wigner crystals.

The second part of the thesis investigates infinite nuclear matter on a cubic lattice. The exact thermal formalism describes nucleons with a Hamiltonian that accommodates on-site and next-neighbor parts of the central, spin-exchange and isospin-exchange interaction. Using auxiliary field Monte Carlo methods, I show that energy and basic saturation properties of nuclear matter can be reproduced. A first order phase transition from an uncorrelated Fermi gas to a clustered system is observed by computing mechanical and thermodynamical quantities such as compressibility, heat capacity, entropy and grand potential. The structure of the clusters is investigated with the help two-body correlations. I compare symmetry energy and first sound velocities with literature and find reasonable agreement. I also calculate the energy of pure neutron matter and search for a similar phase transition, but the survey is restricted by the infamous Monte Carlo sign problem. Also, a regularization scheme to extract potential parameters from scattering lengths and effective ranges is investigated.

Enhanced Thermoelectric Performance at the Superionic Phase Transitions of Mixed Ion-Electron Conducting Materials

The quality of a thermoelectric material is judged by the size of its temperature de- pendent thermoeletric-figure-of-merit (zT ). Superionic materials, particularly Zn₄Sb₃ and Cu₂Se, are of current interest for the high zT and low thermal conductivity of their disordered, superionic phase. In this work it is reported that the super-ionic materials Ag₂Se, Cu₂Se and Cu_1.97Ag_0.03Se show enhanced zT in their ordered, normal ion-conducting phases. The zT of Ag₂Se is increased by 30% in its ordered phase as compared to its disordered phase, as measured just below and above its first order phase transition. The zT ’s of Cu₂Se and Cu_1.97Ag_0.03Se both increase by more than 100% over a 30 K temperatures range just below their super-ionic phase transitions. The peak zT of Cu₂Se is 0.7 at 406 K and of Cu_1.97Ag_0.03Se is 1.0 at 400 K. In all three materials these enhancements are due to anomalous increases in their Seebeck coefficients, beyond that predicted by carrier concentration measurements and band structure modeling. As the Seebeck coefficient is the entropy transported per carrier, this suggests that there is an additional quantity of entropy co-transported with charge carriers. Such co-transport has been previously observed via co-transport of vibrational entropy in bipolaron conductors and spin-state entropy in Na_xCo₂O₄. The correlation of the temperature profile of the increases in each material with the nature of their phase transitions indicates that the entropy is associated with the thermodynamcis of ion-ordering. This suggests a new mechanism by which high thermoelectric performance may be understood and engineered.

Phase transitions from the solid state of monatomic elements

The principle aims of this thesis include the development of models of sublimation and melting from first principles and the application of these models to the rare gases.

A simple physical model is constructed to represent the sublimation of monatomic elements. According to this model, the solid and gas phases are two states of a single physical system. The nature of the phase transition is clearly revealed, and the relations between the vapor pressure, the latent heat, and the transition temperature are derived. The resulting theory is applied to argon, krypton, and xenon, and good agreement with experiment is found.

For the melting transition, the solid is represented by an anharmonic model and the liquid is described by the Percus-Yevick approximation. The behavior of the liquid at high densities is studied on the isotherms kT/∈ = 1.3, 1.8, and 2.0, where k is Boltzmann's constant, T is the temperature, and e is the well depth of the Lennard-Jones 12-6 pair potential. No solutions of the PercusYevick equation were found for ρσ³ above 1.3, where ρ is the particle density and σ is the radial parameter of the Lennard-Jones potential. The liquid structure is found to be very different from the solid structure near the melting line. The liquid pressures are about 50 percent low for experimental melting densities of argon. This discrepancy gives rise to melting pressures up to twice the experimental values.

I. Electronic states of heavily-doped molecular crystals--naphthalene. I. Theoretical. II. Experimental. II. Lowest singlet and triplet excited states of pyrazine

PART I

The energy spectrum of heavily-doped molecular crystals was treated in the Green’s function formulation. The mixed crystal Green’s function was obtained by averaging over all possible impurity distributions. The resulting Green’s function, which takes the form of an infinite perturbation expansion, was further approximated by a closed form suitable for numerical calculations. The density-of-states functions and optical spectra for binary mixtures of normal naphthalene and deuterated naphthalene were calculated using the pure crystal density-of-state functions. The results showed that when the trap depth is large, two separate energy bands persist, but when the trap depth is small only a single band exists. Furthermore, in the former case it was found that the intensities of the outer Davydov bands are enhanced whereas the inner bands are weakened. Comparisons with previous theoretical calculations and experimental results are also made.

PART II

The energy states and optical spectra of heavily-doped mixed crystals are investigated. Studies are made for the following binary systems: (1) naphthalene-h₈ and d₈, (2) naphthalene--h₈ and αd₄, and (3) naphthalene--h₈ and βd₁, corresponding to strong, medium and weak perturbations. In addition to ordinary absorption spectra at 4˚K, band-to-band transitions at both 4˚K and 77˚K are also analyzed with emphasis on their relations to cooperative excitation and overall density-of-states functions for mixed crystals. It is found that the theoretical calculations presented in a previous paper agree generally with experiments except for cluster states observed in system (1) at lower guest concentrations. These features are discussed semi-quantitatively. As to the intermolecular interaction parameters, it is found that experimental results compare favorably with calculations based on experimental density-of-states functions but not with those based on octopole interactions or charge-transfer interactions. Previous experimental results of Sheka and the theoretical model of Broude and Rashba are also compared with present investigations.

PART III

The phosphorescence, fluorescence and absorption spectra of pyrazine-h₄ and d₄ have been obtained at 4˚K in a benzene matrix. For comparison, those of the isotopically mixed crystal pyrazine-h₄ in d₄ were also taken. All these spectra show extremely sharp and well-resolved lines and reveal detailed vibronic structure.

The analysis of the weak fluorescence spectrum resolves the long-disputed question of whether one or two transitions are involved in the near-ultraviolet absorption of pyrazine. The “mirror-image relationship” between absorption and emission shows that the lowest singlet state is an allowed transition, properly designated as ¹B_3u ← ¹A_1g. The forbidden component ¹B_2g, predicted by both “exciton” and MO theories to be below the allowed component, must lie higher. Its exact location still remains uncertain.

The phosphorescence spectrum when compared with the excitation phosphorescence spectra, indicates that the lowest triplet state is also symmetry allowed, showing a strong 0-0 band and a “mirror-image relationship” between absorption and emission. In accordance with previous work, the triplet state is designated as ³B_3u.

The vibronic structure of the phosphorescence spectrum is very complicated. Previous work on the analysis of this spectrum all concluded that a long progression of v_6a exists. Under the high resolution attainable in our work, the supposed v_6a progression proves to have a composite triplet structure, starting from the second member of the progression. Not only is the v_9a hydrogen-bending mode present as shown by the appearance of the C-D bending mode in the d₄ spectrum, but a band of 1207 cm^-1 in the pyrazine in benzene system and 1231 cm^-1 in the mixed crystal system is also observed. This band is assigned as 2v_6b and of a_1g symmetry. Its anonymously strong intensity in the phosphorescence spectrum is interpreted as due to the Fermi resonance with the 2v_6a and v_9a band.

To help resolve the present controversy over the crystal phosphorescence spectrum of pyrazine, detailed vibrational analyses of the emission spectra were made. The fluorescence spectrum has essentially the same vibronic structure as the phosphorescence spectrum.

Nuclear structure effects in internal conversion

Experimental studies of nuclear effects in internal conversion in Ta¹⁸¹ and Lu¹⁷⁵ have been performed. Nuclear structure effects (“penetration” effects), in internal conversion are described in general. Calculation of theoretical conversion coefficients are outlined. Comparisons with the theoretical conversion coefficient tables of Rose and Sliv and Band are made. Discrepancies between our results and those of Rose and Sliv are noted. The theoretical conversion coefficients of Sliv and Band are in substantially better agreement with our results than are those of Rose. The ratio of the M1 penetration matrix element to the M1 gamma-ray matrix element, called λ, is equal to + 175 ± 25 for the 482 keV transition in Ta¹⁸¹ . The results for the 343 keV transition in Lu¹⁷⁵ indicate that λ may be as large as – 8 ± 5. These transitions are discussed in terms of the unified collective model. Precision L subshell measurements in Tm¹⁶⁹ (130keV), W¹⁸² (100 keV), and Ta¹⁸¹ (133 keV) show definite systematic deviations from the theoretical conversion coefficients. The possibility of explaining these deviations by penetration effects is investigated and is shown to be excluded. Other explanations of these anomalies are discussed.

1. Ultrasonic studies of binary liquid structure in the critical region. Theory and experiment for the 2,6-lutidine/water system. 2. Hartree-Fock calculations of electric polarizabilities of some simple atoms and molecules, and their practicality. 3. Calculation of vibrational transition probabilities in collinear atom-diatom and diatom-diatom collisions with Lennard-Jones interaction

Part 1. Many interesting visual and mechanical phenomena occur in the critical region of fluids, both for the gas-liquid and liquid-liquid transitions. The precise thermodynamic and transport behavior here has some broad consequences for the molecular theory of liquids. Previous studies in this laboratory on a liquid-liquid critical mixture via ultrasonics supported a basically classical analysis of fluid behavior by M. Fixman (e. g., the free energy is assumed analytic in intensive variables in the thermodynamics)--at least when the fluid is not too close to critical. A breakdown in classical concepts is evidenced close to critical, in some well-defined ways. We have studied herein a liquid-liquid critical system of complementary nature (possessing a lower critical mixing or consolute temperature) to all previous mixtures, to look for new qualitative critical behavior. We did not find such new behavior in the ultrasonic absorption ascribable to the critical fluctuations, but we did find extra absorption due to chemical processes (yet these are related to the mixing behavior generating the lower consolute point). We rederived, corrected, and extended Fixman's analysis to interpret our experimental results in these more complex circumstances. The entire account of theory and experiment is prefaced by an extensive introduction recounting the general status of liquid state theory. The introduction provides a context for our present work, and also points out problems deserving attention. Interest in these problems was stimulated by this work but also by work in Part 3.

Part 2. Among variational theories of electronic structure, the Hartree-Fock theory has proved particularly valuable for a practical understanding of such properties as chemical binding, electric multipole moments, and X-ray scattering intensity. It also provides the most tractable method of calculating first-order properties under external or internal one-electron perturbations, either developed explicitly in orders of perturbation theory or in the fully self-consistent method. The accuracy and consistency of first-order properties are poorer than those of zero-order properties, but this is most often due to the use of explicit approximations in solving the perturbed equations, or to inadequacy of the variational basis in size or composition. We have calculated the electric polarizabilities of H₂, He, Li, Be, LiH, and N₂ by Hartree-Fock theory, using exact perturbation theory or the fully self-consistent method, as dictated by convenience. By careful studies on total basis set composition, we obtained good approximations to limiting Hartree-Fock values of polarizabilities with bases of reasonable size. The values for all species, and for each direction in the molecular cases, are within 8% of experiment, or of best theoretical values in the absence of the former. Our results support the use of unadorned Hartree-Pock theory for static polarizabilities needed in interpreting electron-molecule scattering data, collision-induced light scattering experiments, and other phenomena involving experimentally inaccessible polarizabilities.

Part 3. Numerical integration of the close-coupled scattering equations has been carried out to obtain vibrational transition probabilities for some models of the electronically adiabatic H₂-H₂ collision. All the models use a Lennard-Jones interaction potential between nearest atoms in the collision partners. We have analyzed the results for some insight into the vibrational excitation process in its dependence on the energy of collision, the nature of the vibrational binding potential, and other factors. We conclude also that replacement of earlier, simpler models of the interaction potential by the Lennard-Jones form adds very little realism for all the complication it introduces. A brief introduction precedes the presentation of our work and places it in the context of attempts to understand the collisional activation process in chemical reactions as well as some other chemical dynamics.