Development of a tight-binding potential for bcc-Zr. Application to the study of vibrational properties

We present a tight-binding potential based on the moment expansion of the density of states, which includes up to the fifth moment. The potential is fitted to bcc and hcp Zr and it is applied to the computation of vibrational properties of bcc-Zr. In particular, we compute the isothermal elastic constants in the temperature range 1200K < T < 2000K by means of standard Monte Carlo simulation techniques. The agreement with experimental results is satisfactory, especially in the case of the stability of the lattice with respect to the shear associated with C'. However, the temperature decrease of the Cauchy pressure is not reproduced. The T=0K phonon frequencies of bcc-Zr are also computed. The potential predicts several instabilities of the bcc structure, and a crossing of the longitudinal and transverse modes in the (001) direction. This is in agreement with recent ab initio calculations in Sc, Ti, Hf, and La.Comment: 14 pages, 6 tables, 4 figures, revtex; the kinetic term of the isothermal elastic constants has been corrected (Eq. (4.1), Table VI and Figure 4

An Effective-Medium Tight-Binding Model for Silicon

A new method for calculating the total energy of Si systems is presented. The method is based on the effective-medium theory concept of a reference system. Instead of calculating the energy of an atom in the system of interest a reference system is introduced where the local surroundings are similar. The energy of the reference system can be calculated selfconsistently once and for all while the energy difference to the reference system can be obtained approximately. We propose to calculate it using the tight-binding LMTO scheme with the Atomic-Sphere Approximation(ASA) for the potential, and by using the ASA with charge-conserving spheres we are able to treat open system without introducing empty spheres. All steps in the calculational method is {\em ab initio} in the sense that all quantities entering are calculated from first principles without any fitting to experiment. A complete and detailed description of the method is given together with test calculations of the energies of phonons, elastic constants, different structures, surfaces and surface reconstructions. We compare the results to calculations using an empirical tight-binding scheme.Comment: 26 pages (11 uuencoded Postscript figures appended), LaTeX, CAMP-090594-

Two-band second moment model and an interatomic potential for caesium

A semi-empirical formalism is presented for deriving interatomic potentials for materials such as caesium or cerium which exhibit volume collapse phase transitions. It is based on the Finnis-Sinclair second moment tight binding approach, but incorporates two independent bands on each atom. The potential is cast in a form suitable for large-scale molecular dynamics, the computational cost being the evaluation of short ranged pair potentials. Parameters for a model potential for caesium are derived and tested

Electronic properties of metal induced gap states at insulator/metal interfaces -- dependence on the alkali halide and the possibility of excitonic mechanism of superconductivity

Motivated from the experimental observation of metal induced gap states (MIGS) at insulator/metal interfaces by Kiguchi {\it et al.} [Phys. Rev. Lett. {\bf 90}, 196803 (2003)], we have theoretically investigated the electronic properties of MIGS at interfaces between various alkali halides and a metal represented by a jellium with the first-principles density functional method. We have found that, on top of the usual evanescent state, MIGS generally have a long tail on halogen sites with a $p_z$-like character, whose penetration depth ($\lambda$) is as large as half the lattice constant of bulk alkali halides. This implies that $\lambda$, while little dependent on the carrier density in the jellium, is dominated by the lattice constant (hence by energy gap) of the alkali halide, where $\lambda_{\rm LiF} < \lambda_{\rm LiCl} < \lambda_{\rm LiI}$. We also propose a possibility of the MIGS working favorably for the exciton-mediated superconductivity.Comment: 7 pages, 9 figure

Ab initio Calculations of Multilayer Relaxations of Stepped Cu Surfaces

