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: When slow ions incident on a surface are neutralized, the excess potential energy is passed on to an electron inside the surface, leading to emission of secondary electrons. The microscopic description of this process, as
well as the calculation of the secondary electron yield, is a challenging problem due to its complexity as well
as its sensitivity to surface properties. One of the first quantitative descriptions was articulated in the 1950s by
Hagstrum, who based his calculation on a parametrization of the density of states of the material. In this paper, we
present a model for calculating the secondary electron yield, derived from Hagstrum’s initial approach. We use
first-principles density functional theory calculations to acquire the necessary input and introduce the concept of
electron cascades to Hagstrum’s model in order to improve the calculated spectra, as well as remove its reliance
on fitting parameters. We apply our model to He+ and Ne+ ions incident on Ge(111) and Si(111) and obtain
yield spectra that match closely to the experimental results of Hagstrum.

Abstract: Lithium-rich layered oxides such as Li₂MnO₃have shown great potential as cathodes in Li-ion batteries, mainly because of their large capacities. However, these materials still suffer from structural degradation as the battery is cycled, reducing the average voltage and capacity of the cell. The voltage fade is believed to be related to the migration of transition metals into the lithium layer, linked to the formation of O-O dimers with a short bond length, which in turn is driven by the presence of oxygen holes due to the participation of oxygen in the redox process. We investigate the formation of O-O dimers for partially charged O1-Li₂MnO₃using a first-principles density functional theory approach by calculating the reaction energy and kinetic barriers for dimer formation. Next, we perform similar calculations for partially charged O1-Li₂IrO₃, a Li-rich material for which the voltage fade was not observed during cycling. When we compare the stability of the oxygen framework, we conclude that the formation of O-O dimers is both thermodynamically and kinetically viable for O1-Li_0.5MnO₃. For O1-Li_0.5IrO₃, we observe that the oxygen lattice is much more stable, either returning to its original state when perturbed, or resulting in a structure with an O-O dimer that is much higher in energy. This can be explained by the mixed redox process for Li₂IrO₃, which is also shown from the calculated magnetic moments. The lack of O-O dimer formation in O1-Li_0.5IrO₃provides valuable insight as to why Li₂IrO₃does not demonstrate a voltage fade as the battery is cycled, which can be used to design Li-rich battery cathodes with an improved cycling performance.

Abstract: The adsorption and diffusion of Na, K, and Ca atoms on MXene/graphene heterostructures of MXene systems Sc2C(OH)(2), Ti2CO2, and V2CO2 are systematically investigated by using first-principles methods. We found that alkali metal intercalation is energetically favorable and thermally stable for Ti2CO2/graphene and V2CO2/graphene heterostructures but not for Sc2C(OH)(2). Diffusion kinetics calculations showed the advantage of MXene/graphene heterostructures over sole MXene systems as the energy barriers are halved for the considered alkali metals. Low energy barriers are found for Na and K ions, which are promising for fast charge/discharge rates. Calculated voltage profiles reveal that estimated high capacities can be fully achieved for Na ion in V2CO2/graphene and Ti2CO2/graphene heterostructures. Our results indicate that Ti2CO2/graphene and V2CO2/graphene electrode materials are very promising for Na ion battery applications. The former could be exploited for low voltage applications while the latter will be more appropriate for higher voltages.

Abstract: Using the non-equilibrium Green's function formalism, we study carrier transport through imperfect two-dimensional (2D) topological insulator (TI) ribbons. In particular, we investigate the effect of vacancy defects on the carrier transport in 2D TI ribbons with hexagonal lattice structure. To account for the random distribution of the vacancy defects, we present a statistical study of varying defect densities by stochastically sampling different defect configurations. We demonstrate that the topological edge states of TI ribbons are fairly robust against a high concentration (up to 2%) of defects. At very high defect densities, we observe an increased inter-edge interaction, mediated by the localisation of the edge states within the bulk region. This effect causes significant back-scattering of the, otherwise protected, edge-states at very high defect concentrations (>2%), resulting in a loss of conduction through the TI ribbon. We discuss how this coherent vacancy scattering can be used to our advantage for the development of TI-based transistors. We find that there is an optimal concentration of vacancies yielding an ON-OFF current ratio of up to two orders of magnitude. Finally, we investigate the importance of spin-orbit coupling on the robustness of the edge states in the TI ribbon and show that increased spin-orbit coupling could further increase the ON-OFF ratio.

Abstract: The classical-map hypernetted-chain (CHNC) scheme, developed for treating fermion fluids at strong coupling and at finite temperatures, is applied to electron-electron and electron-hole double quantum wells. The pair-distribution functions and the local field factors needed in linear-response theory are determined for a range of temperatures, carrier densities, and barrier widths typical for experimental double-quantum-well systems in GaAs-GaAlAs. For electron-hole double quantum wells, a large enhancement in the pair-distribution functions is found for small carrier separations. The CHNC equations for electron-hole systems no longer hold at low densities where bound-state formation occurs.

Abstract: Devising magnetic nanostructures with spatially heterogeneous Dzyaloshinskii-Moriya interaction (DMI) is a promising pathway toward advanced confinement and control of magnetic skyrmions in potential devices. Here we discuss theoretically how a skyrmion interacts with a heterochiral interface using micromagnetic simulations and analytic arguments. We show that a heterochiral interface deflects the trajectory of ferromagnetic (FM) skyrmions, and that the extent of such deflection is tuned by the applied spin-polarized current and the difference in DMI across the interface. Further, we show that this deflection is characteristic of the FM skyrmion, and it is completely absent in the antiferromagnetic (AFM) case. In turn, we reveal that the AFM skyrmion achieves much higher velocities than its FM counterpart, yet experiences far stronger confinement in nanoengineered heterochiral tracks, which reinforces AFM skyrmions as a favorable choice for skyrmion-based devices.

Abstract: Understanding the appearance of local electrical connectivity in liquid filled layered graphene oxide subjected to an external electric field is important to design electrically controlled smart permeable devices and also to gain insight into the physics behind electrical effects on confined water permeation. Motivated by recent experiments [K. G. Zhou et al. Nature (London) 559, 236 (2018)], we introduce a new model with random percolating paths for electrical connectivity in micron thick water filled layered graphene oxide, which mimics parallel resistors connected across the top and bottom electrodes. We find that a strong nonuniform radial electric field of the order similar to 10-50 mV/nm can be induced between layers depending on the current flow through the formed conducting paths. The maxima of the induced fields are not necessarily close to the electrodes and may be localized in the middle region of the layered material. The emergence of electrical connectivity and the associated electrical effects have a strong influence on the surrounding fluid in terms of ionization and wetting which subsequently determines the permeation properties.

Abstract: Reducing the size of a superconductor below its characteristic length scales can either enhance or suppress its critical temperature (T-c). Depending on the bulk value of the electron-phonon coupling strength, electronic and phonon confinement effects will play different roles in the modification of T-c. Experimentally disentangling each contribution has remained a challenge. We have measured both the phonon density of states and T-c of Sn nanowires with diameters of 18, 35, and 100 nm in order to quantify the effects of phonon confinement on superconductivity. We observe a shift of the phonon frequency towards the low-energy region and an increase in the electron-phonon coupling constant that can account for the measured increase in T-c.