Abstract: Bi2Se3 is a three-dimensional topological insulator which has been extensively studied because it has a single Dirac cone on the surface, inside a relatively large bulk band gap. However, the effect of two-dimensional topological insulator Bi bilayers on the properties of Bi2Se3 and vice versa, has not been explored much. Bi bilayers are often present between the quintuple layers of Bi2Se3, since (Bi2)n(Bi2Se3)m form stable ground-state structures. Moreover, Bi2Se3 is a good substrate for growing ultrathin Bi bilayers. By first-principles techniques, we first show that there is no preferable surface termination by either Bi or Se. Next, we investigate the electronic structure of Bi bilayers on top of, or inside a Bi2Se3 slab. If the Bi bilayers are on top, we observe a charge transfer to the quintuple layers that increases the binding energy of the surface Dirac cones. The extra states, originating from the Bi bilayers, were declared to form a topological Dirac cone, but here we show that these are ordinary Rashba-split states. This result, together with the appearance of a new Dirac cone that is localized slightly deeper, might necessitate the reinterpretation of several experimental results. When the Bi bilayers are located inside the Bi2Se3 slab, they tend to split the slab into two topological insulators with clear surface states. Interface states can also be observed, but an energy gap persists because of strong coupling between the neighboring quintuple layers and the Bi bilayers.

Abstract: Based on first-principles calculations, we study the structural stability of the H and T phases of monolayer MoS2 upon Li doping. Our calculations demonstrate that it is possible to stabilize a distorted T phase of MoS2 over the H phase through adsorption of Li atoms on the MoS2 surface. Through molecular dynamics and phonon calculations, we show that the T phase of MoS2 is dynamically unstable and undergoes considerable distortions. The type of distortion depends on the concentration of adsorbed Li atoms and changes from zigzag-like to diamond-like when increasing the Li doping. There exists a substantial energy barrier to transform the stable H phase to the distorted T phases, which is considerably reduced by increasing the concentration of Li atoms. We show that it is necessary that the Li atoms adsorb on both sides of the MoS2 monolayer to reduce the barrier sufficiently. Two processes are examined that allow for such two-sided adsorption, namely, penetration through the MoS2 layer and diffusion over the MoS2 surface. We show that while there is only a small barrier of 0.24 eV for surface diffusion, the amount of energy needed to pass through a pure MoS2 layer is of the order of similar or equal to 2 eV. However, when the MoS2 layer is covered with Li atoms the amount of energy that Li atoms should gain to penetrate the layer is drastically reduced and penetration becomes feasible.

Abstract: We examine the electronic structure of the family of ternary zinc spinel oxides ZnX2O4 (X=Al, Ga and In). The band gap of ZnAl2O4 calculated using density functional theory (DFT) is 4.25 eV and is overestimated compared with the experimental value of 3.83.9 eV. The DFT band gap of ZnGa2O4 is 2.82 eV and is underestimated compared with the experimental value of 4.45.0 eV. Since DFT typically underestimates the band gap in the oxide system, the experimental measurements for ZnAl2O4 probably require a correction. We use two first-principles techniques capable of describing accurately the excited states of semiconductors, namely the GW approximation and the modified BeckeJohnson (MBJ) potential approximation, to calculate the band gap of ZnX2O4. The GW and MBJ band gaps are in good agreement with each other. In the case of ZnAl2O4, the predicted band gap values are >6 eV, i.e. ~2 eV larger than the only reported experimental value. We expect future experimental work to confirm our results. Our calculations of the electron effective masses and the second band gap indicate that these compounds are very good candidates to act as transparent conducting host materials.

Abstract: ZnAl2O4 (gahnite) is a ceramic which is considered a possible transparent conducting oxide (TCO) due to its wide band gap and transparency for UV. Defects play an important role in controlling the conductivity of a TCO material along with the dopant, which is the main source of conductivity in an otherwise insulating oxide. A comprehensive first-principles density functional theory study for point defects in ZnAl2O4 spinel is presented using the Heyd, Scuseria, and Ernzerhof hybrid functional (HSE06) to overcome the band gap problem. We have investigated the formation energies of intrinsic defects which include the Zn, Al, and O vacancy and the antisite defects: Zn at the Al site (ZnAl) and Al at the Zn site (AlZn). The antisite defect AlZn has the lowest formation energy and acts as a shallow donor, indicating possible n-type conductivity in ZnAl2O4 spinel by Al doping.

