Kinetic investigations of heterogeneously catalyzed reactions on the Ru(001) and Pt(110)-(1x2) surfaces

The initial probabilities of activated, dissociative chemisorption of methane and ethane on Pt(110)-(1 x 2) have been measured. The surface temperature was varied from 450 to 900 K with the reactant gas temperature constant at 300 K. Under these conditions, we probe the kinetics of dissociation via trapping-mediated (as opposed to 'direct') mechanism. It was found that the probabilities of dissociation of both methane and ethane were strong functions of the surface temperature with an apparent activation energies of 14.4 kcal/mol for methane and 2.8 kcal/mol for ethane, which implys that the methane and ethane molecules have fully accommodated to the surface temperature. Kinetic isotope effects were observed for both reactions, indicating that the C-H bond cleavage was involved in the rate-limiting step. A mechanistic model based on the trapping-mediated mechanism is used to explain the observed kinetic behavior. The activation energies for C-H bond dissociation of the thermally accommodated methane and ethane on the surface extracted from the model are 18.4 and 10.3 kcal/mol, respectively.

The studies of the catalytic decomposition of formic acid on the Ru(001) surface with thermal desorption mass spectrometry following the adsorption of DCOOH and HCOOH on the surface at 130 and 310 K are described. Formic acid (DCOOH) chemisorbs dissociatively on the surface via both the cleavage of its O-H bond to form a formate and a hydrogen adatom, and the cleavage of its C-O bond to form a carbon monoxide, a deuterium adatom and an hydroxyl (OH). The former is the predominant reaction. The rate of desorption of carbon dioxide is a direct measure of the kinetics of decomposition of the surface formate. It is characterized by a kinetic isotope effect, an increasingly narrow FWHM, and an upward shift in peak temperature with Ɵ_T, the coverage of the dissociatively adsorbed formic acid. The FWHM and the peak temperature change from 18 K and 326 K at Ɵ_T = 0.04 to 8 K and 395 K at Ɵ_T = 0.89. The increase in the apparent activation energy of the C-D bond cleavage is largely a result of self-poisoning by the formate, the presence of which on the surface alters the electronic properties of the surface such that the activation energy of the decomposition of formate is increased. The variation of the activation energy for carbon dioxide formation with Ɵ_T accounts for the observed sharp carbon dioxide peak. The coverage of surface formate can be adjusted over a relatively wide range so that the activation energy for C-D bond cleavage in the case of DCOOH can be adjusted to be below, approximately equal to, or well above the activation energy for the recombinative desorption of the deuterium adatoms. Accordingly, the desorption of deuterium was observed to be governed completely by the desorption kinetics of the deuterium adatoms at low Ɵ_T, jointly by the kinetics of deuterium desorption and C-D bond cleavage at intermediate Ɵ_T, and solely by the kinetics of C-D bond cleavage at high Ɵ_T. The overall branching ratio of the formate to carbon dioxide and carbon monoxide is approximately unity, regardless the initial coverage Ɵ_T, even though the activation energy for the production of carbon dioxide varies with Ɵ_T. The desorption of water, which implies C-O bond cleavage of the formate, appears at approximately the same temperature as that of carbon dioxide. These observations suggest that the cleavage of the C-D bond and that of the C-O bond of two surface formates are coupled, possibly via the formation of a short-lived surface complex that is the precursor to to the decomposition.

The measurement of steady-state rate is demonstrated here to be valuable in determining kinetics associated with short-lived, molecularly adsorbed precursor to further reactions on the surface, by determining the kinetic parameters of the molecular precursor of formaldehyde to its dissociation on the Pt(110)-(1 x 2) surface.

Overlayers of nitrogen adatoms on Ru(001) have been characterized both by thermal desorption mass spectrometry and low-energy electron diffraction, as well as chemically via the postadsorption and desorption of ammonia and carbon monoxide.

