Recent publications:

Zhao, Y. J., S. Picozzi, A. Continenza, W. T. Geng and A. J. Freeman (2002). "Possible impurity-induced ferromagnetism in II-Ge-V-2 chalcopyrite semiconductors." Physical Review B 65, 094415.

            Recently reported room-temperature ferromagnetic (FM) semiconductors Cd1-xMnxGeP2 and Zn1-xMnxGeP2 point to a possible important role of II-IV-V-2 chalcopyrite semiconductors in "spintronic" studies. Here, structural, electronic, and magnetic properties of (i) Mn-doped II-Ge-VI2 (II = Zn or Cd and V = P or As) chalcopyrites and (ii) the role of S as impurity in Cd1-xMnxGeP2 are studied by first- principles density functional calculations. We find that the total energy of the anti ferro magnetic (AFM) state is lower than the corresponding FM state for all systems with Mn composition x = 0.25, 0.50, and 1.0. This prediction is in agreement with a recent experimental finding that Zn1-xMnxGeP2 experiences a FM to AFM transition for T less than 47 K. Furthermore, a possible transition to the half-metallic FM phase is predicted in Cd1-xMnxGeP2 due to the electrons introduced by n-type S doping, which indicates the importance of carriers for FM coupling in magnetic semiconductors. As expected, the total magnetic moment for the FM phase is reduced by one mu(B) with each S substituting P.

Zhao, Y. J., W. T. Geng, A. J. Freeman and B. Delley (2002). "Structural, electronic, and magnetic properties of alpha- and beta-MnAs: LDA and GGA investigations." Physical Review B 65, 113202.

            Zinc-blende (alpha-) and NiAs-type (beta-) MnAs are investigated with a combined first-principles linearized argumented plane wave and DMol(3) study within both the local density approximation (LDA) and the generalized gradient approximation (GGA). First-principles calculations within the GGA predict the lattice volume for b- MnAs much better than LDA (which underestimates it by 15%) compared with experiment. The LDA calculated equilibrium lattice volume of alpha-MnAs is 10% smaller than that of GaAs, which is in contradiction to the well-accepted fact that the lattice volume of Ga1-xMnxAs increases with x. In contrast, the GGA predicts a reasonable lattice volume for alpha-MnAs. The ferromagnetic alpha-MnAs is shown to be a metal at a=5.7 Angstrom, and to undergo a transition to a half-metallic phase when it expands to a>5.8 Angstrom due to the decreased bandwidth. Further, the calculated cohesive energy of beta-MnAs is nearly 0.87 eV greater than that of alpha-MnAs, which provides theoretical support for the instability of alpha-MnAs.

 

Zhao, Y. J. and A. J. Freeman (2002). "First-principles prediction of a new class of ferromagnetic semiconductors." Journal of Magnetism and Magnetic Materials 246, 145.

            Ferromagnetism in I-III-VI2 chalcopyrite semiconductors is predicted to arise from holes provided by Mn doping in the illustrative case of CuGaSe2. The first-principles calculations within the generalized gradient approximation to density functional theory demonstrate that Mn substitutes for Ga sites in CuGaSe2 with a formation energy that is similar to0.25 eV lower than that for Cu sites. Ferromagnetic CuGa1-xMnxSe2 is found to be a half-metallic material with a magnetic moment of 4 mu(B) per Mn for xgreater than or equal to0.25, and its estimated Curie temperature is more than 110 K. We suggest that higher Curie temperatures may be achieved for this new class of ferromagnetic semiconductors based on Mn-doped I-III-VI2 chalcopyrites by employing other materials with a smaller lattice constant such as CuGaSe2 and CuAlSe2. (C) 2002 Elsevier Science B.V. All rights reserved.

Youn, S. J., B. I. Min and A. J. Freeman (2002). "Large anisotropy in the optical conductivity of YNi2B2C." Physical Review B 66, 052512.

            The optical properties of YNi2B2C are studied by using the first-principles full-potential linearized augmented plane-wave method within the local density approximation. Anisotropic behavior is obtained in the optical conductivity, even though the electronic structure shows three-dimensional character. A large peak in sigma(z) is obtained at 2.4 eV. The anisotropic optical properties are analyzed in terms of interband transitions between energy levels, and it is found that the Ni site plays an important role. The electronic energy loss spectroscopy spectra are also calculated to help elucidate the anisotropic properties in this system.

 

Stampfl, C., R. Asahi and A. J. Freeman (2002). "Surface properties of the refractory metal-nitride semiconductor ScN: Screened-exchange LDA-FLAPW investigations." Physical Review B 65, 161204.

            Density-functional theory calculations employing the screened- exchange local-density approximation (SX-LDA) with the full potential linearized augmented plane-wave method have recently shown that the relatively unexplored refractory nitrides ScN, YN, and LaN are semiconductors. For the ScN(001) surface, the present calculations predict that the ideal-relaxed surface has the lowest formation energy for most of the range of the allowed chemical potentials-and is semiconducting-while N- deficient structures, which are predicted to form for Sc-rich conditions, are metallic in nature. Compared to the LDA surface-state band structures, the SX-LDA selectively pushes down the valence bands for the Sc-terminated surface and pushes up the conduction bands for the N-terminated structure.

 Picozzi, S., R. Asahi, C. B. Geller, A. Continenza and A. J. Freeman (2002). "Impact ionization in GaAs: A screened exchange density- functional approach." Physical Review B 65, 113206.

            Results are presented of a fully ab initio calculation of impact ionization rates in GaAs within the density functional theory framework, using a screened-exchange formalism and the highly precise all-electron full-potential linearized augmented plane wave method. The calculated impact ionization rates show a marked orientation dependence in k space, indicating the strong restrictions imposed by the conservation of energy and momentum. This anisotropy diminishes as the impacting electron energy increases. A Keldysh type fit performed on the energy- dependent rate shows a rather soft edge and a threshold energy greater than the direct band gap. The consistency with available Monte Carlo and empirical pseudopotential calculations shows the reliability of our approach and paves the way to ab initio calculations of pair production rates in new and more complex materials.

 

Picozzi, S., A. Continenza, S. Massidda and A. J. Freeman (2002). "Schottky barrier height in GaN/Al junctions: an ab-initio study." Physica Status Solidi a-Applied Research 190, 257.

            We performed ab-initio full-potential linearized augmented plane wave (FLAPW) calculations for [0001]-wurtzite and [111]- zincblende GaN/Al junctions, focusing on the Schottky barrier height. We propose a procedure to evaluate the potential discontinuity in the presence of electric fields, showing that their effect is relatively small (a few tenths of an eV). These calculations assess the rectifying behaviour of the GaN/Al contact, in agreement with experimental values for the barrier.

 

Picozzi, S., G. Profeta, A. Continenza, S. Massidda and A. J. Freeman (2002). "Role of structural relaxations and chemical substitutions on piezoelectric fields and potential lineup in GaN/Al junctions." Physical Review B 65, 165316.

            First-principles full-potential linearized augmented plane wave calculations are performed to clarify the role of the interface geometry on piezoelectric fields and potential lineups in [0001] wurtzite and [111]-zincblende GaN/Al junctions. The electric field (polarity and magnitude) is found to be strongly affected by atomic relaxations in the interface region. A procedure is used to evaluate the Schottky-barrier height in the presence of electric fields, showing that their effect is relatively small (a few tenths of an eV). These calculations assess the rectifying behavior of the GaN/Al contact, in agreement with experimental values for the barrier. We disentangle chemical and structural effects on the relevant properties (such as the potential discontinuity and the electric field) by studying unrelaxed ideal nitride/metal systems. Using simple electronegativity arguments, we outline the leading mechanisms that define the values of the electric field and Schottky barrier in these ideal systems. Finally, the transitivity rule is proved to be well satisfied.

 

