Abstract: In the quantum mechanical simulation of exploratory semiconductor devices, continuum methods based on a k.p/envelope function model have the potential to significantly reduce the computational burden compared to prevalent atomistic methods. However, full-zone k.p/envelope function simulation approaches are scarce and existing implementations are not self-consistent with the calculation of the electrostatic potential due to the lack of a stable procedure and a proper normalization of the multi-band envelope functions. Here, we therefore present a self-consistent procedure based on a full-zone spectral k.p/envelope function band structure model. First, we develop a proper normalization for the multi-band envelope functions in the presence of transmitting boundary conditions. This enables the calculation of the free carrier densities. Next, we construct a procedure to obtain self-consistency of the carrier densities with the electrostatic potential. This procedure is stabilized with an adaptive scheme that relies on the solution of Poisson's equation in the Gummel form, combined with successive underrelaxation. Finally, we apply our procedure to homostructure In0.53Ga0.47As tunnel field-effect transistors (TFETs) and staggered heterostructure GaAs0.5Sb0.5/In0.53Ga0.47As TFETs and show the importance of self-consistency on the device predictions for scaled dimensions. Published by AIP Publishing.

Abstract: Structural and superconducting properties of high-quality niobium nanofilms with different thicknesses are investigated on silicon oxide (SiO2) and sapphire substrates. The role played by the different substrates and the superconducting properties of the Nb films are discussed based on the defectivity of the films and on the presence of an interfacial oxide layer between the Nb film and the substrate. The X-ray absorption spectroscopy is employed to uncover the structure of the interfacial layer. We show that this interfacial layer leads to a strong proximity effect, especially in films deposited on a SiO2 substrate, altering the superconducting properties of the Nb films. Our results establish that the critical temperature is determined by an interplay between quantum-size effects, due to the reduction of the Nb film thicknesses, and proximity effects. The detailed investigation here provides reference characterizations and has direct and important implications for the fabrication of superconducting devices based on Nb nanofilms.

Abstract: We model transport through two-dimensional topological insulator (TI) nanoribbons. To model the quantum transport, we employ the non-equilibrium Green's function approach. With the presented approach, we study the effect of lattice imperfections on the carrier transport. We observe that the topologically protected edge states of TIs are robust against a high percentage (2%) of vacancy defects. We also investigate tunneling of the edge states in two decoupled TI nanoribbons.

Abstract: In this study, we predicted two new stable metallic Re-C based monolayer structures with a rectangular (r-ReC2) and a hexagonal (h-Re2C) crystal symmetry using first-principle calculations based on density functional theory. Our results obtained from mechanical and phonon calculations and high-temperature molecular dynamic simulations clearly proved the stability of these two-dimensional (2D) crystals. Interestingly, Re-C monolayers in common transition metal carbide structures (i.e. MXenes) were found to be unstable, contrary to expectations. We found that the stable structures, i.e. r-ReC2 and h-Re2C, display superior mechanical properties over the well-known 2D materials. The Young's modulus for r-ReC2 and h-Re2C are extremely high and were calculated as 351 (1310) and 617 (804) N/m (GPa), respectively. Both materials have larger Young's modulus values than the most of the well-known 2D materials. We showed that the combination of the short strong directional p-d bonds, the high coordination number of atoms in the unit-cell and high valence electron density result in strong mechanical properties. Due to its crystal structure, the r-ReC2 monolayer has anisotropic mechanical properties and the crystallographic direction parallel to the C-2 dimers is stiffer compared to perpendicular direction due to strong covalent bonding within C-2 dimers. h-Re2C was derived from the corresponding bulk structure for which we determined the critical thickness for the dynamically stable bulk-derived monolayer structures. In addition, we also investigated the electronic of these two stable structures. Both exhibit metallic behavior and Re-5d orbitals dominate the states around the Fermi level. Due to their ultra high mechanical stability and stiffness, these novel Re-C monolayers can be exploited in various engineering applications.