Abstract: We present the quasiparticle band structure of ZnO in its zincblende (ZB) and rocksalt (RS) phases at the Γ point, calculated within the GW approximation. The effect of the pd hybridization on the quasiparticle corrections to the band gap is discussed. We compare three systems, ZB-ZnO which shows strong pd hybridization and has a direct band gap, RS-ZnO which is also hybridized but includes inversion symmetry and therefore has an indirect band gap, and ZB-ZnS which shows a weaker hybridization due to a change of the chemical species from oxygen to sulfur. The quasiparticle corrections are calculated with different numbers of valence electrons in the Zn pseudopotential. We find that the Zn20 + pseudopotential is essential for the adequate treatment of the exchange interaction in the self-energy. The calculated GW band gaps are 2.47 eV and 4.27 eV respectively, for the ZB and RS phases. The ZB-ZnO band gap is underestimated compared to the experimental value of 3.27 by ~ 0.8 eV. The RS-ZnO band gap compares well with the experimental value of 4.5 eV. The underestimation for ZB-ZnO is correlated with the strong pd hybridization. The GW band gap for ZnS is 3.57 eV, compared to the experimental value of 3.8 eV.

Abstract: CdO in the rocksalt structure is an indirect band gap semiconductor. Thus, in order to determine its band gap one needs to calculate the complete band structure. However, in practice, the exact evaluation of the quasiparticle band structure for the large number of k-points which constitute the different symmetry lines in the Brillouin zone can be an extremely demanding task compared to the standard density functional theory (DFT) calculation. In this paper we report the full quasiparticle band structure of CdO using a plane-wave pseudopotential approach. In order to reduce the computational effort and time, we make use of maximally localized Wannier functions (MLWFs). The MLWFs offer a highly accurate method for interpolation of the DFT or GW band structure from a coarse k-point mesh in the irreducible Brillouin zone, resulting in a much reduced computational effort. The present paper discusses the technical details of the scheme along with the results obtained for the quasiparticle band gap and the electron effective mass.

Abstract: We theoretically investigate the ground-state properties of a quantum dot defined on the surface of a strong three-dimensional time-reversal invariant topological insulator. Confinement is realized by ferromagnetic barriers and Coulomb interaction is treated numerically for up to seven electrons in the dot. Experimentally relevant intermediate interaction strengths are considered. The topological origin of the dot has several consequences: (i) spin polarization increases and the ground state exhibits quantum phase transitions at specific angular momenta as a function of interaction strength, (ii) the onset of Wigner correlations takes place mainly in one spin channel, and (iii) the ground state is characterized by a robust persistent current that changes sign as a function of the distance from the center of the dot.

Abstract: We consider an interface between two strong time-reversal invariant topological insulators having surface states with opposite spin chirality, or equivalently, opposite mirror Chern number. We show that such an interface supports gapless modes that are protected by mirror symmetry. The interface states are investigated with a continuum model for the Bi2Se3 class of topological insulators that takes into account terms up to third order in the crystal momentum, which ensures that the model has the correct symmetry. The model parameters are obtained from ab initio calculations. Finally, we consider the effect of rotational mismatch at the interface, which breaks the mirror symmetry and opens a gap in the interface spectrum.