We present trends in the multilayer relaxations of several vicinals of Cu(100) and Cu(111) of varying terrace widths and geometry. The electronic structure calculations are based on density functional theory in the local density approximation with norm-conserving, non-local pseudopotentials in the mixed basis representation. While relaxations continue for several layers, the major effect concentrates near the step and corner atoms. On all surfaces the step atoms contract inwards, in agreement with experimental findings. Additionally, the corner atoms move outwards and the atoms in the adjacent chain undergo large inward relaxation. Correspondingly, the largest contraction (4%) is in the bond length between the step atom and its bulk nearest neighbor (BNN), while that between the corner atom and BNN is somewhat enlarged. The surface atoms also display changes in registry of upto 1.5%. Our results are in general in good agreement with LEED data including the controversial case of Cu(511). Subtle differences are found with results obtained from semi-empirical potentials.Comment: 21 pages and 3 figure

First-principles study of the polar O-terminated ZnO surface in thermodynamic equilibrium with oxygen and hydrogen

Using density-functional theory in combination with a thermodynamic formalism we calculate the relative stability of various structural models of the polar O-terminated (000-1)-O surface of ZnO. Model surfaces with different concentrations of oxygen vacancies and hydrogen adatoms are considered. Assuming that the surfaces are in thermodynamic equilibrium with an O2 and H2 gas phase we determine a phase diagram of the lowest-energy surface structures. For a wide range of temperatures and pressures we find that hydrogen will be adsorbed at the surface, preferentially with a coverage of 1/2 monolayer. At high temperatures and low pressures the hydrogen can be removed and a structure with 1/4 of the surface oxygen atoms missing becomes the most stable one. The clean, defect-free surface can only exist in an oxygen-rich environment with a very low hydrogen partial pressure. However, since we find that the dissociative adsorption of molecular hydrogen and water (if also the Zn-terminated surface is present) is energetically very preferable, it is very unlikely that a clean, defect-free (000-1)-O surface can be observed in experiment.Comment: 10 pages, 4 postscript figures. Uses REVTEX and epsf macro

Close-Packing of Clusters: Application to Al_100

The lowest energy configurations of close-packed clusters up to N=110 atoms with stacking faults are studied using the Monte Carlo method with Metropolis algorithm. Two types of contact interactions, a pair-potential and a many-atom interaction, are used. Enhanced stability is shown for N=12, 26, 38, 50, 59, 61, 68, 75, 79, 86, 100 and 102, of which only the sizes 38, 75, 79, 86, and 102 are pure FCC clusters, the others having stacking faults. A connection between the model potential and density functional calculations is studied in the case of Al_100. The density functional calculations are consistent with the experimental fact that there exist epitaxially grown FCC clusters starting from relatively small cluster sizes. Calculations also show that several other close-packed motifs existwith comparable total energies.Comment: 9 pages, 7 figure

Ab initio and finite-temperature molecular dynamics studies of lattice resistance in tantalum

This manuscript explores the apparent discrepancy between experimental data and theoretical calculations of the lattice resistance of bcc tantalum. We present the first results for the temperature dependence of the Peierls stress in this system and the first ab initio calculation of the zero-temperature Peierls stress to employ periodic boundary conditions, which are those best suited to the study of metallic systems at the electron-structure level. Our ab initio value for the Peierls stress is over five times larger than current extrapolations of experimental lattice resistance to zero-temperature. Although we do find that the common techniques for such extrapolation indeed tend to underestimate the zero-temperature limit, the amount of the underestimation which we observe is only 10-20%, leaving open the possibility that mechanisms other than the simple Peierls stress are important in controlling the process of low temperature slip.Comment: 12 pages and 9 figure

Screened Coulomb interactions in metallic alloys: II Screening beyond the single-site and atomic sphere approximations

A quantitative description of the configurational part of the total energy of metallic alloys with substantial atomic size difference cannot be achieved in the atomic sphere approximation: It needs to be corrected at least for the multipole moment interactions in the Madelung part of the one-electron potential and energy. In the case of a random alloy such interactions can be accounted for only by lifting the atomic sphere and single-site approximations, in order to include the polarization due to local environment effects. Nevertheless a simple parameterization of the screened Coulomb interactions for the ordinary single-site methods, including the generalized perturbation method, is still possible. We obtained such a parameterization for bulk and surface NiPt alloys, which allows one to obtain quantitatively accurate effective interactions in this system.Comment: 24 pages, 2 figure