The nitrogen-adatom overlayer was prepared by decomposing ammonia thermally on the surface at a pressure of 2.8 x 10^(-6) Torr and a temperature of 480 K. The saturated overlayer prepared under these conditions has associated with it a (√247/10 x √247/10)R22.7° LEED pattern, has two peaks in its thermal desorption spectrum, and has a fractional surface coverage of 0.40. Annealing the overlayer to approximately 535 K results in a rather sharp (√3 x √3)R30° LEED pattern with an associated fractional surface coverage of one-third. Annealing the overlayer further to 620 K results in the disappearance of the low-temperature thermal desorption peak and the appearance of a rather fuzzy p(2x2) LEED pattern with an associated fractional surface coverage of approximately one-fourth. In the low coverage limit, the presence of the (√3 x √3)R30° N overlayer alters the surface in such a way that the binding energy of ammonia is increased by 20% relative to the clean surface, whereas that of carbon monoxide is reduced by 15%.

A general methodology for the indirect relative determination of the absolute fractional surface coverages has been developed and was utilized to determine the saturation fractional coverage of hydrogen on Ru(001). Formaldehyde was employed as a bridge to lead us from the known reference point of the saturation fractional coverage of carbon monoxide to unknown reference point of the fractional coverage of hydrogen on Ru(001), which is then used to determine accurately the saturation fractional coverage of hydrogen. We find that ƟSAT/H = 1.02 (±0.05), i.e., the surface stoichiometry is Ru : H = 1 : 1. The relative nature of the method, which cancels systematic errors, together with the utilization of a glass envelope around the mass spectrometer, which reduces spurious contributions in the thermal desorption spectra, results in high accuracy in the determination of absolute fractional coverages.

Quantum Mechanics Studies of Fuel Cell Catalysts and Proton Conducting Ceramics with Validation by Experiment

We carried out quantum mechanics (QM) studies aimed at improving the performance of hydrogen fuel cells. This led to predictions of improved materials, some of which were subsequently validated with experiments by our collaborators.

In part I, the challenge was to find a replacement for the Pt cathode that would lead to improved performance for the Oxygen Reduction Reaction (ORR) while remaining stable under operational conditions and decreasing cost. Our design strategy was to find an alloy with composition Pt3M that would lead to surface segregation such that the top layer would be pure Pt, with the second and subsequent layers richer in M. Under operating conditions we expect the surface to have significant O and/or OH chemisorbed on the surface, and hence we searched for M that would remain segregated under these conditions. Using QM we examined surface segregation for 28 Pt₃M alloys, where M is a transition metal. We found that only Pt₃Os and Pt₃Ir showed significant surface segregation when O and OH are chemisorbed on the catalyst surfaces. This result indicates that Pt₃Os and Pt3Ir favor formation of a Pt-skin surface layer structure that would resist the acidic electrolyte corrosion during fuel cell operation environments. We chose to focus on Os because the phase diagram for Pt-Ir indicated that Pt-Ir could not form a homogeneous alloy at lower temperature. To determine the performance for ORR, we used QM to examine all intermediates, reaction pathways, and reaction barriers involved in the processes for which protons from the anode reactions react with O₂ to form H₂O. These QM calculations used our Poisson-Boltzmann implicit solvation model include the effects of the solvent (water with dielectric constant 78 with pH 7 at 298K). We found that the rate determination step (RDS) was the O_ad hydration reaction (O_ad + H₂O_ad -> OH_ad + OH_ad) in both cases, but that the barrier for pure Pt of 0.50 eV is reduced to 0.48 eV for Pt₃Os, which at 80 degrees C would increase the rate by 218%. We collaborated with the Pu-Wei Wu’s group to carry out experiments, where we found that the dealloying process-treated Pt2Os catalyst showed two-fold higher activity at 25 degrees C than pure Pt and that the alloy had 272% improved stability, validating our theoretical predictions.

We also carried out similar QM studies followed by experimental validation for the Os/Pt core-shell catalyst fabricated by the underpotential deposition (UPD) method. The QM results indicated that the RDS for ORR is a compromise between the OOH formation step (0.37 eV for Pt, 0.23 eV for Pt_2ML/Os core-shell) and H₂O formation steps (0.32 eV for Pt, 0.22 eV for Pt_2ML/Os core-shell). We found that Pt_2ML/Os has the highest activity (compared to pure Pt and to the Pt₃Os alloy) because the 0.37 eV barrier decreases to 0.23 eV. To understand what aspects of the core shell structure lead to this improved performance, we considered the effect on ORR of compressing the alloy slab to the dimensions of pure Pt. However this had little effect, with the same RDS barrier 0.37 eV. This shows that the ligand effect (the electronic structure modification resulting from the Os substrate) plays a more important role than the strain effect, and is responsible for the improved activity of the core- shell catalyst. Experimental materials characterization proves the core-shell feature of our catalyst. The electrochemical experiment for Pt_2ML/Os/C showed 3.5 to 5 times better ORR activity at 0.9V (vs. NHE) in 0.1M HClO₄ solution at 25 degrees C as compared to those of commercially available Pt/C. The excellent correlation between experimental half potential and the OH binding energies and RDS barriers validate the feasibility of predicting catalyst activity using QM calculation and a simple Langmuir–Hinshelwood model.

