Quantum Theory of Advanced Materials

A small amount of catalyst, such as Ti, was found to greatly improve the kinetics of hydrogen reactions in the prototypical hydrogen storage compound sodium alanate (NaAlH(4)). We propose a near-surface alloying mechanism for the rehydrogenation cycle based on a detailed analysis of available experimental data as well as first-principles calculations. The calculated results indicate that the catalyst remains at subsurface sites near the Al surface, reducing the dissociation energy barrier of H(2). The binding between Ti and Al modifies the surface charge distribution, which facilitates hydrogen adsorption and enhances hydrogen mobility on the surface.

It has been demonstrated that replacing Li(2)NH with the mixed imide Li(2)Mg(NH)(2) improves the reaction conditions for the hydrogen-storage system Li(2)NH + H(2) <-> LiNH(2) + LiH, at the expense of reducing the gravitational hydrogen capacity from 6.5% to 5.6%. In this article, we report from first-principles calculations a possible mixed imide Li(6)Mg(NH)(4) that has less Mg concentration and higher gravimetric capacity for hydrogen storage than Li(2)Mg(NH)(2). We find that Li(6)Mg(NH)(4) is thermodynamically more stable than the phase-separated mixture of Li(2)Mg(NH)(2) and Li(2)NH over a large temperature range. The reaction LiH + 1/4Mg(NH(2))(2) + 1/2LiNH(2) <-> 1/4Li(6)Mg(NH)(4) + H(2) can be completed via two steps and releases 6.0 wt % hydrogen in total, at a temperature about 40 degrees C lower than that for the cycling between LiNH(2) and Li(2)NH.

It has been suggested that the diffusion of AlH(3) vacancies plays an essential role in the decomposition of NaAlH(4), a prototypical material for hydrogen storage. We find from first-principles calculations that the AlH(3) vacancy induces several isolated vibrational modes that are highly localized in the vacancy region with frequencies within the phonon gaps of pure NaAlH(4) in both the a and g phases. Thus, the proposed existence of AlH(3) vacancies in the dehydrogenation reaction of NaAlH(4) can be possibly confirmed with the experimental detection of these unique vibrational modes associated with the AlH(3) vacancy.

Lithium imide (Li2NH) has been considered as a promising medium for hydrogen storage with the following reaction: LiNH(2)+LiH <-> Li(2)NH+H(2). All possible phases involved in the reaction need to be fully characterized in order to understand the right pathway connecting the two end compounds LiNH(2) and Li(2)NH and to further improve its reaction condition to meet the requirements of practical applications. We study from first-principles calculations the possible intermediate compounds Li(2-x)NH(1+x) between Li(2)NH and LiNH(2). Based on the energetics results, possible intermediate phases are identified for 0

Using first-principles calculations within density functional theory, we have investigated the electronic properties of H-passivated Si nanowires (SiNWs) oriented along the 112 direction, with the atomic geometries retrieved via global search using genetic algorithm. We show that [112] SiNWs have an indirect band gap in the ultrathin diameter regime, whereas the energy difference between the direct and indirect fundamental band gaps progressively decreases as the wire size increases, indicating that larger [112] SiNWs could have a quasi-direct band gap. We further show that this quasi-direct gap feature can be enhanced when applying uniaxial compressive stress along the wire axis. Moreover, our calculated results also reveal that the electronic band structure is sensitive to the change of the aspect ratio of the cross sections.

The electronic structure of Si/Ge core-shell nanowires along the [110] and [111] directions are studied with first-principles calculations. We identify the near-gap electronic states that are spatially separated within the core or the shell region, making it possible for a dopant to generate carriers in a different region. The confinement energies of these core and shell states provide an operational definition of the "band offset," which is not only size dependent but also component dependent. The optimal doping strategy in Si/Ge core-shell nanowires is proposed based on these energy results.

