No Arabic abstract
An extensive theoretical investigation of the nonpolar (10$bar{1}$0) and (11$bar{2}$0) surfaces as well as the polar zinc terminated (0001)--Zn and oxygen terminated (000$bar{1}$)--O surfaces of ZnO is presented. Particular attention is given to the convergence properties of various parameters such as basis set, k--point mesh, slab thickness, or relaxation constraints within LDA and PBE pseudopotential calculations using both plane wave and mixed basis sets. The pros and cons of different approaches to deal with the stability problem of the polar surfaces are discussed. Reliable results for the structural relaxations and the energetics of these surfaces are presented and compared to previous theoretical and experimental data, which are also concisely reviewed and commented.
The study of zinc oxide, within the homogeneous electron gas approximation, results in overhybridization of zinc $3d$ shell with oxygen $2p$ shell, a problem shown for most transition metal chalcogenides. This problem can be partially overcome by using LDA+$U$ (or, GGA+$U$) methodology. However, in contrast to the zinc $3d$ orbital, Hubbard type correction is typically excluded for the oxygen $2p$ orbital. In this work, we provide results of electronic structure calculations of an oxygen vacancy in ZnO supercell from ab initio perspective, with two Hubbard type corrections, $U_{mathrm{Zn}-3d}$ and $U_{mathrm{O}-2p}$. The results of our numerical simulations clearly reveal that the account of $U_{mathrm{O}-2p}$ has a significant impact on the properties of bulk ZnO, in particular the relaxed lattice constants, effective mass of charge carriers as well as the bandgap. For a set of validated values of $U_{mathrm{Zn}-3d}$ and $U_{mathrm{O}-2p}$ we demonstrate the appearance of a localized state associated with the oxygen vacancy positioned in the bandgap of the ZnO supercell. Our numerical findings suggest that the defect state is characterized by the highest overlap with the conduction band states as obtained in the calculations with no Hubbard-type correction included. We argue that the electronic density of the defect state is primarily determined by Zn atoms closest to the vacancy.
Solid-liquid interfaces are at the heart of many modern-day technologies and provide a challenge to many materials simulation methods. A realistic first-principles computational study of such systems entails the inclusion of solvent effects. In this work we implement an implicit solvation model that has a firm theoretical foundation into the widely used density-functional code VASP. The implicit solvation model follows the framework of joint density functional theory. We describe the framework, our algorithm and implementation, and benchmarks for small molecular systems. We apply the solvation model to study the surface energies of different facets of semiconducting and metallic nanocrystals and the S$_{text{N}} 2$ reaction pathway. We find that solvation reduces the surface energies of the nanocrystals, especially for the semiconducting ones and increases the energy barrier of the S$_{text{N}} 2$ reaction.
In this work we investigate the electronic properties of mercaptocarboxylic acids with several carbon chain lengths adsorbed on ZnO-(10-10) surfaces via density functional theory calculations using semi-local and hybrid exchange-correlation functionals. Amongst the investigated structures, we identify the monodentate adsorption mode to be stable. Moreover, this mode introduces optically active states in the ZnO gap, is further confirmed by the calculation of the dielectric function at PBE0 and TD-PBE0 levels. One interesting finding is that adsorption mode and the dielectric properties of the hybrid system are both rather insensitive to the chain length, since the acceptor molecular state is very localized on the sulphur atom. This indicates that even small molecules can be used to stabilize ZnO surface and to enhance its functionality for opto-electronic applications.
The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, we study the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). We were able to perform large-scale first principles simulations and capture effects of different surface reconstructions by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, we implemented a bulk unfolding scheme by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, we find a good agreement between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. Our combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, e.g., topological qubits utilizing the Majorana zero modes.
The interaction of tungsten hexacarbonyl W(CO)$_6$ precursor molecules with SiO$_2$ substrates is investigated by means of density functional theory calculations with and without inclusion of long range van der Waals interactions. We consider two different surface models, a fully hydroxylated and a partially hydroxylated SiO$_2$ surface, corresponding to substrates under different experimental conditions. For the fully hydroxylated surface we observe only a weak interaction between the precursor molecule and the substrate with physisorption of W(CO)$_6$. Inclusion of van der Waals corrections results in a stabilization of the molecules on this surface, but does not lead to significant changes in the chemical bonding. In contrast, we find a spontaneous dissociation of the precursor molecule on the partially hydroxylated SiO$_2$ surface where chemisorption of a W(CO)$_5$ fragment is observed upon removal of one of the CO ligands from the precursor molecule. Irrespective of the hydroxylation, the precursor molecule prefers binding of more than one of its CO ligands. In the light of these results, implications for the initial growth stage of tungsten nano-deposits on SiO$_2$ in an electron beam induced deposition process are discussed.