In part II, we used QM calculations to study methane stream reforming on a Ni-alloy catalyst surfaces for solid oxide fuel cell (SOFC) application. SOFC has wide fuel adaptability but the coking and sulfur poisoning will reduce its stability. Experimental results suggested that the Ni4Fe alloy improves both its activity and stability compared to pure Ni. To understand the atomistic origin of this, we carried out QM calculations on surface segregation and found that the most stable configuration for Ni₄Fe has a Fe atom distribution of (0%, 50%, 25%, 25%, 0%) starting at the bottom layer. We calculated that the binding of C atoms on the Ni4Fe surface is 142.9 Kcal/mol, which is about 10 Kcal/mol weaker compared to the pure Ni surface. This weaker C binding energy is expected to make coke formation less favorable, explaining why Ni₄Fe has better coking resistance. This result confirms the experimental observation. The reaction energy barriers for CHx decomposition and C binding on various alloy surface, Ni₄X (X=Fe, Co, Mn, and Mo), showed Ni₄Fe, Ni₄Co, and Fe₄Mn all have better coking resistance than pure Ni, but that only Ni₄Fe and Fe₄Mn have (slightly) improved activity compared to pure Ni.

In part III, we used QM to examine the proton transport in doped perovskite-ceramics. Here we used a 2x2x2 supercell of perovskite with composition Ba₈X₇M₁(OH)₁O₂₃ where X=Ce or Zr and M=Y, Gd, or Dy. Thus in each case a 4⁺ X is replace by a 3⁺ M plus a proton on one O. Here we predicted the barriers for proton diffusion allowing both includes intra-octahedron and inter-octahedra proton transfer. Without any restriction, we only observed the inter-octahedra proton transfer with similar energy barrier as previous computational work but 0.2 eV higher than experimental result for Y doped zirconate. For one restriction in our calculations is that the O_donor-O_acceptor atoms were kept at fixed distances, we found that the barrier difference between cerates/zirconates with various dopants are only 0.02~0.03 eV. To fully address performance one would need to examine proton transfer at grain boundaries, which will require larger scale ReaxFF reactive dynamics for systems with millions of atoms. The QM calculations used here will be used to train the ReaxFF force field.

Performance prediction of graphene nanoribbon and carbon nanotube transistors

Carbon nanotubes (CNTs) and graphene nanoribbons (GNRs) field-effect transistor (FET) can be the basis for a quasi-one- dimensional (Q1D) transistor technology. Recent experiments show that the on-off ratio for GNR devices can be improved to level exploration of transistor action is justified. Here we use the tight-binding energy dipersion approximation, to assess the performance of semiconducting CNT and GNR is qualitatively in terms of drain current drive strength, bandgap and density of states for a specified device. By reducing the maximum conductance 4e2/h by half, we observed that our model has a particularly good fit with 50 nm channel single walled carbon nanotube (SWCNT) experimental data. Given the same bandgap, CNTs outperform GNRs due to valley degeneracy. Nevertheless, the variation of the device contacts will decide which transistor will exhibit better conductivity and thus higher ON currents. © 2011 American Institute of Physics.

Electronic structure of a hydrogenic acceptor impurity in semiconductor nano-structures

The electronic structure and binding energy of a hydrogenic acceptor impurity in 2, 1, and 0-dimensional semiconductor nano-structures (i.e. quantum well (QW), quantum well wire (QWW), and quantum dot (QD)) are studied in the framework of effective-mass envelope-function theory. The results show that (1) the energy levels monotonically decrease as the quantum confinement sizes increase; (2) the impurity energy levels decrease more slowly for QWWs and QDs as their sizes increase than for QWs; (3) the changes of the acceptor binding energies are very complex as the quantum confinement size increases; (4) the binding energies monotonically decrease as the acceptor moves away from the nano-structures' center; (5) as the symmetry decreases, the degeneracy is lifted, and the first binding energy level in the QD splits into two branches. Our calculated results are useful for the application of semiconductor nano-structures in electronic and photoelectric devices.

