Quantum Theory of Advanced Materials

A p-T phase diagram of graphene nanoribbons (GNRs) terminated by hydrogen atoms is established based on first-principles calculations, where the stable phase at standard conditions (25 degrees C and 1 bar) is found to be a zigzag GNR (zzGNR). The stability of this new GNR is understood based on an electron-counting model, which predicts semiconducting nonmagnetic zzGNRs. Quantum confinement of Dirac fermions in the stable zzGNRs is found to be qualitatively different from that in ordinary semiconductors. Bifurcation of the band gap is predicted to take place, leading to the formation of polymorphs with distinct band gaps but equal thermodynamic stability. A tight-binding model analysis reveals the role of edge symmetry on the band-gap bifurcation.

Accurate multideterminant ground-state energies of circular quantum dots containing N <= 13 electrons as a function of interaction strength have been evaluated by the diffusion quantum Monte Carlo method. Two unique features are found for these confined two-dimensional systems: (1) as the electron density decreases, the quantum dots favor states with zero orbital angular momentum (L = 0); and (2) for some values of N, the ground state cannot be fully spin-polarized because of a symmetry constraint.

The physical and chemical properties of thin metal films show damped oscillations as a function of film thickness (one-dimensional shell effects). While the oscillation period, determined by subband crossings of the Fermi level, is the same for all properties, the phases can be different. Specifically, oscillations in the work function and surface energy are offset by 1/4 of a period. For Pb(111) films, this offset is similar to 0.18 monolayers, a seemingly very small effect. However, aliasing caused by the discrete atomic layer structure leads to striking out-of-phase beating patterns displayed by these two quantities.

The phonon dispersions of monolayer and few-layer graphene (AB bilayer, and ABA and ABC trilayers) are investigated using the density-functional perturbation theory. Compared with the monolayer, the optical phonon E(2g) mode at Gamma splits into two and three doubly degenerate branches for bilayer and trilayer graphene, respectively, due to the weak interlayer coupling. These modes are of various symmetries and exhibit different sensitivities to either Raman or infrared measurements (or both). The splitting is found to be 5 cm(-1) for bilayer and 2-5 cm(-1) for trilayer graphene. The interlayer coupling is estimated to be about 2 cm(-1). We found that the highest optical modes at K move up by about 12 cm(-1) for bilayer and 18 cm(-1) for trilayer relative to monolayer graphene. The atomic displacements of these optical eigenmodes are analyzed.

By using a low temperature scanning tunneling microscope we have probed the superconducting energy gap of epitaxially grown Pb films as a function of the layer thickness in an ultrathin regime (5-18 ML). The layer-dependent energy gap and transition temperature (T-c) show persistent quantum oscillations down to the lowest thickness without any sign of suppression. Moreover, by comparison with the quantum-well states measured above T-c and the theoretical calculations, we found that the T-c oscillation correlates directly with the density of states oscillation at E-F. The oscillation is manifested by the phase matching of the Fermi wavelength and the layer thickness, resulting in a bilayer periodicity modulated by a longer wavelength quantum beat.

Using classical mechanical and quantum Monte Carlo methods we trace the ground-state behavior with an applied magnetic field of localized electron pair states in a quantum dot. By developing a method to treat nonconserved paramagnetic interactions using variational and diffusion quantum Monte Carlo techniques we find (i) a single-triplet transition at very small magnetic field strengths, (ii) enhanced localization of the two electrons with increasing magnetic field, and (iii) a mechanism for pair breakup that is different from that proposed recently by Wan et al. [Phys. Rev. Lett. 75, 2879 (1995)]. [S0163-1829(98)04016-8].

The solid-solution phase of hydrogen in hexagonal close-packed yttrium (a-YH(x)) is studied using the pseudopotential method within the local-density-functional approximation with a plane-wave basis. The binding energies associated with different interstitial sites are evaluated for several ordered structures: YH0.5, YH0.25, and YH0.167. It is found that the occupation of the tetrahedral site is always energetically favorable. The hydrogen potential-energy curves around the tetrahedral sites along the c axis and along the path connecting the adjacent octahedral sites are also calculated for YH0.25. In particular, the local vibrational mode along the c axis is estimated to be 100 meV, in excellent agreement with that measured in neutron-scattering experiments. Finally, the intriguing pairing phenomenon is investigated by calculating the total energy for various pairing configurations. The possibility of pairing between nearest-neighbor tetrahedral sites is excluded due to the high energy. It is found that the pairing of hydrogen across a metal atom is indeed energetically favorable compared with other kinds of pairs considered and also with isolated tetrahedral hydrogen atoms. The connection with the electronic structure of the system is also examined.

We present the calculation of the full phonon spectrum for silicon and germanium with the pseudopotential method and the local-density approximation without using linear-response theory. The interplanar-force constants for three high-symmetry orientations [(100), (110), and (111)] are evaluated by supercell calculations using the Hellmann-Feynman theorem. By considering the symmetry of the crystal, three-dimensional interatomic-force-constant matrices are determined by a least-squares fit. Interactions up to the eighth nearest neighbors are included. The dynamical matrix, which is the Fourier transform of the force constant matrix, is hence constructed and diagonalized for any arbitrary wave vector in the Brillouin zone, yielding the phonon dispersion. In this paper we will present the calculation details and discuss various aspects of convergence. Phonon dispersions of Si and Ge calculated are in excellent agreement with experiments.

A pseudopotential local-density calculation is performed for YH3 to study the unusual hydrogen displacements previously found in neutron diffraction. These displacements are identified as Peierls distortions associated with (hydrogen) lattice instability in this 3D system. The wave vector of these displacements is close to the vector connecting the electron and hole pockets in the undistorted system. With other electron and hole pockets at GAMMA that still overlap after distortion, the possibility of the existence of an excitonic insulator phase will be discussed.

The structural properties of hexagonal-close-packed yttrium are studied by using the plane-wave basis within the pseudopotential method and local-density-functional approximation. By employing a "soft" pseudopotential proposed by Troullier and Martins, satisfactory convergence is achieved with a plane-wave energy cutoff of 30-40 Ry for this early-transition-metal element. The overall results for the structural properties are in good agreement with experiment. It is found that the charge overlap between core and valence electrons has a substantial effect on the accuracy of the calculated structural properties. Two different calculations are performed with and without the outer-core 4p orbital included as a valence state. In addition, as found in some other local-density calculations, the uncertainty in the results due to different exchange-correlation energy functionals may not be negligible in transition metals.