Picozzi, S., A. Continenza and A. J. Freeman (2002). "Co2MnX (X=Si, Ge, Sn) Heusler compounds: An ab initio study of their structural, electronic, and magnetic properties at zero and elevated pressure." Physical Review B 66, 094421.

            The structural, electronic, and magnetic properties of Co2MnX (X=Si, Ge, Sn) Heusler compounds have been determined by means of all-electron full-potential linearized augmented plane wave (FLAPW) calculations. We focus on the effects on the electronic and magnetic properties induced by: (i) substitution of the X atom, (ii) applied pressure, and (iii) the use of the local spin density approximation (LSDA) vs the generalized gradient approximation (GGA) in density functional theory. A comparison between LSDA and GGA for the exchange-correlation functional shows that GGA is essential for an accurate description of the equilibrium volumes and of the electronic and magnetic properties of these systems. We find that both the energy gap and the spin gap increase as the X atomic number decreases. As a result of the semiconducting (metallic) character found in the minority (majority) spin band structure, the Si and Ge based alloys are predicted to be half-metallic. In contrast, Co2MnSn is found to be a "nearly half-metallic" compound, since the minority valence band maximum crosses the Fermi level. The calculated total magnetization of 5 mu(B) is in excellent agreement with recent experiments. By including a fully self- consistent treatment of spin-orbit coupling, the GGA calculated orbital moments are shown to be very small (about 0.008 mu(B) for Mn and about 0.02 mu(B) for Co), showing that the quenching of the orbital magnetic moment is nearly complete. The calculated hyperfine fields, both at zero and elevated pressure, are compared with available experimental data, and show general agreement, except for Mn. Finally, the calculated Mn 2p exchange splittings, found to be in good agreement with experiment, are proportional to the Mn magnetic moments, suggesting a localized nature of ferromagnetism in these Heusler compounds.

 

Picozzi, S., A. Continenza and A. J. Freeman (2002). "The Co2MnGe Heusler compound: A first principles study of the bulk phase and of the interface with GaAs." Ieee Transactions on Magnetics 38, 2895.

            Within density functional theory, full-potential linearized augmented plane wave calculations have been performed to study the structural, electronic and magnetic properties of the Co2MnGe Heusler compound. As a result of the semiconducting (metallic) character found in the minority (majority) band structure, this material is predicted to be half-metallic, with an integer magnetic moment of 5 mu(B). This peculiar property, highly desired for spintronics applications, stimulated the study of the CO2MnGe/GaAs interface. Our calculations show that the predicted band offset is not suitable for spin-injection and that, due to interface states, the half-metallicity is lost in proximity to the junction.

 

Picozzi, S., R. Asahi, C. B. Geller and A. J. Freeman (2002). "Accurate first-principles detailed-balance determination of Auger recombination and impact ionization rates in semiconductors." Physical Review Letters 89, 197601.

            The technologically important prediction of Auger recombination lifetimes in semiconductors is addressed by means of a fully first-principles formalism, based on precise energy bands and wave functions provided by the full-potential linearized augmented plane wave code. The minority carrier Auger lifetime is determined by two related approaches: (i) a direct evaluation within Fermi's golden rule, and (ii) an indirect evaluation, based on a detailed balance formulation combining Auger recombination and its inverse process, impact ionization, in a unified framework. Lifetimes determined with the direct and indirect methods show excellent consistency between them (i) for n-doped GaAs and (ii) with measured values for GaAs and InGaAs. This indicates the computational formalism as a new sensitive tool for use in materials performance optimization.

 

Picozzi, S., A. Continenza, Y. J. Zhao, W. T. Geng and A. J. Freeman (2002). "Structural, electronic and magnetic properties of chalcopyrite magnetic semiconductors: A first-principles study." Journal of Vacuum Science & Technology a-Vacuum Surfaces and Films 20, 2023.

            Stimulated by recent experimental observations of room- temperature ferromagnetism of MnxCd1-xGeP2 and MnxZn1-xGeP2, we investigate the structural, electronic, and magnetic properties of this class of systems (II-Ge-V-2, II=Zn, Cd, and V=As, P) as a function of Mn concentration and chemical constituents by means of first-principles density-functional-theory-based codes. Our calculations indicate that, for Mn substituting the II element, the anti ferromagnetic alignment is the most stable ordering for all the systems studied. For Zn- and Cd-rich systems, the total magnetic moments per Mn atom of the ferromagnetic phase is very close to the ideal value of 5 muB, since the Mn 3d states in the minority spin channel are nearly empty; on the other hand, for Mn rich compounds, the stronger p-d hybridization lowers the total magnetic moment to about 4.4 muB. (C) 2002 American Vacuum Society.

 

Nakamura, K., A. J. Freeman, D. S. Wang, L. P. Zhong and J. Fernandez-de-Castro (2002). "Magnetic structures at the ferromagnetic NiFe and antiferromagnetic NiMn interface in exchange-biased films: Role of noncollinear magnetism and roughness." Physical Review B 65, 012402.

            The magnetic structures at the compensated ferromagnetic/antiferromagnetic (FM/AFM) interface of exchange bias FM NiFe/AFM NiMn films were determined with first- principles full-potential linearized augmented plane-wave calculations including noncollinear magnetism. The results predict that the magnetic moments of the FM NiFe layer lie perpendicular to those of the AFM NiMn layer. The intra-atomic noncollinear magnetism that arises near the interface is found to play an important role in stabilizing the perpendicular coupling that leads to a large biquadratic exchange energy (BEE). The BEE is large enough to require formation of a magnetic domain wall (with an estimated thickness similar to370 Angstrom) in the AFM NiMn layers, which may account for the observed large coercivity and exchange bias. We also discuss magnetic structures at a rough FM/AFM interface-as simulated in model calculations with the inclusion of line step defects- which may contribute to a unidirectional magnetic anisotropy.

 

Nakamura, K. and A. J. Freeman (2002). "Canted ferromagnetism in RuSr2GdCu2O8." Physical Review B 66, 140405.

            First principles calculations using the full-potential linearized augmented plane wave method including intra-atomic noncollinear magnetism have been performed to determine the magnetic structures of RuSr2GdCu2O8. The magnetism clearly arises from the RuO6 octahedra where the moments on neighboring Ru sites order antiferromagnetically but cant perpendicular to the antiferromagnetic (AFM) axis-and so induce a weak ferromagnetism. The projected Ru moments along the AFM and FM axes result in magnetic moments of 1.16 and 0.99mu(B), respectively. The results are consistent with the possible coexistence of canted ferromagnetism and superconductivity in the RuSr2GdCu2O8-inferred from experiments.

 