Quantum-dot growth simulation on periodic stress of substrate

InAs quantum dots (QDs) are grown on the cleaved edge of an InxGa1-xAs/GaAs supperlattice experimentally and a good linear alignment of these QDs on the surface of an InxGa1-xAs layer has been realized. The modulation effects of periodic strain on the substrate are investigated theoretically using a kinetic Monte Carlo method. Our results show that a good alignment of QDs can be achieved when the strain energy reaches 2% of the atomic binding energy. The simulation results are in excellent qualitative agreement with our experiments. (C) 2005 American Institute of Physics.

Photoluminescence degradation in GaN induced by light enhanced surface oxidation

The exponential degradation of the photoluminescence (PL) intensity at the near-band-gap was observed in heavily doped or low-quality GaN with pristine surface under continuous helium-cadmium laser excitation. In doped GaN samples, the degradation speed increased with doping concentration. The oxidation of the surface with laser irradiation was confirmed by x-ray photoemission spectroscopy measurements. The oxidation process introduced many oxygen impurities and made an increase of the surface energy band bending implied by the shift of Ga 3d binding energy. The reason for PL degradation may lie in that these defect states act as nonradiative centers and/or the increase of the surface barrier height reduces the probability of radiative recombination.

Electronic states of a hydrogenic donor impurity in semiconductor nano-structures

We propose a method for uniformly calculating the electronic states of a hydrogenic donor impurity in low-dimensional semiconductor nano-structures in the framework of effective-mass envelope-function theory, and we study the electronic structures of this systems. Compared to previous methods, our method has the following merits: (a) It can be widely applied in the calculation of the electronic states of hydrogenic donor impurities in nano-structures of various shapes; (b) It can easily be extended to study the effects of external fields and other complex cases; (c) The excited states are more easily calculated than with the variational method; (d) It is convenient to calculate the change of the electronic states with the position of a hydrogenic donor impurity in nano-structures; (e) The binding energy can be calculated explicitly. (c) 2007 Elsevier B.V. All rights reserved.

Hydrogenic impurity in double quantum dots

The ground state binding energy and the average interparticle distances for a hydrogenic impurity in double quantum dots with Gaussian confinement potential are studied by the variational method. The probability density of the electron is calculated, too. The dependence of the binding energy on the impurity position is investigated for GaAs quantum dots. The result shows that the binding energy has a minimum as a function of the distance between the two quantum dots when the impurity is located at the center of one quantum dot or at the center of the edge of one quantum dot. When the impurity is located at the center of the two dots, the binding energy decreases monotonically. (c) 2006 Elsevier B.V. All rights reserved.

The Stark effect on excitons in a bilayer system

The Stark effect on excitons in a bilayer system is investigated theoretically within the framework of the effective-mass approximation. The calculations indicate that the energy of the excitons decreases as the value of the in-plane electric field F increases at a fixed value of the distance d between the layers. However, the energy of the excitons increases with d at a fixed value of F. In particular, it increases linearly at small values of d but increases as 1/d at large values. Therefore, it can be concluded that excitons in a bilayer system have a small binding energy equal to the absolute value of the excitonic energy at large d or small F. In addition, the radiative lifetime of heavy-hole excitons in this system is calculated and is found to be short at small values of both F and d. The radiative lifetime of heavy-hole excitons in a bilayer system can be increased by two orders by an in-plane electric field of 2 kV/cm when d is twice the excitonic Rydberg. (c) 2006 American Institute of Physics.

Evolution of the amount of InAs in wetting layers in an InAs/GaAs quantum-dot system studied by reflectance difference spectroscopy

