Quantum Theory of Advanced Materials

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.

We have used the combination of the coupling-constant integration procedure and the variational quantum Monte Carlo method to study the exchange-correlation (XC) interaction in small molecules: Si-2, C2H2, C2H4, and C2H6. In this paper we report the calculated XC energy density, a central quantity in density functional theory, as deduced from the interaction between the electron and its XC hole integrated over the interaction strength. Comparing these "exact" XC energy densities with results using the local-density approximation (LDA), one can analyze the errors in this widely used approximation. Since the XC energy is an integrated quantity, error cancellation among the XC energy density in different regions is possible. Indeed we find a general error cancellation between the high-density and low-density regions. Moreover, the error distribution of the exchange contribution is out of phase with the error distribution of the correlation contribution. Similar to what is found for bulk silicon and an isolated silicon atom, the spatial variation of the errors of the LDA XC energy density in these molecules largely follows the sign and shape of the Laplacian of the electron density. Some noticeable deviations are found in Si-2 in which the Laplacian peaks between the atoms, while the LDA error peaks in the regions "behind" atoms where a good portion of the charge density originates from an occupied 1 sigma(u) antibonding orbital. Our results indicate that, although the functional form could be quite complex, an XC energy functional containing the Laplacian of the energy is a promising possibility for improving LDA.

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.

We present a detailed study of the exchange-correlation hole and exchange-correlation energy per particle in the Si crystal as calculated by the variational Monte Carlo method and predicted by various density-functional models. Nonlocal density-averaging methods prove to be successful in correcting severe errors in the local-density approximation (LDA) at low densities where the density changes dramatically over the correlation length of the LDA hole. but fail to provide systematic improvements at higher densities where the effects of density inhomogeneity are more subtle. Exchange and correlation considered separately show a sensitivity to the nonlocal semiconductor-crystal environment, particularly within the Si bond. which is not predicted by the nonlocal approaches based on density averaging. The exchange hole is well described by a bonding-orbital picture, while the correlation hole has a significant component due to the polarization of the nearby bonds, which partially screens out the anisotropy in the exchange hole.

We have studied the pair-correlation function, the exchange-correlation hole, and the exchange-correlation energy density of the valence electrons in the Si atom using the Coulomb-coupling constant integration technique with the variational quantum Monte Carlo method. These quantities are compared to those derived from various approximate models within the Kohn-Sham density functional theory. We find that the charge density prefactor in the expression for the exchange-correlation hole dominates the errors found in the local spin density approximation (LSDA), that the generalized gradient approximation improves energy calculations by improving the LSDA at long ranges, and that the weighted spin density approximation, which uses the correct charge density prefactor, gives the lowest root mean square error for the exchange-correlation energy density.

The charge-density-wave transition in TiSe2, which results in a commensurate (2 X 2 X 2) superlattice at temperatures below similar to 200 K, presumably involves softening of a zone-boundary phonon mode. For the first time, this phonon-softening behavior has been examined over a wide temperature range by synchroton x-ray thermal diffuse scattering.

When disilane (Si2H6) is used in the homoepitaxial growth of Si by chemical vapor deposition (CVD), the fragment SiH2 is believed to be the basic unit adsorbed on the surface. The bonding site of SiH2 on Si(100) has been proposed in the literature to be either on top of a dimer (the on-dimer site) or between two dimers in the same row (the intrarow site). Since the pathway of SiH2 combination is dependent on the adsorption site, a first-principles calculation will shed light on the underlying process. We have performed self-consistent pseudopotential density-functional calculations within the local-density approximation. On the bare Si(100) surface, the on-dimer site is found to be more stable than the intrarow site, even though the former has unfavorable Si-Si bond angles. This is ascribed to the extra dangling bond created in the latter geometry when the weak dimer a bonds are broken. However, the presence of hydrogen adatoms eliminates this difference and makes the intrarow site more favorable than the on-dimer site. It is therefore revealed in this theoretical study that hydrogen, an impurity unavoidable in the CVD process, plays an important role in determining the stable configuration of adsorbed SiH2 on Si(100) and hence affects the growth mechanism. [S0163-1829(98)52544-1].

We have evaluated the surface free energies of hydrogen-covered (100), (111), and (110) surfaces of diamond and silicon as a function of the hydrogen chemical potential using first-principles methods. The change in surface-energy anisotropy and equilibrium crystal shape due to hydrogen adsorption is examined. The three low-index facets are affected differently by the presence of hydrogen and unexpected differences are found between diamond and silicon. Taking into account possible formation of local facets on the hydrogen-covered (100) surfaces, we find that the dihydride phase is not stable on both C(100) and Si(100). nor is the 3x1 phase on C(100).

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].

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.

The bare and hydrogen-covered diamond (100) surfaces were investigated through pseudopotential density-functional calculations within the local-density approximation. Different hydrogen coverages, ranging from one to two, were considered. These corresponded to different structures (1x1, 2x1, and 3x1) and different hydrogen-carbon arrangements (monohydride, dihydride, and configurations in between). Assuming the system was in equilibrium with a hydrogen reservoir, the formation energy of each phase was expressed as a function of hydrogen chemical potential. As the chemical potential increased, the stable phase successively changed from bare 2x1 to (2x1):H, to (3x1):1.33H, and finally to the canted (1x1):2H. Setting the chemical potential at the energy per hydrogen in H-2 and in a free atom gave the (3x1):1.33H and the canted (1x1):2H phase as the most stable one, respectively. However, after comparing with the formation energy of CH4, only the (2x1):H and (3x1):1.33H phases were stable against spontaneous formation of CH4. The former existed over a chemical potential range ten times wider than the latter, which may explain why the latter, despite having a low energy, has not been observed so far. Finally, the vibrational energies of the C-H stretch mode were calculated for the (2x1):H phase.