Activated O2 dissociation and formation of oxide islands on the Be(0001) surface: Another atomistic model for metal oxidation

By simulating the dissociation of O2 molecules on the Be(0001) surface using the first-principles molecular dynamics approach, we propose a new atomistic model for the surface oxidation of sp metals. In our model, only the dissociation of the first oxygen molecule needs to overcome an energy barrier, while the subsequent oxygen molecules dissociate barrierlessly around the adsorption area. Consequently, oxide islands form on the metal surface, and grow up in a lateral way. We also discover that the firstly dissociated oxygen atoms are not so mobile on the Be(0001) surface, as on the Al(111) surface. Our atomistic model enlarges the knowledge on metal surface oxidations by perfectly explaining the initial stage during the surface oxidation of Be, and might be applicable to some other sp metal surfaces.Comment: 5 pages, 4 figure

First-principles calculations for the adsorption of water molecules on the Cu(100) surface

First-principles density-functional theory and supercell models are employed to calculate the adsorption of water molecules on the Cu(100) surface. In agreement with the experimental observations, the calculations show that a H2O molecule prefers to bond at a one-fold on-top (T1) surface site with a tilted geometry. At low temperatures, rotational diffusion of the molecular axis of the water molecules around the surface normal is predicted to occur at much higher rates than lateral diffusion of the molecules. In addition, the calculated binding energy of an adsorbed water molecule on the surfaces is significantly smaller than the water sublimation energy, indicating a tendency for the formation of water clusters on the Cu(100) surface.Comment: 5 pages, 3 figures, submitted to Phys. Rev.

Metal-to-metal charge transfer between dopant and host ions: Photoconductivity of Yb-doped CaF2 and SrF2 crystals

