
Atom/Electron Papers bibTeX EntriesSee also our main Papers and Publications page, as well as the publications pages of the Fink, Johnson, and Soljacic groups. @article{JoannopoulosCo72, author = {J. D. Joannopoulos and Marvin L. Cohen}, title = {Comparison of the electronic structure of amorphous and crystalline polytopes of {Ge}}, journal = {Solid State Communications}, year = 1972, volume = 11, number = 4, pages = {549553}, month = {August}, doi = {doi:10.1016/00381098(72)901238}, abstract = {We have calculated the density of states and the imaginary part of the dielectric constant as a function of energy for four polytypes of Ge using the EPM and the tight binding model developed by Weaire. The increasing complexity of the crystal structures indicates that short range disorder is able to account well for the optical properties and density of states of amorphous Ge. Furthermore, we predict the form of the optical properties and density of states for Ge III and Ge IV.} } @article{JoannopoulosCo72PLA, author = {J. D. Joannopoulos and Marvin L. Cohen}, title = {Electronic structure of crystalline polytypes and amorphous {Si}}, journal = {Physics Letters~A}, year = 1992, volume = 41, number = 1, pages = {7172}, month = {August}, doi = {doi:10.1016/03759601(72)90638X}, abstract = {We have calculated the density of states and the imaginary part of the energy dependent dielectric constant for four polytypes of Si using pseudopotential and tight binding models. Comparisons are made with the experimental results for amorphous Si.} } @article{PollakLe72, author = {R. A. Pollak and L. Ley and S. Kowalczyk and D. A. Shirley and J. D. Joannopoulos and D. J. Chadi and Marvin L. Cohen}, title = {Xray photoemission valenceband spectra and theoretical valenceband densities of states for {Ge}, {GaAs}, and {ZnSe}}, journal = {Physical Review Letters}, year = 1972, volume = 29, pages = {11031105}, month = {October}, url = {http://link.aps.org/abstract/PRL/v29/p1103}, abstract = {The first highresolution valenceband xray photoemission spectra for Ge, GaAs, and ZnSe are compared with pseudopotential densityofstates calculations. The general agreement between the experimental and theoretical results is quite good. For ZnSe the photoemission spectrum shows the Zn $3d$ states to be higher in energy than the lowest valenceband s state. In order to obtain this ordering of states in the theoretical calculation, a pseudopotential with an explicit energy dependence is required.} } @article{JoannopoulosCo73I, author = {J. D. Joannopoulos and Marvin L. Cohen}, title = {Electronic properties of complex crystalline and amorphous phases of {Ge} and {Si}. {I}. Density of states and band structures}, journal = {Physical Review~B}, year = 1973, volume = 7, pages = {26442657}, month = {March}, url = {http://link.aps.org/abstract/PRB/v7/p2644}, abstract = {We present calculations of the band structures and densities of states of Ge and Si in the diamond, wurzite, SiIII (BC8), and GeIII (ST12) structures using the empiricalpseudopotential method and the tightbinding model used recently by Weaire. The increasing complexity of the crystal structures indicates that shortrange disorder is able to account well for the density of states and optical properties of amorphous Ge and Si. This calculation also provides a method for explaining various features in the amorphous density of states and shows what structural aspects of the amorphous state are responsible for these features.} } @article{JoannopoulosCo73II, author = {J. D. Joannopoulos and Marvin L. Cohen}, title = {Electronic properties of complex crystalline and amorphous phases of {Ge} and {Si}. {II}. Band structure and optical properties}, journal = {Physical Review~B}, year = 1973, volume = 8, pages = {27332755}, month = {September}, url = {http://link.aps.org/abstract/PRB/v8/p2733}, abstract = {We present calculations of the band structures and the imaginary part of the dielectric function ε2 as a function of energy for Ge and Si in the diamond, wurtzite, SiIII (BC8) and GeIII (ST12) structures using the empirical pseudopotential method. In particular we have obtained the symmetries of wave functions along important symmetry directions and identified the major contributions to the optical structure. A further study is made into the optical properties of amorphous Ge and Si using our shortrangedisorder model. We find that, unlike longrangedisorder models, shortrange disorder can explain \emph{both} the amorphous density of states and the amorphous $\varepsilon_2$. In particular we find that the $\varepsilon_2$ spectrum has the same form as an averaged matrix element as a function of frequency.} } @article{JoannopoulosCo73, author = {J. D. Joannopoulos and Marin L. Cohen}, title = {New insight into the optical properties of amorphous {Ge} and {Si}}, journal = {Solid State Communications}, year = 1973, volume = 13, number = 8, pages = {11151118}, month = {October}, doi = {doi:10.1016/00381098(73)905450}, abstract = {A short range disorder model, unlike present long range disorder theories, has been able to account well for both the density of states and the optical properties of amorphous Ge and Si. Our results indicate that the imaginary part of the dielectric function for amorphous Ge and Si has the same form as an averaged gradient matrix element as a function of energy. This conclusion should be valid for all tetrahedrally bonded amorphous solids.} } @article{JoannopoulosCo74, author = {J. D. Joannopoulos and Marvin L. Cohen}, title = {Electronic density of states of amorphous {III}{V} semiconductors}, journal = {Solid State Communications}, year = 1974, volume = 15, number = 2, pages = {105108}, month = {July}, doi = {doi:10.1016/00381098(74)907182 }, abstract = {A shortrange disorder model is used to predict the predominant features of the density of states of amorphous IIIV semiconductors. Confirmation of this model would help establish the existence of IIIIII and VV bonds in these materials.} } @article{JoannopoulosCo74disorder, author = {J. D. Joannopoulos and Marvin L. Cohen}, title = {Effects of disorder on the electronic density of states of {III}{V} compounds }, journal = {Physical Review~B}, year = 1974, volume = 10, pages = {15451559}, month = {August}, url = {http://link.aps.org/abstract/PRB/v10/p1545}, abstract = {We investigate the effects of two types of disorder on the electronic density of states of IIIV semiconductors using simple tightbinding models and the empirical pseudopotential method. For the first type of disorder we consider a stoichiometric system with fourfold coordination, all bonds satisfied, variations in the bond lengths and angles, and only unlikeatom bonds. The second type of disorder includes the properties of the first with the addition of likeatom bonds. These two types of disorder are studied explicitly by taking GaAs as a prototype and making various GaAs structures using the atomic positions of certain crystal structures with shortrange disorder. These structures are crystals; however, they have atoms in the primitive cells arranged in varying fashions. A comparison of the trends observed in the densities of states with the inclusion of different types of disorder reveals valuable information concerning the relationship of the structural nature of an amorphous system to its density of states. We present a model of the density of states of our amorphous prototype GaAs, for each type of disorder, which we believe would be consistent with some of the present experimental radialdistributionfunction data. The effects of these types of disorder are discussed in general, and hopefully they will be useful in identifying specific types of disorder in amorphous samples.} } @article{JoannopoulosCo74surface, author = {J. D. Joannopoulos and Marvin L. Cohen}, title = {Intrinsic surface states of (110) surfaces of group {IV} and {III}{V} semiconductors }, journal = {Physical Review~B}, year = 1974, volume = 10, pages = {50755081}, month = {December}, url = {http://link.aps.org/abstract/PRB/v10/p5075}, abstract = {Electronic local and total density of states calculations have been performed using tight binding models on the (110) surface of a group IV and IIIV semiconductor. Ge and GaAs are taken as prototypes and the surface is assumed to be unrelaxed. Several new surface states are obtained near the bottom of the valence bands. The origin, localization, and character of the surface states are examined.} } @article{YndurainJo74, author = {F. Yndurain and J. D. Joannopoulos and Marvin L. Cohen and L. M. Falicov}, title = {New theoretical method to study densities of states of tetrahedrally coordinated solids}, journal = {Solid State Communications}, year = 1974, volume = 15, number = 3, pages = {617620}, month = {August}, doi = {doi:10.1016/00381098(74)911569}, abstract = {A new simple method is proposed to calculate local densities of states of arbitrary tetrahedrally coordinated solids. It involves the selection of a finite cluster of atoms connected to an infinite Bethe lattice of coordination four. The method is accurate, is easily handled numerically, and converges fast. Loworder approximations yield sufficient information which is susceptible to consistent physical interpretation. This has been made for the diamond, BC8 and ST12 structures in terms of the ring topology around a given atom. Comparison with exact calculations is very good and the ring interpretation is physical and conceptually appealing.} } @article{JoannopoulosYn74, author = {J. D. Joannopoulos and F. Yndurain}, title = {``Cluster{Bethe}lattice'' method: Electronic density of states of amorphous and crystalline homopolar solids}, journal = {Physical Review~B}, year = 1974, volume = 10, pages = {51645174}, month = {December}, url = {http://link.aps.org/abstract/PRB/v10/p5164}, abstract = {A new method is developed to study the electronic density of states of infinite networks of atoms. The method involves treating part of the system exactly as a cluster and simulating the effects of the rest of the environment by connecting a Bethe lattice (Cayley tree) to the surface of the cluster. Calculations show that the local ringlike topologies of each atom are of primary importance in determining structure in the electronic density of states. The densities of states of the diamond, BC8, and ST12 structures are studied in detail using this method. These calculations are in excellent agreement with the exact results. Because of this, the method is used to obtain the density of states of the Polk and Connell randomnetwork models. These models give the same radial distribution functions but exhibit striking differences in their densities of states which are interpreted in terms of topology.} } @article{SchlueterJo74, author = {M. Schl{\"{u}}ter and J. D. Joannopoulos and Marvin L. Cohen and L. Ley and S. P. Kowalczyk and R. A. Pollak and D. A. Shirley}, title = {Structural nature of amorphous {Se} and {Te}}, journal = {Solid State Communications}, year = 1974, volume = 15, number = 6, pages = {10071010}, month = {September}, doi = {doi:10.1016/00381098(74)905195}, abstract = {We present new photoemission measurements on amorphous and trigonal Te. These results and other recent photoemission experiments on amorphous and trigonal Se are interpreted using density of states and charge density calculations. These studies result in new insights into the possible structural nature of the amorphous phase.} } @article{SchueterJo74xray, author = {M. Schl{\"{u}}ter and J. D. Joannopoulos and Marvin L. Cohen}, title = {New interpretation of the soft{X}ray absorption spectra of several alkali halides}, journal = {Physical Review Letters}, year = 1974, volume = 33, pages = {8991}, month = {July}, url = {http://link.aps.org/abstract/PRL/v33/p89}, abstract = {New pseudopotential calculations of the densities of states of trigonal Se and Te (in excellent agreement with recent photoemission measurements) are used to show that two easily resolved peaks in the photoemission data are directly related to interchain and intrachain bonding. This identification is accomplished by calculating electronic charge densities as a function of energy for different energy regions. Finally we introduce a new method for determining bonding charges by extracting shortwavelength components of the electronic charge density.} } @article{SchueterJo74trigonal, author = {M. Schl{\"{u}}ter and J. D. Joannopoulos and Marvin L. Cohen}, title = {New interpretation of photoemission measurements of trigonal {Se} and {Te}}, journal = {Physical Review Letters}, year = 1974, volume = 33, pages = 337, month = {July}, url = {http://link.aps.org/abstract/PRL/v33/p337/s2} } @article{JoannopoulosSc75, author = {J. D. Joannopoulos and M. chl{\"{u}}ter and Marvin L. Cohen}, title = {Electronic structure of trigonal and amorphous {Se} and {Te}}, journal = {Physical Review~B}, year = 1975, volume = 11, pages = {21862199}, month = {March}, url = {http://link.aps.org/abstract/PRB/v11/p2186}, abstract = {The electronic structure of trigonal and amorphous Se and Te is investigated using the empirical pseudopotential method (EPM), chargedensity calculations, and simple tightbinding models. Band structures and electronic densities of states are obtained which are in excellent agreement with recent photoemission measurements. The tightbinding models are used to obtain analytic expressions for the energy bands and to interpret the EPM band structures in terms of realspace orbitalorbital interactions. Chargedensity calculations obtained as a function of energy and evaluated within specific energy intervals are used to interpret various structure in the density of states. Specifically certain easily resolvable peaks in the experimental photoemission spectra are associated with intrachain and interchain localized states, respectively. By taking only shortwavelength components of the charge density, a bonding charge can be defined which gives an estimate of the intrachain vs interchain bonding strengths. The trigonal results along with model calculations to investigate the effects of bondangle variations on chains and the presence of eight and sixfold rings of bonds are used to interpret the changes observed in the experimental spectra of amorphous Se and Te. A new model of amorphous Se is proposed.} } @article{YndurainJo75, author = {F. Yndurain and J. D. Joannopoulos}, title = {``Cluster{Bethe}lattice'' method: The electronic density of states of heteropoloar systems}, journal = {Physical Review~B}, year = 1975, volume = 11, pages = {29572964}, month = {April}, url = {http://link.aps.org/abstract/PRB/v11/p2957}, abstract = {The ``clusterBethelattice'' method is extended to the study of heteropolar systems. The Bethe lattice is solved for binary compounds of arbitrary coordination using simple tightbinding models. In particular systems with tetrahedral coordination, such as the zincblende, BC8, and randomnetwork structures are examined in detail. The results are compared with recent experimental photoemission data on the amorphous phases of binary compounds and interpreted in terms of topology.} } @article{JoannopoulosKa75, author = {J. D. Joannopoulos and M. Kastner}, title = {Evidence for weak $s$$p$ hybridization in chalcogens}, journal = {Solid State Communications}, year = 1975, volume = 17, number = 2, pages = {221224}, month = {July}, doi = {doi:10.1016/00381098(75)900472}, abstract = {Recent arguments supporting strong $s$$p$ hybridization in crystalline S and Se based on ultraviolet photoemission spectroscopy and molecular orbital models are inconsistent with the results of Xray spectroscopy. However, all this data can be explained by weak $s$$p$ hybridization and the inclusion of a small admixture of $d$ states.} } @article{JoannopoulosYn75, author = {J. D. Joannopoulos and F. Yndurain}, title = {Moments and averages of the electronic density of states of amorphous and crystalline homopolar solids}, journal = {Physics Letters~A}, year = 1975, volume = 51, number = 2, pages = {7980}, month = {February}, doi = {doi:10.1016/03759601(75)902327}, abstract = {We have calculated the first fifteen moments of the density of states of the BC8 and ST12 structures and of the Polk random network model. Using these moments we calculate exact averages of the density of states of these structures.} } @article{YndurainJo76, author = {Felix Yndurain and J. D. Joannopoulos}, title = {Study of the electronic local density of states using the cluster{Bethe}lattice method: Application to amorphous {III}{V} semiconductors }, journal = {Physical Review~B}, year = 1976, volume = 14, pages = {35693577}, month = {October}, url = {http://link.aps.org/abstract/PRB/v14/p3569}, abstract = {The clusterBethelattice method is extended to study the WeaireThorpe $sp^3$ Hamiltonian. We derive a transformation which makes it possible to calculate the \emph{local} density of states of a system described with a fourorbital Hamiltonian by using the local Green's function of a oneorbital Hamiltonian. We study the effects on the density of states of the presence of likeatom bonds in amorphous IIIV semiconductors. We analyze two models: (a) likeatom bond in tetrahedral configuration, and (b) likeatom bond in triangular pyramidal configuration. Our results are in agreement with the experimental data.} } @article{JoannopoulosPo76, author = {J. D. Joannopoulos and W. B. Pollard}, title = {Electrons and phonons in amorphous pyramidally bonded solids}, journal = {Solid State Communications}, year = 1976, volume = 20, number = 10, pages = {947950}, month = {December}, doi = {doi:10.1016/00381098(76)904798}, abstract = {Electronic and phonon state densities of amorphous pyramidally bonded solids are calculated with As as a prototype. Using a Bethelattice, the GreavesDavis random network model and the ClusterBethelattice method, the results show: (1) the phonon states are isomorphic to, and \emph{obtained} from, the plike electron states, (2) the photoemission and infrared measurements are easily interpreted, and (3) the electronic and phonon state densities are \emph{extremely} sensitive to the local topology.} } @article{MeleJo77, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {{Gaussian} relaxation method. {I}. Homopolar tetrahedral solids}, journal = {Physical Review~B}, year = 1977, volume = 15, pages = {901908}, month = {January}, url = {http://link.aps.org/abstract/PRB/v15/p901}, abstract = {We present a method for extending an empirical tightbinding theory to the calculation of densities of states of very large but finite systems. We also derive a set of bondcentered analytic wave functions which can be used to obtain realistic energydependent charge densities from a tightbinding calculation. As an example, we apply the method to a model homopolar tetrahedral solid representative of Ge.} } @article{MeleJo77surface, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Electronic states near the band gap for the {GaAs}(110) surface}, journal = {Surface Science}, year = 1977, volume = 66, number = 1, pages = {3844}, month = {August}, doi = {doi:10.1016/00396028(77)903983}, abstract = {The recent theoretical controversy regarding empty surface states at the GaAs cleavage plane is examined by the exact solution of several tight binding models applied to a semiinfinite crystal. We correlate the existence of an empty surface gap state with the anion p/cation p character of this state. This work isolates the specific properties of a Hamiltonian for GaAs which yield band gap states on a relaxed surface and clarifies the essential differences among current theoretical treatments.} } @article{JoannopoulosSt77, author = {J. D. Joannopoulos and Th. Starkloff and Marc Kastner}, title = {Theory of pressure dependence of the density of states and reflectivity of selenium}, journal = {Physical Review Letters}, year = 1977, volume = 38, pages = {660663}, month = {March}, url = {http://link.aps.org/abstract/PRL/v38/p660}, abstract = {The pressure dependence of the density of states and the reflectivity of trigonal Se is calculated using a new selfconsistent pseudopotential technique. The results show that pressuredependent reflectivity experiments may be interpreted without the need for microscopic fields.} } @article{Joannopoulos77, author = {J. D. Joannopoulos}, title = {Theory of fluctuations and localized states in amorphous tetrahedrally bounded solids}, journal = {Physical Review~B}, year = 1977, volume = 16, pages = {27642774}, month = {September}, url = {http://link.aps.org/abstract/PRB/v16/p2764}, abstract = {We present a general solution of a Bethe lattice with arbitrary coordination for a Hamiltonian with an arbitrary number of degrees of freedom per site and an arbitrary number of interaction integrals. This solution is used in conjunction with a realistic tightbinding Hamiltonian to study the effects of rings and bondangle fluctuations on the plike region and gap region of the electronic density of states of amorphous tetrahedrally bonded solids. It is shown that even for completely connected networks with no dangling bonds, bondangle fluctuations create welldefined localized states which lie predominantly at the top of the valence band. These fluctuations also account for the steepening of the valenceband edge with disorder as observed experimentally in photoemission measurements. It is shown that rings do not play a direct role in this effect.} } @article{LaughlinJo77, author = {R. B. Laughlin and J. D. Joannopoulos}, title = {Phonons in amorphous silica}, journal = {Physical Review~B}, year = 1977, volume = 16, pages = {29422952}, month = {September}, url = {http://link.aps.org/abstract/PRB/v16/p2942}, abstract = {A new theory of lattice vibrations in amorphous silicon dioxide is presented in which the randomness of the solid is treated separately from its chemistry. The theory attributes all measurable properties of phonons in silica to the nearly crystalline nearestneighbor geometry of the lattice and to the disruptive effects of bondangle disorder. Neutron, infrared, and Raman spectra are calculated and compared with experiment. The theory is an application of the recently developed clusterBethelattice approach to studying amorphous solids.} } @article{StarkloffJo77, author = {Th. Starkloff and J. D. Joannopoulos}, title = {Local pseudopotential theory for transition metals }, journal = {Physical Review~B}, year = 1977, volume = 16, pages = {52125215}, month = {December}, url = {http://link.aps.org/abstract/PRB/v16/p5212}, abstract = {A simple localpseudopotential formalism capable of describing dband metals is introduced. The theory involves a straightforward method for generating a local pseudopotential for both valence and outermostcoreshell electrons. Niobium is used as a prototype to demonstrate the applicability of this method. Excellent agreement with nonlocal selfconsistent pseudopotential and augmentedplanewave band structures is obtained using a modest number of plane waves. The electronic density of states and core $4p$ absorption spectra are calculated. Two regions of distinctly different electronic character in the density of states could be demonstrated.} } @article{MeleJo78, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Intrinsic surface states and {Fermi}level pinning at metalsemiconductor interfaces}, journal = {Journal of Vacuum Science and Technology}, year = 1978, volume = 15, number = 4, pages = {13701373}, month = {July}, doi = {doi:10.1116/1.569789}, abstract = {We find that the wellknown transition in the Fermilevel pinning behavior of metalsemiconductor contacts can be quantitatively characterized in a fundamental way. Our model makes use of the clean semiconductor surface state energies, predicts the onset of the covalentionic transition, and explains strong Fermilevel pinning at metal contacts to covalent materials with large optical gaps and low bond polarizabilities.} } @article{JoannopoulosMe78, author = {J. D. Joannopoulos and Eugene J. Mele}, title = {Extrinsic surface states for oxygen chemisorbed on the {GaAs} (110) surface}, journal = {Journal of Vacuum Science and Technology}, year = 1978, volume = 15, number = 4, pages = {12871289}, month = {July}, doi = {doi:10.1116/1.569754}, abstract = {The electronic structure of various configurations of oxygen chemisorbed on GaAs (110) is studied using a localized orbital approach. Recent conflicting interpretations of experimental measurements, concerning the site of chemisorption, are resolved. It is concluded that for low averages, oxygen prefers to chemisorb to the surface arsenics and chemisorbs as an O$_2$ molecule.} } @article{MeleJo78chemisorption, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Site of oxygen chemisorption on the {GaAs}(110) surface}, journal = {Physical Review Letters}, year = 1978, volume = 40, pages = {341344}, month = {January}, url = {http://link.aps.org/abstract/PRL/v40/p341}, abstract = {We reconcile energyloss spectroscopy and chemicalshift studies of the oxidation of the GaAs(110) surface, which previously have led to contradictory conclusions about the oxygen bonding site. We have calculated densities of states and 100eV ultravioletphotoelectronspectroscopy (UPS) valenceband spectra to examine the site and molecular species of the chemisorbed oxygen. We conclude that for low coverages oxygen prefers to chemisorb to the surface arsenic atoms and chemisorbs as an O$_2$ molecule.} } @article{LaughlinJo78, author = {R. B. Laughlin and J. D. Joannopoulos and C. A. Murray and K. J. Hartnett and T. J. Greytak}, title = {Intrinsic surface phonons in porous glass}, journal = {Physical Review Letters}, year = 1978, volume = 40, pages = {461465}, month = {February}, url = {http://link.aps.org/abstract/PRL/v40/p461}, abstract = {Raman scattering and infrared reflectivity measurements have been performed on samples of porous Vycor glass. Three features in the spectra are due to intrinsic surface phonons. This identification is supported by a simple new theory for the surface.} } @article{MeleJo78metal, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Theory of metalsemiconductor interfaces}, journal = {Physical Review~B}, year = 1978, volume = 17, pages = {15281539}, month = {February}, url = {http://link.aps.org/abstract/PRB/v17/p1528}, abstract = {We present a microscopic model for the formation of Schottky barriers at metalsemiconductor contacts. The theory proposes that Schottky barriers are determined by ``metalinduced gap states'' at the semiconductor surface, which are danglingbond derived resonances. The danglingbond character of these states implies that the energies of surface states at clean semiconductor surfaces are very important in Schottkybarrier formation. This work introduces an ionicity parameter, obtained from atomicterm values, to quantitatively characterize the wellknown transition from covalent to ionic behavior at metalsemiconductor interfaces. This model directly associates this transition with a truly fundamental change in the electronic structure of the semiconductor substrate and provides a natural interpretation of strong Fermilevel pinning at metal contacts to covalent materials with large optical gaps and low bond polarizabilities.} } @article{PollardJo78I, author = {W. B. Pollard and J. D. Joannopoulos}, title = {Excitations in amorphous pyramidally bonded solids. {I}. Electrons}, journal = {Physical Review~B}, year = 1978, volume = 17, pages = {17701777}, month = {February}, url = {http://link.aps.org/abstract/PRB/v17/p1770}, abstract = {The electronic and phonon state densities of amorphous pyramidally bonded solids are calculated. These solids are shown to exhibit a unique isomorphism between the electronic and vibrational state densities. This atypical correspondence between the state densities suggests that a study of the electronic states can be used to understand the vibrational states. Arsenic, which epitomizes pyramidally bonded solids, is chosen as a prototype. Our results reproduce the trends founds in the xray photoemission data (XPS) for rhombohedral and amorphous As. The effects of bondangle variations and topological disorder upon the s and p like regions of the XPS spectra are investigated. A prescription is given by which the behavior of the plike electrons in pyramidally bonded solids can be understood in terms of the behavior of the slike electrons.} } @article{PollardJo78II, author = {W. B. Pollard and J. D. Joannopoulos}, title = {Excitations in amorphous pyramidally bonded solids. {II}. Phonons}, journal = {Physical Review~B}, year = 1978, volume = 17, pages = {17781784}, month = {February}, url = {http://link.aps.org/abstract/PRB/v17/p1778}, abstract = {Phonon state densities are calculated for crystalline and amorphous As. The results indicate a strong analogy between the electronic p states and the vibrational states. This isomorphism is used to derive relationships between the electronic and phonon coupling constants. The effects of topology upon the vibrational excitations in pyramidal solids are studied. A new model for infrared absorption is presented. The resulting spectra agree well with experimental measurements. Finally, inelastic neutronstructure factors are calculated and are used to interpret recent experimental results.} } @article{MeleJo78relax, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Electronic states at unrelaxed and relaxed {GaAs} (110) surfaces}, journal = {Physical Review~B}, year = 1978, volume = 17, pages = {18161827}, month = {February}, url = {http://link.aps.org/abstract/PRB/v17/p1816}, abstract = {We have shown that, using a general class of Hamiltonians, the transfermatrix technique may be used to obtain exact solutions for the electronic states at any crystal surface bounded by semiinfinite bulk. This result is formally generalized as a theorem and is used to study the electronic states at a clean GaAs (110) surface. The calculation employs an empirical tightbinding Hamiltonian which realistically models the GaAs surface and allows meaningful comparison with both experiments and selfconsistent pseudopotential calculations. Surface states are calculated for the clean (110) surface, and a variety of structural relaxations are studied.} } @article{LaughlinJo78vibrations, author = {R. B. Laughlin and J. D. Joannopoulos}, title = {Effect of secondnearestneighbor forces on the vibrations of amorphous {SiO}$_2$}, journal = {Physical Review~B}, year = 1978, volume = 17, pages = {27902792}, month = {March}, url = {http://link.aps.org/abstract/PRB/v17/p2790}, abstract = {A method is presented of solving for the vibrational Green's function of a silicondioxide Bethe lattice when it contains secondnearestneighbor interactions. The method is used to find the density of states of a Bethe lattice constructed with a Keating Hamiltonian. The primary effect of the secondnearestneighbor interactions is a softening of the rocking bands below 450 cm$^{1}$.} } @article{LaughlinJo78phonons, author = {R. B. Laughlin and J. D. Joannopoulos}, title = {Theory of surface phonons in amorphous silica}, journal = {Physical Review~B}, year = 1978, volume = 17, pages = {49224930}, month = {June}, url = {http://link.aps.org/abstract/PRB/v17/p4922}, abstract = {A new theoretical approach to studying surface vibrations in amorphous solids is applied to silicon dioxide. The method entails modeling the surface as a Bethe lattice with a dangling bond, and treating the twodimensional nature of the surface and the surface topography as small perturbations. The theory successfully describes the nature and origin of the surface states, as well as their relative intensities as observed in recent infraredreflectivity and Ramanscattering experiments. Both intrinsic surface effects and those caused by the presence of adsorbates are discussed.} } @article{MeleJo78chemisorption2, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Theory of oxygen chemisorption on {GaAs}(110)}, journal = {Physical Review~B}, year = 1978, volume = 18, pages = {69997010}, month = {December}, url = {http://link.aps.org/abstract/PRB/v18/p6999}, abstract = {Using a localizedorbital representation, we have calculated the electronic structure of oxygen chemisorbed on the GaAs (110) surface in various bonding configurations. Two principal results which emerge are (i) that the effects of chemisorption on the surface species are strongly affected by allowing the surface to relax, and (ii) that molecularlike correlations persist and are important for localized nonbonding orbitals in the adsorbate. Our results provide detailed interpretations of xray photoemission spectroscopy, ultraviolet photoemission spectroscopy, and electronloss spectroscopy studies of the oxidation of GaAs(110), and a comparison with all the experimental data strongly indicates that in the initial stages of oxidation O$_2$ chemisorbs principally on the surface As, with the surface maintaining its relaxed configuration.} } @article{NemanichLu78, author = {R. J. Nemanich and G. Lucovsky and W. Pollard and J. D. Joannopoulos}, title = {Spectroscopic evidence for bonding coordination defects in amorphous {As}}, journal = {Solid State Communications}, year = 1978, volume = 26, number = 3, pages = {137139}, month = {April}, doi = {doi:10.1016/00381098(78)910554}, abstract = {Spectroscopic evidence for fourfold and possibly, twofold coordinated defects in amorphous As is presented. Structure in Raman scattering and infrared absorption measurements is interpreted using forceconstant models and the clusterBethelattice method. The positions and intensities of certain peaks suggest that these defects may be the analogues of valence alternation pairs in amorphous chalcogenide glasses.} } @article{StarkloffJo79, author = {Th. Starkloff and J. D. Joannopoulos}, title = {Theory of the pressure dependence of the electronic and optical properties of trigonal {Te}}, journal = {Physical Review~B}, year = 1979, volume = 19, pages = {10771088}, month = {January}, url = {http://link.aps.org/abstract/PRB/v19/p1077}, abstract = {Using a selfconsistent pseudopotential formalism based on realistic pseudopotentials of the free ion and atom, the influence of hydrostatic pressure on the electronic structure and the optical spectra of trigonal Te is investigated. The pressure derivative of the fundamental gap is found to be in good agreement with experimental results. The joint density of states, $\varepsilon_1$ and $\varepsilon_2$, energyloss function, reflectivity, and absorption coefficient are calculated from 0 to 13 eV. Very good agreement is obtained with experimental spectra. The very strong pressureinduced enhancement of the reflectivity can be reproduced. The anisotropy anomaly in the reflectivity at high energies is explained and traced to distinct transitions in the Brillouin zone. Predictions are made about the pressure dependence of the optical spectra above 4 eV. Finally, selfconsistent charge densities and the electronic density of states are also presented.} } @article{MeleJo79dangling, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Doubledanglingbond defects and band bending at the {GaAs} (110) surface}, journal = {Physical Review~B}, year = 1979, volume = 19, pages = {29282932}, month = {March}, url = {http://link.aps.org/abstract/PRB/v19/p2928}, abstract = {We note that contactpotential measurements which have shown band bending at the nonpolar GaAs surface can be interpreted with a simple model for cleavage defects at the surface. We find that defects with two dangling bonds at the surface produce both filled and empty defect states in the gap. This result then supports experimental work indicating that band bending at GaAs(110) is defect related.} } @article{PollardJo79, author = {W. B. Pollard and J. D. Joannopoulos}, title = {Vibrational excitations at defect sites in amorphous tetrahedral and pnictide semiconductors}, journal = {Journal of NonCrystalline Solids}, year = 1979, volume = {3536}, pages = {11791184}, month = {January}, doi = {doi:10.1016/00223093(80)903579}, abstract = {The vibrational properties of defects in amorphous As, amorphous Si, and fluorinated amorphous Si are investigated. Infrared absorption, Raman scattering, and densities of states spectra are calculated for a variety of bonding coordination defects within the clusterBethelattice formalism. A new technique for solving arbitrary Bethe lattices is also briefly discussed.} } @article{PollardJo79arsenic, author = {W. B. Pollard and J. D. Joannopoulos}, title = {Electronic structure of defects in amorphous arsenic}, journal = {Physical Review~B}, year = 1979, volume = 19, pages = {42174223}, month = {April}, url = {http://link.aps.org/abstract/PRB/v19/p4217}, abstract = {The clusterBethelattice method, selfconsistent pseudopotentials, and tightbinding Hamiltonians are used to study the electronic structure of various possible defects in amorphous As. The electronic densities of states of isolated twofold coordinated As atoms, fourfold coordinated As atoms, ``intimate valence alternation pairs,'' and vacancies are calculated. The origin, character, and localization of the resulting gap states are discussed. It is found that isolated fourfold coordinated As atoms give rise to states in the gap that are nondegenerate, fairly localized, and predominantly plike in character. Moreover, our investigation indicates that the character and localization of the gap states arising from the ``valence alternation pair'' are consistent with the properties of optically induced paramagnetic states observed in amorphous As, while the electronic gap states of the vacancy defect may be consistent with the recently observed thermally generated paramagnetic states in this material.} } @article{Joannopoulos79, author = {J. D. Joannopoulos}, title = {Theoretical calculations of defect states in amorphous semiconductors}, journal = {Journal of NonCrystalline Solids}, year = 1979, volume = {3536}, pages = {781792}, month = {January}, doi = {doi:10.1016/00223093(80)902987}, abstract = {A review of realistic theoretical calculations of the electronic structure of defects in amorphous As, amorphous Se, amorphous As$_2$Se$_3$ and hydrogenated amorphous Si is presented.} } @article{Joannopoulos79cayley, author = {J. D. Joannopoulos}, title = {Use of {Cayley} trees to study excitations in disordered solids}, journal = {Journal of NonCrystalline Solids}, year = 1979, volume = 32, pages = {241255}, month = {February}, doi = {doi:10.1016/00223093(79)900759}, abstract = {It is shown how, and why, Cayley trees (Bethelattices) are useful models for studying elementary excitations in disordered solids. Special attention is focused on the effects of disorder on the electronic states of various amorphous bonded solids. Specific examples are drawn from problems involving: (1) bondangle distortions and localized states in amorphous Ge, (2) hydrogen configurations in hydrogenated amorphous Si, (3) bondingcoordination defects in amorphous As and (4) amorphous SiSiO$_2$ interfaces.} } @article{LaughlinJo79, author = {R. B. Laughlin and J. D. Joannopoulos and D. J. Chadi}, title = {Use of the cluster{Bethe}lattice method in surface studies}, journal = {Journal of Vacuum Science and Technology}, year = 1979, volume = 16, number = 5, pages = {13271330}, month = {September}, doi = {doi:10.1116/1.570151}, abstract = {A general review of the clusterBethelattice method is presented which emphasizes its use as a tool in surface studies. The basic elements of the method are summarized, with simple models used to illustrate the most important ideas. The strengths and weaknesses of the method are contrasted, and the types of questions the method can best deal with ar identified. The method is then illustrated with an example. Recent results are presented which show that both bondangle and structural disorder can produce localized states at the SiSiO$_2$ interface which have no analog in either bulk silicon or bulk SiO$_2$. Inhomogeneities at the interface, particularly small ''peninsulas'' of silicon protruding into the oxide, are found to produce strong interface resonances approximately 2 eV below the valence band maximum. Distortions of the SiOSi angle near the interface are found to cause tailing into the gap in the oxide near 10 eV below the silicon valence band edge. Finally, distortions of the tetrahedral angles on a silicon atom at the interface are investigated and found to have effects similar to those observed in bulk elemental silicon.} } @article{LaughlinJo79silica, author = {R. B. Laughlin and J. D. Joannopoulos and D. J. Chadi}, title = {Bulk electronic structure of {SiO}$_2$}, journal = {Physical Review~B}, year = 1979, volume = 20, pages = {52285237}, month = {December}, url = {http://link.aps.org/abstract/PRB/v20/p5228}, abstract = {The electronic structures of crystalline and amorphous SiO$_2$ are examined via the tightbinding method. A new tightbinding Hamiltonian, fit to experiment and to the pseudopotential band structure of α quartz, is used to calculate densities of states for both α quartz and an SiO$_2$ Bethe lattice. These are shown to compare favorably with xray photoemission spectra of $\alpha$ quartz and amorphous SiO$_2$. The computational results are analyzed qualitatively using the bondorbital approach. For both crystalline and amorphous SiO$_2$ it is suggested that oxygen 2s character in the lower conduction bands may be necessary to account for the large gap. Local symmetries of the lonepairlike bands of possible relevance to the optical properties are discussed.} } @article{MeleJo79, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Electronic structure of {Al} chemisorbed on {GaAs}(110)}, journal = {Journal of Vacuum Science and Technology}, year = 1979, volume = 16, number = 5, pages = {11541158}, month = {September}, doi = {doi:10.1116/1.570180}, abstract = {We report results of tight binding calculation for ordered half monolayer coverages of Al on GaAs(110) in a variety of chemisorption configurations. We find that differences between the chemisorption geometries and relaxations of these configurations lead to a variety of electrical characteristics of the surfaces. The relation of this result to the barrier formed at higher coverages is discussed.} } @article{MeleJo79A1, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Surfacebarrier formation for {A1} chemisorbed on {GaAs}(110)}, journal = {Physical Review Letters}, year = 1979, volume = 42, pages = {10941097}, month = {April}, url = {http://link.aps.org/abstract/PRL/v42/p1094}, abstract = {Using a localized orbital theory, we have studied the electronic and electrical properties of A1 chemisorbed on GaAs(110). We find that both the microscopic data and the macroscopic induced barrier are fully explained by the electronic structure of a new AlAsGa complex on the surface. This complex results from an exchange reaction in which Al replaces the surface Ga and an unexpected structural relaxation induced by the chemisorption of the metal.} } @article{DoerreJo79, author = {G. Doerre and J. D. Joannopoulos}, title = {Electronic states of {Te} above the highpressure phase transition}, journal = {Physical Review Letters}, year = 1979, volume = {10401042}, number = 43, pages = {10401042}, month = {October}, url = {http://link.aps.org/abstract/PRL/v43/p1040}, abstract = {A selfconsistent pseudopotential bandstructure calculation for the recently determined highpressure (above 38 kbar) monoclinic tellurium structure is presented. The density of states and charge densities for representative regions within the unit cell are calculated. Our results show that highpressure monoclinic tellurium is metallic, and that a crystal bondlength asymmetry is caused by an electronically driven distortion. The results also predict a highly anisotropic conductivity.} } @article{VanderbiltJo79, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Calculation of defect states in amorphous selenium}, journal = {Physical Review Letters}, year = 1979, volume = 42, pages = {10121015}, month = {April}, url = {http://link.aps.org/abstract/PRL/v42/p1012}, abstract = {Selfconsistent pseudopotential and tightbinding techniques are used to study glassy Se. Both onefold and threefold corrdinated defects give rise to nondegenerate, nonhydrogenic gap states. The former state arises from an interaction, unique to the chalcogenides, between a dangling bond and a neighboring lone pair. The latter state is surprisingly delocalized, involving interactions between many orbitals. It is shown that defect electronic configurations, and repulsive interatomic terms, are crucial to the calculation of total energies.} } @article{WangJo80, author = {Y. Wang and J. D. Joannopoulos}, title = {The {Ga} core exciton at unrelaxed {GaAs} (110)}, journal = {Journal of Vacuum Science and Technology}, year = 1980, volume = 17, number = 5, pages = {9971000}, month = {September}, doi = {doi:10.1116/1.570656}, abstract = {ore excitonic effects at an unrelaxed GaAs (110) surface are studied using a localized orbital approach. A realistic Hamiltonian for the conduction bands is obtained with inclusion of $5s$ orbitals and third nearest neighbor interactions. The binding energies of core excitons both in the bulk and at the surface are calculated. The results are in good agreement with experimental measurements, with a contact potential of $\angle$1.0eV. In addition it is shown that the binding energy has a maximum value and the limiting case can be thought of as an ideal vacancy.} } @article{VanderbiltJo80, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Structural excitation energies in {Selenium}}, journal = {Solid State Communications}, year = 1980, volume = 35, number = 7, pages = {535538}, month = {August}, doi = {doi:10.1016/00381098(80)908911}, abstract = {A realistic selfconsistent pseudopotential approach, capable of generating accurate structural ground state and excitation energies, is used to settle the controversy regarding neutral bond coordination defects in glassy Se. Many surprising results emerge. The threefold defect is not even metastable, but decays immediately into a onefold defect. The latter is the lowerst energy simple defect, with a remarkably low formation energy of reverse similar, equals $\simeq$0.5eV. Relaxations are found to be unexpectedly small.} } @article{LeeJo80, author = {D. H. Lee and J. D. Joannopoulos}, title = {Surface states at unrelaxed {ZnO}($10\bar{1}0$)}, journal = {Journal of Vacuum Science and Technology}, year = 1980, volume = 17, pages = {987988}, month = {September}, doi = {doi:10.1116/1.570653}, abstract = {The electronic structure of an ideal ZnO($10\bar{1}0$) surface has been studied using a tightbinding approach, with the surface modeled in a semiinfinite geometry. It is found that the surface states follow the band edge very closely due to the ionic character of the surface. Moreover, the close spacing between Zn $3d$ and O $2p$ levels produces a small but not negligible amount of $d$ character in the surface states.} } @article{LeeJo81, author = {D. H. Lee and J. D. Joannopoulos}, title = {A new theory of electronic surface states}, journal = {Journal of Vacuum Science and Technology}, year = 1981, volume = 19, number = 3, pages = {355359}, month = {September}, doi = {doi:10.1116/1.571063}, abstract = {Presently, stateoftheart theoretical studies of elementary excitations at surfaces involve large, time consuming, and often cumbersome calculations, which are not very practical for surfaces with low symmetry. A new method has just been developed that directly yields more information and is computationally faster (by an order of magnitude) than any current method available. This method will allow calculations of complicated surface systems not previously possible. The ingredients of the theory are any semiinfinite system of atomic layers, any set of basis functions spanning each layer, and one single matrix which depends only on the interaction matrix elements among the atomic layers. What is remarkable is that the eigenanalysis of this matrix gives a complete quantum mechanical description of the elementary excitations of the system. Therefore with just a diagonalization of this matrix, the knowledge of its eigenvalues and eigenvectors give directly (1) the projected bandstructure (without having to diagonalize the bulk Hamiltonian for a large number of $K_\perp$'s); (2) the bona fide surface bands with high accuracy (without the need of a bulk Green function first); (3) the surface Green function as a function of energy $E+i\delta$, in the correct limit that $\delta\to 0$! (without the need to solve slowly convergent iterative equations for \emph{complex energies}); and (4) the precise decay length, orbital character, and symmetry of the wave function. The method, although versatile, is also conceptually simple. All the properties mentioned above are obtained from a simple set of rules. A tutorial discussion of the formalism will be presented along with several instructive examples illustrating its use.} } @article{IhmJo82, author = {J. Ihm and J. D. Joannopoulos}, title = {Firstprinciples determination of the structure of the {Al}/{GaAs}(110) surface}, journal = {Journal of Vacuum Science and Technology}, year = 1982, volume = 21, number = 2, pages = {340343}, month = {July}, doi = {doi:10.1116/1.571777 }, abstract = {The structure of Al atoms deposited on the GaAs (110) surface has been studied using a first principles pseudopotential energy minimization calculation. The lowest energy configurations are determined by minimizing the total energy of the system with respect to its structural degrees of freedom. The most stable configuration is such that Al atoms replace the second (or deeper) layer Ga atoms. At temperatures where this reaction cannot be activated, two important processes are found to exist. In the low coverage limit, atoms favor twofold sites connecting an As atom and a Ga atom (in the next zigzag chain) on the surface, forming strong chemical bonds with the substrate. At higher coverages, on the other hand, Al atoms tend to cluster and make new bonds among themselves. The chemisorption energy map over the surface has been obtained and the possible path of the migration of Al atoms for clustering is investigated.} } @article{LeeJoannopoulos82, author = {D. H. Lee and J. D. Joannopoulos}, title = {Ideal and relaxed surfaces of {SiC}}, journal = {Journal of Vacuum Science and Technology}, year = 1982, volume = 21, number = 2, pages = {351357}, month = {July}, doi = {doi:10.1116/1.571779}, abstract = {During the past few years there has been considerable interest in the relaxation of GaAs(110) and Si(100)($2\times 1$) surfaces. The surfaces of SiC, however, provide an intermediate system between these heteropolar and homopolar structures. It is interesting therefore to investigate the kinds of relaxation that might occur at the zincblende and wurtzite surfaces of SiC. We perform calculations using Chadi's energy minimization scheme, with the coefficients of linear, quadratic, and cubic energy correction terms fitted to the bulk lattice constant, bulk modulus, and thermal expansion coefficient. To check the validity of this model, we calculate six experimentally known phonon frequencies. The agreement between theory and measured values is quite good. With this model, we investigate the relaxation at the SiC(110) surface in the zincblende structure and SiC($10\bar{1}0$) and ($11\bar{2}0$) surfaces in the wurtzite structure. The results show a combination of downward movement and buckling for all three surfaces. The reduction in total energy is about 0.21 eV/atom for the (110) surface, and 0.24 eV/atom, 0.18 eV/atom for ($10\bar{1}0$) and ($11\bar{2}0$), respectively. Further results include the determination of ideal and relaxed electronic structure and optimum relaxed geometries. Finally, a new theory which extends Chadi's scheme to do first principles phonon calculations is presented and discussed.} } @article{ChadiIh83, author = {D. J. Chadi and J. Ihm and C. Tanner and J. D. Joannopoulos}, title = {Theoretical study of the {As}(100) surface reconstruction of {GaAs}}, journal = {Physica~B+C}, year = 1983, volume = {117118}, pages = {798800}, month = {March}, doi = {doi:10.1016/03784363(83)906575}, abstract = {The reconstruction of the Asterminated GaAs(100) surface was examined using both the tightbinding and the selfconsistent pseudopotential methods. Dimer models with unit cells as large as the c$4\times 4$ one seen in electrondiffraction experiments were tested. Unlike the Si(100) surface, symmetric as well as asymmetric dimers are predicted to occur in equal densities at the surface.} } @article{VanderbiltJo83, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Total energies of structural defects in glassy {Se}}, journal = {Journal of NonCrystalline Solids}, year = 1983, volume = {5960}, pages = {937944}, month = {December}, doi = {doi:10.1016/00223093(83)90323X}, abstract = {The effective Hubbard $U$ is calculated for the interconversion of 1fold and 3fold bonding coordination defects in glassy Se. This is accomplished by applying local density total energy calculations directly to charged defects in a superlattice configuration. It is found that the defect remains 1fold coordinated in the $D^$ and $D^0$ charge states, but spontaneously forms a 3fold center in the $D^+$ charge state. The structural relaxation energy involved in the bond switching gives rise to a sizable negative contribution to $U$, but a still larger Coulomb repulsion gives rise to an overall positive $U$. This result is not, however, inconsistent with a negative $U$ in the compound chalcogenide glasses.} } @article{IhmLe83, author = {J. Ihm and D. H. Lee and J. D. Joannopoulos and A. N. Berker}, title = {Study of high order reconstructions of the {Si}(100) surface}, journal = {Journal of Vacuum Science and Technology~B}, year = 1983, volume = 1, number = 3, pages = {705708}, month = {July}, doi = {doi:10.1116/1.582584}, abstract = {The structure of the reconstructed Si(100) surface is studied using the tightbinding energy minimization method combined with a renormalizationgroup calculation. ($2\times 1$), $c(2\times 2)$, $p(2\times 2)$, $(4\times 1)$, $c(4\times 2)$, and $p(4\times 2)$ structures are studied and compared. They are classified into two groups which are separated by a large activation barrier. In one group, $(2\times 1)$, $p(2\times 2)$, and $c(4\times 2)$ structures are almost degenerate in energy, which can result in diffuse $(2\times 1)$ LEED patterns. In the other group, a secondorder phase transition can occur from the $c(2\times 2)$ to the disordered structure at $\sim 800$ K. Domains of different structures may coexist on the actual surface.} } @article{TarnowAn85, author = {Eugene Tarnow and Alex Antonelli and J. D. Joannopoulos}, title = {Ground state properties of crystalline {As}$_2${Se}$_3$}, journal = {Journal of NonCrystalline Solids}, year = 1985, volume = {7778}, pages = {9598}, month = {December}, doi = {doi:10.1016/00223093(85)906192}, abstract = {An abinitio pseudopotential calculation in the local density approximation is performed on crystalline As$_2$Se$_3$. Total energies, charge densities, and density of states are presented and compared with experiment.} } @article{BarYamJo85, author = {Y. BarYam and J. D. Joannopoulos}, title = {Correlation energy of deep level traps in {aSi:H}}, journal = {Journal of NonCrystalline Solids}, year = 1985, volume = {7778}, pages = {99102}, month = {December}, doi = {doi:10.1016/00223093(85)906209}, abstract = {We report the first realistic calculation of the effective electronic correlation energy of the dangling bond in amorphous hydrogenated silicon. Total energies obtained using abinitio pseudopotential density functional theory yield a full configuration coordinate diagram. When host distortions of the amorphous material are included by use of a simple model, both the effective electronic correlation energy and the position of the electronic states in the gap are affected. We obtain a distribution of effective electronic correlation energies, and a distortion of electronic transition energies in the gap.} } @article{BarYamJo87, author = {Y. BarYam and J. D. Joannopoulos}, title = {Theories of defects in amorphous semiconductors}, journal = {Journal of NonCrystalline Solids}, year = 1987, volume = {9798}, pages = {467474}, month = {December}, doi = {doi:10.1016/00223093(87)901104}, abstract = {The existence of defects in amorphous semiconductors is itself an problem of fundamental interest. Conventional approaches accept their existence only phenomenologically. We described a new approach to defining and understanding defects. The existence of defects is a natural outgrowth of our thermodynamical ensemble theory of electronic states in disordered systems. Specific properties are also predicted through properties of thermodynamic ensembles and the relationship between defect electronic properties and structural energies. Thus we obtain directly electronic properties without a need for detailed microscopic information.} } @article{AlerhandKa89, author = {O. L. Alerhand and Efthimios Kaxiras and J. D. Joannopoulos and G. W. Turner}, title = {Model of epitaxial growth of {GaAs} on {Si}(100): Nucleation at surface steps}, journal = {Journal of Vacuum Science and Technology~B}, year = 1989, volume = 7, number = 4, pages = {695699}, month = {July}, doi = {doi:10.1116/1.584627}, abstract = {A microscopic model of epitaxial growth of GaAs on Si(100) is proposed, which explains how doublelayer steps on the Si(100) surface act as nucleation sites in the initial stages of growth. On the flat terraces of Si(100), there is a strong tendency towards the formation of mixed layers of Ga and As, which inhibit growth of zincblende GaAs. The bonding topology of the surface steps supresses this mixing, and drives the growth of ideally bonded GaAs in a threedimensional mode with islands forming along the steps edges. Totalenergy calculations indicate that the proposed model for nucleation of GaAs at doublelayer steps on the surface is energetically stable.} } @article{VanderbiltAl89, author = {David Vanderbilt and O. L. Alerhand and Robert D. Meade and J. D. Joannopoulos}, title = {Elastic stress domains of the {Si}(100) surface}, journal = {Journal of Vacuum Science and Technology~B}, year = 1989, volume = 7, number = 4, pages = {10131016}, month = {July}, doi = {doi:10.1116/1.584792}, abstract = {It is proposed that surfaces of crystals that reconstruct with broken orientational symmetry and exhibit an anisotropic stress tensor are unstable to the formation of elastic stress domains. The ground state of such a surface is not uniform, but corresponds to a mosaic pattern of the different domains. This result is applied together with microscopic calculations of surface stress to study the Si(100) surface. A structural transition from singlelayer to doublelayer steps for vicinal Si(100) surfaces is also studied.} } @article{PayneRo89, author = {M. C. Payne and N. Roberts and R. J. Needs and M. Needels and J. D. Joannopoulos}, title = {Total energy and stress of metal and semiconductor surfaces}, journal = {Surface Science}, year = 1989, volume = {211212}, pages = {120}, month = {April}, doi = {doi:10.1016/00396028(89)907486}, abstract = {Energies and stresses of reconstructed surfaces of metals and semiconductors are calculated using pseudopotential methods. Surface stresses of aluminium, iridium, platinum, gold and silicon surfaces are calculated and the implications for reconstruction in these systems are discussed. The energies of dimer reconstructions on the (001) surfaces of silicon and germanium have been calculated and the stability of buckled dimers is investigated. The mechanism responsible for stabilising the Takayanagi reconstruction on the (111) surfaces of silicon and strained germanium is deduced from total energy calculations for substructures of the reconstructed surface.} } @article{KaxirasJo89, author = {Efthimios Kaxiras and J. D. Joannopoulos}, title = {On the possibility of twodimensional growth of {GaAs} on atomically flat {Si}(100) surfaces}, journal = {Surface Science}, year = 1989, volume = 224, pages = {515524}, month = {December}, doi = {doi:10.1016/00396028(89)909291}, abstract = {The possibility of growing GaAs on atomically flat Si(100) surfaces is investigated at monolayer and bilayer coverages. Based on abinitio totalenergy calculations at zerotemperature we compare the stability of separate domains of Ga or As versus mixed monolayers, the stability of pure GaAs bilayers versus mixed bilayers and the formation of GaAs bilayer islands versus wetting of the surface. The results predict that twodimensional growth of bulk GaAs on atomically flat regions (terraces) of Si(100) is severely inhibited at the initial stages, due to chemical and rehybridization reactions of the Ga and As atoms on the Si surface.} } @article{KaxirasAl90, author = {Efthimios Kaxiras and O. L. Alerhand and J. D. Joannopoulos and G. W. Turner}, title = {Thermodynamic and kinetic aspects of {GaAs} growth on {Si}(100)}, journal = {Superlattices and Microstructures}, year = 1990, volume = 8, number = 2, pages = {229232}, doi = {doi:10.1016/07496036(90)90098R}, abstract = {A microscopic model of GaAs growth on Si(100) is discussed, which takes into account the role of kinetics and the morphology of the substrate. The model predicts that growth at the initial stages is dominated by GaAs islands nucleated at steps of the substrate, which begin forming upon deposition concurrently with layered growth on the terraces. This mode of growth does not fall within any of the traditional thermodynamic models of heteroepitaxial growth, and allows for different growth channels depending on deposition conditions. The predicted shape and structure of the GaAs seeds nucleated at Si steps are in agreement with recent experimental observations.} } @article{AlerhandWa91, author = {O. L. Alerhand and J. Wang and J. D. Joannopoulos}, title = {Growth of {As} overlayers on vicinal {Si}(100) surfaces}, journal = {Journal of Vacuum Science and Technology~B}, year = 1991, volume = 9, number = 4, pages = {24232426}, month = {July}, doi = {doi:10.1116/1.585714 }, abstract = {Firstprinciples totalenergy calculations are used to study the structure of As overlayers deposited on a stepped Si(100) surface. It is predicted that deposition of As at low substrate temperatures allows the As to grow directly on top of the Si surface, while growth at high temperatures results in a rearrangement of the surface as if the As replaced the original top Si layer. This result explains the sublattice orientation dilemma in GaAsonSi epitaxy.} } @article{RappeJo92, author = {Andrew M. Rappe and J. D. Joannopoulos and P. A. Bash}, title = {Test of utility of plane waves for the study of molecules from first principles}, journal = {Journal of the American Chemical Society}, year = 1992, volume = 114, number = 16, pages = 6466, month = {July}, abstract = {This paper studies the applicability of a planewave basis set for density functinal calculations of the properties of molecules from first principles. The main features of the planewave method are described, including pseudopotentials and supercells. The results for a number of small molecules are reported. The close agreement with experiment and with a standard method of quantum chemistry calculation indicates the promise which this method holds for chemical and biochemical firstprinciples computations.} } @article{IretaGa98, author = {Joel Ireta and Marcelo Galv{\'{a}}n and Kyeongjae Cho and John D. Joannopoulos}, title = {Local reactivity of charybdotoxin, a $K^+$ channel blocker}, journal = {Journal of the American Chemical Society}, year = 1998, volume = 120, pages = {97719778}, month = {September}, abstract = {Charybdotoxin (ChTX) is a 37residue polypeptide that has been extensively used in sitedirected mutagenesis experiments as a template to deduce models for the external pore appearance of $K^+$ channels. The microscopic details of the ChTXchannel interaction, however, remain as a challenge for experimental and theoretical approaches. In this work, regional chargetransfer abilities, measured by chemical softness, $s(\mathbf{r})$, are used as companion properties of the electrostatic potential, $V(\mathbf{r})$, in the search for a qualitative structurefunction relationship in the ChTX$K^+$ channel interaction. Both quantities were obtained with an ab initio methodology in massively parallel computers. In the analysis of $s(\mathbf{r})$ and $V(\mathbf{r})$, regions of the size of amino acids were considered because this is the appropriate scale to correlate with sitedirected mutagenesis experiments. The correspondence between experimentally identified crucial amino acids sites and regional softnesses indicates that charge transfer to ChTX could be one of the main stabilization effects in the ChTXchannel complex. Also, it provides an explanation for the strong dependence of the dissociation constant of the complex on mutations of crucial amino acids. In addition, it is shown to be feasible to find structurefunction relationships by combining local reactivity parameters and experimental data involving site directed mutagenesis.} } @article{GalanDalPino93, author = {Marcelo Galvan and Dal Pino, Jr., Arnaldo and Jing Wang and J. D. Joannopoulos}, title = {Local softness, scanning tunneling microscopy, and surface reactivity}, journal = {Journal of Physical Chemistry}, year = 1993, volume = 97, pages = {783785}, month = {January}, abstract = {Within the finite temperature extension of density functional theory, it is shown that local softness measures the ability of a chemical system to gain or donate charge as the external potential changes. This feature allows us to establish a formal connection between local softness and scanning tunneling microscopy (STM) images. We show that, under appropriate conditions, STM images can be used to measure local softness for surfaces. Finally, a potential application of those images as reactivity criteria within the context of the hard and soft acids and bases (HSAB) principle is discussed.} } @article{KaxirasAl92, author = {Efthimios Kaxiras and O. L. Alerhand and J. Wang and J. D. Joannopoulos}, title = {Theoretical modeling of heteroepitaxial growth initiation}, journal = {Materials Science & Engineering~B}, year = 1992, volume = 14, number = 3, pages = {245253}, month = {November}, doi = {doi:10.1016/09215107(92)90306T}, abstract = {We review recent theoretical investigations on the initiation of heteroepitaxial growth. The thermodynamic premises on which the various growth models are based are considered. A qualitatively new model which applies to polaronnonpolar heteroepitaxial growth is discussed. Specific consequences of the new growth model are explored and shown to be consistent with experimental observations of GaAs growth on Si(100) vicinal surfaces. Firstprinciples calculations of the total energy of competing structures provide support for the proposed model. The totalenergy comparisons afford a resolution to a longstanding puzzle concerning the effect of temperature on the orientation of GaAs overlayers on vicinal silicon substrates.} } @article{AriasJo92, author = {T. A. Arias and J. D. Joannopoulos}, title = {Ab initio prediction of dopant segregation at elemental semiconductor grain boundaries without coordination defects}, journal = {Physical Review Letters}, year = 1992, volume = 69, number = 23, pages = {33303333}, month = {December}, url = {http://link.aps.org/abstract/PRL/v69/p3330}, abstract = {We present results of an extensive stateoftheart ab initio study of the segregation of dopants to the $\Sigma=5$ tilt (310) grain boundary in germanium. Despite the lack of reactive bond coordination defects in this boundary, we predict a significant tendency for $n$type dopants to segregate to the boundary. Our results lead to a general theory of dopant segregation at tetrahedrally bonded semiconductor defects, where lattice relaxation effects are irrelevant and segregation is driven by the mixing of localized carrier states of the defect with the hydrogenic states of the dopant.} } @article{BrommerNe92, author = {Karl D. Brommer and M. Needels and B. E. Larson and J. D. Joannopoulos}, title = {Ab initio theory of the {Si}(111)($7\times 7$) surface reconstruction: A challenge for massively parallel computation}, journal = {Physical Review Letters}, year = 1992, volume = 68, number = 9, pages = {13551358}, month = {March}, url = {http://link.aps.org/abstract/PRL/v68/p1355}, abstract = {An ab initio investigation of the Si(111)($7\times 7$) surface reconstruction is undertaken using the state of the art in massively parallel computation. Calculations of the total energy of an $\sim 700$ effectiveatom supercell are performed to determine (1) the fully relaxed atomic geometry, (2) the scanningtunnelingmicroscope images as a function of bias voltage, and (3) the energy difference between the ($7 \times 7$) and the ($2 \times 1$) reconstructions. The ($7 \times 7$) reconstruction is found to be energetically favorable to the ($2 \times 1$) surface by 60 meV per ($1 \times 1$) unit cell.} } @article{AriasPa92, author = {T. A. Arias and M. C. Payne and J. D. Joannopoulos}, title = {Ab initio molecular dynamics: Analytically continued energy functionals and insights into iterative solutions}, journal = {Physical Review Letters}, year = 1992, volume = 69, number = 7, pages = {10771080}, month = {August}, url = {http://link.aps.org/abstract/PRL/v69/p1077}, abstract = {We present a new method for performing finitetemperature ab initio totalenergy calculations at long length scales, which we demonstrate with a dynamics calculation of 50 $\AA$long phonon modes in silicon. The method involves \emph{both} a prescription for the analytic continuation of tradional fermionic energy functionals into the space of nonorthonormal singleparticle orbitals (speeding convergence to the minimum) \emph{and} insights into the common computational physics problem of solving by iterative refinement for the state of a complex system as a function of a continuous external parameter.} } @article{GalvanPi93, author = {Marcelo Galv{\'{a}}n and Dal Pino, Jr., Arnaldo and John D. Joannopoulos}, title = {Hardness and softness in the ab initio study of polyatomic systems }, journal = {Physical Review Letters}, year = 1993, volume = 70, pages = {2124}, month = {January}, url = {http://link.aps.org/abstract/PRL/v70/p21}, abstract = {A combination of the maximum hardness (MH) principle, and the local version of the hard and soft acids and bases (HSAB) principle is tested as a tool to describe the stability and reactivity of polyatomic systems. The local HSAB principle describes regional differences in reactivity whereas the MH principle selects the most stable configuration. This approach is tested using ab initio density functional theory techniques. Totalenergy calculations fully confirm the validity of both principles.} } @article{ChoJo93, author = {K. Cho and J. D. Joannopoulos and Leonard Kleinman}, title = {Constanttemperature molecular dynamics with momentum conservation}, journal = {Physical Review~E}, year = 1993, volume = 47, pages = {31453151}, month = {May}, url = {http://link.aps.org/abstract/PRE/v47/p3145}, abstract = {The Nos{\'{e}} theorem of the extendedsystem method of the constanttemperature molecular dynamics is generalized by including the conservation of the total virtual momentum. It is proved that a canonical ensemble of an ($N1$)particle system is generated from an extended system of an $N$particle system only if the total virtual momentum is zero. It is also shown that the resulting ($N1$)particle system has a slightly different mass spectrum than that of the original $N$particle system. The consequences of this new mass spectrum are relevant in the calculation of dynamical properties and the relaxation times of the system, but irrelevant to thermodynamic averages. For practical considerations, numerical simulations are performed and tested against this theorem. The differences in application of the Nos{\'{e}} theorem and the generalized No{\'{e}} theorem are discussed.} } @article{ChoJo93ergodicity, author = {K. Cho and J. D. Joannopoulos}, title = {Ergodicity and dynamical properties of constanttemperature molecular dynamics}, journal = {Physical Review~A}, year = 1992, volume = 45, pages = {70897097}, month = {May}, url = {http://link.aps.org/abstract/PRA/v45/p7089}, abstract = {The assumption of ergodicity in Nos{\'{e}}'s original formulation of the constanttemperature molecular dynamics is tested for a LennardJonespotential system. With the performance of very long simulations, it is shown that the extended system of a LennardJonespotential system is ergodic for all values of thermostat parameters tested. It is also shown, however, that the rate of convergence to the canonical ensemble strongly depends on the value of thermostat effective mass $Q$. The dynamical properties of the extended system are also studied using the velocity autocorrelation function and the power spectral density. From the analysis of the simulations, it is found that the dynamical properties are not correctly represented for arbitrary values of thermostat parameters. A prescription and a set of quantitative criteria are introduced to generate physically meaningful dynamics. Thus the results of this work show that with a special choice of thermostat parameters it is possible to obtain both the correct canonical ensemble and physically meaningful dynamical behavior of the physical system.} } @article{ChoJo93tip, author = {K. Cho and J. D. Joannopoulos}, title = {Tipsurface interactions in scanning tunneling microscopy}, journal = {Physical Review Letters}, year = 1993, volume = 71, pages = {13871390}, month = {August}, url = {http://link.aps.org/abstract/PRL/v71/p1387}, abstract = {The tipsurface interactions in scanning tunneling microscopy (STM) of the Si(100) surface are investigated with ab initio total energy pseudopotential calculations. The results of the calculations lead to a new understanding of the microscopic STM measurement process. It is found that under typical conditions the influence of the tip is large enough to effectively flip a dimer on this surface. This leads to a reinterpretation of the ``symmetric'' dimer STM image as an asymmetric dimer configuration that flips as it follows the motion of the scanning tip.} } @article{DevenyiCh93, author = {A. Devenyi and K. Cho and T. A. Arias and J. D. Joannopoulos}, title = {Adaptive {Riemannian} metric for allelectron calculations }, journal = {Physical Review~B}, year = 1994, volume = 49, pages = {1337313376}, month = {May}, url = {http://link.aps.org/abstract/PRB/v49/p13373}, abstract = {We present two techniques that make feasible the application of the adaptive Riemannian metric technique to allelectron localdensityfunctional calculations. The first overcomes both the real ($r=0$) and Fourier ($G=0$) space divergences of the nuclear Coulomb potential by computing the electronion energy as the smooth periodic electrostatic potential due to the electrons measured at the positions of the ions. The second overcomes the problem of slow convergence of the extreme metrics which the $r=0$ Coulomb divergence necessitates by giving an explicit prescription for a suitable metric for arbitrary ionic configurations. Allelectrondiamond calculations then serve as a proving ground for these ideas and demonstrate the viability of adaptive Riemannian methods for bypassing the pseudopotential approximation in solidstate calculations.} } @article{DalPinoGa93, author = {Dal Pino, Jr., A. and M. Galv{'{a}}n and T. A. Arias and J. D. Joannopoulos}, title = {Chemical softness and impurity segregation at grain boundaries}, journal = {Journal of Chemical Physics}, year = 1993, volume = 98, pages = {16061610}, month = {January}, note = {See erratum: ibid, vol. 98, p. 10106.}, url = {http://link.aip.org/link/?JCPSA6/98/1606/1}, abstract = {We analyze the process of impurity segregation at grain boundaries as a chemical reaction between the impurity and the interface. From this point of view, we test the ability of the concepts of local softness and hardness to predict the most probable sites for impurity accumulation. As a test, an ab initio investigation of the $\Sigma=5$ tilt [310] grain boundary in germanium is presented and the energetics of the accumulation of dopant atoms in this interface are studied. Our results support the utility of an analysis in terms of softness for impurity segregation problems.} } @article{BrommerLa93, author = {Karl D. Brommer and B. E. Larson and M. Needels and J. D. Joannopoulos}, title = {Modeling large surface reconstructions on the {Connection Machine}}, journal = {Japanese Journal of Applied Physics, Part 1}, year = 1993, volume = 32, number = {3B}, pages = {13601367}, month = {March}, url = {http://jjap.ipap.jp/link?JJAP/32/1360/}, abstract = {Using a massively parallel computer, we undertake an ab initio investigation of the Si(111)($7 \times 7$) surface reconstruction. Calculation of the total energy of an $\sim 700$ effective atom supercell at an 8 Ry plane wave cutoff allows us to determine (1) the energy difference between the ($7 \times 7$) and ($2 \times 1$) reconstructions, (2) the relaxed atomic geometry, (3) the scanning tunneling microscope topographs as a function of bias voltage and (4) the occupied and unoccupied electronic states.} } @article{ChoAr93, author = {K. Cho and T. A. Arias and J. D. Joannopoulos and Pui K. Lam}, title = {Wavelets in electronic structure calculations}, journal = {Physical Review Letters}, year = 1993, volume = 71, pages = {18081811}, month = {September}, url = {http://link.aps.org/abstract/PRL/v71/p1808}, abstract = {A threedimensional wavelet analysis is employed to develop a new formalism for electronic structure calculations. The wavelet formalism provides a systematically improvable and tractable description of electronic wave functions and overcomes limitations of conventional basis expansions. The potential power of the wavelet formalism for ab initio electronic structure calculations is demonstrated by a calculation of $1s$ states for all the naturally occurring nuclei on the periodic table and the interaction energies of the hydrogen molecule ion.} } @article{AriasJo94, author = {T. A. Arias and J. D. Joannopoulos}, title = {Ab initio theory of dislocation interactions: from closerange spontaneous annihilation to the longrange continuum limit}, journal = {Physical Review Letters}, year = 1994, volume = 73, pages = {680683}, month = {August}, url = {http://link.aps.org/abstract/PRL/v73/p680}, abstract = {Parallel supercomputer technology now permits \emph{ab initio} studies of systems of sufficient size to explore the interactions among dislocations in a solid. This study shows that the silicon shuffleset (110) screw dislocation is stable against spontaneous dissociation, provides an \emph{ab initio} value for the dislocation core energy, demonstrates a dislocationantidislocation interaction approaching the classical limit within a few tens of angstroms, and reveals a pathway for the spontaneous mutual annihilation of a dislocation dipole of the type that occurs when a FrankRead source emits a dislocation loop.} } @article{BrommerGa94, author = {Karl D. Brommer and M. Galv{\'{a}}n and Dal Pino, Jr., A. and J. D. Joannopoulos}, title = {Theory of adsorption of atoms and molecules on {Si}(111)($7\times 7$)}, journal = {Surface Science}, year = 1994, volume = 314, pages = {5770}, month = {July}, doi = {doi:10.1016/00396028(94)902127}, abstract = {Ab initio density functional theory calculations are performed to analyze a regional reactivity index based on chemical concepts of local softness and electronegativity. This leads to the prediction of a general qualitative pattern for the interaction of the Si(111)($7\times 7$) surface with atoms and molecules that preserve the reconstruction. We successfully compare this general pattern with experimental information available. Based on this approach, we predict the reactivity behavior of the reconstruction with a wide range of chemical reactants.} } @article{CapazLi95, author = {R. B. Capaz and H. Lim and J. D. Joannopoulos}, title = {Ab initio studies of {GaN} epitaxial growth on {SiC}}, journal = {Physical Review~B}, year = 1995, volume = 51, pages = {1775517757}, month = {June}, url = {http://link.aps.org/abstract/PRB/v51/p17755}, abstract = {Ab initio methods were used to investigate the initial stages of GaN epitaxial growth on (0001) 6HSiC. Total energies of four types of interfaces were calculated. Polarity matching at the interface plays a fundamental role in determining the lowerenergy structures, yielding strong binding for SiN and CGa interfaces and very weak binding for SiGa and CN. We therefore predict that Siterminated substrates will produce ideally Gaterminated films, whereas Cterminated substrates will produce ideally Nterminated films. This prediction suggests reinterpretation of recent experiments.} } @article{CapazCh95, author = {R. B. Capaz and K. Cho and J. D. Joannopoulos}, title = {Signatures of bulk and surface arsenic antisite defects in {GaAs}(110)}, journal = {Physical Review Letters}, year = 1995, volume = 75, pages = {18111814}, month = {August}, url = {http://link.aps.org/abstract/PRL/v75/p1811}, abstract = {Scanning tunneling microscopy (STM) has recently been used in the study of bulk arsenic antisite defects in GaAs. In this work, we report extensive theoretical calculations of such defects in the vicinity of a GaAs(110) surface, which provide essential information for the interpretation of experiments. Defects display remarkably distinct properties depending on whether they are fourfold or threefold coordinated. The nature of ``satellite peaks'' observed in experiment is elucidated. We predict the conditions under which STMmeasured properties will be faithful to the properties of the bulk defect.} } @article{ChoJo95, author = {K. Cho and J. D. Joannopoulos}, title = {Mechanical hysteresis on an atomic scale}, journal = {Surface Science}, year = 1995, volume = 328, number = 3, pages = {320324}, month = {May}, doi = {doi:10.1016/00396028(95)000437}, abstract = {A micromechanical hysteresis associated with intimate atomic force microscopy (AFM) on the Si(100) surface is discovered using ab initio total energy pseudopotential density functional calculations. It is predicted that it is possible to cycle repeatedly between two buckled configurations of a surface dimer along two different paths, one of which involves a discontinuous change of the equilibrium dimerangle, by varying appropriately the vertical movement of the AFM tip. This microscopic mechanical hysteresis effect should be detectable experimentally in low temperature AFM measurements.} } @article{ChoJo95tip, author = {K. Cho and J. D. Joannopoulos}, title = {Tipinduced modifications in {STM} and {AFM}}, journal = {Scanning Microscopy}, year = 1995, volume = 9, pages = {381} } @article{LimCh95, author = {H. Lim and K. Cho and I. Park and J. D. Joannopoulos and Efthimios Kaxiras}, title = {Ab initio study of hydrogen adsorption on the {Si}(111)($7\times 7$) surface}, journal = {Physical Review~B}, year = 1995, volume = 52, pages = {1723117237}, month = {December}, url = {http://link.aps.org/abstract/PRB/v52/p17231}, abstract = {Firstprinciples totalenergy pseudopotential calculations are performed to investigate the adsorption interaction of a H atom with dangling bonds on the Si(111)($7\times 7$) surface. he binding energies for adsorption of H at the adatom, rest atom, and corner hole sites are calculated to be 2.9, 3.2, and 3.5 eV, respectively. Spectral analysis of the electronic states shows that nonlocal changes of chemical reactivity are induced by \emph{charge transfer} upon H adsorption. It is found that H adsorption on the adatoms or rest atoms induces a charge transfer onto the SiH bond and a shift in energy of the remaining danglingbond states. Adsorption on the corner hole, however, does not involve any charge transfer. The relationship between charge transfer and binding energies is discussed.}, note = {Erratum: ibid., vol. 54, p. 5179 (1996).} } @article{LeeLe96, author = {Sik Lee and Sang Joo Lee and Jin Yong Lee and Jongseob Kim and Kwang S. Kim and Ickjin Park and K. Cho and J. D. Joannopoulos}, title = {Ab initio study of water hexamer anions}, journal = {Chemical Physics Letters}, year = 1996, volume = 254, pages = {128134}, month = {May}, doi = {doi:10.1016/00092614(96)003004}, abstract = {The wet electronan electron interacting with a small cluster of water moleculesis a simple yet fundamental system for understanding the behavior of electrons in complex molecular systems. A comprehensive post HartreeFock ab initio study is performed on the wet electron in various water hexamer clusters including the lowlying energy conformers of the neutral species. The predicted geometries, total energies, photoemission ionization energies, electronic structure and orbital character of the excess electron in ground and excited states are discussed. To understand the behavior of the excess electron in the clusters, the sorbitallike character of the HOMOs and the porbitallike character of the LUMOs are investigated.} } @article{ChoJo96, author = {K. Cho and J. D. Joannopoulos}, title = {Intrinsic surface atom manipulations in {STM} and {AFM}}, journal = {Applied Surface Science}, year = 1996, volume = {104105}, pages = {286290}, month = {September}, doi = {doi:10.1016/S01694332(96)001596}, abstract = {Stateoftheart firstprinciples simulations are performed to investigate the possibility of using a tungsten tip in atomic force microscopy (AFM) as a mechanical tool to manipulate the surface atoms of Si(100). Calculations of total energy and electronic structure are used to study the energetics and bonding properties associated with the tip and surface for a variety of atomic configurations. The results predict that under certain protocols the tip can be used effectively to flip dimers on the surface, from one buckled configuration to another, reversibly, and without inducing damage to either the intrinsic surface or the tip. This leads directly to the exciting possibility of using the intrinsic (100) surface of silicon as an ultrahigh memory storage device approaching the ultimate limit of one bit of data per atom.} } @article{MirbtMo99, author = {S. Mirbt and N. Moll and A. Kley and J. D. Joannopoulos}, title = {A general rule for surface reconstructions of {III}{V} semiconductors}, journal = {Surface Science}, year = 1999, volume = 422, number = {13}, pages = {L177L182}, month = {February}, doi = {doi:10.1016/S00396028(98)008723}, abstract = {First principles total energy calculations are performed for a large number (70) of IIIV semiconductor surfaces in order to establish a database from which a general rule is extracted to help isolate and predict the lowest energy atomic surface geometries for these complex systems. The general rule involves minimizing a single, material and geometryindependent, parameter, whose value depends only on a weighted sum of specific surface atom and bond structural units.} } @article{MollWa99, author = {N. Moll and Tairan Wang and Kyeongjae Cho and J. D. Joannopoulos}, title = {Semiconductor alloys for monolithic integration with {Si} microelectronics}, journal = {Materials Science & Engineering~B}, year = 1999, volume = 67, pages = {1722}, month = {December}, doi = {doi:10.1016/S09215107(99)002044}, abstract = {We introduce a novel class of semiconductors which is optically active and solves the problems with latticeconstant mismatches and polarity mismatches that are normally an issue in heteroepitaxial growth of IIIV alloys on silicon substrates. Using ab initio total energy calculations and quasiparticle GW calculations we examine physical properties of various configurations of these novel materials and identify a particular suitable configuration. This configuration is latticeconstant matched to Si and has a direct band gap very close to the operating wavelength of optical fibers. Therefore, this novel class could lead to monolithic integration of optical materials and Si circuits.} } @article{SolerGa01, author = {J. M. Soler and I. L. Garz{\'{o}}n and J. D. Joannopoulos}, title = {Structural patterns of unsupported gold clusters}, journal = {Solid State Communications}, year = 2001, volume = 117, number = 10, pages = {621625}, month = {February}, doi = {doi:10.1016/S00381098(00)004932}, abstract = {The structure of metal clusters is essential to predict many of their physical and chemical properties. Using first principles density functional calculations it was recently found that even ``magic'' cluster sizes, for which very compact and symmetric structures exist, have lowerenergy structures. The origin of these structures was shown to lie in the nonpairwise metallic interactions; while the compact ordered geometries are very stable for pair potentials, they are destabilized by the tendency of metallic bonds to contract at the surface. Here we identify important patterns of the resulting ``amorphous'' structures, showing why they are optimal for the metallic potential, and how they can be used to predict structures for other cluster sizes.} } @article{ParkCh01, author = {Ickjin Park and Kyeongjae Cho and Sik Lee and Kwang S. Kim and J. D. Joannopoulos}, title = {Ab initio atomistic dynamical study of an excess electron in water}, journal = {Computational Materials Science}, year = 2001, volume = 21, number = 3, pages = {291300}, month = {July}, doi = {doi:10.1016/S09270256(01)001525}, abstract = {The microscopic transport processes of an excess electron in bulk water are studied using hybrid ab initio molecular dynamics calculations. In contrast to the typical cavity obtained with solvated anions, the electron cavity structure is found to be much more variable, with water molecules easily exchanging at the surface of the cavity. The microscopic mechanism of electron transport involves a novel sequence of opportunistic electron redistributions driven by a positive feedback between thermal fluctuations and the attraction of the electron to hydrogen atoms that are not saturated in hydrogen bonding.} } @article{SkorobogatiyPa05, author = {M. Skorobogatiy and I. J. Park and J. D. Joannopoulos}, title = {The nature of a floating electron}, journal = {Computational Materials Science}, year = 2005, volume = 32, number = 1, pages = {96106}, month = {January}, doi = {doi:10.1016/j.commatsci.2004.04.007}, abstract = {Dynamics of electron solvation at the water surface is studied using extensive ab initio simulations. Calculations have been performed on semiclassical water at 300 K temperature with an excess quantum electron on its surface. It is found that after a very fast 3050 fs initial localization, there follow fast 5070 fs rotationally mediated transitions of the excess electron between surface trap states with a lifetime of 150400 fs. In less then 2 ps the excess electron gets trapped in an ordered ``floating'' electron state on the surface of water with a lifetime of more than 7 ps. The excess electron diffusion coefficient and spectrum of its velocity autocorrelation function change over time to reflect a transition from a very mobile phase (first 2 ps) to a trapped ``floating'' electron phase (27 ps).} } @article{WolfeNe92, author = {R. H. Wolfe and M. Needels and T. Arias and J. D. Joannopoulos}, title = {Visual revelations from silicon ab initio calculations}, journal = {IEEE Computer Graphics and Applications}, year = 1992, volume = 12, number = 4, pages = {4553}, month = {July}, doi = {doi:10.1109/38.144826 }, abstract = {The use of volumetric rendering for interpreting the computed electronic structure of solids is discussed. The structure, manifested as an atomicscale 3D electronic charge density, is computed from a pure mathematical model. The mathematical technique represents a class of ab initio computations that predict the atomic structure by solving the equations of quantum mechanics, wherein all configuration parameters except the atomic number are under the analyst's control. Some physics aspects of electronic structure of solids are reviewed. The use of this visualization technique to interpret the computations for silicon with an oxygen impurity and for roomtemperature silicon is described.} } @article{AllanJo80, author = {Douglas C. Allan and J. D. Joannopoulos}, title = {Electronic Structure of Hydrogenated Amorphous Silicon }, journal = {Physical Review Letters}, year = 1980, volume = 44, pages = {4347}, month = {January}, url = {http://link.aps.org/abstract/PRL/v44/p43}, abstract = {alculations of the electronic states and total energies of various bonding conformations in hydrogenated amorphous Si are presented. Various surprising results emerge, including identification of peaks in photoemission spectra as signatures of nearestneighbor SiH configurations, a gap increasing with increasing H content while the conduction band remains essentially unchanged, and localized states in the gap arising from various defects whose energies are in sharp contrast to recently proposed simple model estimates.} } @article{PollardJo80, author = {W. B. Pollard and J. D. Joannopoulos}, title = {Vibrational excitations of defect sites in amorphous group{V} semiconductors}, journal = {Physical Review~B}, year = 1980, volume = 21, pages = {760766}, month = {January}, url = {http://link.aps.org/abstract/PRB/v21/p760}, abstract = {A study of the vibrational excitations of structural defects in amorphous pnictide semiconductors is presented. Amorphous As is chosen as a prototype for this study. The local phonon densities of states, infrared absorption, and Ramanscattering spectra of twofold coordinated As atoms, fourfold coordinated As atoms, intimate valence alternation pairs, and vacancy defects are presented. The effects of relaxations, disorders, and distortions upon the vibrational states of the defects are discussed. The results suggest that only the relaxed twofold coordinated As atoms (or vacancylike defects) give rise to localized vibrational modes with optical activities consistent with sharp weak features observed in the experimental infrared and Raman spectra of bulk amorphous As. The implications of our results for amorphous P and Sb, as well as the relationships between spectroscopic measurements and other optical measurements are also discussed.} } @article{LaughlinJo80, author = {R. B. Laughlin and J. D. Joannopoulos and D. J. Chadi}, title = {Theory of the electronic structure of the {Si}{SiO}$_2$ interface}, journal = {Physical Review~B}, year = 1980, volume = 21, pages = {57335744}, month = {June}, url = {http://link.aps.org/abstract/PRB/v21/p5733}, abstract = {A theory of the SiSiO$_2$ interface based on recent experimental findings for silicon surfaces and their oxidation is presented. It is proposed that a simple localorbital picture can simultaneously describe silicon, its oxidation, and the SiSiO$_2$ interface and that two dimensionality is not essential to physically meaningful calculations of interface local densities of states. Calculations are performed which show that interface states do not arise simply from the presence of a boundary. It is argued that band tailing at the interface, like that in amorphous silicon, is due primarily to strain rather than to charged centers, and that dangling bonds at the interface should give rise to an inhomogeneously broadened discrete level at midgap.} } @article{VanderbiltJo80glassy, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Theory of defect states in glassy selenium}, journal = {Physical Review~B}, year = 1980, volume = 22, pages = {29272939}, month = {September}, url = {http://link.aps.org/abstract/PRB/v22/p2927}, abstract = {A realistic approach to the electronic theory of bondcoordination defects in chalcogenides, based on selfconsistent pseudopotential calculations, is used to study glassy Se. The results of the pseudopotential calculations are interpreted in terms of simpler tightbinding models. The onefold and threefold coordination defects are both found to give rise to nondegenerate, nonhydrogenic gap states, whose properties are unique to the chalcogenides in several respects. The existence of $\pi$ interactions between nonbonding orbitals at defect sites is found to be crucial to an understanding of the electronic structure. These interactions are responsible for large charge transfers between atoms and consequently large energy shifts of atomic valence orbitals, which make these defects quite unlike those in other semiconductors. The origin, character, energy location, and localization of the defect states associated with bondcoordination defects, and with defect pairs and certain relaxed defects, are discussed.} } @article{VanderbiltJo81, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Theory of defect states in glassy {As}$_2${Se}$_3$}, journal = {Physical Review~B}, year = 1981, volume = 23, pages = {25962606}, month = {March}, url = {http://link.aps.org/abstract/PRB/v23/p2596}, abstract = {Structural defects in glassy As$_2$Se$_3$ are classified and labeled according to the constituent likeatom bonds and malcoordinated atoms. A selection rule is formulated to reduce the number of allowed Fermilevel pinning reactions. An elementary Bethelattice model is introduced as a starting point for a discussion of the electronic structure of simple defects. The defect states are found to be very different from those in Se; deep gap states arise in As$_2$Se$_3$ because of unique bond orbitals, whereas they occur in Se due to unique $\pi$ interactions between orbitals. Malcoordinated Se atoms are expected to give rise to \emph{hydrogenic} levels in As$_2$Se$_3$, in contrast to Se. Surprisingly, the undercoordinated pnictide defects are \emph{positively} charged in this model. Finally, defect creation energies and densities at $T_g$ are discussed.} } @article{LeeJo81I, author = {D. H. Lee and J. D. Joannopoulos}, title = {Simple scheme for surfaceband calculations. {I}}, journal = {Physical Review~B}, year = 1981, volume = 23, pages = {49884996}, month = {May}, url = {http://link.aps.org/abstract/PRB/v23/p4988}, abstract = {A simple, efficient scheme for calculating the electronic structure of a surface is presented. The scheme is applicable to any general Hamiltonian that can be described within a localizedorbital basis. The method is much faster than the current techniques available. The basic concept is that of wavefunction matching through a transfer matrix. The eigensolutions of this matrix then provide all the information concerning the projected band structure, surfacestate energies, orbital character, and decay lengths. A rather detailed discussion of the formalism is presented for a general surface system. A comprehensive and illustrative example is also presented for readers who are interested in using the scheme but not in the details of the theory.}, note = {Erratum: Phys. Rev. Lett., vol. 52, p. 1054 (1984).} } @article{LeeJo81II, author = {D. H. Lee and J. D. Joannopoulos}, title = {Simple scheme for surfaceband calculations. {II}. The {Green}'s function}, journal = {Physical Review~B}, year = 1981, volume = 23, pages = {49975004}, month = {May}, url = {http://link.aps.org/abstract/PRB/v23/p4997}, abstract = {We present a very simple scheme for calculating the Green's function of a semiinfinite surface system described within a localized orbital basis. By generating a series of matching conditions for the Green's function we can calculate its matrix elements much faster than any method currently available. We present the formalism for a specific class of systems and include a simple example to illustrate the use of the technique.} } @article{PollardJo81, author = {W. B. Pollard and J. D. Joannopoulos}, title = {Vibrational properties of amorphous silicon alloys}, journal = {Physical Review~B}, year = 1981, volume = 23, pages = {52635268}, month = {May}, url = {http://link.aps.org/abstract/PRB/v23/p5263}, abstract = {A theoretical model of the vibrational states of recently reported amorphous siliconfluorine alloys is presented. The infrared absorption, polarized Raman scattering, and phonon state density spectra are calculated for different local siliconfluorine bonding arrangements. Changes in the infrared and Raman spectra associated with the incorporation of oxygen and hydrogen in the amorphous alloy is also discussed. Our results suggest that unlike the situation in amorphous siliconhydrogen alloys, the presence of SiF bonds in the fluorinated materials significantly disturb the bulk Si vibrational bands. It is found that the signatures of isolated SiF$_2$ and SiF$_3$ complexes are stretching modes in the region of 800 to 1000 cm$^{1}$ accompanied by one (SiF$_2$) or two (SiF$_3$) sharp resonances in the region of 400 to 500 cm$^{1}$, while the SiF unit is distinguished by a single stretching mode near 800 cm$^{1}$. Finally, our results also indicate that SiH bonds in close proximity with SiF bonds alter the fluorine vibrational modes. SiO bonds, on the other hand, do not affect the SiF stretching mode near 800 cm$^{1}$.} } @article{IhmJo81, author = {J. Ihm and J. D. Joannopoulos}, title = {Structural Energies of {A1} Deposited on the {GaAs}(110) Surface}, journal = {Physical Review Letters}, year = 1981, volume = 47, pages = {679682}, month = {August}, url = {http://link.aps.org/abstract/PRL/v47/p679}, abstract = {Major controversies over the structure of A1 atoms deposited on the GaAs(110) surface have been resolved by using a firstprinciples energyminimization method. The most stable configuration is that of A1 atoms replacing the second layer Ga atoms beneath the surface. At temperatures where this reaction cannot be activated, two important processes are found to exist. In the lowcoverage limit, A1 atoms favor twofold sites to form strong bonds with the substrate atoms. At higher coverages, A1 atoms tend to cluster and make new bonds among themselves.} } @article{MeleJo81, author = {Eugene J. Mele and J. D. Joannopoulos}, title = {Electronic structure of the zincblende and rocksalt phases of {InSb}}, journal = {Physical Review~B}, year = 1981, volume = 24, pages = {31453154}, month = {September}, url = {http://link.aps.org/abstract/PRB/v24/p3145}, abstract = {The electronic structure of InSb in the common zincblendecrystal phase and in a rocksaltcrystal phase (which is metastable at standard temperature and pressure) are investigated using a selfconsistent pseudopotential formalism including relativistic effects. For the zincblende structure we find that a local $s$$p$ potential for the valence electrons yields, in a selfconsistent calculation for the solid, a charge density in excellent agreement with previous calculations employing empirical nonlocal potentials. Relativistic effects are found to be very important in order to obtain a good description of the band gap and overall bandwidth. For the rocksalt phase we obtain a metallic solid, in agreement with experiment, and observe (in comparison with the zincblende results) substantial changes in the valenceband density of states. These results are in very good agreement with the experimental xrayphotoemissionspectroscopy studies of these two phases. Unlike the situation for the covalently bonded zincblende crystal, we obtain very large charge transfer from the cations to the anions (estimated to be 0.9$e^$) in the metallic rocksalt phase, which we speculate helps to stabilize the solid. Bandstructure, densitiesofstates, chargedensity, and Fermisurface results are presented.} } @article{StoneJo81, author = {A. Douglas Stone and J. D. Joannopoulos}, title = {Probability distribution and new scaling law for the resistance of a onedimensional {Anderson} model}, journal = {Physical Review~B}, year = 1981, volume = 24, pages = {35923595}, month = {September}, url = {http://link.aps.org/abstract/PRB/v24/p3592}, abstract = {The exact probability distributions of the resistance $\rho$, the conductance, and $\ln(1+\rho)$ are calculated for the 1D Anderson model with purely offdiagonal disorder at $E=0$. Analysis of the distribution yields the surprising results that $\langle \rho \rangle$ grows exponentially with length, despite previous studies indicating that the state at $E=0$ is extended, and the typical resistance $\tilde{\rho}$ increases as $\exp(L^{1/2})$. The relationship between this behavior and the temperature dependence of the resistance of thin wires is discussed.} } @article{HoJo82, author = {KaiMing Ho and J. Ihm and J. D. Joannopoulos}, title = {Dielectric matrix scheme for fast convergence in selfconsistent electronicstructure calculations}, journal = {Physical Review~B}, year = 1982, volume = 25, pages = {42604262}, month = {March}, url = {http://link.aps.org/abstract/PRB/v25/p4260}, abstract = {A scheme is devised which drastically reduces the number of iterations required to reach selfconsistency in electronicstructure calculations. This scheme is particularly helpful in calculations for systems with large unit cells.} } @article{LeeJo82, author = {D. H. Lee and J. D. Joannopoulos}, title = {Simple Scheme for Deriving Atomic Force Constants: Application to {SiC}}, journal = {Physical Review Letters}, year = 1982, volume = 48, pages = {18461849}, month = {June}, url = {http://link.aps.org/abstract/PRL/v48/p1846}, abstract = {A scheme is presented for deriving atomic force constants that is fairly accurate, yet simple enough to allow determination of the dynamical matrix of bulk solids, relaxed and reconstructed surfaces, defects, etc. For illustrative purposes, the method is applied to bulk SiC to predict the phonon dispersion curves throughout the entire Brillouin zone, and to a carbon vacancy in SiC to study the change of atomic force constants due to atomic relaxations.} } @article{VanderbiltJo82, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Bonding Coordination Defect in {g}{Se}: A ``Positive{$U$}'' System}, journal = {Physical Review Letters}, year = 1982, volume = 49, pages = {823826}, month = {September}, url = {http://link.aps.org/abstract/PRL/v49/p823}, abstract = {The effective Hubbard $U$ for the bonding coordination defect in glassy Se is investigated. This is accomplished by applying localdensity totalenergy calculations directly to charged defects. The existence of a sizable negative contribution to $U$ from structural relaxation is confirmed. However, a still larger Coulomb repulsion gives rise to an overall positive $U$. This result is not inconsistent with a negative $U$ in the compound chalcogenide glasses.} } @article{VanderbiltJo82occupation, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Offdiagonal occupation numbers in localdensity theory }, journal = {Physical Review~B}, year = 1982, volume = 26, pages = {32033210}, month = {September}, url = {http://link.aps.org/abstract/PRB/v26/p3203}, abstract = {We introduce a new method for specifying the occupations of states in localdensity theory by allowing ``offdiagonal occupation numbers'' (i.e., a generalized density matrix). We show that this technique has important applications to superlattice and slab geometries, allowing bondingtype interactions between defects, molecules, or surfaces to be effectively eliminated. Moreover, previously inaccessible charge states may now be studied. The method is shown to have the full sanction of localdensity theory.} } @article{IhmJo82PRB, author = {J. Ihm and J. D. Joannopoulos}, title = {Structure of the {Al}{GaAs}(110) interface from an energyminimization approach}, journal = {Physical Review~B}, year = 1982, volume = 26, pages = {44294435}, month = {October}, url = {http://link.aps.org/abstract/PRB/v26/p4429}, abstract = {The structure of the AlGaAs(110) interface at the initial stages of Al deposition has been studied. In this system, the adatomsubstrate and the adatomadatom interactions are not necessarily weak, and the interface region may be strongly perturbed. The pseudopotentialenergy minimization method is ideally suited for determining the stable structure of such systems involving significant energy change, and this method is successfully applied to the AlGaAs system in the present study. By the calculation of the chemisorption energy of the Al atom at different sites on the entire GaAs surface, an energy contour map is obtained for the first time. Favorable channels of migration of Al atoms on this map are identified, which suggests preferential surface diffusion along Ga theAs bonding chain direction. The calculated surface hopping diffusion is sufficiently high even at room temperature so that Al atoms frequently encounter and interact with other Al atoms to form Al clusters. The Al clusters are found to have lower free energy than chemisorbed Al atoms. Al atoms eventually replace substrate Ga atoms and Al formAs bonds in the interface region to minimize the energy of the system.} } @article{StoneAl83, author = {A. Douglas Stone and Douglas C. Allan and J. D. Joannopoulos}, title = {Phase randomness in the onedimensional {Anderson} model}, journal = {Physical Review~B}, year = 1983, volume = 27, pages = {836843}, month = {January}, url = {http://link.aps.org/abstract/PRB/v27/p836}, abstract = {The probability distribution $P(\nu)$ of the phase $\nu$ which determines the scaling properties of the dimensionless resistance $R/T$ is analyzed for the onedimensional Anderson model with diagonal disorder. For weak disorder ($W/V \ll 1$) and $E \neq 0$ we find that $P(\nu)$ is initially sharply peaked for short chains but becomes uniform over $2\pi$ as $N$ increases. Moreover, the length at which $P(\nu)$ becomes uniform is found to be independent of $W/V$ (as long as $W/V \ll 1$); instead it appears to depend only on the distance in energy from band center, i.e., on the extent to which the electron's wavelength is incommensurate with the lattice. Thus, unlike the disordered KronigPenny model, there is a phase randomness length, shorter than the localization length but longer than the lattice spacing. $P(\nu)$ is also studied in the limit of large disorder and analytic expressions for the mean and variance of the inverse localization length are derived for weak and strong disorder. Our results suggest that uniform phase randomness is indeed a property of most disordered 10 systems in the limit of weak disorder.} } @article{IhmCh83, author = {J. Ihm and D. J. Chadi and J. D. Joannopoulos}, title = {Study of the reconstructed {GaAs}(100) surface}, journal = {Physical Review~B}, year = 1983, volume = 27, pages = {51195121}, month = {April}, url = {http://link.aps.org/abstract/PRB/v27/p5119}, abstract = {The reconstruction of an Asterminated GaAs(100) surface has been studied using the selfconsistent pseudopotential method. Total energies for the $(1 \times 1)$ ideal surface and the $c(2 \times 2)$ and $p(2\times 2)$ econstructed surfaces within the dimer model are compared. Unlike the Si(100) surface, at least two inequivalent dimers are required to produce the semiconducting surface and stabilize the system. Important features in the density of states and the valencecharge distribution are presented.} } @article{VanderbiltJo83I, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Total energies in {Se}. {I}. The trigonal crystal}, journal = {Physical Review~B}, year = 1983, volume = 27, pages = {62966301}, month = {May}, url = {http://link.aps.org/abstract/PRB/v27/p6296}, abstract = {The equilibrium crystal structure (lattice constants and chain radius) of trigonal selenium, a molecular crystal, are determined by minimizing the total energy, as calculated in the localdensity and frozencore approximations. The cohesive energy and $\Gamma_1(A_1)$ phonon frequency are also computed. Comparison with experiment shows excellent agreement for intrachain properties, and satisfactory agreement for interchain properties. This indicates that ab initio localdensity totalenergy calculations are viable for the case of molecular crystals.} } @article{VanderbiltJo83II, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Total energies in {Se}. {II}. Vacancy in the crystal}, journal = {Physical Review~B}, year = 1983, volume = 27, pages = {63026310}, month = {May}, url = {http://link.aps.org/abstract/PRB/v27/p6302}, abstract = {Ab initio totalenergy calculations are used to determine theoretically the structural configuration of the vacancy in trigonal Se. The method consists of calculating the forces, as well as the total energies, within the localdensity and frozencore approximations, for a superlattice structure containing a vacancy. In this way, relaxations at the vacancy can be fully taken into account, including a possible selfhealing of the vacancy. A slightly relaxed symmetric version of the ideal vacancy is the lowestenergy structure found; neither asymmetric relaxation nor valence alternation appears to occur. A simple Hubbard Hamiltonian is used to analyze the spin configuration of the lowestenergy structure.} } @article{VanderbiltJo83III, author = {David Vanderbilt and J. D. Joannopoulos}, title = {Total energies in {Se}. {II}. Defects in the glass}, journal = {Physical Review~B}, year = 1983, volume = 27, pages = {63116321}, month = {May}, url = {http://link.aps.org/abstract/PRB/v27/p6311}, abstract = {The effective Hubbard $U$ for the bonding coordination defect in glassy Se is investigated. This is accomplished by applying localdensity totalenergy calculations directly to charged defects in a superlattice configuration. The existence of a large negative contribution to $U$, arising from interconversion between danglingbond and threefoldcoordinated structures, is confirmed. However, a still larger Coulomb repulsion gives rise to an overall \emph{positive} $U$.} } @article{IhmLe83phase, author = {J. Ihm and D. H. Lee and J. D. Joannopoulos and J. J. Xiong}, title = {Structural Phase Diagrams for the Surface of a Solid: A TotalEnergy, RenormalizationGroup Approach }, journal = {Physical Review Letters}, year = 1983, volume = 51, pages = {18721875}, month = {November}, url = {http://link.aps.org/abstract/PRL/v51/p1872}, abstract = {Totalenergy calculations based on microscopic electronic structure are combined with positionspace renormalizationgroup calculations to predict the structural phase transitions of the Si(100) surface as a function of temperature. It is found that two distinct families of reconstructed geometries can exist on the surface, with independent phase transitions occurring within each. Two critical temperatures representing orderdisorder transitions are calculated.} } @article{LeeJo84, author = {D. H. Lee and J. D. Joannopoulos}, title = {Magnetic properties of the {Si}(111) unreconstructed surface}, journal = {Physical Review~B}, year = 1984, volume = 29, pages = {14721473}, month = {February}, url = {http://link.aps.org/abstract/PRB/v29/p1472}, abstract = {The solution of the twodimensional Hubbard Hamiltonian for describing the magnetic properties of the Si(111) unreconstructed surface is discussed.} } @article{LeeJo84xy, author = {D. H. Lee and J. D. Joannopoulos and J. W. Negele and D. P. Landau}, title = {DiscreteSymmetry Breaking and Novel Critical Phenomena in an Antiferromagnetic Planar ({$XY$}) Model in Two Dimensions}, journal = {Physical Review Letters}, year = 1984, volume = 42, pages = {433436}, month = {February}, url = {http://link.aps.org/abstract/PRL/v52/p433}, abstract = {LandauGinzburgWilson symmetry analyses and Monte Carlo calculations for the classical antiferromagnetic planar ($XY$) model on a triangular lattice reveal a wealth of interesting critical phenomena. From this simple model arise a zerofield transition to a state of longrange order, a new mechanism for spin disordering, and a critical point associated with a possible new universality class.} } @article{LeeCa84, author = {D. H. Lee and R. G. Caflisch and J. D. Joannopoulos and F. Y. Wu}, title = {Antiferromagnetic classical {$XY$} model: A meanfield analysis}, journal = {Physical Review~B}, year = 1984, volume = 29, pages = {26802684}, month = {March}, url = {http://link.aps.org/abstract/PRB/v29/p2680}, abstract = {The antiferromagnetic classical $XY$ (planarrotator) model is analyzed under the meanfield approximation. Phase diagrams are obtained and found to be strongly dependent on the underlying lattice geometry. For bipartite lattices, there exists a secondorder transition across a unique phase boundary. For tripartite lattices, there exist two phase boundaries, separating an intermediate ``nonhelical'' phase from a lowtemperature ``helical'' phase and the hightemperature paramagnetic phase. The two phase boundaries merge into a single critical point at finite temperature and zero magnetic field.} } @article{BarYamJo84, author = {Y. BarYam and J. D. Joannopoulos}, title = {Barrier to Migration of the Silicon SelfInterstitial }, journal = {Physical Review Letters}, year = 1984, volume = 52, pages = {11291132}, month = {March}, url = {http://link.aps.org/abstract/PRL/v52/p1129}, abstract = {The first totalenergy calculations of barriers to interstitial migration have been used to study silicon selfinterstitial migration. Migration occurs through the lowelectrondensity path. Relaxation was found to be important in determining the barrier for both Si$^{(0)}$ and Si$^{(++)}$. Electronassisted migration has been demonstrated. Si$^{(++)}$ was found to have lower energy at the tetrahedral site while Si$^{(0)}$ has lower energy at the hexagonal site.} } @article{LeeJo84monte, author = {D. H. Lee and J. D. Joannopoulos and J. W. Negele}, title = {{Monte} {Carlo} solution of antiferromagnetic quantum {Heisenberg} spin systems }, journal = {Physical Review~B}, year = 1984, volume = 30, pages = {15991602}, month = {August}, url = {http://link.aps.org/abstract/PRB/v30/p1599}, abstract = {A Monte Carlo method is introduced that overcomes the problem of alternating signs in Handscomb's method of simulating antiferromagnetic quantum Heisenberg systems. The scheme is applied to both bipartite and frustrated lattices. Results of internal energy, specific heat, and uniform and staggered susceptibilities are presented suggesting that quantum antiferromagnets may now be studied as extensively as classical spin systems using conventional Monte Carlo techniques.} } @article{BarYamJo84PRB, author = {Y. BarYam and J. D. Joannopoulos}, title = {Electronic structure and totalenergy migration barriers of silicon selfinterstitials }, journal = {Physical Review~B}, year = 1984, volume = 30, pages = {18441852}, month = {August}, url = {http://link.aps.org/abstract/PRB/v30/p1844}, abstract = {Pseudopotential densityfunctional calculations have been performed on silicon selfinterstitial supercell model geometries to yield extensive information on electronic states. We present band structures, charge densities, and densities of states to identify and characterize electronic states associated with silicon selfinterstitials in the geometries studied. Total energies obtained yield migration barriers for both Si$^{(0)}$ and Si$^{(2+)}$ interstitials. We also present the results of preliminary totalenergy relaxation studies and show their effects on electronic states and totalenergy calculations, demonstrating the importance of relaxation in determining migration barriers. Electronassisted migration is shown to occur, thus solving the mystery of the disappearing selfinterstitial and providing an initial understanding of migration in lowtemperature irradiated silicon.} } @article{BarYamJo84multiple, author = {Y. BarYam and J. D. Joannopoulos}, title = {Silicon selfinterstitial migration: Multiple paths and charge states}, journal = {Physical Review~B}, year = 1984, volume = 30, pages = {22162218}, month = {August}, url = {http://link.aps.org/abstract/PRB/v30/p2216}, abstract = {Total energy calculations performed to study the migration of silicon selfinterstitials reveal a great complexity of migration paths. Si$^{(++)}$ and Si$^{(0)}$ each migrate along more than one path and different paths from each other. Some paths involve interchange with bulk atoms. Electronassisted transport occurs for the Si$^{(++)}$ stable tetrahedral site in all directions. Si$^{(0)}$ has several almostdegenerate lowestenergy configurations. The lowest is the "exchange" configuration which also has the possibility of being negative (Si$^{()}$).} } @article{RabeJo85, author = {K. M. Rabe and J. D. Joannopoulos}, title = {Ab initio relativistic pseudopotential study of the zerotemperature structural properties of {SnTe} and {PbTe}}, journal = {Physical Review~B}, year = 1985, volume = 32, pages = {23022314}, month = {August}, url = {http://link.aps.org/abstract/PRB/v32/p2302}, abstract = {The equilibrium lattice parameters, bulk moduli, cohesive energies, $k=0$ TOphonon frequencies and elastic constant $C_{44}$ of SnTe and PbTe are calculated completely from first principles with use of the localdensityfunctional pseudopotential totalenergy method including relativistic effects. Good agreement with experiment is obtained. The approximation of neglecting spinorbit coupling is found to be adequate for most purposes. Band structures and valence charge densities are presented. The latter are seen to differ significantly from previous empiricalpseudopotentialmethod calculations. The form of $\omega_\mathrm{TO(P)}$ in PbTe is found to be linear for small pressures $P$. Some possibilities for further applications of this method to these materials are discussed.} } @article{LeeJo86, author = {D. H. Lee and J. D. Joannopoulos and J. W. Negele and D. P. Landau}, title = {Symmetry analysis and {Monte} {Carlo} study of a frustrated antiferromagnetic planar ($XY$) model in two dimensions }, journal = {Physical Review~B}, year = 1986, volume = 33, pages = {450475}, month = {January}, url = {http://link.aps.org/abstract/PRB/v33/p450}, abstract = {The antiferromagnetic planar ($XY$) model on a triangular lattice is investigated using grouptheoretical symmetry arguments combined with extensive Monte Carlo simulations and a detailed finitesizescaling analysis. This approach allows for a systematic exploration of all possible symmetrybreaking transitions and their associated critical behavior. The entire magneticfieldversustemperature phase diagram is deduced. A rich class of possible new critical phenomena is obtained, including the introduction of a new multicritical point describing the confluence of the Ising and KosterlitzThouless universality classes. The existence of the critical point is associated with a domainwallinduced vortexantivortexunbinding transition.} } @article{AntonelliTa86, author = {A. Antonelli and Eugen Tarnow and J. D. Joannopoulos}, title = {New insight into the electronic structure of a {As}$_2${Se}$_3$}, journal = {Physical Review~B}, year = 1986, volume = 33, pages = {29682971}, month = {February}, url = {http://link.aps.org/abstract/PRB/v33/p2968}, abstract = {Results of a realistic ab initio investigation of the geometric structure and electronic properties of the compound chalcogenide prototype As$_2$Se$_3$ are presented. Structural parameters are accurate to within a few percent of their experimental values. The nature of bonding within and between layers is revealed and the character of states near the band edges is obtained, including predictions for the electron and hole effective masses.} } @article{KaxirasBarYam86PRB, author = {E. Kaxiras and Y. BarYam and J. D. Joannopoulos and K. C. Pandev}, title = {$(2\times 2)$ reconstructions of the $\{111\}$ polar surfaces of {GaAs}}, journal = {Physical Review~B}, year = 1986, volume = 33, pages = {44064409}, month = {March}, url = {http://link.aps.org/abstract/PRB/v33/p4406}, abstract = {Ab initio totalenergy calculations were used to examine $(2 \times 2)$ reconstruction models for the (111) and $(\bar{1}\bar{1}\bar{1})$ surfaces of GaAs. For the (111) surface the lowestenergy Gavacancy geometry is determined; several mechanisms for Gavacancy formation are examined and other reconstructions are discussed. For the $(\bar{1}\bar{1}\bar{1})$ surface it is shown that the Asvacancy model is unlikely and other geometries are considered.} } @article{BarYamJo86, author = {Y. BarYam and J. D. Joannopoulos}, title = {Dangling Bond in {a}{Si}:{ZZH}}, journal = {Physical Review Letters}, year = 1986, volume = 56, pages = {22032206}, month = {May}, url = {http://link.aps.org/abstract/PRL/v56/p2203}, abstract = {The configurationcoordinate diagram for a model danglingbond geometry is obtained by ab initio totalenergy calculations. The results yield equilibrium geometries for different charge states, thermodynamic and optical transition energies, large Stokes shifts (0.50.6 eV), and a \emph{negative} effective correlation energy ($0.2 \pm 0.2$ eV). Implications of these results for the current experimental and theoretical understanding are discussed.} } @article{KaxirasPa86, author = {E. Kaxiras and K. C. Pandey and Y. BarYam and J. D. Joannopoulos}, title = {Role of chemical potentials in surface reconstruction: A new model and phase transition of {GaAs}(111)$2 \times 2$}, journal = {Physical Review Letters}, year = 1986, volume = 56, pages = {28192822}, month = {June}, url = {http://link.aps.org/abstract/PRL/v56/p2819}, abstract = {The role of chemical potentials in surface reconstructions is examined and shown to be crucial for binary semiconductor srufaces such as GaAs(111)$2 \times 2$. We predict that under Asrich conditions a new model, the As triangle, is the lowestenergy geometry, whereas the Gavacancy model is appropriate for Garich conditions. A change in the relative chemical potential of Ga and As should produce a phase transition between the two structures.} } @article{KaxirasBarYam86, author = {E. Kaxiras and Y. BarYam and J. D. Joannopoulos and K. C. Pandey}, title = {Variable stoichiometry surface reconstructions: New models for {GaAs}$(\bar{1}\bar{1}\bar{1})(2\times 2)$ and $(\sqrt{19}\times\sqrt{19})$}, journal = {Physical Review Letters}, year = 1986, volume = 57, pages = {106109}, month = {July}, url = {http://link.aps.org/abstract/PRL/v57/p106}, abstract = {The $(\bar{1}\bar{1}\bar{1})$ surface of GaAs exhibits three stable reconstructions. Two are $(2\times 2)$, As stabilized and Ga stabilized respectively, and the third is $(\sqrt{19}\times\sqrt{19})$. Transitions between these structures are obtained by variation of the experimental conditions. We propose new models for all of the above reconstructions, based on ab initio totalenergy calculations and experimental information regarding surface composition.} } @article{BarYamAd86, author = {Y. BarYam and D. Adler and J. D. Joannopoulos}, title = {Structure and Electronic States in Disordered Systems }, journal = {Physical Review Letters}, year = 1986, volume = 57, pages = {467470}, month = {July}, url = {http://link.aps.org/abstract/PRL/v57/p467}, abstract = {We adopt the view that the structure of disordered systems is itself determined by the formation energy of structural deviations from an ideally bonded network. By developing an equilibrium statistical mechanical model for disorder we predict properties of electronic band tails, defects, defect densities of states, and impurity doping.} } @article{TarnowAn86electronic, author = {Eugen Tarnow and A. Antonelli and J. D. Joannopoulos}, title = {Crystalline {As}$_2${Se}$_3$: Electronic and geometric structure}, journal = {Physical Review~B}, year = 1986, volume = 34, pages = {40594073}, month = {September}, url = {http://link.aps.org/abstract/PRB/v34/p4059}, abstract = {Realistic ab initio totalenergy calculations are performed to investigate the electronic and geometric structure of crystalline As$_2$Se$_3$. Results include the following: total energies for various distortions, the equilibrium geometric structure, intralayer and interlayer cohesive energies, charge densities, the density of states, the band structure, rigidlayer phonon frequencies, and effective masses. Specific attention is focused on studying the nature of bonding within and between layers, as well as elucidating the nature of the electron and hole wave functions in the vicinity of the fundamental gap. In particular it is found that the interactions between the layers are crucial in determining the properties of the wave functions at the band edges. Moreover, the electron and hole masses are predicted to exhibit a large and unusual anisotropy.} } @article{TarnowAn86optical, author = {Eugen Tarnow and A. Antonelli and J. D. Joannopoulos}, title = {Crystalline {As}$_2${Se}$_3$: Optical properties}, journal = {Physical Review~B}, year = 1986, volume = 34, pages = {87188727}, month = {December}, url = {http://link.aps.org/abstract/PRB/v34/p8718}, abstract = {Realistic ab initio totalenergy calculations are performed to investigate the optical properties of crystalline As$_2$Se$_3$. The complete dielectric tensor is calculated using interband transitions. This includes the offdiagonal element, which has not previously been considered in the literature. The presence of the offdiagonal element leads to optical axes which are not aligned with the crystal axes and whose directions change with frequency. Comparisons with experimental results on diagonal elements are made and a prediction of the offdiagonal reflectivity is given. The origin of the offdiagonal element is investigated, and the element is found to arise solely from interlayer interactions. Accordingly, experimental measurements of the offdiagonal element may reveal unique and interesting features of the environment in between the layers.} } @article{PayneBr87, author = {M. C. Payne and P. D. Bristowe and J. D. Joannopoulos}, title = {Ab initio determination of the structure of a grain boundary by simulated quenching}, journal = {Physical Review Letters}, year = 1987, volume = 58, pages = {13481351}, month = {March}, url = {http://link.aps.org/abstract/PRL/v58/p1348}, abstract = {Results of the first completely ab initio investigation of the microscopic structure of a grain boundary in a semiconductor are presented. By use of the moleculardynamicssimulated annealing method for performing totalenergy calculations within the localdensityfunctional and pseudopotential approximation, the $\Sigma=5$ (001) twist boundary in germanium is studied. A number of rotationandtranslation states are investigated leading to a prediction for the structure of this geometry. Evidence for the possible presence of novel defects and glasslike tunneling mode states at grain boundaries is presented.} } @article{NeedelsPa87, author = {M. Needels and M. C. Payne and J. D. Joannopoulos}, title = {Ab initio molecular dynamics on the {Ge}(100) surface}, journal = {Physical Review Letters}, year = 1987, volume = 58, pages = {17651768}, month = {April}, url = {http://link.aps.org/abstract/PRL/v58/p1765}, abstract = {An ab initio totalenergy approach to molecular dynamics and simulated quenching is applied for the first time and used to study the structure of the Ge(100) surface. The geometry of the lowestenergy $c(4\times2)$symmetry dimer model is predicted. A soft energy surface for displacement of dimers in the $[1\bar{1}0]$ direction is discovered which could lead to a new type of symmetry breaking.} } @article{KaxirasBarYam87I, author = {E. Kaxiras and Y. BarYam and J. D. Joannopoulos and K. C. Pandey}, title = {Ab initio theory of polar semiconductor surfaces. I. Methodology and the $(2\times2)$ reconstructions of {GaAs}(111)}, journal = {Physical Review~B}, year = 1987, volume = 35, pages = {96259635}, month = {June}, url = {http://link.aps.org/abstract/PRB/v35/p9625}, abstract = {A methodology is developed for the theoretical study of the polar surfaces of compound semiconductors. It is based on the calculation of the total energy in the context of densityfunctional theory in the pseudopotential approximation. The method is used to investigate the $(2\times2)$ reconstructions of GaAs(111). Emphasis is given to the relative chemical potential, which plays a crucial role in determining the lowestenergy geometry for surfaces with different stoichiometries. The totalenergy versus chemicalpotential curves indicate that there are at least two stable reconstructions. We predict one to be an Astriangle geometry and the other the Ga vacancy.} } @article{KaxirasBarYam87II, author = {E. Kaxiras and Y. BarYam and J. D. Joannopoulos and K. C. Pandey}, title = {Ab initio theory of polar semiconductor surfaces. II. $(2\times2)$ reconstructions and related phase transitions of {GaAs}$(\bar{1}\bar{1}\bar{1})$}, journal = {Physical Review~B}, year = 1987, volume = 35, pages = {96369643}, month = {June}, url = {http://link.aps.org/abstract/PRB/v35/p9636}, abstract = {The $(2\times2)$ reconstructions of GaAs$(\bar{1}\bar{1}\bar{1})$ are studied with use of a theoretical approach based on the calculation of the total energy in the context of densityfunctional theory and the pseudopotential approximation. New models are proposed for the Asrich and Garich reconstructions. The relative chemical potential plays a crucial role in determining the lowestenergy configuration. The totalenergy versus chemicalpotential curves indicate the possibility of phase transitions between different configurations. One such transition concerning the experimentally observed $(sqrt{19}\times\sqrt{19})$ reconstruction can be explained as an intermediate phase between the proposed lowenergy $(2\times2)$ reconstructions.} } @article{RabeJo87, author = {K. M. Rabe and J. D. Joannopoulos}, title = {Ab initio determination of a structural phase transition temperature}, journal = {Physical Review Letters}, year = 1987, volume = 59, pages = {570573}, month = {August}, url = {http://link.aps.org/abstract/PRL/v59/p570}, abstract = {The temperature of the structural phase transition in GeTe is calculated complete ab initio from a model latticedynamical Hamiltonian that is constructed from microscopic quantummechanical totalenergy calculations. The phase transition in the model system is studied through a renormalizationgrouptheory approach, leading to the prediction of a fluctuationdriven firstorder transition at $657 \pm 100$ K.} } @article{RabeJo87T0, author = {K. M. Rabe and J. D. Joannopoulos}, title = {Structural properties of {GeTe} at {$T=0$}}, journal = {Physical Review~B}, year = 1987, volume = 36, pages = {33193324}, month = {August}, url = {http://link.aps.org/abstract/PRB/v36/p3319}, abstract = {The equilibrium lattice parameters, bulk modulus, shear elastic constant $C_{44}$, and cohesive energy of GeTe are calculated with use of the ab initio scalar relativistic pseudopotential in the localdensity approximation. The convergence of the structural properties in basisset size and Brillouinzone averaging is discussed. For the rocksalt structure, valence charge densities and the band structure are presented. For the rhombohedral structure, the band structure is presented and the shear deformation potential is calculated. Good agreement of calculated quantities with available experimental results is obtained.} } @article{RabeJo87PRB, author = {K. M. Rabe and J. D. Joannopoulos}, title = {Theory of the structural phase transition of {GeTe}}, journal = {Physical Review~B}, year = 1987, volume = 36, pages = {66316639}, month = {October}, url = {http://link.aps.org/abstract/PRB/v36/p6631}, abstract = {A completely ab initio theoretical investigation of the rocksaltrhombohedral structural phase transition of GeTe is described. Starting from an anharmonic lattice Hamiltonian, a model Hamiltonian that includes coupling of the order parameter to longwavelength strain is constructed. The parameters appearing in the model are calculated using the selfconsistent ab initio pseudopotential totalenergy method. The phase transition in the model system is studied through a momentumspace renormalizationgrouptheory approach, leading to the prediction of a fluctuation driven firstorder transition at $657 \pm 100$ K. The strain coupling is found to be crucial in determining the firstorder character of the transition. A discussion of approximations made in this approach and considerations relevant to its application to structural transitions in other systems is included.} } @article{GomezSantosJo87, author = {G. GomezSantos and J. D. Joannopoulos}, title = {Application of spinwave theory to the ground state of {$XY$} quantum Hamiltonians}, journal = {Physical Review~B}, year = 1987, volume = 36, pages = {87078711}, month = {December}, url = {http://link.aps.org/abstract/PRB/v36/p8707}, abstract = {Spinwave theory is successfully applied to study the ground state of the quantum $XY$ Hamiltonian. It is found that a judicious choice of the quantized spin axis removes difficulties of previous applications. Results for the simplest one, two, and threedimensional lattices are shown to compare favorably with either exact or other approximate calculations.} } @article{PayneNe88, author = {M. C. Payne and M. Needels and J. D. Joannopoulos}, title = {Symmetry breaking in the moleculardynamics method for ab initio totalenergy calculations }, journal = {Physical Review~B}, year = 1988, volume = 37, pages = {81388144}, month = {May}, url = {http://link.aps.org/abstract/PRB/v37/p8138}, abstract = {An analysis of the symmetry of the electron wave functions as they evolve under the moleculardynamics equation of motion is presented. It is found that symmetry in the initial conditions, corresponding to wave functions of two or more bands being partner functions, can prevent wave functions from converging to eigenfunctions of the Hamiltonian. This artificial symmetry is broken by the GramSchmidt orthogonalization procedure so any convenient initial conditions will lead to convergence to eigenstates of the Hamiltonian.} } @article{KaxirasJo88, author = {E. Kaxiras and J. D. Joannopoulos}, title = {Hydrogenation of semiconductor surfaces: Si and Ge (111)}, journal = {Physical Review~B}, year = 1988, volume = 37, pages = {88428848}, month = {May}, url = {http://link.aps.org/abstract/PRB/v37/p8842}, abstract = {The relaxations of hydrogenated Si and Ge (111) surfaces are determined using ab initio selfconsistent calculations in a slab configuration. The SiH and GeH bonds are found to be considerably larger than the sum of covalent radii. The substrate relaxations are small and their physical origin can be explained in terms of electronic charge transfer which eliminates the surface dipole moment. The calculated frequencies of the hydrogen vibrational modes are in excellent agreement with experiment. A surfaceatom vibrational mode is compared to similar modes in the amorphous hydrogenated materials. The comparison predicts that internal surfaces (microvoids) in the amorphous network are locally much softer than the corresponding crystalline surface configuration.} } @article{WangBarYam88, author = {Xiaomei Wang and Y. BarYam and D. Adler and J. D. Joannopoulos}, title = {{dc} conductivity and the {Meyer}{Neldel} rule in {a}{Si}:{H}}, journal = {Physical Review~B}, year = 1988, volume = 38, pages = {16011604}, month = {July}, url = {http://link.aps.org/abstract/PRB/v38/p1601}, abstract = {Disorder in amorphous semiconductors results in unusual properties of dc conductivity. We demonstrate a quantitative description of the temperature dependence of conductivity in aSi:H. The universal activation energy dependence of the conductivity prefactor (the MeyerNeldel rule) is reproduced. Excellent agreement with experimental results is obtained by describing disorder and defects using the general thermodynamic ensemble theory for the structure of disordered systems.} } @article{NeedelsPa88, author = {M. Needels and M. C. Payne and J. D. Joannopoulos}, title = {Highorder reconstructions of the {Ge}(100) surface}, journal = {Physical Review~B}, year = 1988, volume = 38, pages = {55435546}, month = {September}, url = {http://link.aps.org/abstract/PRB/v38/p5543}, abstract = {We have performed ab initio densityfunctional calculations of total energies of the Ge(100) surface to compare the ground states of $(2\times1)$, $c(4\times2)$, $p(2\times2)$, and $p(4\times1)$ symmetry dimer reconstructions. We find that $p(2\times2)$ and $c(4\times2)$ are the lowestenergy reconstructions and are nearly degenerate in energy. From these ab initio total energies, we compute the coupling constants for a model Hamiltonian for the surface and predict the phasetransition temperature from either an ordered $c(4\times2)$ or $p(2\times2)$ state to a disordered buckled $b(2\times1)$ state.} } @article{TarnowPa88, author = {Eugen Tarnow and M. C. Payne and J. D. Joannopoulos}, title = {Pairing of Electrons by a Point Defect in $c${As}$_2${Se}$_3$}, journal = {Physical Review Letters}, year = 1988, volume = 61, pages = {17721775}, month = {October}, url = {http://link.aps.org/abstract/PRL/v61/p1772}, abstract = {By use of ab initio densityfunctional totalenergy calculations a point defect with a negative Hubbard $U$ is identified, for the first time, in a chalcogenide system. It is found that in $c$As$_2$Se$_3$ two noninteracting neutral Se antisite defects are unstable toward the formation of a pair of noninteracting oppositely charged Se antisites. The correlation energy gained is 0.3 eV per electron pair.} } @article{AlerhandVa88, author = {O. L. Alerhand and David Vanderbilt and Robert D. Meade and J. D. Joannopoulos}, title = {Spontaneous Formation of Stress Domains on Crystal Surfaces}, journal = {Physical Review Letters}, year = 1988, volume = 61, pages = {19731976}, month = {October}, url = {http://link.aps.org/abstract/PRL/v61/p1973}, abstract = {It is proposed that surfaces of crystals which reconstruct with broken orientational symmetry and exhibit an anisotropic intrinsic stress tensor are unstable to the formation of elasticstress domains. Thus the ground state of such a surface is not uniform, contrary to previous expectations. Evidence of this for the Si(001) surface is discussed.} } @article{BarYamPa89, author = {Y. BarYam and S. T. Pantelides and J. D. Joannopoulos}, title = {Ab initio pseudopotential solidstate calculations of highly electronegative firstrow elements}, journal = {Physical Review~B}, year = 1989, volume = 39, pages = {33963399}, month = {February}, url = {http://link.aps.org/abstract/PRB/v39/p3396}, abstract = {Densityfunctional pseudopotential totalenergy calculations have been a powerful tool in studies of the electronic and structural properties of semiconductors and simple metals. The inclusion of highly electronegative firstrow elements such as oxygen and fluorine has not, in general, been tractable. In this paper we address the inherent difficulties and demonstrate that they can be overcome. Focusing on oxygensilicon systems the precision of our method is demonstrated particularly through the study of angular stability of SiOSi bridge bonds in silica. The results predict bending to an equilibrium angle of $140^\circ \pm 5^\circ$, as observed, with an energy gain of only $\sim 0.15$ eV.} } @article{AriasJo89, author = {T. A. Arias and J. D. Joannopoulos}, title = {Reexamination of magnetic effects in the {Bose} gas}, journal = {Physical Review~B}, year = 1989, volume = 39, pages = {40714078}, month = {March}, url = {http://link.aps.org/abstract/PRB/v39/p4071}, abstract = {The breakdown of the MeissnerOchsenfeld effect in the threedimensional Bose gas as the applied field passes through its critical value is an entropydriven weakly firstorder transition, rather than the secondorder transition usually ascribed to the system. The transition is second order at the usual Bose condensation temperature $T_c$ as well as at $T=0$, with a line of firstorder transitions connecting these critical points. The firstorder transitions make the Bose gas resemble familiar superconductors, and a LandauGinzburg analysis indicates that the Bose gas is always a typeI superconductor.} } @article{GomezSantosJo89, author = {G. GomezSantos and J. D. Joannopoulos and J. W. Negele}, title = {{Monte} {Carlo} study of the quantum spin1/2 {Heisenberg} antiferromagnet on the square lattice}, journal = {Physical Review~B}, year = 1989, volume = 39, pages = {44354443}, month = {March}, url = {http://link.aps.org/abstract/PRB/v39/p4435}, abstract = {A new Monte Carlo sampling technique is used to study the spin1/2 Heisenberg antiferromagnet on square lattices of up to 256 spins. The energy, specific heat, uniform and staggered susceptibilities, and correlation function are calculated. The correlation length as a function of temperature is deduced from correlation functions and is shown to be consistent with a quantumrenormalized classical approximation. Using the coupling constant provided by Raman scattering experiments, good agreement is found between the calculated and experimental correlation lengths in La$_2$CuO$_4$.} } @article{TarnoJo89, author = {Eugen Tarnow and J. D. Joannopoulos and M. C. Payne}, title = {Antisites, antistructures, and bondswitching reactions in layered chalcogenides}, journal = {Physical Review~B}, year = 1989, volume = 39, pages = {60176024}, month = {March}, url = {http://link.aps.org/abstract/PRB/v39/p6017}, abstract = {State of the art ab initio totalenergy calculations are used to investigate the properties of a large number of antisite and antistructure defects in $c$As$_2$Se$_3$. Calculations of their formation energies are made, and an examination of their electronic states and relaxed geometries is performed. Surprisingly it is found that a select set of the antistructures undergoes bondswitching reactions that give rise to a new peculiar class of defects. These bondswitched antistructure defects are found to have lower formation energies than any of the antistructures. One of these is of particular interest as it forms a bridging bond between the layers of $c$As$_2$Se$_3$. This crosslinking defect (CLD) is examined in more detail. It is found that the associated bondswitching reaction proceeds with small or no barrier for most charge states. This is surprising as the atomic movements associated with the formation of the CLD are quite large.} } @article{KaxirasAl89, author = {Efthimios Kaxiras and O. L. Alerhand and J. D. Joannopoulos and G. W. Turner}, title = {Microscopic model of heteroepitaxy of {GaAs} on {Si}(100)}, journal = {Physical Review Letters}, year = 1989, volume = 62, pages = {24842486}, month = {May}, url = {http://link.aps.org/abstract/PRL/v62/p2484}, abstract = {A new microscopic model of heteroepitaxial growth is introduced using GaAs on Si(100) as a prototype. This model takes into account specific features of surface topology, predicts that in the prototype system conventional twodimensional epitaxy should be inhibited, and provides a fundamental explanation for threedimensional nature of the initial stages of growth. The ingredients of the model, which are supported by totalenergy calculations, include new structural geometries for each state of growth and the chemical and rehybridization reactions linking these stages.} } @article{AlerhandJo89, author = {O. L. Alerhand and J. D. Joannopoulos and E. J. Mele}, title = {Thermal amplitudes of surface atoms on {Si}(111) $2\times1$ and {Si}(001) $2\times1$}, journal = {Physical Review~B}, year = 1989, volume = 39, pages = {1262212629}, month = {June}, url = {http://link.aps.org/abstract/PRB/v39/p12622}, abstract = {Atomic displacementdisplacement correlation functions for the surface atoms of the Si(111) $2\times1$ and Si(001) $2\times1$ surfaces are calculated using phonon energies and eigenvectors obtained from a tightbinding theory of the lattice dynamics of Si. The temperature dependence and the possible anisotropies of the amplitudes of vibration of the surface atoms are investigated. For the $\pi$bonded chain model of Si(111) $2\times1$ it is found that the surface atoms have larger vibrational amplitudes perpendicular to the surface than parallel to it, and this anisotropy (1.8 in the meansquared amplitudes at 270 K) increases rapidly with increasing temperature. This is attributed to the existence of a smallwavelength acoustic phonon with an unusually low energy. For the surface atoms on Si(001) $2\times1$ the anisotropy between perpendicular and parallel vibrations is not as marked. The interatomic vibrational correlations between surface atoms reveal the strength of the new bonds that are formed in the reconstruction of these surfaces.} } @article{ShirleyAl89, author = {Eric L. Shirley and Douglas C. Allan and Richard M. Martin and J. D. Joannopoulos}, title = {Extended normconserving pseudopotentials}, journal = {Physical Review~B}, year = 1989, volume = 40, pages = {36523660}, month = {August}, url = {http://link.aps.org/abstract/PRB/v40/p3652}, abstract = {Atomic pseudopotentials simplify electronic calculations by eliminating atomic core levels and the potentials that bind them. Outside some core radius, normconserving pseudopotentials produce the same scattering properties (radial logarithmic derivatives of wave functions for angular momenta of interest) as fullatomic potentials to zeroth and first order in energy about valencelevel eigenvalues. We extend the correctness of the radial logarithmic derivative one order further in energy and present analytic and numerical results showing that this extension improves higherorder energy derivatives as well. We also show how our change improves predictions of excited singleparticle eigenvalues in a wide variety of atoms, as well as highenergy scattering properties, with effects visible in a bandstructure calculation. Our potentials converge nearly as quickly in reciprocal space as the Vanderbilt (modified HamannSchl{\"{u}}terChiang) potentials from which they are derived, and are easily generated.} } @article{RappeRa90, author = {Andrew M. Rappe and Karin M. Rabe and Efthimios Kaxiras and J. D. Joannopoulos}, title = {Optimized pseudopotentials}, journal = {Physical Review~B}, year = 1990, volume = 41, pages = {12271230}, month = {January}, url = {http://link.aps.org/abstract/PRB/v41/p1227}, abstract = {A planewave basis has great advantages for many calculations in the physics solids. To apply this basis to a wider class of materials, the atomic characteristic of a pseudopotential is identified which leads to rapid convergence in the solid, and a new method for generating pseudopotentials optimized according to this criterion is shown. As a test case, an ab initio planewave basis determination of the structural properties of fcc copper is performed. The results indicate that these optimized pseudopotentials will facilitate study of transition metals and firstrow nonmetals.} } @article{AlerhandBe90, author = {O. L. Alerhand and A. Nihat Berker and J. D. Joannopoulos and David Vanderbilt and R. J. Hamers and J. E. Demuth}, title = {Finitetemperature phase diagram of vicinal {Si}(100) surfaces }, journal = {Physical Review Letters}, year = 1990, volume = 64, pages = {24062409}, month = {May}, url = {http://link.aps.org/abstract/PRL/v64/p2406}, abstract = {The phase diagram of vicinal Si(100) as a function of misorientation angle and temperature is calculated. Contrary to previous suggestions that only doublelayer steps should appear on the equilibrium surface, it is predicted that the singlelayer stepped surface is at equilibrium for small misorientation angles. This structure is stabilized by strain relaxation and by the thermal roughening of the steps. For annealed surfaces the critical angle at which the transition between the single and doublelayer stepped surface occurs is calculated to be $\theta_c \approx 2^\circ$.} } @article{TarnowDa90, author = {E. Tarnow and P. Dallot and P. D. Bristowe and J. D. Joannopoulos and G. P. Francis and M. C. Payne}, title = {Structural complexity in grain boundaries with covalent bonding}, journal = {Physical Review~B}, year = 1990, volume = 42, pages = {36443657}, month = {August}, url = {http://link.aps.org/abstract/PRB/v42/p3644}, abstract = {The structural properties of two shortperiod twist boundaries in germanium are explored using a stateoftheart totalenergy calculation. The structures of these boundaries are found to be very complex, with boundary bonds that are distorted and weak. These systems are found to exhibit a large degeneracy in the number of local energy minima. Thus the boundaries have difficulty in arriving at a locally ordered state. The situation may be unique to the semiconducting twist grain boundaries due to the inherent frustration present between the tendency to form directional bonds and the imposed twist geometry which makes the bond formation improbable. This study focuses on the energy, coordination, volume change, and electronic states characteristic of the local minima. A trend towards dimerization is found especially in the highestangle twist boundary.} } @article{NeedelsJo91, author = {M. Needels and J. D. Joannopoulos and Y. BarYam and S. T. Pantelides}, title = {Oxygen complexes in silicon}, journal = {Physical Review~B}, year = 1991, volume = 43, pages = {42084215}, month = {February}, url = {http://link.aps.org/abstract/PRB/v43/p4208}, abstract = {Ab initio totalenergy calculations are performed to investigate the initial stages of oxygen aggregation in silicon. Since the volume per SiO$_2$ unit in silicon dioxide is approximately twice the volume per silicon atom in crystalline silicon, it is generally believed that oxygen clustering in silicon will require the creation of selfinterstitials. In a similar vein, it is generally believed that neighboring interstitial oxygen configurations are unbound in a weakly distorted silicon lattice. In this work, several interesting phenomena that are associated with aggregation and do not require the creation of silicon selfinterstitials are predicted. These predictions are (1) oxygen atoms can cluster with large binding energies of about 1 eV; (2) this clustering can only occur for very specific geometries; and (3) the geometries are chainlike arrangements of oxygen bridging configurations. Indeed, it is the quasionedimensional nature of the proposed oxygen clusters that allows the complexes to form without the accompanying large lattice stress that can lead to the ejection of silicon interstitials. The thermal treatment necessary for aggregation is identified, and predictions of the observable effects are made.} } @article{AlerhandWa91stepped, author = {O. L. Alerhand and Jing Wang and J. D. Joannopoulos and Efthimios Kaxiras and R. S. Becker}, title = {Adsorption of {As} on stepped {Si}(100): Resolution of the sublatticeorientation dilemma }, journal = {Physical Review~B}, year = 1991, volume = 44, pages = {65346537}, month = {September}, url = {http://link.aps.org/abstract/PRB/v44/p6534}, abstract = {Firstprinciples calculations are used to investigate the energetics of an As overlayer adsorbed on a stepped Si(110) surface. We show that the growth of As directly on top of the Si surface produces a metastable structure, while the replacement of the original top Si layer by As leads to a lowerenergy configuration. In the latter case, the rearrangement of the surface is driven by the relaxation of stress by surface steps. This result explains the sublatticeorientation dilemma in GaAsonSi heteroepitaxy.} } @article{AriasPa92largescale, author = {T. A. Arias and M. C. Payne and J. D. Joannopoulos}, title = {Ab initio moleculardynamics techniques extended to largelengthscale systems}, journal = {Physical Review~B}, year = 1992, volume = 45, pages = {15381549}, month = {January}, url = {http://link.aps.org/abstract/PRB/v45/p1538}, abstract = {The BornOppenheimer approximation divides the problem of quantum molecular dynamics into two familiar problems: (1) solution for the electronic wave functions for a given instantaneous arrangement of ions and (2) the motion of the atomic cores under the influence of those wave functions. A combination of conjugategradient methods to solve (1) with standard molecular dynamics to solve (2) results in a scheme that is at least two orders of magnitude more accurate than previously possible, thus allowing accurate calculation of dynamic correlation functions while maintaining tolerable energy conservation for microcanonical averages of those correlation functions over picosecond time scales. By employing conjugategradient techniques, this method is used to extend the applicability of finitetemperature ab initio techniques to systems with large length scales.} } @article{DalPinoNe92, author = {Dal Pino, Jr., A. and M. Needels and J. D. Joannopoulos}, title = {Oxygeninduced brokenbond defect in silicon}, journal = {Physical Review~B}, year = 1992, volume = 45, pages = {33043308}, month = {February}, url = {http://link.aps.org/abstract/PRB/v45/p3304}, abstract = {The electronic and geometric structure of a recently proposed oxygenrelated defect in silicon is studied using a precise Iab initioP quantummoleculardynamics method based on localdensityfunctional theory, nonlocal pseudopotentials, and the supercell approximation. This defect is a metastable configuration of interstitial oxygen in Si (Si:O$_i$), and it is closely related to a vacancyinterstitial pair. This defect has a relatively low formation energy despite the presence of a silicon broken bond. A simple and straightforward procedure is used to identify the static barrier for the formation of an intimate vacancyinterstitial pair in silicon and the static barrier for the recombination of the broken bond in metastable Si:O$_i$. The similarities in geometric configurations for both cases are explored to introduce a mechanism for the creation of intimate vacancyinterstitiallike defects in silicon.} } @article{RappeDalPino92, author = {Andrew M. Rappe and al Pino, Jr., A. and Mark Needels and J. D. Joannopoulos}, title = {Mixedbasis pseudopotential method applied to iterative diagonalization techniques}, journal = {Physical Review~B}, year = 1992, volume = 46, pages = {73537357}, month = {September}, url = {http://link.aps.org/abstract/PRB/v46/p7353}, abstract = {We apply a mixedbasis formulation of pseudopotential totalenergy electronic structure calculations in the context of iterative diagonalization techniques. The formulation combines a small set of auxiliary functions to describe the localized part of the wave function with a planewave basis set. The method is tested on lowsymmetry configurations of interstitial oxygen in silicon. This method provides accuracy comparable to planewave basisset calculations and requires less computational effort. Therefore, it appears to be a promising tool for the description of the electronic structure of systems with localized valence electrons.} } @article{PayneTe92, author = {M. C. Payne and M. P. Teter and D. C. Allan and T. A. Arias and J. D. Joannopoulos}, title = {Iterative minimization techniques for ab initio totalenergy calculations: molecular dynamics and conjugate gradients}, journal = {Reviews of Modern Physics}, year = 1992, volume = 64, pages = {10451097}, month = {October}, url = {http://link.aps.org/abstract/RMP/v64/p1045}, abstract = {This article describes recent technical developments that have made the totalenergy pseudopotential the most powerful ab initio quantummechanical modeling method presently available. In addition to presenting technical details of the pseudopotential method, the article aims to heighten awareness of the capabilities of the method in order to stimulate its application to as wide a range of problems in as many scientific disciplines as possible.} } @article{NeedelsRa92, author = {M. Needels and A. M. Rappe and P. D. Bristowe and J. D. Joannopoulos}, title = {Ab initio study of a grain boundary in gold }, journal = {Physical Review~B}, year = 1992, volume = 46, pages = {97689771}, month = {October}, url = {http://link.aps.org/abstract/PRB/v46/p9768}, abstract = {The total energy of a grain boundary in a transition metal is calculated in an ab initio manner. Such calculations are now feasible using a planewave basis set provided that an optimally convergent pseudopotential is used. The investigation focuses on resolving two competing atomic models for the $\Sigma=5$ [001] twist boundary in gold. It is found that the model of lower energy, which involves small atomic displacements, corresponds to a structure determined both experimentally, using quantitative xraydiffraction techniques, and theoretically, using the embeddedatom method.} } @article{GalvanDalPino93, author = {Marcelo Galv{\'{a}}n and Dal Pino, Jr., Arnaldo and John D. Joannopoulos}, title = {Hardness and softness in the ab initio study of polyatomic systems}, journal = {Physical Review Letters}, year = 1993, volume = 70, pages = {2124}, month = {January}, url = {http://link.aps.org/abstract/PRL/v70/p21}, abstract = {A combination of the maximum hardness (MH) principle, and the local version of the hard and soft acids and bases (HSAB) principle is tested as a tool to describe the stability and reactivity of polyatomic systems. The local HSAB principle describes regional differences in reactivity whereas the MH principle selects the most stable configuration. This approach is tested using ab initio density functional theory techniques. Totalenergy calculations fully confirm the validity of both principles.} } @article{WangAr93, author = {Jing Wang and T. A. Arias and J. D. Joannopoulos and G. W. Turner and O. L. Alerhand}, title = {Scanningtunnelingmicroscopy signatures and chemical identifications of the (110) surface of Sidoped GaAs}, journal = {Physical Review~B}, year = 1993, volume = 47, pages = {1032610334}, month = {April}, url = {http://link.aps.org/abstract/PRB/v47/p10326}, abstract = {The electronic properties of the (110) surface of both pure and Sidoped bulk GaAs are studied using firstprinciples totalenergy calculations within the localdensity functional and pseudopotential approximations. The wave functions of the relaxed configurations are used to generate theoretical scanningtunnelingmicroscopy (STM) images. For the clean surface, the buckling angle of the surface GaAs bond is found to be $26^\circ$ and the theoretically generated STM images are in good agreement with those obtained from experiment. For the Sidoped GaAs(110) surface, the extra electron of the Si substitutional at a Ga site on the surface is found to be welllocalized around the Si atom. In addition, danglingbond states of surface As atoms bordering the Si substitutional are found to be altered due to the distinctively different chemical property of the Si substitutionals. These features should act as a signature for the location of the substitutional surface Si atoms via either voltagedependent STM imaging or currentvoltage ($I$$V$) measurements of various positions on the surface.} } @article{WangAr93dimer, author = {Jing Wang and T. A. Arias and J. D. Joannopoulos}, title = {Dimer vacancies and dimervacancy complexes on the Si(100) surface}, journal = {Physical Review~B}, year = 1993, volume = 47, pages = {1049710508}, month = {April}, url = {http://link.aps.org/abstract/PRB/v47/p10497}, abstract = {Ab initio totalenergy calculations are performed using a conjugategradients minimization technique to calculate the properties of various dimer vacancies and dimervacancy complexes on the Si(100) surface. A dimer vacancy is found to be a much more probable intrinsic defect for this surface than the recently proposed dimer interstitial. Two mechanisms are identified that contribute to the low formation energy and stability of dimer vacancies on the surface: (i) the need to eliminate dangling bonds in the second layer and (ii) the need for atoms to readjust in order to relax the strain. The results of the calculations provide a good quantitative explanation of the major features of scanningtunnelmicroscopy images obtained from the Si(100) surface. These include the abundance of dimer vacancies, the tendency of dimer vacancies to cluster, and the details of the nonthermal population distribution of distinctive vacancy complexes.} } @article{DalPinoRa93, author = {Dal Pino, Jr., A. and Andrew M. Rappe and J. D. Joannopoulos}, title = {Ab initio investigation of carbonrelated defects in silicon}, journal = {Physical Review~B}, year = 1993, volume = 47, pages = {1255412557}, month = {May}, url = {http://link.aps.org/abstract/PRB/v47/p12554}, abstract = {Ab initio totalenergy calculations based on the localdensityfunctional, pseudopotential, and supercell approximations are performed to investigate carbon defects in silicon. The geometry and the formation energy of substitutional and impurityvacancy defects are studied including the relaxation of nearest and nextnearest neighbors. Results for substitutional carbon appear to be consistent with a recently suggested reinterpretation of the available experimental formation energy data. Results for the interaction energy between a carbon atom and a silicon vacancy predict a small binding energy of 0.19 eV.} } @article{AriasJo94grain, author = {T. A. Arias and J. D. Joannopoulos}, title = {Electron trapping and impurity segregation without defects: Ab initio study of perfectly rebonded grain boundaries}, journal = {Physical Review~B}, year = 1994, volume = 49, pages = {45254531}, month = {February}, url = {http://link.aps.org/abstract/PRB/v49/p4525}, abstract = {We present the results of an extensive ab initio study of the $\Sigma=5$ [310] grain boundary in germanium. We find that the boundary reliably reconstructs to the tetrahedrally bonded network observed in highresolution electron microscopy experiments without the proliferation of false local minima observed in similar twist boundaries. The reduced density of bonds crossing the grainboundary plane leads us to conjecture that the boundary may be a preferred fracture interface. Though there are no dangling bonds or miscoordinated sites in the reconstruction, the boundary presents electrontrap states just below the conduction band. Further, we show that lattice relaxation effects are irrelevant to the segregation of impurities to tetrahedrally reconstructed defects and that the interfacial electrontrap states give rise to an electronic frustration mechanism that selectively drives the segregation of only ntype dopants to the boundary.} } @article{CapazDalPino94, author = {R. B. Capaz and Dal Pino, Jr., A. and J. D. Joannopoulos}, title = {Identification of the migration path of interstitial carbon in silicon}, journal = {Physical Review~B}, year = 1994, volume = 50, pages = {74397442}, month = {September}, url = {http://link.aps.org/abstract/PRB/v50/p7439}, abstract = {We have performed ab initio totalenergy calculations of groundstate properties and migration paths of interstitial carbon in silicon. The ground state involves threefoldcoordinated carbon and silicon atoms and its geometry suggests primarily $p$ and $sp$ bonding for carbon, rather than $sp^2$ one would naively expect. Examination of possible migration paths reveals that only three correspond to small ``jumps'' involving a single ``bond breaking.'' Of these, we predict that only one has a barrier of considerably lower energy ($\sim 0.5$ eV) and involves an intermediate ``saddlepoint'' configuration of $C_2$ symmetry} } @article{ChoJo96vicinal, author = {K. Cho and J. D. Joannopoulos and A. Nihat Berker}, title = {Vicinal {Si}(100) surfaces under external strain }, journal = {Physical Review~B}, year = 1996, volume = 53, pages = {10021005}, month = {January}, url = {http://link.aps.org/abstract/PRB/v53/p1002}, abstract = {The phase diagram of vicinal Si(100) surfaces is calculated as a function of misorientation angle, temperature, and applied external strain. It is shown that a change of the applied external strain can lead to a phase transition between the singlelayer step surface phase and the doublelayer step surface phase. The effect of temperature is shown to be negligible up to 300 K. The order parameter of the surface phase transition is also calculated to make contact with experimental measurements.} } @article{KimPa96, author = {Kwang S. Kim and Ickjin Park and Sik Lee and K. Cho and Jin Yong Lee and Jongseob Kim and J. D. Joannopoulos}, title = {The Nature of a Wet Electron}, journal = {Physical Review Letters}, year = 1996, volume = 76, pages = {956959}, month = {February}, url = {http://link.aps.org/abstract/PRL/v76/p956}, abstract = {A comprehensive stateoftheart ab initio study is performed on the wet electronn electron interacting with a small cluster of molecules waterin the water hexamer system. Predictions include two previously unknown distinctive geometries which bind the excess electron as internal and external states, photoemission ionization energies in agreement with experiment, identification of generic electrophilic sites involving dangling hydrogen atoms, and the tendency of all hydrogen atoms to be saturated in hydrogen bonding or in interaction with the excess electron. An emerging insight is the capability of electrophilic sites to be actuators of electron transport pathways in biomolecular systems.} } @article{ChoJo96flipping, author = {K. Cho and J. D. Joannopoulos}, title = {Flipping silicon dimers on {Si}(100) using scanning tip microscopy: A theoretical investigation}, journal = {Physical Review~B}, year = 1996, volume = 53, pages = {45534556}, month = {February}, url = {http://link.aps.org/abstract/PRB/v53/p4553}, abstract = {Stateoftheart firstprinciples calculations are performed to investigate the possibility of using a tungsten tip in scanning tunneling microscopy or atomic force microscopy as a tool to manipulate the surface atoms of Si(100). Calculations of total energy and electronic structure are used to study the energetics and bonding properties associated with the tip and surface for a variety of atomic configurations. The results predict that under certain protocols the tip can be used to flip dimers on the surface, from one buckled configuration to another, reversibly, and without inducing damage to either the intrinsic surface or the tip. The implications of these results for using the intrinsic (100) surface of silicon as an ultrahigh density memory storage medium approaching one bit of data per dimer are discussed.} } @article{LimCh96, author = {H. Lim and K. Cho and R. B. Capaz and J. D. Joannopoulos and K. D. Brommer and B. E. Larson}, title = {Ab initio studies of adatom vacancies on the {Si}(111)$(7\times7)$ surface}, journal = {Physical Review~B}, year = 1996, volume = 53, pages = {1542115424}, month = {June}, url = {http://link.aps.org/abstract/PRB/v53/p15421}, abstract = {Ab initio totalenergy calculations are used to investigate adatom vacancies on the Si(111)$(7\times7)$ surface. In striking contrast to recent experimental estimates, vacancy formation energies are calculated to be 0.9 eV on average, with $\sim 0.1$eV variations depending on the type of adatom. We find that faulted or corner adatoms can be removed more easily than unfaulted or edge adatoms, respectively. Structural relaxations induce large changes in the electronic structure of the surface states. Calculation of scanning tunneling microscopy (STM) images show that the predicted variations should be readily observed in differential STM measurements.} } @article{CapazJo96, author = {R. B. Capaz and J. D. Joannopoulos}, title = {Unified approach for the calculation of force constants and accelerated convergence of atomic coordinates }, journal = {Physical Review~B}, year = 1996, volume = 54, pages = {1340213405}, month = {November}, url = {http://link.aps.org/abstract/PRB/v54/p13402}, abstract = {A method which simultaneously computes the forceconstant matrix (FCM) and relaxes the atomic coordinates to the ground state of a generic system is introduced. The method is based on the optimization of the FCM as the atomic coordinates evolve to their groundstate configuration, using the variations in atomic positions and forces from previous iterations as inputs. Striking accelerated convergence as compared to traditional methods is obtained, as well as accurate vibrational frequencies with no additional cost.} } @article{AydinolKo97, author = {M. K. Aydinol and A. F. Kohan and G. Ceder and K. Cho and J. Joannopoulos}, title = {Ab initio study of lithium intercalation in metal oxides and metal dichalcogenides}, journal = {Physical Review~B}, year = 1997, volume = 56, pages = {13541365}, month = {July}, url = {http://link.aps.org/abstract/PRB/v56/p1354}, abstract = {A study of the average voltage to intercalate lithium in various metal oxides is presented. By combining the ab initio pseudopotential method with basic thermodynamics the average intercalation voltage can be predicted without the need for experimental data. This procedure is used to systematically study the effect of metal chemistry, anion chemistry, and structure. It is found that Li is fully ionized in the intercalated compounds with its charge transferred to the anion and to the metal. The substantial charge transfer to the anion is responsible for the large voltage difference between oxides, sulfides, and selenides. Ionic relaxation, as a result of Li intercalation, causes nonrigidband effects in the density of states of these materials. Suggestions for compounds that may have a substantially larger voltage than currently used materials are also presented.} } @article{ChoKa97, author = {Kyeongjae Cho and Efthimios Kaxiras and J. D. Joannopoulos}, title = {Theory of Adsorption and Desorption of H$_2$ Molecules on the {Si}(111)$(7\times7)$ surface}, journal = {Physical Review Letters}, year = 1997, volume = 79, pages = {50785081}, month = {December}, url = {http://link.aps.org/abstract/PRL/v79/p5078}, abstract = {The physics of adsorption and desorption of H$_2$ molecules on the Si(111)$(7\times7)$ surface is investigated through firstprinciples density functional theory calculations. The calculated adsorption and desorption energy barriers are 0.8 and 2.4 eV, respectively. The process of adatom backbond breaking is identified as the fundamental microscopic mechanism determining the adsorption energy barrier and the kinetic energy of the desorbed H$_2$ molecule. These results shed light on controversial experimental findings for this classic surfacemolecule system.} } @article{CapazDalPino98, author = {R. B. Capaz and Dal Pino, Jr., A. and J. D. Joannopoulos}, title = {Theory of carboncarbon pairs in silicon}, journal = {Physical Review~B}, year = 1998, volume = 58, pages = {98459850}, month = {October}, url = {http://link.aps.org/abstract/PRB/v58/p9845}, abstract = {Interstitialsubstitutional carbon pairs (C$_i$C$_s$) in silicon display interesting metastable behavior associated with two different structural configurations. In this work, we perform extensive ab initio calculations on this system. Our results show the following. (i) The metastable configuration for the neutral charge state displays $C_{1h}$ symmetry and it is reminiscent of the isolated interstitial carbon configuration, i.e., a split interstitial CSi pair with the substitutional carbon bonded to the silicon interstitial. (ii) The groundstate configuration also has $C_{1h}$ symmetry, but it consists of a single silicon interstitial twofold coordinated in an unusual bridge configuration between two substitutional carbon atoms. With an activation energy of 0.07 eV, this configuration becomes a motionalaveraged state with $C_{3v}$ symmetry. (iii) The ground state is lower in energy by 0.11 eV with respect to the metastable state. The jump from one configuration to the other corresponds to bondswitching mechanism with a calculated energy barrier of 0.13 eV. (iv) Both configurations have two electronic states in the gap, with gapstate wave functions consistent with the local bonding of the defect complex in each case. (v) Analysis of localmode vibrations on the groundstate configuration indicates a stronger component in one of the carbon atoms, which explains the experimentally observed isotope splittings. Vibrational frequencies for the metastable configuration are also predicted. All of these results are in satisfactory agreement with experiments. } } @article{CapazAs99, author = {R. B. Capaz and L. V. C. Assali and L. C. Kimerling and K. Cho and J. D. Joannopoulos}, title = {Mechanism for hydrogenenhanced oxygen diffusion in silicon}, journal = {Physical Review~B}, year = 1999, volume = 59, pages = {48984900}, month = {February}, url = {http://link.aps.org/abstract/PRB/v59/p4898}, abstract = {Oxygen diffuses in silicon with an activation energy of 2.532.56 eV. In hydrogenated samples, this activation energy is found to decrease to 1.62.0 eV. In this paper, a microscopic mechanism for hydrogenenhanced oxygen diffusion in pdoped silicon is proposed. A path for joint diffusion of O and H is obtained from an ab initio moleculardynamics simulation in which the O atom is ``kicked'' away from its equilibrium position with a given initial kinetic energy. After reaching a maximum potential energy of 1.46 eV above the ground state, the system relaxes to a metastable state on which a SiSi bond is broken and the H atom saturates one of the dangling bonds. With an extra 0.16 eV, the SiH bond is broken and the system relaxes to an equivalent groundstate configuration. Therefore, the migration pathway is an intriguing twostep mechanism. This path represents a 0.54eV reduction in the static barrier when compared with the diffusion of isolated O in Si, in excellent agreement with experiments. This mechanism elucidates the role played by the H atom in the process: it not only serves to ``open up'' a SiSi bond to be attacked by the oxygen, but it also helps in reducing the energy of an important intermediate state in the diffusion pathway. } } @article{WangMo99, author = {Tairan Wang and N. Moll and Kyeongjae Cho and J. D. Joannopoulos}, title = {Deliberately Designed Materials for Optoelectronics Applications}, journal = {Physical Review Letters}, year = 1999, volume = 82, pages = {33043307}, month = {April}, url = {http://link.aps.org/abstract/PRL/v82/p3304}, abstract = {A novel class of semiconductors is introduced, based on computational design, to solve the longstanding problem of lattice and polarity mismatch in heteroepitaxial growth of IIIV alloys on silicon substrates. Ab initio totalenergy calculations and quasiparticle GW calculations are used to investigate the physical properties of these new semiconductors. One particular configuration is designed to match lattice constant and polarity with the Si(100) surface and to possess a direct band gap of 1.59 $\mu$m, which is close to the canonical frequency used by the optoelectronics industry. These results could pave the way for eventual monolithic integration of optical materials on silicon.} } @article{SkorobogatiyJo99, author = {M. Skorobogatiy and J. D. Joannopoulos}, title = {Nonzerotemperature pathintegral method for fermions and bosons: A grand canonical approach }, journal = {Physical Review~B}, year = 1999, volume = 60, pages = {14331436}, month = {July}, url = {http://link.aps.org/abstract/PRB/v60/p1433}, abstract = {The calculation of the density matrix for fermions and bosons in the grand canonical ensemble allows an efficient way for the inclusion of fermionic and bosonic statistics at all temperatures. It is shown that in a path integral formulation the oneparticle density matrix can be expressed via an integration over a novel representation of the universal temperaturedependent functional. In this paper we discuss a representation for the universal functional in terms of Hankel functions which is convenient for computational applications. Temperature scaling for the universal functional and its derivatives is also introduced thus allowing an efficient rescaling rather then recalculation of the functional at different temperatures. We expect that our method will give rise to a numerically efficient pathintegral approach for calculation of a density matrix in the grand canonical ensemble.} } @article{MirbtMo99cation, author = {S. Mirbt and N. Moll}, title = {Cationrich (100) surface reconstructions of InP and GaP}, journal = {Physical Review~B}, year = 1999, volume = 60, pages = {1328313286}, month = {November}, url = {http://link.aps.org/abstract/PRB/v60/p13283}, abstract = {The trimer reconstruction of the (100) InP surface which has been discovered experimentally is confirmed by firstprinciple calculations. The charge density of atomic configuration, which has the lowest surface energy is in perfect agreement with experimental scanning tunneling microscopy images. We predict the same trimer reconstruction also to be observable on GaP (100) surfaces and discuss how local stress makes this reconstruction energetically unfavorable for GaAs.} } @article{AppelbaumWa04, author = {Ian Appelbaum and Tairan Wang and J. D. Joannopoulos and V. Narayanamurti}, title = {Ballistic hotelectron transport in nanoscale semiconductor heterostructures: Exact selfenergy of a threedimensional periodic tightbinding {Hamiltonian}}, journal = {Physical Review~B}, year = 2004, volume = 69, pages = 165301, month = {April}, url = {http://link.aps.org/abstract/PRB/v69/e165301}, abstract = {As the length scale for semiconductor heterostructures approaches the regime of the lattice constant, our current theory for calculating ballistic hotelectron transport becomes inapplicable. In this case, a method such as the Green's function formalism should be used to calculate ballistic electron transmission functions from the exact, periodic lattice potential. We present a method for directly calculating the exact surface Green's function for threedimensional periodic leads which is necessary for such a scheme. Except in cases of high crystal symmetry, the method is limited by the difficulty to solve a nonsymmetric matrix Riccati equation.} } @article{ReedFr03, author = {Evan J. Reed and Laurence E. Fried and J. D. Joannopoulos}, title = {A Method for Tractable Dynamical Studies of Single and Double Shock Compression }, journal = {Physical Review Letters}, year = 2003, volume = 90, pages = 235503, month = {June}, url = {http://link.aps.org/abstract/PRL/v90/e235503}, abstract = {A new multiscale simulation method is formulated for the study of shocked materials. The method combines molecular dynamics and the Euler equations for compressible flow. Treatment of the difficult problem of the spontaneous formation of multiple shock waves due to material instabilities is enabled with this approach. The method allows the molecular dynamics simulation of the system under dynamical shock conditions for orders of magnitude longer time periods than is possible using the popular nonequilibrium molecular dynamics approach. An example calculation is given for a model potential for silicon in which a computational speedup of $10^5$ is demonstrated. Results of these simulations are consistent with the recent experimental observation of an anomalously large elastic precursor on the nanosecond time scale.} } @article{AppelbaumJo02, author = {Ian Appelbaum and J. D. Joannopoulos and V. Narayanamurti}, title = {Alternative paradigm for physical computing }, journal = {Physical Review~E}, year = 2002, volume = 66, pages = 066612, month = {December}, url = {http://link.aps.org/abstract/PRE/v66/e066612}, abstract = {We identify a different class of physical systems that are able to form universal logic gates. By analogy with Si(100) surface dimers, we present a model to analyze the trajectories of the fixed points (interpreted as logic states) under variation of the basic parameters. Using the perspective of catastrophe theory, we show that information processing is the result of cycling the parameters of such systems through a path containing a cusptype catastrophe. We apply this analysis to the construction of an example based on magnetic memory.} } @article{WangMo00, author = {Tairan Wang and N. Moll and Kyeongjae Cho and J. D. Joannopoulos}, title = {Computational design of compounds for monolithic integration in optoelectronics}, journal = {Physical Review~B}, year = 2000, volume = 63, pages = 035306, month = {January}, url = {http://link.aps.org/abstract/PRB/v63/e035306}, abstract = {A class of semiconductors is introduced and their physical properties are examined using both ab initio totalenergy calculations and quasiparticle GW calculations. These compounds are designed to address problems of latticeconstant mismatch and polarity mismatch that are common issues in heteroepitaxial growth of IIIV alloys on silicon substrates. A variety of configurations of these materials is explored. It is found that their lattice constants and band gaps fall into a region of phase space different from that of conventional semiconductors, making them potential candidates for the basis of optical devicesinfrared emitters and detectors. A particular suitable configuration is identified that is latticeconstant matched to Si and has a direct band gap of 0.8 eV. This gap corresponds to the canonical wavelength of 1.5 $\mu$m in optoelectronics. Thus this material could ultimately enable tractable monolithic integration of optics with electronics. The characteristics of this particular configuration are examined in depth, including its temperature dependence, its bulk energetics, and its growth energetics. The results of these analyses indicate that fabrication of these compounds using heteroepitaxial growth techniques should be feasible.} } @article{ReedJo00, author = {Evan J. Reed and J. D. Joannopoulos and Laurence E. Fried}, title = {Electronic excitations in shocked nitromethane }, journal = {Physical Review~B}, year = 2000, volume = 62, pages = {1650016509}, month = {December}, url = {http://link.aps.org/abstract/PRB/v62/p16500}, abstract = {The nature of electronic excitations in crystalline solid nitromethane under conditions of shock loading and static compression are examined. Densityfunctional theory calculations are used to determine the crystal bandgap under hydrostatic stress, uniaxial strain, and shear strain. Bandgap lowering under uniaxial strain due to molecular defects and vacancies is considered. Ab initio moleculardynamics simulations are done of all possible nearestneighbor collisions at a shock front, and of crystal shearing along a sterically hindered slip plane. In all cases, the bandgap is not lowered enough to produce a significant population of excited states in the crystal. The nearly free rotation of the nitromethane methyl group and localized nature of the highest occupied molecular orbital and lowest unoccupied molecular orbital states play a role in this result. Dynamical effects have a more significant effect on the bandgap than static effects, but relative molecule velocities in excess of 6 km/s are required to produce a significant thermal population of excited states.} } @article{HuangWa05PRB, author = {Kerwyn Casey Huang and Tairan Wang and John D. Joannopoulos}, title = {Nanoscale properties of melting at the surface of semiconductors}, journal = {Physical Review~B}, year = 2005, volume = 72, pages = 195314, month = {November}, url = {http://link.aps.org/abstract/PRB/v72/e195314}, abstract = {In this study, we use densityfunctional simulations to investigate how nanoscale surface coatings can alter the melting dynamics of semiconductor surfaces. We demonstrate that a singlemonolayer coating of GaAs can dramatically reduce the diffusive motion of the surface atoms of a Ge(110) crystal and cause superheating of the bulk at temperatures well above the Ge melting point on the 10 ps time scale. In direct contrast, a singlemonolayer coating of Ge will induce surface melting of a GaAs(110) structure 300 K below the GaAs melting point. We also identify a metallization of the band structure and bond alterations in the charge density near the melting transition. In addition, we suggest that the Ge monolayer causes the GaAs(110) surface to melt through transient penetration of the Ge atoms into the bulk, which locally initiates the collective diffusive motion of large groups of Ga and As atoms. These studies on the effect of coatings in semiconductors clearly point to the surface material as the dominant determinant of the melting characteristics of a hybrid structure. These simulations have important implications for hightemperature materials design while simultaneously probing the fundamental features of the melting transition.} } @article{EvansZh05, author = {M. H. Evans and X.G. Zhang and J. D. Joannopoulos and S. T. Pantelides}, title = {FirstPrinciples Mobility Calculations and AtomicScale Interface Roughness in Nanoscale Structures}, journal = {Physical Review Letters}, year = 2005, volume = 95, pages = 106802, month = {September}, url = {http://link.aps.org/abstract/PRL/v95/e106802}, abstract = {Calculations of mobilities have so far been carried out using approximate methods that suppress atomicscale detail. Such approaches break down in nanoscale structures. Here we report the development of a method to calculate mobilities using atomicscale models of the structures and density functional theory at various levels of sophistication and accuracy. The method is used to calculate the effect of atomicscale roughness on electron mobilities in ultrathin doublegate silicononinsulator structures. The results elucidate the origin of the significant reduction in mobility observed in ultrathin structures at low electron densities.} } @article{EvansJo05, author = {M. H. Evans and J. D. Joannopoulos and S. T. Pantelides}, title = {Electronic and mechanical properties of planar and tubular boron structures}, journal = {Physical Review~B}, year = 2005, volume = 72, pages = 045434, month = {July}, url = {http://link.aps.org/abstract/PRB/v72/e045434}, abstract = {We report the results of firstprinciples calculations showing that boron can form a wide variety of metastable planar and tubular forms with unusual electronic and mechanical properties. The preferred planar structure is a buckled triangular lattice that breaks the threefold ground state degeneracy of the flat triangular plane. When the plane is rolled into a tube, the ground state degeneracy leads to a strong chirality dependence of the binding energy and elastic response, an unusual property that is not found in carbon nanotubes. The achiral ($n$,0) tubes derive their structure from the flat triangular plane. The achiral ($n$,$n$) boron nanotubes arise from the buckled plane, and have large cohesive energies and different structures as a result. ($n$,$n$) boron nanotubes have an internal relaxation mechanism that results in a very low Poisson ratio. The strong variation in elastic properties of boron nanotubes makes them the mechanical analogue of carbon nanotubes, and may make them ideal candidates for applications in composite materials and nanoelectromechanical systems.} } @article{HuangWa05, author = {Kerwyn Casey Huang and Tairan Wang and John D. Joannopoulos}, title = {Superheating and Induced Melting at Semiconductor Interfaces}, journal = {Physical Review Letters}, year = 2005, volume = 94, pages = 175702, month = {May}, url = {http://link.aps.org/abstract/PRL/v94/e175702}, abstract = {We present ab initio densityfunctional simulations of the state of several semiconductor surfaces at temperatures near the bulk melting temperatures. We find that the solidliquid phasetransition temperature at the surface can be altered via a microscopic (singlemonolayer) coating with a different latticematched semiconducting material. Our results show that a singlemonolayer GaAs coating on a Ge(110) surface above the Ge melting temperature can dramatically reduce the diffusion coefficient of the germanium atoms, going so far as to prevent melting of the bulk layers on the 10 ps time scale. In contrast, a singlemonolayer coating of Ge on a GaAs(110) surface introduces defects into the bulk and induces melting of the top layer of GaAs atoms 300 K below the GaAs melting point. To our knowledge, these calculations represent the first ab initio investigation of the superheating and induced melting phenomena.} } @article{WangMo99JVSTB, author = {Tairan Wang and N. Moll and Kyeongjae Cho and J. D. Joannopoulos}, title = {Deliberately designed interfaces for monolithic integration in optoelectronics}, journal = {Journal of Vacuum Science and Technology~B}, year = 1999, volume = 17, number = 4, pages = {16121616}, month = {March}, doi = {doi:10.1116/1.590799}, abstract = {A novel class of semiconductors is introduced, based on computational design, to solve the longstanding problem of lattice and polarity mismatch in heteroepitaxial growth of IIIV alloys on silicon substrates. Ab initio totalenergy calculations and quasiparticle GW calculations are used to investigate the physical properties of these new semiconductors. One particular configuration is designed to match lattice constant and polarity with the Si(100) surface and to possess a direct band gap of 1.59 $\mu$m, which is close to the canonical frequency used by the optoelectronics industry. These results could pave the way for eventual monolithic integration of optical materials on silicon.} } @article{ChoJo96reversible, author = {K. Cho and J. D. Joannopoulos}, title = {Reversible Tipinduced Structural Modifications in Scanning Tip Microscopy}, journal = {Japanese Journal of Applied Physics}, year = 1996, volume = 35, pages = {37143718} } @article{TepeschKo96, author = {Patrick D. Tepesch and Adrian F. Kohan and Gerardo D. Garbulsky and Gerbrand Ceder and Crystal Coley and Harold T. Stokes and Larry L. Boyer and Michael J. Mehl and Benjamin P. Burton and Kyeongjae Cho and John Joannopoulos}, title = {A Model to Compute Phase Diagrams in Oxides with Empirical or FirstPrinciples Energy Methods and Application to the Solubility Limits in the {CaO}{MgO} System}, journal = {Journal of the American Ceramic Society}, year = 1996, volume = 79, number = 8, pages = {20332040}, month = {August}, doi = {doi:10.1111/j.11512916.1996.tb08934.x}, abstract = {The CaOMgO system is used as a prototype system to evaluate the accuracy of several energy and entropy approximations for predicting solidstate phase diagrams in ionic materials. Configurational disorder between the cations is parameterized with the cluster expansion technique. The vibrational contribution to the free energy is incorporated with a harmonic model that accounts for the dependence of the vibrational density of states on the cation configuration. The CaOMgO hase diagram can be predicted very accurately with quantum mechanical energy methods, without the use of any adjustable parameters. Published empirical potential parameters for the CaOMgO system reproduce the qualitative features of the phase diagram but significantly underestimate the solubility limits. Parameters that reasonably reproduce the quantum mechanical results are presented.} } This file was generated by bibtex2html 1.90. 