Kim, I. G., J. I. Lee and A. J. Freeman (2002). "Effects of d-band hole doping on electronic structure, magnetism, and superconductivity in Fe- and Co-doped MgCNi3." Physical Review B 65, 064525.

            Electronic structures and magnetism of the nonoxide perovskite MgCNi3-xTx compounds (x=0, 1, 2, and 3; T=Co and Fe), based on recently found superconducting MgCNi3, have been investigated with first-principles all-electron full-potential linearized augmented plane-wave calculations within the local-(spin) density approximation. From the results of total-energy calculations, it is expected that (i) the suppression of superconductivity occurs faster for the Fe-doped case than for the Co-doped case with increase of x, and (ii) ferromagnetic transitions occur when xgreater than or equal to2 for the Co- doped cases, while the Fe-doped cases become ferromagnetic before x=2. From the calculated density of states, it is found that the superconductor MgCNi3 becomes paramagnetic and then ferromagnetic, as we increase the number of minority spin d- band holes via Co and Fe doping at Ni sites. A hypothetical system, MgC(FeCoNi), has also been calculated and found to be ferromagnetic with magnetic moments of 0.97, 0.24, and 0.03 mu(B) for Fe, Co, and Ni, respectively, which are roughly proportional to the number of d-band holes of minority spin. These features are well explained by d-band spin splitting, as in conventional metallic ferromagnets.

 Freeman, A. J. (2002). "Materials by design and the exciting role of quantum computation/simulation." Journal of Computational and Applied Mathematics 149, 27.

            It is now well-recognized that we are witnessing a golden age of innovation with novel materials, with discoveries important for both basic science and device applications-some of which will be treated at this Workshop. In this talk, we discuss the role of computation and simulation in the dramatic advances of the past and those we are witnessing today. We will also describe the growing acceptance and impact of computational materials science as a major component of materials research and its import for the future. In the process, we will demonstrate how the well-recognized goal driving computational physics/computational materials science-simulations of ever- increasing complexity on more and more realistic models-has been brought into greater focus with the introduction of greater computing power that is readily available to run sophisticated and powerful software codes like our highly precise full-potential linearized augmented plane wave (FLAPW) method, now also running on massively parallel computer platforms. We will then describe some specific advances we are witnessing today, and computation and simulation as a major component of quantum materials design and its import for the future, with the goal-to synthesize materials with desired properties in a controlled way via materials engineering on the atomic scale. The theory continues to develop along with computing power. With the universality and applicability of these methods to essentially all materials and properties, these simulations are starting to fill the increasingly urgent demands of material scientists and engineers. (C) 2002 Elsevier Science B.V. All rights reserved.

Asahi, R., A. Wang, J. R. Babcock, N. L. Edleman, A. W. Metz, M. A. Lane, V. P. Dravid, C. R. Kannewurf, A. J. Freeman and T. J. Marks (2002). "First-principles calculations for understanding high conductivity and optical transparency in InxCd1-x films." Thin Solid Films 411, 101.

            We investigate InxCd1-xO materials, where x = 0.0, 0.031, 0.063 and 0.125, to understand their high electrical conductivity and optical transparency windows, using the full-potential linearized augmented plane wave (FLAPW) method. In addition, we employ the screened exchange LDA (sX-LDA) method to evaluate accurate band structures including band gap that is underestimated by the LDA calculations. The results show a dramatic Burstein-Moss shift of the absorption edge by the In doping, reflecting the small effective mass of the Cd 5s conduction band. The calculated direct band gaps, 2.36 eV for x=0.0 and 3.17 eV for x= 0.063, show excellent agreement with experiment. The effective mass of the conduction band of CdO is calculated to be 0.24 in, (in the Delta direction), in good agreement with an experimental value of 0.27m(e), explaining its high electrical conductivity. The hybridization between the Cd 5s and the In 5s states yields complex many-body effects in the conduction bands: a hybridization gap in the conduction bands and a band-gap narrowing which cancels the further Burstein-Moss shift for higher In doping. (C) 2002 Elsevier Science B.V. All rights reserved.

 

Zhao, Y. J., W. T. Geng, A. J. Freeman and T. Oguchi (2001). "Magnetism of chalcopyrite semiconductors: Cd1-xMnxGeP2." Physical Review B 63, 201202.

            The recently reported room-temperature ferromagnetism in Cd1- xMnxGeP2 was investigated for x = 1.0, 0.5, and 0.25 by the local density first-principles full-potential linearized augmented plane wave (FLAPW) and DMOL3 methods within both local-density approximation (LDA) and generalized gradient approximation (GGA). We find that the total energy of the antiferromagnetic (AFM) state is lower than the corresponding ferromagnetic (FM) state fur all x studied. The GGA gives a better description of magnetic properties than LDA mainly due to its better prediction of structure, particularly for high Mn concentrations. The total spin moment of Cd1-xMnxGeP2 is similar to5.0 mu (B) per Mn atom. The FM alignment between Mn and P increases the total energy of the Mn-Mn FM coupling and makes the AFM ordering preferable.

 

Zhao, Y. J., W. T. Geng, K. T. Park and A. J. Freeman (2001). "Electronic and magnetic properties of Ga1-xMnxAs: Role of Mn defect bands." Physical Review B 64, 035207.

            Ga1-xMnxAs and related semiconductors are under intense investigation for the purpose of understanding the ferromagnetism in these materials, pursuing higher T-c, and, finally, for realizing semiconductor electronic devices that use both charge and spin. In this work, the electronic and magnetic structures of Ga1-xMnxAs (x=3.125%, 6.25%, 12.5%, 25.0%, 50.0%) are studied by first-principles full-potential linearized augmented plane wave calculations with the generalized-gradient approximation. The ferromagnetic state is lower in energy than the paramagnetic and antiferromagnetic states. It is confirmed that Mn atoms stay magnetic with well localized magnetic moments. The calculated band structure shows that Mn doping also forms defect bands, and makes (Ga,Mn)As p- type conducting by providing holes. Furthermore, an s-d population inversion is found in the Mn electronic configuration, which results from the strong Mn p-d mixing. The induced As moments are substantial (about -0.15 mu (B) per Mn atom, and almost independent of x)-in accord with a recent observed negative As magnetic circular dichroism signal.

 

Youn, S. J. and A. J. Freeman (2001). "First-principles electronic structure and its relation to thermoelectric properties of Bi2Te3." Physical Review B 63, 085112.

            The electronic structure of Bi2Te3. which is a major constituent of the best thermoelectric material operating at room temperature, was calculated by using the first principles full-potential linearized-augmented-planewave method with spin- orbit interaction included by a second variational method. A search of the whole Brillouin zone shows that the band edges are located off the symmetry lines, with locations that are in accord with the phenomenological six-band model. In doped Bi2Te3. Fermi surfaces near the band edges show a nonparabolic behavior. At a high doping concentration, the Fermi surfaces display elongated features, i.e., a knifelike Fermi surface for the valence band and spoonlike Fermi surfaces for the conduction band. which can be attributed to the layered structure of Bi2Te3. The effect of the anisotropic electronic structure combined with a low lattice thermal conductivity of Bi2Te3 gives a large figure of merit.

 

Youn, S. J., W. Mannstadt and A. J. Freeman (2001). "Analytic spin-orbit coupling matrix element formulae in FLAPW calculations." Journal of Computational Physics 172, 387.

            The full-potential, linearized augmented, plane wave (FLAPW) method is used widely for accurate electronic structure calculations. For the electronic structure of solids with heavy elements, it is necessary to include spin-orbit coupling interactions. We present simple analytical formulae for calculating FLAPW spin-orbit matrix elements. These wills erve to simplify the calculations and save computational time. (C) Academic Press.

 

Wang, A., J. R. Babcock, N. L. Edleman, A. W. Metz, M. A. Lane, R. Asahi, V. P. Dravid, C. R. Kannewurf, A. J. Freeman and T. J. Marks (2001). "Indium-cadmium-oxide films having exceptional electrical conductivity and optical transparency: Clues for optimizing transparent conductors." Proceedings of the National Academy of Sciences of the United States of America 98, 7113.

            Materials with high electrical conductivity and optical transparency are needed for future flat panel display, solar energy, and other opto-electronic technologies. InxCd1-xO films having a simple cubic microstructure have been grown on amorphous glass substrates by a straightforward chemical vapor deposition process. The x = 0.05 film conductivity of 17.000 S/cm, carrier mobility of 70 cm(2)/Vs, and visible region optical transparency window considerably exceed the corresponding parameters for commercial indium-tin oxide. Ab initio electronic structure calculations reveal small conduction electron effective masses, a dramatic shift of the CdO band gap with doping, and a conduction band hybridization gap caused by extensive Cd 5s + In 5s mixing.

 

