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.

Abstract: The magnetic-field dependence of the energy spectrum, wave function, binding energy, and oscillator strength of exciton states confined in a circular graphene quantum dot (CGQD) is obtained within the configuration interaction method. We predict that (i) excitonic effects are very significant in the CGQD as a consequence of a combination of geometric confinement, magnetic confinement, and reduced screening; (ii) two types of excitons (intravalley and intervalley) are present in the CGQD because of the valley degree of freedom in graphene; (iii) the intravalley and intervalley exciton states display different magnetic-field dependencies due to the different electron-hole symmetries of the single-particle energy spectra; (iv) with increasing magnetic field, the exciton ground state in the CGQD undergoes an intravalley to intervalley transition accompanied by a change of angular momentum; (v) the exciton binding energy does not increase monotonically with the magnetic field due to the competition between geometric and magnetic confinements; and (vi) the optical transitions of the intervalley and intravalley excitons can be tuned by the magnetic field, and valley-dependent excitonic transitions can be realized in a CGQD.

Abstract: The report on square ice sandwiched between two graphene layers by Algara-Siller et al. [Nature (London) 519, 443 (2015)] has generated a large interest in this system. By applying high lateral pressure on nanoconfined water, we found that monolayer ice is transformed to bilayer ice when the two graphene layers are separated by H = 6,7 angstrom. It was also found that three layers of a denser phase of ice with smaller lattice constant are formed if we start from bilayer ice and apply a lateral pressure of about 0.7 GPa with H = 8,9 angstrom. The lattice constant (2.5-2.6 angstrom) in both transitions is found to be smaller than those typical for the known phases of ice and water, i.e., 2.8 angstrom. We validate these results using ab initio calculations and find good agreement between ab initio O-O distance and those obtained from classical molecular dynamics simulations. The reversibility of the mentioned transitions is confirmed by decompressing the systems.

Abstract: We report on a first-principles study of the structural and electronic properties of a Sigma3 (112) grain boundary model in CuInSe2. The study focuses on a coherent, stoichiometry preserving, cation–Se terminated grain boundary, addressing the properties of the grain boundary as such, as well as the effect
of well known defects in CuInSe2. We show that in spite of its apparent simplicity, such a grain boundary exhibits a very rich phenomenology, providing an explanation for several of the experimentally observed properties of grain boundaries in CuInSe2 thin films. In particular, we show that the combined effect of Cu vacancies and cation antisites can result in the observed Cu depletion with no In enrichment at the grain boundaries. Furthermore, Cu vacancies are unlikely to produce a hole barrier at the grain boundaries, but Na may indeed have such an effect. We find that Na-on-Cu defects will tend to form abundantly at
the grain boundaries, and can provide a mechanism for the carrier depletion and/or type inversion experimentally reported.

Abstract: The effect of B and P dopants on the band structure of Si nanowires is studied using electronic structure calculations based on density functional theory. At low concentrations a dispersionless band is formed, clearly distinguishable from the valence and conduction bands. Although this band is evidently induced by the dopant impurity, it turns out to have purely Si character. These results can be rigorously analyzed in the framework of effective mass theory. In the process we resolve some common misconceptions about the physics of hydrogenic shallow impurities, which can be more clearly elucidated in the case of nanowires than would be possible for bulk Si. We also show the importance of correctly describing the effect of dielectric confinement, which is not included in traditional electronic structure calculations, by comparing the obtained results with those of G(0)W(0) calculations.

Abstract: Two-dimensional (2D) carbon nitride materials play an important role in energy-harvesting, energy-storage and environmental applications. Recently, a new carbon nitride, 2D polyaniline (C3N) was proposed [PNAS 113 (2016) 7414-7419]. Based on the structure model of this C3N monolayer, we propose two new carbon nitride monolayers, named dumbbell (DB) C4N-I and C4N-II. Using first-principles calculations, we systematically study the structure, stability, and band structure of these two materials. In contrast to other carbon nitride monolayers, the orbital hybridization of the C/N atoms in the DB C4N monolayers is sp(3). Remarkably, the band structures of the two DB C4N monolayers have a Dirac cone at the K point and their Fermi velocities (2.6/2.4 x 10(5) m/s) are comparable to that of graphene. This makes them promising materials for applications in high-speed electronic devices. Using a tight-binding model, we explain the origin of the Dirac cone. (C) 2017 Elsevier Ltd. All rights reserved.