No Arabic abstract
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.
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 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.
We propose a method to decompose the total energy of a supercell containing defects into contributions of individual atoms, using the energy density formalism within density functional theory. The spatial energy density is unique up to a gauge transformation, and we show that unique atomic energies can be calculated by integrating over Bader and charge-neutral volumes for each atom. Numerically, we implement the energy density method in the framework of the Vienna ab initio simulation package (VASP) for both norm-conserving and ultrasoft pseudopotentials and the projector augmented wave method, and use a weighted integration algorithm to integrate the volumes. The surface energies and point defect energies can be calculated by integrating the energy density over the surface region and the defect region, respectively. We compute energies for several surfaces and defects: the (110) surface energy of GaAs, the mono-vacancy formation energies of Si, the (100) surface energy of Au, and the interstitial formation energy of O in the hexagonal close-packed Ti crystal. The surface and defect energies calculated using our method agree with size-converged calculations of the difference between the total energies of the system with and without the defect. Moreover, the convergence of the defect energies with size can be found from a single calculation.
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.
We demonstrate the accuracy of the hybrid functional HSE06 for computing band offsets of semiconductor alloy heterostructures. The highlight of this study is the computation of conduction band offsets with a reliability that has eluded standard density functional theory. A high-quality special quasirandom structure models an infinite random pseudobinary alloy for constructing heterostructures along the (001) growth direction. Our excellent results for a variety of heterostructures establish HSE06s relevance to band engineering of high-performance electrical and optoelectronic devices.