Stampfl, C., W. Mannstadt, R. Asahi and A. J. Freeman (2001). "Electronic structure and physical properties of early transition metal mononitrides: Density-functional theory LDA, GGA, and screened-exchange LDA FLAPW calculations." Physical Review B 63, 155106.

            The desirable physical properties of hardness, high temperature stability, and conductivity make the early transition metal nitrides important materials for various technological applications. To learn more about the nature of these materials, first-principles density-functional theory calculations using the full-potential linearized augmented plane wave (FLAPW) method within the local-density approximation (LDA) and with the generalized gradient approximation (GGA) have been performed. We investigate the bulk electronic and physical properties of a series of early transition metal mononitrides, namely, those formed with 3d metals (ScN, TiN, VN), 4d metals (YN, ZrN, NbN), and 5d metals (LaN, HfN, TaN) in the rocksalt structure. In particular, lattice constants, bulk moduli, heats of formation, and cohesive energies as well as bulk band structures and densities of states are reported, and trends discussed. We find that the GGA yields 1%-2% larger lattice constants, 10%-20% smaller bulk moduli, and 10%-30% lower heats of formation compared to the LDA. The GGA slightly overestimates lattice constants, but leads overall to an improved agreement with experiment compared to the LDA, for which the values rue too small. These materials are metallic with the exception of ScN, YN, and LaN, which appear to be semimetals within the LDA (and GCA), but are in fact semiconductors with indirect band gaps of 1.58, 0.85, and 0.75 eV, respectively, as revealed by calculations performed using the screened-exchange LDA approach. These last, relatively unexplored. refractory III-V nitrides may therefore have potential use in device applications; in particular, ScN is well lattice matched to GaN, a wide-band-gap semiconductor that is of great current interest in relation to optoelectronic devices, and high temperature and high power electronic applications.

 

Shishidou, T., A. J. Freeman and R. Asahi (2001). "Effect of GGA on the half-metallicity of the itinerant ferromagnet CoS2." Physical Review B 64, 180401.

            The half-metallicity of CoS2, important for its possible use in spintronic applications. is investigated by means of density functional full-potential linearized augmented plane wave calculations within both the local spin-density approximation (LSDA) and the generalized gradient approximation (GGA). In contrast with earlier and our own predictions employing the LSDA, we find that CoS2 is clearly half-metallic with use of the GGA for the exchange-correlation functional. While the GGA induced modifications of the band structure are significant, they are limited to only the e(g down arrow) state which is responsible for producing the half-metallicity, while leaving all other features the same as those obtained by LSDA. We address the recently measured reflectivity spectra and rectify the assignments given for their underlying optical transitions.

 

Profeta, G., A. Continenza and A. J. Freeman (2001). "Energetics and bonding properties of the Ni/beta-SiC (001) interface: An ab initio study." Physical Review B 64, 045303.

            We investigate the adsorption of a Ni monolayer on the beta - SiC (001) surface by means of highly precise first-principles all-electron full-potential linearized augmented plane wave calculations. Total-energy calculations for the Si- and C- terminated surfaces reveal high Ni-SiC adsorption energies, with respect to other metals, confirmining the strong reactivity and the stability of the transition metal/SiC interface. These high binding energies, about 7.3-7.4 eV, are shown to be related to strong p-d hybridization, common to both surface terminations and different adsorption sites and despite the large mismatch, can stabilize overlayer growth. A detailed analysis of the bonding mechanism, hybridization of the surface states, charge transfer, and surface core level shifts reveals the strong covalent character of the bonding. After a proper accounting of the Madelung term. the core-level shift is shown to follow the charge-transfer trend.

 

Profeta, G., A. Blasetti, S. Picozzi, A. Continenza and A. J. Freeman (2001). "Termination effects at metal/ceramic junctions: Schottky barrier heights and interface properties of the beta- SiC(001)/Ni systems." Physical Review B 64, 235312.

            First-principles full-potential linearized augmented plane wave calculations for the beta -SiC[001]/Ni interface are presented, focused on the effects of different termination,, on the structural and electronic properties We find a strong reactivity of the interface, as confirmed by the high adhesion energies that are larger for the C-terminated junction than for the Si-terminated junction, in agreement with that previously found fur Ti and Al. The metal-induced gap state, are efficiently screened (Angstrom)in b(o)th termination,,, resulting in a decay length of about 1 . The calculated dependence of the Schottky barrie(r) height on different termination,, is not very strong and we investigate the observed differences between Si- and C-terminated junctions in terms of Born effective charges, electronegativity arguments, interface geometries, and screening effects. The agreement with available experimental data is excellent. thus confirming the strong rectifying behavior of this metal/ceramic contact.

 

Nakamura, K., K. T. Park, A. J. Freeman and J. D. Jorgensen (2001). "Magnetic and electronic structures of superconducting RuSr2GdCu2O8." Physical Review B 63, 024507.

            The coexistence of ferromagnetism and superconductivity in RuSr2GdCu2O8 was reported both from experiments (by Tallon et al.) and first-principles calculations (by Pickett et al.). Here we report that our first principles full-potential linearized augmented plane wave calculations, employing the precise crystal structure with structural distortions (i.e., RuO6 rotations) determined by neutron diffraction, demonstrate that antiferromagnetic ordering of the Ru moments is energetically favored over the previously proposed ferromagnetic ordering. Our results are consistent with recently performed magnetic neutron diffraction experiments (Lynn et al.). Ru t(2g) states, which are responsible for the magnetism, have only a very small interaction with Cu e (g) states, which results in a small exchange splitting of these states. The Fermi surface, characterized by strongly hybridized dp sigma orbitals, has nesting features similar to those in the two-dimensional high-T-c cuprate super conductors.

 

Kim, I. G., J. I. Lee, B. I. Min and A. J. Freeman (2001). "Surface electronic structures of superconducting thin film MgB2(0001)." Physical Review B 64, 020508.

            The surface electronic structures of superconducting MgB2 were investigated using the all-electron full-potential linearized augmented plane-wave method. Hexagonal (0001) surfaces with both B terminated (B-Term) and Mg terminated (Mg-Term) were considered. Due to the nearly-free-electron nature of the Mg surface layer, the vacuum screening range of Mg-Term is shorter than that of B-Term, which shows the covalent bonding nature of the B surface layer. Considerably enhanced densities of states near the Fermi level are found at the surface layers especially for B-Term, which is expected to yield an enhanced superconductivity in the surface of thin film MgB2 over that in bulk-assuming no large changes in the electron-ion matrix elements and phonon frequency contributions. While this expectation is contrary to the weakened superconductivity observed in surface-oriented experiments, we attribute this discrepancy to extrinsic surface effects.

 

Ingram, B. J., T. O. Mason, R. Asahi, K. T. Park and A. J. Freeman (2001). "Electronic structure and small polaron hole transport of copper aluminate." Physical Review B 64, 155114.

            High-temperature electrical property measurements (electrical conductivity, thermoelectric coefficient) on polycrystalline CuAlO2 exhibited characteristic small polaron features, i.e., low mobilities (0.1-0.4 cm(2)/V s) and an activation energy of similar to0.14 eV. The thermopower was p type (similar to 440 muV/K) and roughly temperature independent. The local-density full-potential linearized augmented-plane-wave method was used to calculate the band structure, densities of states, and optical properties (within the electric-dipole approximation) in order to account for the unique electronic and optical properties of this potential transparent conducting oxide.

 

