Local Charge Excesses in Metallic Alloys: a Local Field Coherent Potential Approximation Theory

Electronic structure calculations performed on very large supercells have shown that the local charge excesses in metallic alloys are related through simple linear relations to the local electrostatic field resulting from distribution of charges in the whole crystal. By including local external fields in the single site Coherent Potential Approximation theory, we develop a novel theoretical scheme in which the local charge excesses for random alloys can be obtained as the responses to local external fields. Our model maintains all the computational advantages of a single site theory but allows for full charge relaxation at the impurity sites. Through applications to CuPd and CuZn alloys, we find that, as a general rule, non linear charge rearrangements occur at the impurity site as a consequence of the complex phenomena related with the electronic screening of the external potential. This nothwithstanding, we observe that linear relations hold between charge excesses and external potentials, in quantitative agreement with the mentioned supercell calculations, and well beyond the limits of linearity for any other site property.Comment: 11 pages, 1 table, 7 figure

First principles studies of band offsets at heterojunctions and of surface reconstruction using Gaussian dual-space density functional theory

The use of localized Gaussian basis functions for large scale first principles density functional calculations with periodic boundary conditions (PBC) in 2 dimensions and 3 dimensions has been made possible by using a dual space approach. This new method is applied to the study of electronic properties of II–VI (II=Zn, Cd, Hg; VI=S, Se, Te, Po) and III–V (III=Al, Ga; V=As, N) semiconductors. Valence band offsets of heterojunctions are calculated including both bulk contributions and interfacial contributions. The results agree very well with available experimental data. The p(2 × 1) cation terminated surface reconstructions of CdTe and HgTe (100) are calculated using the local density approximation (LDA) with two-dimensional PBC and also using the ab initio Hartree–Fock (HF) method with a finite cluster. The LDA and HF results do not agree very well

Systematic treatment of displacements, strains and electric fields in density-functional perturbation theory

The methods of density-functional perturbation theory may be used to calculate various physical response properties of insulating crystals including elastic, dielectric, Born charge, and piezoelectric tensors. These and other important tensors may be defined as second derivatives of the total energy with respect to atomic-displacement, electric-field, or strain perturbations, or as mixed derivatives with respect to two of these perturbations. The resulting tensor quantities tend to be coupled in complex ways in polar crystals, giving rise to a variety of variant definitions. For example, it is generally necessary to distinguish between elastic tensors defined under different electrostatic boundary conditions, and between dielectric tensors defined under different elastic boundary conditions. Here, we describe an approach for computing all of these various response tensors in a unified and systematic fashion. Applications are presented for two materials, wurtzite ZnO and rhombohedral BaTiO3, at zero temperature.Comment: 14 pages. Uses REVTEX macros. Also available at http://www.physics.rutgers.edu/~dhv/preprints/xfw_sys/index.htm

Onset of magnetism in B2 transition metals aluminides

Ab initio calculation results for the electronic structure of disordered bcc Fe(x)Al(1-x) (0.4<x<0.75), Co(x)Al(1-x) and Ni(x)Al(1-x) (x=0.4; 0.5; 0.6) alloys near the 1:1 stoichiometry, as well as of the ordered B2 (FeAl, CoAl, NiAl) phases with point defects are presented. The calculations were performed using the coherent potential approximation within the Korringa-Kohn-Rostoker method (KKR-CPA) for the disordered case and the tight-binding linear muffin-tin orbital (TB-LMTO) method for the intermetallic compounds. We studied in particular the onset of magnetism in Fe-Al and Co-Al systems as a function of the defect structure. We found the appearance of large local magnetic moments associated with the transition metal (TM) antisite defect in FeAl and CoAl compounds, in agreement with the experimental findings. Moreover, we found that any vacancies on both sublattices enhance the magnetic moments via reducing the charge transfer to a TM atom. Disordered Fe-Al alloys are ferromagnetically ordered for the whole range of composition studied, whereas Co-Al becomes magnetic only for Co concentration >0.5.Comment: 11 pages with 9 embedded postscript figures, to be published in Phys.Rev.

LATTICE DISTORTION NEAR VACANCIES IN DIAMOND AND SILICON .1.

A dynamical relaxation procedure, coupled with a valence force potential, has been used to calculate the distortion around point defects in a diamond-type crystal. The method has been applied to the vacancy in diamond and silicon. The response of the lattice to symmetrized forces on the nearest neighbours to the vacancy was calculated. The results can be used in estimates of point defect properties which depend on lattice distortion, including the jahn-teller effect, and formation energies. The ratios of the atomic displacements under uniform external stresses for the perfect lattice and for the lattice with a vacancy are also determined

Vibrational and thermoelastic properties of bcc iron from selected EAM potentials

A comprehensive, critical study of the vibrational, thermodynamic and thermoelastic properties of bcc iron is presented, using well established semi-empirical embedded-atom method potentials available in the literature. Classical molecular dynamics simulations are used to address temperature effects, where dynamical matrices are constructed as a time average of the second moment of the atomic displacements. The $C_{11}, C_{44}, C'$ elastic constants are then obtained from the sound velocities along high symmetry directions in reciprocal space. Results are compared to ultrasonic measurements and highlight the limitations of the potentials considered here in describing thermoelastic properties.Comment: 22 pages, 10 figures, 1 tabl

Origin of the Verwey transition in magnetite: Group theory, electronic structure, and lattice dynamics study

The Verwey phase transition in magnetite has been analyzed using the group theory methods. It is found that two order parameters with the symmetries $X_3$ and $\Delta_5$ induce the structural transformation from the high-temperature cubic to the low-temperature monoclinic phase. The coupling between the order parameters is described by the Landau free energy functional. The electronic and crystal structure for the cubic and monoclinic phases were optimized using the {\it ab initio} density functional method. The electronic structure calculations were performed within the generalized gradient approximation including the on-site interactions between 3d electrons at iron ions -- the Coulomb element $U$ and Hund's exchange $J$. Only when these local interactions are taken into account, the phonon dispersion curves, obtained by the direct method for the cubic phase, reproduce the experimental data. It is shown that the interplay of local electron interations and the coupling to the lattice drives the phonon order parameters and is responsible for the opening of the gap at the Fermi energy. Thus, it is found that the metal-insulator transition in magnetite is promoted by local electron interactions, which significantly amplify the electron-phonon interaction and stabilize weak charge order coexisting with orbital order of the occupied $t_{2g}$ states at Fe ions. This provides a scenario to understand the fundamental problem of the origin of the Verwey transition in magnetite.Comment: 17 pages, 5 figures, 8 tables. Accepted version to be published in Phys. Rev.