Abstract: Computational fluid dynamics (CFD) is increasingly applied in the respiratory domain. The ability to simulate the flow through a bifurcating tubular system has increased the insight into the internal flow dynamics and the particular characteristics of respiratory flows such as secondary motions and inertial effects. The next step in the evolution is to apply the technique to patient-specific cases, in order to provide more information about pathological airways. This study presents a patient-specific approach where both the geometry and the boundary conditions (BC) are based on individual imaging methods using computed tomography (CT). The internal flow distribution of a 73-year-old female suffering from chronic obstructive pulmonary disease (COPD) is assessed. The validation is performed through the comparison of lung ventilation with gamma scintigraphy. The results show that in order to obtain agreement within the accuracy limits of the gamma scintigraphy scan, both the patient-specific geometry and the BC (driving pressure) play a crucial role. A minimal invasive test (CT scan) supplied enough information to perform an accurate CFD analysis. In the end it was possible to capture the pathological features of the respiratory system using the imaging and computational fluid dynamics techniques. This brings the introduction of this new technique in the clinical practice one step closer.

Abstract: Doped Cr2O3 has been shown to be a p-type transparent conducting oxide (TCO). Its conductivity, however, is low. As for most p-type TCOs, the main problem is the high effective hole mass due to flat valence bands. We use first-principles methods to investigate whether one can increase the valence band dispersion (i.e. reduce the hole mass) by anion alloying with sulfur, while keeping the band gap large enough for transparency. The alloying concentrations considered are given by Cr(4)SxO(6-x), with x = 1-5. To be able to describe the electronic properties of these materials accurately, we first study Cr2O3, examining critically the accuracy of different density functionals and methods, including PBE, PBE+U, HSE06, as well as perturbative approaches within the GW approximation. Our results demonstrate that Cr4S2O4 has an optical band gap of 3.08 eV and an effective hole mass of 1.8 m(e). This suggests Cr4S2O4 as a new p-type TCO host candidate.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.

Abstract: We modeled the electron and hole states in Si/SiO2 quantum wells within a basis of standing waves using the 30-band k . p theory. The hard-wall confinement potential is assumed, and the influence of the peculiar band structure of bulk silicon on the quantum-well sub-bands is explored. Numerous spurious solutions in the conduction-band and valence-band energy spectra are found and are identified to be of two types: (1) spurious states which have large contributions of the bulk solutions with large wave vectors (the high-k spurious solutions) and (2) states which originate mainly from the spurious valley outside the Brillouin zone (the extravalley spurious solutions). An algorithm to remove all those nonphysical solutions from the electron and hole energy spectra is proposed. Furthermore, slow and oscillatory convergence of the hole energy levels with the number of basis functions is found and is explained by the peculiar band mixing and the confinement in the considered quantum well. We discovered that assuming the hard-wall potential leads to numerical instability of the hole states computation. Nonetheless, allowing the envelope functions to exponentially decay in a barrier of finite height is found to improve the accuracy of the computed hole states.

Abstract: The valley splitting of the conduction-band states in a thin silicon-on-insulator layer is investigated using the 30-band k . p theory. The system composed of a few nm thick Si layer embedded within thick SiO2 layers is analyzed. The valley split states are found to cross periodically with increasing quantum well width, and therefore the energy splitting is an oscillatory function of the quantum well width, with period determined by the wave vector K-0 of the conduction band minimum. Because the valley split states are classified by parity, the optical transition between the ground hole state and one of those valley split conduction band states is forbidden. The oscillations in the valley splitting energy decrease with electric field and with smoothing of the composition profile between the well and the barrier by diffusion of oxygen from the SiO2 layers to the Si quantum well. Such a smoothing also leads to a decrease of the interband transition matrix elements. The obtained results are well parametrized by the effective two-valley model, but are found to disagree from previous 30-band calculations. This discrepancy could be traced back to the fact that the basis for the numerical solution of the eigenproblem must be restricted to the first Brillouin zone in order to obtain quantitatively correct results for the valley splitting.

Abstract: The DFT-1/2 method is a band-gap correction with GW precision at a density functional theory (DFT) computational cost. The method was also extended to correct the gap between defect levels, allowing for the calculation of optical transitions. However, this method fails when the atomic character of the occupied and unoccupied defect levels is similar as we illustrate by two examples, the tetrahedral hydrogen interstitial and the negatively charged vacancy in diamond. We solve this problem by decoupling the effect of the occupied and unoccupied defect levels and call this the decoupled DFT-1/2 method for defects.