Geng, W. T., A. J. Freeman and R. Q. Wu (2001). "Magnetism at high-index transition-metal surfaces and the effect of metalloid impurities: Ni(210)." Physical Review B 63, 064427.

            Structural, electronic, and magnetic properties near the Ni(210) surface and the effect of Li, B, P, and Ca impurities are determined by means of the all-electron total energy/atomic force full-potential linearized augmented plane wave method with the generalized gradient approximation. For the Ni(210) clean surface, simulated with an 11-layer slab, multilayer relaxation is found to be confined to the top three layers. The magnetic moment of the surface Ni layer, 0.79/mu (B), is enhanced by 27%, compared with its bulk value. This result confirms the well-accepted understanding that the reduced coordination number at a clean transition-metal surface with a high index leads to enhanced magnetic moments. Boron and P strongly alter the atomic structure at the Ni(210) surface, whereas Li and Ca have only a slight influence, due to the weak, chemical bonding with the nickel substrate. it is found that all four selected elements exert detrimental effects on the Ni(210) surface magnetism. The effects of B and P are stronger than those of Li and Ca, mainly due to their stronger hybridization with the nickel d stares. An analysis of the results for B and P suggests that it is the stronger magnetization of its free standing monolayer that makes P more detrimental than B on the nickel surface magnetism.

 

Geng, W. T., A. J. Freeman and G. B. Olson (2001). "Influence of alloying additions on grain boundary cohesion of transition metals: First-principles determination and its phenomenological extension." Physical Review B 63, 165415.

            The toughness and ductility of ultrahigh-strength alloys is often limited by intergranular embrittlement, particularly under conditions of unfavorable environmental interactions such as hydrogen embrittlement and stress corrosion cracking. Here we investigated the mechanism by which the segregated substitutional additions cause intergranular embrittlement. An electronic level phenomenological theory is proposed to predict unambiguously the effect of a substitutional alloying addition on grain boundary cohesion of metallic alloys, based on first- principles full-potential linearized augmented plane-wave method (FLAPW) calculations on the strengthening and embrittling effects of the metals Mo and Pd on the Fe grain boundary cohesion. With the bulk properties of substitutional alloying addition A and the matrix element M as inputs, the strengthening or embrittling effect of A at the grain boundary of M can be predicted without carrying out first-principles calculations once the atomic structure of the corresponding clean grain boundary is determined. Predictions of the embrittlement potency of a large number of metals, including the 3d, 4d, and 5d transition metals, are presented for the Fe Sigma3 (111) and the Ni Sigma5 (210) grain boundaries. Rigorous FLAPW calculations on the effect of Co, Ru, W, and Re on the Fe Sigma3 (111) grain boundary and Ca on the Ni Sigma5 (210) grain boundary cohesion confirm the predictions of our model. This model is expected to be applicable to other high-angle boundaries in general and instructive in the quantum design of ultrahigh-strength alloys with resistance to intergranular fracture.

 

Geng, W. T., M. Kim and A. J. Freeman (2001). "Multilayer relaxation and magnetism of a high-index transition metal surface: Fe(310)." Physical Review B 63, 245401.

            Structural, electronic, and magnetic properties of the Fe(310) surface are studied using first-principles full-potential linearized augmented plane wave method within the generalized gradient approximation. Sizable multilayer relaxation is found to extend to the seventh layer from the surface. While low- energy electron diffraction (LEED) and first-principles calculations on multilayer relaxations generally agree for low- index surfaces, there is a disagreement for this high-index surface. We conjecture that this disagreement might be due to the small data set and variational freedom in the LEED analysis. The spin magnetic moment of the Fe(310) surface and subsurface atoms is enhanced to 2.85 mu (B) and 2.65 mu (B), from a bulk value of 2.23 mu (B). The surface layer enhancement is smaller than that in Fe(100) and larger than that in Fe(111), although-all three surfaces have the same coordination number. Subsurface layers are found to play an important role in the magnetization of the surface atoms in the case of an open surface, where the vacuum affects more atomic layers.

 

Geng, W. T., A. J. Freeman and G. B. Olson (2001). "Influence of alloying additions on the impurity induced grain boundary embrittlement." Solid State Communications 119, 585.

            We have investigated the influence of Mo as a substitutional alloying addition on the grain boundary embrittlement induced by interstitial impurity P by using the first-principles full- potential linearized augmented plane wave total energy/atomic force method within the generalized gradient approximation. With Mo segregated to a clean Fe Sigma3(111) grain boundary, the embrittling potency of P increases from +0.19 to +0.88 eV, suggesting a strong detrimental ternary effect in P-Mo couple. Since Mo has a direct strengthening effect of -0.90 eV, the combined effect of a P-Mo couple is -0.02 eV, which confirms the experimental evidence that Mo overcompensates the P embrittlement. Since Mo has a larger atomic size than Fe, it reduces the volume available for P at the grain boundary and hence an increased elastic energy. On the other hand, the strong chemical bonding between Mo and the Fe surface makes the top Fe layers more saturated and hence a weakened vertical P-Fe bonding. Together, they induce a strong embrittling ternary effect in a P-Mo couple. This understanding can also explain the embrittling ternary effect of Mn on P embrittlement and is expected to be applicable to more general cases and instructive in quantum design of ultrahigh-strength alloys. (C) 2001 Elsevier Science Ltd. All rights reserved.

 

Geng, W. T. and A. J. Freeman (2001). "Multilayer relaxation of Cu(331)." Physical Review B 64, 115401.

            Both quantitative low-energy electron diffraction (LEED) analysis and theoretical predictions based on the embedded-atom method suggested that Cu(331) has a sequence of interlayer spacing contractions and expansions dissimilar to that expected for a terrace with three atomic layers. To understand the anomalous behavior of Cu(331), we have performed a comparative study of Cu(331) and Cu(211) by using the first-principles full-potential linearized augmented plane-wave method within the generalized gradient approximation. For both surfaces, our calculations reproduce the multilayer relaxation sequences given by LEED measurements. The electronic structure analysis indicates that the seemingly anomalous multilayer relaxation in Cu(331) can also be understood in the charge smoothing picture. However. the trend of charge smoothing in such transition metal surfaces depends not only on how many atomic layers a terrace has, but also on the bond-length-bond-order interplays.

 

Geller, C. B., W. Wolf, S. Picozzi, A. Continenza, R. Asahi, W. Mannstadt, A. J. Freeman and E. Wimmer (2001). "Computational band-structure engineering of III-V semiconductor alloys." Applied Physics Letters 79, 368.

            Accurate band structures of binary semiconductors AB (A=Al, Ga, In and B=P, As, Sb) and selected ternary III-V semiconductors were calculated using an all-electron screened exchange approach within the full potential linearized augmented plane- wave method. Fundamental band gaps and Gamma -L and Gamma -X separations in higher-lying conduction bands are predicted with an accuracy of a few tenths of 1 eV. Screened exchange also performs better than the local density approximation for calculating conduction-band effective masses. Highly n-doped InPAs materials with compositions near InP0.2As0.8 offer lower effective masses, greater optical band-gap shifts, and potentially higher electron mobility than n-doped InGaAs materials with comparable band gaps. (C) 2001 American Institute of Physics.

 

Continenza, A., S. Picozzi, W. T. Geng and A. J. Freeman (2001). "Coordination and chemical effects on the structural, electronic, and magnetic properties in Mn pnictides." Physical Review B 64, 085204.

            Simple structures of MnX binary compounds, namely hexagonal NiAs and zincblende, are studied as a function of the anion (X=Sb,As, P) by means of the full-potential linearized augmented plane wave (FLAPW) method within the local spin density and generalized gradient approximations. An accurate analysis of the structural, electronic, and magnetic properties reveals that the cubic structure greatly favors the magnetic alignment in these compounds leading to high magnetic moments and nearly half-metallic behavior for MnSb and MnAs. The effect of the anion chemical species is related to both its size and the possible hybridization with the Mn d states; both contributions are seen to hinder the magnitude of the magnetic moment for small and light anions. Our results are in very good agreement with experiment, where available, and show that the generalized gradient approximation is essential to correctly recover both the equilibrium volume and magnetic moment.

 