Aluminum hydride (alane) AlH(3) is an important material in hydrogen storage applications. It is known that AlH(3) exists in multiply forms of polymorphs, where alpha-AlH(3) is found to be the most stable with a hexagonal structure. Recent experimental studies on gamma-AlH(3) reported an orthorhombic structure with a unique double-bridge bond between certain Al and H atoms. This was not found in alpha-AlH(3) or other polymorphs. Using density functional theory, we have investigated the energetics, and the structural, electronic, and phonon vibrational properties for the newly reported gamma-AlH(3) structure. The current calculation concludes that gamma-AlH(3) is less stable than alpha-AlH(3) by 1.2 KJ/mol, with the zero-point energy included. Interesting binding features associated with the unique geometry of gamma-AlH(3) are discussed from the calculated electronic properties and phonon vibrational modes. The binding of H-s with higher energy Al-p,d orbitals is enhanced within the double-bridge arrangement, giving rise to a higher electronic energy for the system. Distinguishable new features in the vibrational spectrum of gamma-AlH(3) were attributed to the double-bridge and hexagonal-ring structures.

We report first-principles calculations of Pb (100) films up to 22 monolayers to study variations in the surface energy and work function as a function of film thickness. An even-odd oscillation is found in these two quantities, while a jelliumlike model for this s-p metal predicts a periodicity of about three monolayers. This unexpected result is explained by considering a coherent superposition of contributions from quantum-well states centered at both the Gamma and M points in the two-dimensional Brillouin zone, demonstrating the importance of crystal band structure in studying the quantum size effect in metal thin films.

Recently it was discovered that a total of 5.6 wt. % H-2 could be released from the 1:2 mixture of lithium amide and magnesium hydride at temperatures as low as 150 degrees C. With a reaction enthalpy of 44 KJ/mol H-2, this system has high potential for on-board hydrogen storage applications. The fully desorbed product is believed to be a mixed lithium and magnesium imide Li2Mg(NH)(2). In this work, the crystal structure of this mixed imide is studied from total-energy density-functional calculations. Based on a recent experimentally established space group, possible ordered configurations are examined. Important local orderings are identified for the experimentally observed disordered phase at room temperature. These unique local arrangements are also connected with the observed structural transitions above room temperature. In addition, the local ordering in Mg(NH2)(2) is analyzed. The similarity and difference of local arrangements among hydrogen, cations, and vacancies are discussed for the three amide (imide) systems: LiNH2, Mg(NH2)(2), and Li2Mg(NH)(2). The identification of the cation and hydrogen local orderings are expected to facilitate the design of new mixed imides and amides as hydrogen storage materials with desired physical properties.

Hydrogenation-induced metal-semiconductor transitions usually occur in simple systems based on rare earths and/or magnesium, accompanied by major reconstructions of the metal host (atom shifts >2 Angstrom). We report on the first such transition in a quaternary system based on a transition element. Metallic LaMg2Ni absorbs hydrogen near ambient conditions, forming the nonmetallic hydride LaMg2NiH7 which has a nearly unchanged metal host structure (atom shifts <0.7 Angstrom). The transition is induced by a charge transfer of conduction electrons into tetrahedral [NiH4](4-) complexes having closed-shell electron configurations.

We investigate the structural, electronic, and optical properties of hydrogen-passivated silicon nanowires along [110] and [111] directions with diameter d up to 4.2 nm from first principles. The size and orientation dependence of the band gap is investigated and the local-density gap is corrected with the GW approximation. Quantum confinement becomes significant for d<2.2 nm, where the dielectric function exhibits strong anisotropy and new low-energy absorption peaks start to appear in the imaginary part of the dielectric function for polarization along the wire axis.

Atomically uniform Pb films are successfully prepared on Si(111), despite a large lattice mismatch. Angle-resolved photoemission measurements of the electronic structure show layer-resolved quantum well states which can be correlated with dramatic variations in thermal stability. The odd film thicknesses N=5, 7, and 9 monolayers show sharp quantum well states. The even film thicknesses N=6 and 8 do not, but are much more stable than the odd film thicknesses. This correlation is discussed in terms of a total energy calculation and Friedel-like oscillations in properties.

We have studied the properties of quantum well states in supported Ag(100) films on the Fe substrate by first-principles density-functional calculations. The energies of these quantum well states as a function of thickness N are examined in terms of the characteristic phase shift of the electronic wave function at the interface. These energy-dependent phase shifts are determined numerically for both the film-substrate and film-vacuum interfaces. It is also found that the substrate has a major effect on film stability, enhancing the stability of the N=5 film and reversing that of the N=2 film.