Abstract: Superconductivity is one of the most important discoveries of the last century. With many applications in physics, engineering, and technology, superconductors are crucial to our way of living. Several material and engineering issues however prevent their widespread usage in everyday life. Comprehensive studies are being directed at these materials and their properties to come up with new technologies that will address these challenges and enhance their superconductive capabilities. In this context, numerical modeling plays an important role in the search of new solutions to existing material and engineering issues. The time-dependent Ginzburg-Landau (TDGL) theory is a powerful predictive tool for modeling the macroscopic behavior of superconductors. However most of the numerical algorithms developed so far are incapable of describing many basic properties of real superconducting devices, and are too slow on current hardware for large-scale numerical simulations necessary for their accurate description. Therefore, the purpose of this thesis is to develop high-performing numerical solutions that can correctly describe material features to be used as modeling tools of laboratory experiments. Some important innovations introduced in this work include the numerical modeling of nonrectangular geometrical shapes with complex electrical and insulating components, the inclusion of dynamic heating of the material, and the description of different types of material inhomogeneities. These encompass the principal features necessary for a complete description of the superconductive physics in real material samples. In this thesis a numerical solution is developed for modeling superconducting thin films and used to study the superconductive properties of three experimental configurations: the dynamics of vortex matter in a Corbino disk, the motion of ultrafast vortices in an hourglass-shaped microbridge, and the photon detection process in a meander-patterned nanowire. Moreover, a numerical solution is developed for modeling three-dimensional superconductors which are studied here for the first time in the type-I superconducting regime. These numerical algorithms are optimized to exploit the computational horsepower of graphics processing units (GPUs) and multicore central-processing unit (CPU) clusters such that they can achieve high-performance and be used to model large-scale problems previously impossible on conventional machines. Several computational tools are also designed to assist with the modeling of superconducting devices. These include a numerical library of the TDGL equations, a novel mechanism for the generation of complex geometries, a closed-form solver to conduct numerical simulations, and a graphics user interface (GUI) to visualize the dynamic behavior of superconductors. The contributions in this thesis ultimately push the boundaries on what is possible in state-of-the-art numerical modeling of superconductivity.

Abstract: Moore’s Law has driven the electronics industry for the past half century. However, the doubling of transistors about every two years is beginning to break down, owing to fundamental limits that arise as they approach the atomic length. As a result, the search for new pathways for electronics has become crucial. Among potential candidates, the discovery of magnetic textures known as skyrmions has attracted considerable interest and attention in spintronic technology, which relies on both the electron charge and its spin. The unusual topological and particle-like behavior launched skyrmions into the spotlight of scientific research. Topological protected stability, nanoscale size, and low driving currents needed to move them make skyrmions promising candidates for future consumer nanoelectronics. Recent advances in the field have provided all of the basic functions needed for carrying and processing information. In this thesis, we procure to advance the current understanding of skyrmion physics, and explore their potential to replace conventional electronics technology. First, the fundamental properties and lifetimes of racetrack skyrmions at room temperature are investigated. We discover that skyrmions can easily collapse at the boundary in laterally finite systems, and propose ways to improve their stability for constrained geometries. Then, pinning of single skyrmions on atomic defects of distinct origins are studied. We reveal that the preferred pinning positions depend on the skyrmion size and type of defect being considered, and discuss applications where control of skyrmions by defects is of particular interest. Next, we explore other magnetic configurations that can compete with skyrmions when considering new materials, and describe a previously unseen mechanism for collapse of skyrmions into cycloidal spin backgrounds. Finally, switching and interactions between skyrmions with distinct topologies are reported. We find that skyrmions transition to higher or lower topologies by absorbing a unit spin texture. The interactions between skyrmions of different topological charges can be attractive or repulsive, leading to the formation of arranged clusters. We conclude with a numerical library for simulating magnetic skyrmions in various scenarios.

Abstract: The newest members of a two-dimensional material family, involving transition metal carbides and nitrides (called MXenes), have garnered increasing attention due to their tunable electronic and thermal properties depending on the chemical composition and functionalization. This flexibility can be exploited to fabricate efficient electrochemical energy storage (batteries) and energy conversion (thermoelectric) devices. In this study, we calculated the Seebeck coefficients and lattice thermal conductivity values of oxygen terminated M2CO2 (where M = Ti, Zr, Hf, Sc) monolayer MXene crystals in two different functionalization configurations (model-II (MD-II) and model-III (MD-III)), using density functional theory and Boltzmann transport theory. We estimated the thermoelectric figure-of-merit, zT, of these materials by two different approaches, as well. First of all, we found that the structural model (i.e. adsorption site of oxygen atom on the surface of MXene) has a paramount impact on the electronic and thermoelectric properties of MXene crystals, which can be exploited to engineer the thermoelectric properties of these materials. The lattice thermal conductivity kappa(l), Seebeck coefficient and zT values may vary by 40% depending on the structural model. The MD-III configuration always has the larger band gap, Seebeck coefficient and zT, and smaller kappa(l) as compared to the MD-II structure due to a larger band gap, highly flat valence band and reduced crystal symmetry in the former. The MD-III configuration of Ti2CO2 and Zr2CO2 has the lowest kappa(l) as compared to the same configuration of Hf2CO2 and Sc2CO2. Among all the considered structures, the MD-II configuration of Hf2CO2 has the highest kappa(l), and Ti2CO2 and Zr2CO2 in the MD-III configuration have the lowest kappa(l). For instance, while the band gap of the MD-II configuration of Ti2CO2 is 0.26 eV, it becomes 0.69 eV in MD-III. The zT(max) value may reach up to 1.1 depending on the structural model of MXene.

