Stability of Electrode-Electrolyte Interfaces during Charging in Lithium Batteries

In this thesis we study the growth of a Li electrode-electrolyte interface in the presence of an elastic prestress. In particular, we focus our interest on Li-air batteries with a solid electrolyte, LIPON, which is a new type of secondary or rechargeable battery. Theoretical studies and experimental evidence show that during the process of charging the battery the replated lithium adds unevenly to the electrode surface. This phenomenon eventually leads to dendrite formation as the battery is charged and discharged numerous times. In order to suppress or alleviate this deleterious effect of dendrite growth, we put forth a study based on a linear stability analysis. Taking into account all the mechanisms of mass transport and interfacial kinetics, we model the evolution of the interface. We find that, in the absence of stress, the stability of a planar interface depends on interfacial diffusion properties and interfacial energy. Specifically, if Herring-Mullins capillarity-driven interfacial diffusion is accounted for, interfaces are unstable against all perturbations of wavenumber larger than a critical value. We find that the effect of an elastic prestress is always to stabilize planar interfacial growth by increasing the critical wavenumber for instability. A parametric study results in quantifying the extent of the prestress stabilization in a manner that can potentially be used in the design of Li-air batteries. Moreover, employing the theory of finite differences we numerically solve the equation that describes the evolution of the surface profile and present visualization results of the surface evolution by time. Lastly, numerical simulations performed in a commercial finite element software validate the theoretical formulation of the interfacial elastic energy change with respect to the planar interface.

The organometallic chemistry of the Ru(001) surface

The organometallic chemistry of the hexagonally close-packed Ru(001) surface has been studied using electron energy loss spectroscopy and thermal desorption mass spectrometry. The molecules that have been studied are acetylene, formamide and ammonia. The chemistry of acetylene and formamide has also been investigated in the presence of coadsorbed hydrogen and oxygen adatoms.

Acetylene is adsorbed molecularly on Ru(001) below approximately 230 K, with rehybridization of the molecule to nearly sp^3 occurring. The principal decomposition products at higher temperatures are ethylidyne (CCH_3) and acetylide (CCH) between 230 and 350 K, and methylidyne (CH) and surface carbon at higher temperatures. Some methylidyne is stable to approximately 700 K. The preadsorption of hydrogen does not alter the decomposition products of acetylene, but reduces the saturation coverage and also leads to the formation of a small amount of ethylene (via an η^2-CHCH_2 species) which desorbs molecularly near 175 K. Preadsorbed oxygen also reduces the saturation coverage of acetylene but has virtually no effect on the nature of the molecularly chemisorbed acetylene. It does, however, lead to the formation of an sp^2-hybridized vinylidene (CCH_2) species in the decomposition of acetylene, in addition to the decomposition products that are formed on the clean surface. There is no molecular desorption of chemisorbed acetylene from clean Ru(001), hydrogen-presaturated Ru(001), or oxygen-presaturated Ru(001).

The adsorption and decomposition of formamide has been studied on clean Ru(001), hydrogen-presaturated Ru(001), and Ru(001)-p(1x2)-O (oxygen adatom coverage = 0.5). On clean Ru(001), the adsorption of low coverages of formamide at 80 K results in CH bond cleavage and rehybridization of the carbonyl double bond to produce an η^2 (C,O)-NH_2CO species. This species is stable to approximately 250 K at which point it decomposes to yield a mixture of coadsorbed carbon monoxide, ammonia, an NH species and hydrogen adatoms. The decomposition of NH to hydrogen and nitrogen adatoms occurs between 350 and 400 K, and the thermal desorption products are NH_3 (-315 K), H_2 (-420 K), CO (-480 K) and N_2 (-770 K). At higher formamide coverages, some formamide is adsorbed molecularly at 80 K, leading both to molecular desorption and to the formation of a new surface intermediate between 300 and 375 K that is identified tentatively as η^1(N)-NCHO. On Ru(001)- p(1x2)-O and hydrogen-presaturated Ru(001), formamide adsorbs molecularly at 80 K in an η^1(O)- NH_2CHO configuration. On the oxygen-precovered surface, the molecularly adsorbed formamide undergoes competing desorption and decomposition, resulting in the formation of an η^2(N,O)-NHCHO species (analogous to a bidentate formate) at approximately 265 K. This species decomposes near 420 K with the evolution of CO and H_2 into the gas phase. On the hydrogen precovered surface, the Η^1(O)-NH_2CHO converts below 200 K to η^2(C,O)-NH_2CHO and η^2(C,O)-NH^2CO, with some molecular desorption occurring also at high coverage. The η^2(C,O)-bonded species decompose in a manner similar to the decomposition of η^2(C,O)-NH_2CO on the clean surface, although the formation of ammonia is not detected.