Real-time in situ x-ray studies of continuous Pb deposition on Si(111)-(7x7) at 180 K reveal an unusual growth behavior. A wetting layer forms first to cover the entire surface. Then islands of a fairly uniform height of about five monolayers form on top of the wetting layer and grow to fill the surface. The growth then switches to a layer-by-layer mode upon further deposition. This behavior of alternating layer and island growth can be attributed to spontaneous quantum phase separation based on a first-principles calculation of the system energy.

The work function of atomically uniform Ag films grown on Fe(100) is measured as a function of film thickness. It shows layer-resolved variations as a result of quantum confinement of the valence electrons. A first-principles calculation reproduces the observed variations except for very thin films (one and two monolayers), and the differences can be attributed, in part, to strain effects caused by the lattice mismatch between Ag and Fe. These results illustrate the close interaction between interface effects and surface properties.

We have carried out first-principles calculations of Pb (111) films up to 25 monolayers to study the oscillatory quantum size effects exhibited in the surface energy and work function. These oscillations are correlated with the thickness dependence of the energies of confined electrons, which can be properly modeled by an energy-dependent phase shift of the electronic wave function upon reflection at the interface. It is found that a quantitative description of these quantum size effects requires a full consideration of the crystal band structure.

A combination of the coupling constant integration technique and the quantum Monte Carlo method is used to investigate the most relevant quantities in Kohn-Sham density-functional theory. Variational quantum Monte Carlo is used to construct realistic many-body wave functions for diamond-structure silicon at different values of the Coulomb coupling constant. The exchange-correlation energy density along with the coupling constant dependence and the coupling-constant-integrated form of the pair-correlation function, the exchange-correlation hole, and the exchange-correlation energy are presented. Comparisons of these functions an mode with results obtained from the local-density approximation, the average density approximation, the weighted density approximation, and the generalized gradient approximation. We discuss reasons for the success of the local-density approximation. The insights provided by this approach will make it possible to carry out stringent tests of the effectiveness of exchange-correlation functionals and in the long term aid in the search for better functionals. [S0163-1829(98)02115-8].

A model interaction is introduced for quantum many-body simulations of Coulomb systems using periodic I boundary conditions. The interaction gives much smaller finite size effects than the standard Ewald interaction and is also much faster to compute. Variational quantum Monte Carlo simulations of diamond-structure silicon with up to 1000 electrons demonstrate the effectiveness of our method.

Realistic many-body wave functions for diamond-structure silicon are constructed for different values of the Coulomb coupling constant. The coupling-constant-integrated pair correlation function, the exchange-correlation hole, and the exchange-correlation energy density are calculated and compared with those obtained from the local density and average density approximations. We draw conclusions about the reasons for the success of the local density approximation and suggest a method for testing the effectiveness of exchange-correlation functionals.

We studied total energies of various ordered structures of PdHx (in which hydrogen occupies the octahedral sites within the fee Pd lattice) using the pseudopotential method and a plane-wave basis within the local-density-functional approximation. The structures considered include the (420)-plane ordering of hydrogen atoms at different concentrations. For x greater than or equal to 1/2 we found that the NiMo- and Ni4Mo (D1(a))-type structures at x=1/2 and x=4/5, respectively, were energetically favored phases, in agreement with the superlattice reflections found in previous neutron-scattering measurements. For the intermediate concentrations, linear variation of the formation energy as a function of x in several (420)-ordered structures explained the observed short-range order. In contrast to an earlier proposal, we did not find the Fermi surface imaging effect responsible in this case. The overall energy variation in different phases indicates the importance of going beyond pairwise interactions between interstitial hydrogen atoms in this system.

We report the first implementation of orthonormal wavelet bases in self-consistent electronic structure calculations within the local-density approximation. These local bases of different scales efficiently describe localized orbitals of interest. As an example, we studied two molecules, H-2 and O-2, using pseudopotentials and supercells. Considerably fewer bases are needed compared with conventional plane-wave approaches, yet calculated binding properties are similar. Our implementation employs fast wavelet and Fourier transforms, avoiding evaluating any three-dimensional integral numerically.