Dopant-to-host electron transfer is calculated using ab initio wavefunction-based embedded cluster methods for Yb/Ca pairs in CaF2 and Yb/Sr pairs in SrF2 crystals to investigate the mechanism of photoconductivity. The results show that, in these crystals, dopant-to-host electron transfer is a two-photon process mediated by the 4fN-15d excited states of Y b2+: these are reached by the first photon excitation; then, they absorb the second photon, which provokes the Y b2+ + Ca2+ (Sr2+) → Y b3+ + Ca+ (Sr+) electron phototransfer. This mechanism applies to all the observed Y b2+ 4f-5d absorption bands with the exception of the first one: Electron transfer cannot occur at the first band wavelengths in CaF2:Y b2+ because the Y b3+-Ca+ states are not reached by the two-photon absorption. In contrast, Yb-to-host electron transfer is possible in SrF2:Y b2+ at the wavelengths of the first 4f-5d absorption band, but the mechanism is different from that described above: first, the two-photon excitation process occurs within the Y b2+ active center, then, non-radiative Yb-to-Sr electron transfer can occur. All of these features allow to interpret consistently available photoconductivity experiments in these materials, including the modulation of the photoconductivity by the absorption spectrum, the differences in photoconductivity thresholds observed in both hosts, and the peculiar photosensitivity observed in the SrF2 host, associated with the lowest 4f-5d bandThis work was partly supported by a grant from Ministerio de Economía y Competitividad, Spain (Dirección General de Investigación y Gestión del Plan Nacional de I+D+i, Grant Nos. MAT2011-24586 and MAT2014-54395-P

Transferable Pair Potentials for CdS and ZnS Crystals

A set of interatomic pair potentials is developed for CdS and ZnS crystals. We show that a simple energy function, which has been used to describe the properties of CdSe [J. Chem. Phys. 116, 258 (2002)], can be parametrized to accurately describe the lattice and elastic constants, and phonon dispersion relations of bulk CdS and ZnS in the wurtzite and rocksalt crystal structures. The predicted coexistence pressure of the wurtzite and rocksalt structures, as well as the equation of state are in good agreement with experimental observations. These new pair potentials enable the study of a wide range of processes in bulk and nanocrystalline II-VI semiconductor materials

Spin-Orbit Coupling in Iridium-Based 5d Compounds Probed by X-ray Absorption Spectroscopy

We have performed x-ray absorption spectroscopy (XAS) measurements on a series of Ir-based 5d transition metal compounds, including Ir, IrCl3, IrO2, Na2IrO3, Sr2IrO4, and Y2Ir2O7. By comparing the intensity of the "white-line" features observed at the Ir L2 and L3 absorption edges, it is possible to extract valuable information about the strength of the spin-orbit coupling in these systems. We observe remarkably large, non-statistical branching ratios in all Ir compounds studied, with little or no dependence on chemical composition, crystal structure, or electronic state. This result confirms the presence of strong spin-orbit coupling effects in novel iridates such as Sr2IrO4, Na2IrO3, and Y2Ir2O7, and suggests that even simple Ir-based compounds such as IrO2 and IrCl3 may warrant further study. In contrast, XAS measurements on Re-based 5d compounds, such as Re, ReO2, ReO3, and Ba2FeReO6, reveal statistical branching ratios and negligible spin-orbit coupling effects.Comment: 9 pages, 4 figure

First Principles Studies on 3-Dimentional Strong Topological Insulators: Bi2Te3, Bi2Se3 and Sb2Te3

Bi2Se3, Bi2Te3 and Sb2Te3 compounds are recently predicted to be 3-dimentional (3D) strong topological insulators. In this paper, based on ab-initio calculations, we study in detail the topological nature and the surface states of this family compounds. The penetration depth and the spin-resolved Fermi surfaces of the surface states will be analyzed. We will also present an procedure, from which highly accurate effective Hamiltonian can be constructed, based on projected atomic Wannier functions (which keep the symmetries of the systems). Such Hamiltonian can be used to study the semi-infinite systems or slab type supercells efficiently. Finally, we discuss the 3D topological phase transition in Sb2(Te1-xSex)3 alloy system.Comment: 8 pages,17 figure

Phase Decomposition and Chemical Inhomogeneity in Nd2-xCexCuO4

Extensive X-ray and neutron scattering experiments and additional transmission electron microscopy results reveal the partial decomposition of Nd2-xCexCuO4 (NCCO) in a low-oxygen-fugacity environment such as that typically realized during the annealing process required to create a superconducting state. Unlike a typical situation in which a disordered secondary phase results in diffuse powder scattering, a serendipitous match between the in-plane lattice constant of NCCO and the lattice constant of one of the decomposition products, (Nd,Ce)2O3, causes the secondary phase to form an oriented, quasi-two-dimensional epitaxial structure. Consequently, diffraction peaks from the secondary phase appear at rational positions (H,K,0) in the reciprocal space of NCCO. Additionally, because of neodymium paramagnetism, the application of a magnetic field increases the low-temperature intensity observed at these positions via neutron scattering. Such effects may mimic the formation of a structural superlattice or the strengthening of antiferromagnetic order of NCCO, but the intrinsic mechanism may be identified through careful and systematic experimentation. For typical reduction conditions, the (Nd,Ce)2O3 volume fraction is ~1%, and the secondary-phase layers exhibit long-range order parallel to the NCCO CuO2 sheets and are 50-100 angstromsthick. The presence of the secondary phase should also be taken into account in the analysis of other experiments on NCCO, such as transport measurements.Comment: 15 pages, 17 figures, submitted to Phys. Rev.

Path integral Monte Carlo simulations of silicates

We investigate the thermal expansion of crystalline SiO$_2$ in the $\beta$-- cristobalite and the $\beta$-quartz structure with path integral Monte Carlo (PIMC) techniques. This simulation method allows to treat low-temperature quantum effects properly. At temperatures below the Debye temperature, thermal properties obtained with PIMC agree better with experimental results than those obtained with classical Monte Carlo methods.Comment: 27 pages, 10 figures, Phys. Rev. B (in press

Phonon dispersion and electron-phonon coupling in MgB_2 and AlB_2

We present a first principles investigation of the lattice dynamics and electron-phonon coupling of the superconductor MgB_2 and the isostructural AlB_2 within the framework of density functional perturbation theory using a mixed-basis pseudopotential method. Complete phonon dispersion curves and Eliashberg functions \alpha^2F are calculated for both systems. We also report on Raman measurements, which support the theoretical findings. The calculated generalized density-of-states for MgB_2 is in excellent agreement with recent neutron-scattering experiments. The main differences in the calculated phonon spectra and \alpha^2F are related to high frequency in-plane boron vibrations. As compared to AlB_2, they are strongly softened in MgB_2 and exhibit an exceptionally strong coupling to electronic states at the Fermi energy. The total coupling constants are \lambda_{MgB_2}=0.73 and \lambda_{AlB_2}=0.43. Implications for the superconducting transition temperature are briefly discussed.Comment: 10 pages, 4 figures, to appear in Phys. Rev. Let

Assessing the Performance of Recent Density Functionals for Bulk Solids

We assess the performance of recent density functionals for the exchange-correlation energy of a nonmolecular solid, by applying accurate calculations with the GAUSSIAN, BAND, and VASP codes to a test set of 24 solid metals and non-metals. The functionals tested are the modified Perdew-Burke-Ernzerhof generalized gradient approximation (PBEsol GGA), the second-order GGA (SOGGA), and the Armiento-Mattsson 2005 (AM05) GGA. For completeness, we also test more-standard functionals: the local density approximation, the original PBE GGA, and the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA. We find that the recent density functionals for solids reach a high accuracy for bulk properties (lattice constant and bulk modulus). For the cohesive energy, PBE is better than PBEsol overall, as expected, but PBEsol is actually better for the alkali metals and alkali halides. For fair comparison of calculated and experimental results, we consider the zero-point phonon and finite-temperature effects ignored by many workers. We show how Gaussian basis sets and inaccurate experimental reference data may affect the rating of the quality of the functionals. The results show that PBEsol and AM05 perform somewhat differently from each other for alkali metal, alkaline earth metal and alkali halide crystals (where the maximum value of the reduced density gradient is about 2), but perform very similarly for most of the other solids (where it is often about 1). Our explanation for this is consistent with the importance of exchange-correlation nonlocality in regions of core-valence overlap.Comment: 32 pages, single pdf fil