Ammonia adsorbs reversibly on Ru(001) at 80 K, with negligible dissociation occurring as the surface is annealed The EEL spectra of ammonia on Ru(001) are very similar to those of ammonia on other metal surfaces. Off-specular EEL spectra of chemisorbed ammonia allow the v(Ru-NH_3) and ρ(NH_3) vibrational loss features to be resolved near 340 and 625 cm^(-1), respectively. The intense δ_g (NH_3) loss feature shifts downward in frequency with increasing ammonia coverage, from approximately 1160 cm^(-1) in the low coverage limit to 1070 cm^(-1) at saturation. In coordination compounds of ammonia, the frequency of this mode shifts downward with decreasing charge on the metal atom, and its downshift on Ru(001) can be correlated with the large work function decrease that the surface has previously been shown to undergo when ammonia is adsorbed. The EELS data are consistent with ammonia adsorption in on-top sites. Second-layer and multilayer ammonia on Ru(001) have also been characterized vibrationally, and the results are similar to those obtained for other metal surfaces.

Geological evolution of two crustal scale shear zones. Part I. The Rand thrust complex, northwestern Mojave Desert, California. Part II. The Magdalena metamorphic core complex, north central Sonora, Mexico

The geology and structure of two crustal scale shear zones were studied to understand the partitioning of strain within intracontinental orogenic belts. Movement histories and regional tectonic implications are deduced from observational data. The two widely separated study areas bear the imprint of intense Late Mesozoic through Middle Cenozoic tectonic activity. A regional transition from Late Cretaceous-Early Tertiary plutonism, metamorphism, and shortening strain to Middle Tertiary extension and magmatism is preserved in each area, with contrasting environments and mechanisms. Compressional phases of this tectonic history are better displayed in the Rand Mountains, whereas younger extensional structures dominate rock fabrics in the Magdalena area.

In the northwestern Mojave desert, the Rand Thrust Complex reveals a stack of four distinctive tectonic plates offset along the Garlock Fault. The lowermost plate, Rand Schist, is composed of greenschist facies metagraywacke, metachert, and metabasalt. Rand Schist is structurally overlain by Johannesburg Gneiss (= garnet-amphibolite grade orthogneisses, marbles and quartzites), which in turn is overlain by a Late Cretaceous hornblende-biotite granodiorite. Biotite granite forms the fourth and highest plate. Initial assembly of the tectonic stack involved a Late Cretaceous? south or southwest vergent overthrusting event in which Johannesburg Gneiss was imbricated and attenuated between Rand Schist and hornblende-biotite granodiorite. Thrusting postdated metamorphism and deformation of the lower two plates in separate environments. A post-kinematic stock, the Late Cretaceous Randsburg Granodiorite, intrudes deep levels of the complex and contains xenoliths of both Rand Schist and mylonitized Johannesburg? gneiss. Minimum shortening implied by the map patterns is 20 kilometers.

Some low angle faults of the Rand Thrust Complex formed or were reactivated between Late Cretaceous and Early Miocene time. South-southwest directed mylonites derived from Johannesburg Gneiss are commonly overprinted by less penetrative north-northeast vergent structures. Available kinematic information at shallower structural levels indicates that late disturbance(s) culminated in northward transport of the uppermost plate. Persistence of brittle fabrics along certain structural horizons suggests a possible association of late movement(s) with regionally known detachment faults. The four plates were juxtaposed and significant intraplate movements had ceased prior to Early Miocene emplacement of rhyolite porphyry dikes.