Cho, S. L., S. J. Youn, Y. Kim, A. DiVenere, G. K. L. Wong, A. J. Freeman and J. B. Ketterson (2001). "Polarity inversion in polar-nonpolar-polar heterostructures." Physical Review Letters 87, 126403.

            We have observed an epilayer-thickness-dependent polarity inversion for the growth of CdTe on Sb(Bi)/CdTe(111)B. For films with Sb(Bi) thicknesses of less than 40 Angstrom (15 Angstrom), the CdTe layer shows a B (Te-terminated) face, but it switches to an A (Cd-terminated) face for thicker layers. On the other hand, a CdTe layer grown on Bi(Sb)/Me(111)A always shows the A face regardless of Sb or Bi layer thicknesses. In order to address the observations we have performed ab initio calculations, which suggest that the polarity of a polar material on a nonpolar one results from the binding energy difference between the two possible surface configurations.

Zhong, L. P., R. Q. Wu, A. J. Freeman and G. B. Olson (2000). "Charge transfer mechanism of hydrogen-induced intergranular embrittlement of iron." Physical Review B 62, 13938.

            Impurity-induced reduction of intergranular cohesion-a major factor limiting the usable strength level of ultrahigh-strength steels-is particularly severe when aggravated by mobile hydrogen through environmental interaction, as in the case of hydrogen stress corrosion cracking. As an aid in establishing an understanding on the electronic level, the influence of hydrogen on the cohesion of an iron grain boundary was determined using the full-potential linearized augmented plane wave (FLAPW) method with the generalized gradient approximation. Through precise calculations on both grain boundary and free surface environments, we found that hydrogen is a strong embrittler. Analysis of the results in terms of structural relaxation, bonding character, and magnetic interactions shows that the hydrogen-iron chemical bond is stronger on the free surface and a charge-transfer mechanism is found to play a dominant role for the hydrogen-induced reduction of cohesion across the iron grain boundary. These results provide a quantitative explanation from first principles for the technologically important phenomenon of hydrogen-induced intergranular embrittlement.

 

Smith, R. W., W. T. Geng, C. B. Geller, R. Wu and A. J. Freeman (2000). "The effect of Li, He and Ca on grain boundary cohesive strength in Ni." Scripta Materialia 43, 957.

           

Profeta, G., A. Continenza, L. Ottaviano, W. Mannstadt and A. J. Freeman (2000). "Structural and electronic properties of the Sn/Si(111)root 3x root 3R30 degrees surface." Physical Review B 62, 1556.

            The structural and electronic properties of the Sn/Si(111)root 3x root 3 surface are determined by means of all-electron local density full-potential augmented plane wave thin film calculations. We find strong similarities with the more extensively studied isolectronic systems (Sn/Ge and Pb/Ge) as far as the relaxed structure and the electronic properties are concerned. In analogy with these systems, we find, within the local density approximation (LDA), two surface states weakly dispersed in the Brillouin zone which are responsible for the hexagonal patterns observed in scanning tunneling microscopy (STM) experiments. When our calculated results, including STM images, are compared with available structural data and STM images, we find that the LDA predictions well reproduce the electron distribution at the surface and the structural properties, leading to a complete description of this system at room temperature.

 

Picozzi, S., A. Continenza, R. Asahi, W. Mannstadt, A. J. Freeman, W. Wolf, E. Wimmer and C. B. Geller (2000). "Volume and composition dependence of direct and indirect band gaps in ordered ternary III-V semiconductor compounds: A screened-exchange LDA study." Physical Review B 61, 4677.

            Quasiparticle electronic band structures for ternary InGaAs and InGaSb systems, modeled by a luzonite structure, have been obtained using the screened-exchange local-density approximation approach. We focus our attention on energy transitions relevant for the electron-hole pair lifetime and, in particular, on the direct (E<(Gamma)over bar>-<(Gamma)over bar>) and lowest indirect (E<(Gamma)over bar>-(L) over bar) band gaps as a function of volume, common anionic species, and composition. All these degrees of freedom can be used to tune the transitions considered: as a first result, a large range of band gaps (almost 1 eV) can be obtained by varying the lattice constants within a few percent with respect to their equilibrium values. Moreover, the lowest indirect transition can be strongly reduced by replacing As with Sb. Finally, the composition dependence of these transitions shows that, due to symmetry properties of the potential in ternary luzonites, In- rich systems have the smallest indirect band gaps of all compounds investigated in this work.

 

Picozzi, S., A. Continenza, G. Satta, S. Massidda and A. J. Freeman (2000). "Metal-induced gap states and Schottky barrier heights at nonreactive GaN/noble-metal interfaces." Physical Review B 61, 16736.

            We present ab initio local density FLAPW calculations on nonreactive N-terminated [001] ordered GaN/Ag and GaN/Au interfaces and compare the results (such as metal induced gap states and Schottky barrier heights) with those obtained for GaN/Al, in order to understand the dependence of the relevant electronic properties on the deposited metal. Our results show that the density-of-gap states is appreciable only in the first semiconductor laver close to the interface. The decay length of the gap states in the semiconductor side is about 2.0+/-0.1 Angstrom and is independent of the deposited metal, therefore being to a good extent a bulk property of GaN. Our calculated values of the Schottky barrier heights are Phi(Bp) (GaN/Ag) = 0.87 eV and Phi(Bp) (GaN/Au) = 1.08 eV; both values are smaller than the GaN/Al Value [Phi(Bp) (GaN/Al) = 1.51 eV] and this quite large spread of values excludes the possibility of a Fermi-level pinning within the GaN band gap. Because of the low screening in GaN, the potential barrier at the junction is strongly affected by the structural arrangement of the first metal layer at the interface. This leads to quite large variations of the Schottky barrier height as a function of the metal. in contrast with the behavior of GaAs/metal interfaces.

 

Nakamura, K., M. Kim, A. J. Freeman, L. Zhong and J. Fernandez-de-Castro (2000). "Magnetism and magnetic anisotropy in FM NiFe, AFM NiMn, and their interface." IEEE Transactions on Magnetics 36, 3269.

            First-principles local density full-potential linearized augmented plane wave (FLAPW) calculations have been performed to investigate the magnetism and magnetic anisotropy in ferromagnetic NiFe, antiferromagnetic NiMn and their interface. We found that the magnetic moments of the Fe atom in NiFe and that of the Mn atom in NiMn are enhanced at their surface compared to their bulk states, The magneto-crystalline anisotropy (MCA) is also drastically influenced by the environment. Bulk NiFe and NiMn show slightly perpendicular and in-plane easy magnetic axes, respectively, However, the effects of the surface change their MCA direction and strength. For the interface between NiFe and NiMn, where the magneto-crystalline anisotropy is dominated by the NiMn layer, the interface effects may be important for discussing the exchange bias.

 

Lee, J. I., S. C. Hong, W. Mannstadt and A. J. Freeman (2000). "Rippled surface structure and electronic and magnetic properties of Ni3Al(001)." Physical Review B 62, 6982.

            Structural, electronic, and magnetic properties of Ni3Al(001) are investigated by the all-electron thin him full-potential linearized augmented plane wave method based on the local density approximation. A stable rippled surface atomic geometry is determined by atomic force and total energy calculations. The surface Ni atoms contract down to the bulk region by 4.2% of the bulk interlayer spacing while all of the other atoms including the surface AI atoms remain close to their bulk positions. The amount of rippling found by calculation (0.06 Angstrom) is almost within experimental error (0.02+/-0.03 Angstrom). Charge densities, calculated work functions and densities of states for the relaxed rippled and unrelaxed surfaces are reported. The spin polarized calculation predicts that the Ni3Al(001) surface is "magnetically dead," unlike the bulk region.

 Kim, M., W. T. Geng, A. J. Freeman, L. P. Zhong and J. Fernandez-de-Castro (2000). "First-principles calculations for the structural and magnetic properties of ordered NiFe(001) thin films with and without a Ta overlayer." Journal of Applied Physics 87, 5735.

            The first-principles calculations within the local density approximation using the full potential linearized augmented plane wave (FLAPW) method were performed to investigate the structural and magnetic properties of the Ta/NiFe interface for both clean NixFe1-x (001) thin films and with a Ta overlayer. To study the composition dependence, we adopted x = 0.5 (L1(0) structure with either Fe or Ni layers on the surface) and 0.75 (L1(2) structure with mixed Fe-Ni layers or with Ni on the surface). The equilibrium overlayer/substrate distance and the preferred site position of Ta were obtained by tructural optimization employing atomic-force calculations and total energy comparisons for several possible adsorption sites of Ta. By comparing with results for the clean surface of five-layer NixFe1-x (001) films, we found that Ta has a significant detrimental effect on the magnetic properties of NiFe with its induced magnetic moment coupled ferro- or antiferro- magnetically with the substrate depending sensitively on the surface layer. (C) 2000 American Institute of Physics. [S0021- 8979(00)73908-9].

 

