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.

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.