In the Magdalena region of north central Sonora, components of a pre-Middle Cretaceous stratigraphy are used as strain markers in tracking the evolution of a long lived orogenic belt. Important elements of the tectonic history include: (1) Compression during the Late Cretaceous and Early Tertiary, accompanied by plutonism, metamorphism, and ductile strain at depth, and thrust driven? syntectonic sedimentation at the surface. (2) Middle Tertiary transition to crustal extension, initially recorded by intrusion of leucogranites, inflation of the previously shortened middle and upper crustal section, and surface volcanism. (3) Gravity induced development of a normal sense ductile shear zone at mid crustal levels, with eventual detachment and southwestward displacement of the upper crustal stratigraphy by Early Miocene time.

Elucidation of the metamorphic core complex evolution just described was facilitated by fortuitous preservation of a unique assemblage of rocks and structures. The "type" stratigraphy utilized for regional correlation and strain analysis includes a Jurassic volcanic arc assemblage overlain by an Upper Jurassic-Lower Cretaceous quartz pebble conglomerate, in turn overlain by marine strata with fossiliferous Aptian-Albian limestones. The Jurassic strata, comprised of (a) rhyolite porphyries interstratified with quartz arenites, (b) rhyolite cobble conglomerate, and (c) intrusive granite porphyries, are known to rest on Precambrian basement north and east of the study area. The quartz pebble conglomerate is correlated with the Glance Conglomerate of southeastern Arizona and northeastern Sonora. The marine sequence represents part of an isolated arm? of the Bisbee Basin.

Crosscutting structural relationships between the pre-Middle Cretaceous supracrustal section, younger plutons, and deformational fabrics allow the tectonic sequence to be determined. Earliest phases of a Late Cretaceous-Early Tertiary orogeny are marked by emplacement of the 78 ± 3 Ma Guacomea Granodiorite (U/Pb zircon, Anderson et al., 1980) as a sill into deep levels of the layered Jurassic series. Subsequent regional metamorphism and ductile strain is recorded by a penetrative schistosity and lineation, and east-west trending folds. These fabrics are intruded by post-kinematic Early Tertiary? two mica granites. At shallower crustal levels, the orogeny is represented by north directed thrust faulting, formation of a large intermontane basin, and development of a pronounced unconformity. A second important phase of ductile strain followed Middle Tertiary? emplacement of leucogranites as sills and northwest trending dikes into intermediate levels of the deformed section (surficial volcanism was also active during this transitional period to regional extension). Gravitational instabilities resulting from crustal swelling via intrusion and thermal expansion led to development of a ductile shear zone within the stratigraphic horizon occupied by a laterally extensive leucogranite sill. With continued extension, upper crustal brittle normal faults (detachment faults) enhanced the uplift and tectonic denudation of this mylonite zone, ultimately resulting in southwestward displacement of the upper crustal stratigraphy.

Strains associated with the two ductile deformation events have been successfully partitioned through a multifaceted analysis. R_f/Ø measurements on various markers from the "type" stratigraphy allow a gradient representing cumulative strain since Middle Cretaceous time to be determined. From this gradient, noncoaxial strains accrued since emplacement of the leucogranites may be removed. Irrotational components of the postleucogranite strain are measured from quartz grain shapes in deformed granites; rotational components (shear strains) are determined from S-C fabrics and from restoration of rotated dike and vein networks. Structural observations and strain data are compatable with a deformation path of: (1) coaxial strain (pure shear?), followed by (2) injection of leucogranites as dikes (perpendicular to the minimum principle stress) and sills (parallel to the minimum principle stress), then (3) southwest directed simple shear. Modeling the late strain gradient as a simple shear zone permits a minimum displacement of 10 kilometers on the Magdalena mylonite zone/detachment fault system. Removal of the Middle Tertiary noncoaxial strains yields a residual (or pre-existing) strain gradient representative of the Late Cretaceous-Early Tertiary deformation. Several partially destrained cross sections, restored to the time of leucogranite emplacement, illustrate the idea that the upper plate of the core complex bas been detached from a region of significant topographic relief. 50% to 100% bulk extension across a 50 kilometer wide corridor is demonstrated.

Late Cenozoic tectonics of the Magdalena region are dominated by Basin and Range style faulting. Northeast and north-northwest trending high angle normal faults have interacted to extend the crust in an east-west direction. Net extension for this period is minor (10% to 15%) in comparison to the Middle Tertiary detachment related extensional episode.