Geng, W. T., A. J. Freeman, R. Wu and G. B. Olson (2000). "Effect of Mo and Pd on the grain-boundary cohesion of Fe." Physical Review B 62, 6208.

            The effects of Mo and Pd segregation on the cohesion of the Fe Sigma 3(111) grain boundary an investigated by using the first- principles full-potential linearized augmented-plane-wave total-energy-atomic-force method with the generalized gradient approximation. Based on the Rice-Wang model, our total energy calculations show that Mo has a significant beneficial effect on the Fe grain-boundary cohesion, while Pd behaves as a weak embrittler. An analysis of the geometry optimization indicates that Mo has a moderate atomic size to fit well in the grain- boundary hole, whereas Pd introduces a larger perturbation on the atomic structure near the grain boundary. The elastic energy associated with the Mo and Pd segregation is estimated with a rigid environment approximation. It is found that both Mo and Pd introduce a beneficial volume effect. Studies of the electronic structures show that its strong bonding capability makes Mo a cohesion enhancer (-0.90 eV) for the Fe Sigma 3(111) grain-boundary. By comparison, its weak bonding capability leads Pd to be a weak embrittler (+0.08 eV). Our first- principles quantum-mechanical results support the main idea of the atomistic theories in that the elemental cohesive energy difference between the substitutional element and the host element plays an important role in determining its effect on the grain-boundary cohesion. However, the numerical results for Pd, which has a similar elemental cohesive energy to that of Fe, point to the importance of the role played by the volume effect. It is expected that in a lower-angle Fe grain boundary which has a larger grain-boundary volume expansion, Pd can possibly become a cohesion enhancer.

           

Asahi, R., Y. Taga, W. Mannstadt and A. J. Freeman (2000). "Electronic and optical properties of anatase TiO2." Physical Review B 61, 7459.

            First-principles calculations using the full-potential linearized augmented plane-wave method have been performed to investigate detailed electronic and optical properties of TiO2 in the anatase structure. The fully optimized structure, obtained by minimizing the total energy and atomic fords, are in good agreement with experiment. Stabilization of the structure by the trade off between a favorable coordination in the sp(2) hybridization and the Coulomb repulsion among oxygen atoms is also demonstrated. We calculate band structure, densities of states and charge densities, and interpret their features in terms of the bonding structure in the molecular orbital picture. The optical properties, calculated within the dipole approximation, are found to agree with recent experiments on single crystals of anatase TiO2. Near the absorption edge, the results show a significant optical anisotropy in the components parallel and perpendicular to the c axis. We demonstrate that this large dichroism results from the existence of nonbonding d(xy) orbitals located at the bottom of the conduction bands, which allows direct dipole transitions dominantly for the perpendicular component.

 

Asahi, R., W. Mannstadt and A. J. Freeman (2000). "Screened-exchange LDA methods for films and superlattices with applications to the Si(100)2X1 surface and InAs/InSb superlattices." Physical Review B 62, 2552.

            We have developed efficient screened-exchange local-density approximation (sX-LDA) methods for films and superlattices (FLM/SL) with which to calculate self-consistent electronic structures for both occupied and unoccupied states. Considering nonuniform charge densities and local-field effects in the z direction for FLM/SL, we have employed nonlocal Thomas-Fermi wave vectors to define the screened-exchange interaction. Three methods, for bulk, superlattice, and film, have been implemented in the full-potential linearized augmented plane wave method. The film sX-LDA method was then applied to the Si(100)2x1 surface. The calculated occupied surface states show very good agreement with experiment. On the other hand, an underestimation of the correction to the unoccupied surface states, by about 0.2 eV, was estimated in comparison with available GW calculations. The ionization energy of Si was evaluated with the film geometry to be 5.35 eV by virtue of the quasiparticle corrections, showing good agreement with the experimental value of 5.15 +/- 0.08 eV. Rie also present an application of the superlattice sX-LDA method to [001] ordered InAs/InSb heterojunctions and superlattices. Band gaps and band offsets under strained conditions were directly calculated by sX-LDA without any experimental data as input. Slightly larger valence-band offsets than the LDA results. by about 0.08 eV, agree with the consequence of the GW calculations, indicating an increase of the potential negativity in the InAs region. This potential change along with the charge redistribution at the interface is found to be crucial to evaluate accurate band gaps of the superlattices.

 

Wu, R. Q. and A. J. Freeman (1999). "Spin-orbit induced magnetic phenomena in bulk metals and their surfaces and interfaces." Journal of Magnetism and Magnetic Materials 200, 498.

            First-principles electronic structure studies based on local spin density functional theory and performed on extremely complex simulations of ever increasingly realistic systems, play a very important role in explaining and predicting surface and interface magnetism. This review deals with what is a major issue for first-principles theory, namely the theoretical/computational treatment of the weak spin-orbit coupling in magnetic transition metals and their alloys and its important physical consequences: magneto-crystalline anisotropy, magnetostriction, magneto-optical Kerr effects and X-ray magnetic circular dichroism. As is demonstrated, extensive first-principles calculations and model analyses now provide simple physical insights and guidelines to search for new magnetic recording and sensor materials. (C) 1999 Elsevier Science B.V. All rights reserved.

 

Shick, A. B., Y. N. Gornostyrev and A. J. Freeman (1999). "Magnetoelastic mechanism of spin-reorientation transitions at step edges." Physical Review B 60, 3029.

            The symmetry-induced magnetic anisotropy due to monoatomic steps at strained Ni films is determined using results of first-principles relativistic full-potential linearized augmented-plane-wave calculations and an analogy with the Neel model. We show that there is a magnetoelastic anisotropy contribution to the uniaxial magnetic anisotropy energy in the vicinal plane of a stepped surface. In addition to the known spin-direction reorientation transition at a flat Ni/Cu(001) surface, we propose a spin-direction reorientation transition in the vicinal plane for a stepped Ni/Cu surface due to the magnetoelastic anisotropy. We show that with an increase of Ni- film thickness, the magnetization in the vicinal plane turns perpendicular to the step edge at a critical thickness calculated to be in the range of 16-24 Ni layers for the Ni/Cu(1,1,13) stepped surface. [S0163-1829(99)13629-4].

 

Shick, A. B., J. B. Ketterson, D. L. Novikov and A. J. Freeman (1999). "Electronic structure, phase stability, and semimetal- semiconductor transitions in Bi." Physical Review B 60, 15484.

            The structural stability of bulk Bi is studied using the local- density full-potential linear muffin-tin orbital method. The effect of both the trigonal shear angle and internal displacement on the electronic structure is determined. It is shown that the internal displacement changes the Bi electronic structure from a metal to a semimetal, in qualitative agreement with a Jones-Peierls-type transition. The total energy is calculated to have a double-well dependence on the internal displacement, and to provide a stabilization of the trigonal phase. We show that an increase of trigonal shear angle (towards the cubic value of 60 degrees) leads to a semimetal- semiconductor transition in Bi. Using coherency strain arising from a film/substrate lattice constant mismatch, this may provide a way to induce semiconducting behavior in Bi films, and with it to control their thermoelectric properties. [S0163- 1829(99)02947-1].

 