The wetting layer (WL) in InAs/GaAs quantum-dot systems has been studied by reflectance difference spectroscopy (RDS). Two structures related to the heavy-hole (HH) and light-hole (LH) related transitions in the WL have been observed. On the basis of a calculation model that takes into account the segregation effect and exciton binding energies, the amount of InAs in the WL (t(WL)) and its segregation coefficient ( R) have been determined from the HH and LH transition energies. The evolutions of tWL and R exhibit a close relation to the growth modes. Before the formation of InAs dots, t(WL) increases linearly from similar to 1 to similar to 1.6 monolayer (ML), while R increases almost linearly from similar to 0.8 to similar to 0.85. After the onset of dot formation, t(WL) is saturated at similar to 1.6 ML and R decreases slightly from 0.85 to 0.825. The variation of tWL can be interpreted by using an equilibrium model. Different variations of in-plane optical anisotropy before and after dot formation have been observed.

Lifetime study of N impurity states in GaAs1-xNx (x=0.1%) under hydrostatic pressure

The lifetimes of a series of N-related photoluminescence lines (A(2)-A(6)) in GaAs1-xNx (x=0.1%) were studied under hydrostatic pressures at similar to 30 K. The lifetimes of A(5) and A(6) were found to increase rapidly with increasing pressure: from 2.1 ns at 0 GPa to more than 20 ns at 0.92 GPa for A(5) and from 3.2 ns at 0.63 GPa to 10.8 ns at 0.92 GPa for A(6). The lifetime is found to be closely correlated with the binding energy of the N impurity states, which is shown either in the pressure dependence for a given emission line or in the lifetime variation from A(2) to A(6). (c) 2006 American Institute of Physics.

Photoluminescence from ZnS1-xTex alloys under hydrostatic pressure

ZnS1-xTex (0.02less than or equal toxless than or equal to0.3) alloys are studied by photoluminescence under hydrostatic pressure at room temperature. Only a wide emission band is observed for each sample. Its peak energy is much lower than the corresponding band gap of alloys. These bands are ascribed to the radiative annihilation of excitons bound at Te-n(ngreater than or equal to2) isoelectronic centers. The pressure coefficients of the emission bands are smaller than those of alloy band gaps from 48% to 7%. The difference of the pressure coefficient of the emission band and the band gap increases when the binding energy of Te-n centers decreases. It seems contrary to our expectation and needs further analysis. The integrated intensities of emission bands decrease with increasing pressure due to the decreasing of the absorption coefficient associated with the Te-n centers under pressure. According to this model the Stokes shifts between the emission and absorption bands of the Te-n centers are calculated, which decrease with the increasing Te composition in alloys.

Type I-type II transition of self-assembled In0.55Al0.45As/Al0.5Ga0.5As quantum dots

The photoluminescence (PL) of In0.55Al0.45As/Al0.5Ga0.5As self-assembled quantum dots has been measured at 15 and 80 K under hydrostatic pressure. The lateral size of the dots ranges from 7 to 62 nm. The emissions from the dots with 26, 52 and 62 nm size have a blue shift under pressure, indicating that these quantum dots have the normal type-I structure with lowest conduction band at the Gamma -valley. However, the PL peak of dots with 7 nm diameter moves to lower energy with increasing pressure. It is a typical character for the X-related transition. Then these small dots have a type-II structure with the X-valley as the lowest conduction level. An envelope-function calculation confirms that the Gamma -like exciton transition energy will rise above the X-like transition energy in the In0.55Al0.45As/Al0.5Ga0.5As structure if the dot size is small enough.

GAAS/GAALAS GRADED INDEX SEPARATE CONFINEMENT SINGLE QUANTUM-WELL SINGLE-MODE WAVE-GUIDE ELECTROABSORPTION LIGHT-MODULATOR

A GaAs/GaAlAs graded-index separate confinement single quantum well heterostructure single-mode ridge waveguide electroabsorption modulator was fabricated and investigated. For the modulator with a quantum well width of 100 angstrom and device length of 700-mu-m, an on/off ratio of 29.7 dB and estimated absorption insertion loss of 3 dB were obtained for TE polarised light with wavelength 8650 angstrom, and for TM polarisation the on/off ratio was 28.5 dB. With a switching voltage of 1 V, an on/off ratio of 15 dB was achieved. Photocurrent spectra exhibited a red shift of 600 angstrom of the absorption edge when the voltage applied to the PIN diode was varied from 0.5 to -7 V. The corresponding shift of the room temperature exciton peak energy was 96 meV.