Abstract: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom-thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.

Abstract: We investigate the role of transition metal atoms of group V-b (V, Nb, Ta) and VI-b (Cr, Mo, W) as n- or p-type dopants in anatase TiO$_2$ using thermodynamic
principles and density functional theory with the Heyd-Scuseria-Ernzerhof HSE06 hybrid functional. The HSE06 functional provides a realistic value for the band gap, which ensures a correct classification of dopants as shallow or deep donors or acceptors. Defect formation energies and thermodynamic transition levels are calculated taking into account the constraints imposed by the stability of TiO$_2$ and the solubility limit of the impurities.
Nb, Ta, W and Mo are identified as shallow donors. Although W provides two electrons, Nb and Ta show a considerable lower formation energy, in particular under O-poor conditions. Mo donates in principle one electron, but under specific conditions can turn into a double donor. V impurities are deep donors and Cr
shows up as an amphoteric defect, thereby acting as an electron trapping center in n-type TiO$_2$ especially under O-rich conditions. A comparison with the available experimental data yields excellent agreement.

Abstract: Recent studies showed that positron annihilation methods can provide key insights into the nanostructure and electronic structure of thin film solar cells. In this study, positron annihilation lifetime spectroscopy (PALS) is applied to investigate CdSe quantum dot (QD) light absorbing layers, providing evidence of positron trapping at the surfaces of the QDs. This enables one to monitor their surface composition and electronic structure. Further, 2D-Angular Correlation of Annihilation Radiation (2D-ACAR) is used to investigate the nanostructure of divacancies in photovoltaic-high-quality a-Si:H films. The collected momentum distributions were converted by Fourier transformation to the direct space representation of the electron-positron autocorrelation function. The evolution of the size of the divacancies as a function of hydrogen dilution during deposition of a-Si:H thin films was examined. Finally, we present a first positron Doppler Broadening of Annihilation Radiation (DBAR) study of the emerging class of highly efficient thin film solar cells based on perovskites.

Abstract: A modeling approach, based on an analytical solution of the semiclassical multi-subband Boltzmann transport equation, is presented to study resistivity scaling in metallic thin films and nanowires due to grain boundary and surface roughness scattering. While taking into account the detailed statistical properties of grains, roughness and barrier material as well as the metallic band structure and quantum mechanical aspects of scattering and confinement, the model does not rely on phenomenological fitting parameters. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: The Aharonov-Bohm (AB) effect in square phosphorene quantum rings, with armchair and zigzag edges, is investigated using the tight-binding method. The energy spectra and wave functions of such rings, obtained as a function of the magnetic flux Phi threading the ring, are strongly influenced by the ringwidthW, an in-plane electric field E-p, and a side-gating potential V-g. Compared to a square dot, the ring shows an enhanced confinement due to its inner edges and an interedge coupling along the zigzag direction, both of which strongly affect the energy spectrum and the wave functions. The energy spectrum that is gapped consists of a regular part, of conduction (valence) band states, that shows the usual AB oscillations in the higher-(lower-) energy region, and of edge states, in the gap, that exhibit no AB oscillations. As the width W decreases, the AB oscillations become more distinct and regular and their period is close to Phi(0)/2, where the flux quantum Phi(0) = h/e is the period of an ideal circular ring (W -> 0). Both the electric field E-p and the side-gating potential V-g reduce the amplitude of the AB oscillations. The amplitude can be effectively tuned by E-p or V-g and exhibits an anisotropic behavior for different field directions or side-gating configurations.

Abstract: The Dzyaloshinskii-Moriya interaction (DMI) is a chiral interaction that favors formation of domain walls. Recent experiments and ab initio calculations show that there are multiple ways to modify the strength of the interfacially induced DMI in thin ferromagnetic films with perpendicular magnetic anisotropy. In this paper we reveal theoretically the effects of spatially varied DMI on the magnetic state in thin films. In such heterochiral 2D structures we report several emergent phenomena, ranging from the equilibrium spin canting at the interface between regions with different DMI, over particularly strong confinement of domain walls and skyrmions within high-DMI tracks, to advanced applications such as domain tailoring nearly at will, design of magnonic waveguides, and much improved skyrmion racetrack memory.