Picozzi, S., A. Continenza and A. J. Freeman (1999). "Surface states and Fermi-level pinning at clean and Al covered GaN surfaces." Physical Review B 59, 1609.

            Ab initio electronic structure calculations have been performed for the clean GaN surface (both N- and Ga-terminated) and for the GaN/Al 1-ML system, obtained from the clean surface by addition of an Al monolayer. In the most stable configuration, Al is found to occupy the bridge site at the clean unreconstructed N-terminated surface. All the systems studied show an appreciable surface charge density due to midgap states strongly localized in the surface region; this localization is particularly enhanced for the N-terminated case, showing the high instability of the unreconstructed clean surface. On the other hand, the Ga-terminated surface is found to be more stable and to behave very similarly to the Al-covered system. The Fermi level E-F is shifted toward the conduction band with respect to the clean N-terminated surface. Additional coverage of Al layers does not significantly affect the initial position of E-F within the GaN band gap, showing that the Schottky barrier height is well established after the first Al layer. [S0163-829(98)02643-5].

 

Lee, J. I., W. Mannstadt and A. J. Freeman (1999). "Multilayer-relaxed structure of the (1 x 2) Pt(110) surface." Physical Review B 59, 1673.

            The multilayer-relaxed structure and electronic properties of the (1 x 2) Pt(110) surface have been investigated by the self- consistent all-electron full-potential linearized augmented plane-wave method. The relaxed geometry, determined by total energy and atomic force calculations, shows large contractions in the first and second interlayer spacings, significant buckling in the third layer, and a lateral displacement in the fourth (center) layer of the slab. In general, our calculated results are consistent with experimental data. The microscopic origin of the relaxed structure is discussed using the calculated electronic structures. The large inward relaxation of the surface atoms is attributed to the more localized nature of their 5d electrons, which weakens the d-d hybridization. [S0163-1829(99)06804-6].

 

Kim, M., A. J. Freeman and R. Q. Wu (1999). "Surface effects and structural dependence of magneto-optical spectra: UItrathin Co films and CoPtn alloys and multilayers." Physical Review B 59, 9432.

            First-principles calculations of the optical conductivities and magneto-optical Kerr effects for bulk and thin films of Co and for CoPtn alloys (n = 1 and 3) and multilayer (n = 3) are carried out using the full potential linearized augmented plane wave method with the local spin density approximation. For fee bulk Co, we obtain a negligible dependence of the polar Kerr rotation on the direction of magnetization. By contrast, the calculated Kerr spectra of CoPtn are found to be very sensitive to structural changes. For Co(001) and Co(111) ultrathin films, the surface effects cause a redshift for the high-energy peaks in the magneto-optic Kerr effect (MOKE) spectra, mainly due to band narrowing. Except for the free standing Co(001) monolayer, the profile of their MOKE spectra, however, is very close to that of bulb Co. [S0163-1829(99)04313-1].

 

Geng, W. T., A. J. Freeman, R. Wu, C. B. Geller and J. E. Raynolds (1999). "Embrittling and strengthening effects of hydrogen, boron, and phosphorus on a Sigma 5 nickel grain boundary." Physical Review B 60, 7149.

            The embrittling and strengthening effects of hydrogen, boron, and phosphorus on a Sigma 5(210) [100] nickel grain boundary are investigated by means of the full-potential linearized augmented plane-wave method with the generalized-gradient approximation formula. Optimized geometries for both the free surface and grain-boundary systems are obtained by atomic-force calculations. The results obtained show that hydrogen and phosphorus are embrittlers and that boron acts as a cohesion enhancer. An analysis of the atomic, electronic, and magnetic structures indicates that atomic size and the bonding behavior of the impurity with the surrounding nickel atoms play important roles in determining its relative embrittling or cohesion enhancing behavior. [S0163-1829(99)02926-4].

 

Freeman, A. J., R. Q. Wu, M. Y. Kim and V. I. Gavrilenko (1999). "Magnetism, magneto-crystalline anisotropy, magnetostriction and MOKE at surfaces and interfaces." Journal of Magnetism and Magnetic Materials 203, 1.

            A major issue for first-principles theory in magnetism is the treatment of the weak spin-orbit coupling (SOC) and its subsequent effects in magnetic transition metal bulk, surfaces and multilayers. Using either a perturbative or a self- consistent approach for SOC, we have recently investigated important phenomena such as magnetic crystalline anisotropy (MCA), magnetostriction, and magneto-optical Kerr effects (MOKE) in various transition metal systems using the highly precise local density full-potential linearized augmented plane wave (FLAPW) method. With the aid of accurate total energy and atomic force approaches based on the LDA and GGA formalisms, the atomic structures of all the surfaces and interfaces are fully relaxed. Excellent agreement with experiment has been achieved for most of the systems investigated in terms of their equilibrium geometries, MCA energies, magnetostrictive coefficients and MOKE spectra. (C) 1999 Elsevier Science B.V. All rights reserved.

 

Charache, G. W., D. M. DePoy, J. E. Raynolds, P. F. Baldasaro, K. E. Miyano, T. Holden, F. H. Pollak, P. R. Sharps, M. L. Timmons, C. B. Geller, W. Mannstadt, R. Asahi, A. J. Freeman and W. Wolf (1999). "Moss-Burstein and plasma reflection characteristics of heavily doped n-type InxGa1-xAs and InPyAs1-y." Journal of Applied Physics 86, 452.

            Degenerately doped (> 10(19) cm(-3)) n-type InxGa1-xAs (x similar to 0.67) and InPyAs1-y (y similar to 0.65) possess a number of intriguing electrical and optical properties relevant to electro-optic devices and thermophotovoltaic devices in particular. Due to the low electron effective mass of these materials (m*< 0.2) and the demonstrated ability to incorporate n-type dopants into the high 10(19) cm(-3) range, both the Moss-Burstein band gap shift and plasma reflection characteristics are particularly dramatic. For InGaAs films with a nominal undoped band gap of 0.6 eV and N=5x10(19) cm(- 3), the fundamental absorption edge increased to 1.27 eV. InPAs films exhibit a shorter plasma wavelength (lambda(p)similar to 5 mu m) in comparison to InGaAs films (lambda(p)similar to 6 mu m) with similar doping concentrations. The behavior of the plasma wavelength and the fundamental absorption edge are investigated in terms of conduction band nonparabolicity and Gamma-L valley separation using detailed band structure measurements and calculations. (C) 1999 American Institute of Physics. [S0021-8979(99)01612-6].

 

Asahi, R., W. Mannstadt and A. J. Freeman (1999). "Optical properties and electronic structures of semiconductors with screened-exchange LDA." Physical Review B 59, 7486.

            Results of first-principles calculations of the optical properties and electronic structure, determined by employing the screened-exchange LDA method to obtain accurate electronic structures as well as self-consistent eigenfunctions, are presented for some cubic semiconductors (Si, Ge, InSb, and GaAs). As implemented in the full-potential linearized augmented-plane-wave method, this approach does not require any adjustable parameters to calculate optical properties. The inadequacy of using the transverse expression of the dielectric function (involving ep matrix elements) is shown numerically for the optical spectra and dielectric constants even when the renormalization factor associated with the LDA eigenvalues is included. This is because the nonlocality of the Hamiltonian and the difference of the self-consistent eigenfunctions from the LDA ones are not negligible. In contrast, the optical properties evaluated using the longitudinal expression (with e(iq.r) matrix elements) yield excellent agreement with experiment considering the neglect of lifetime, local-field, and excitonic effects. [S0163-1829(99)01708-7].



© Quantum Materials Design, Inc. 2002