Abstract: The quantum electrostatic heterostructure (QEH) model allows for efficient computation of the dielectric screening properties of layered van der Waals (vdW)-bonded heterostructures in terms of the dielectric functions of the individual two-dimensional (2D) layers. Here, we extend the QEH model by including (1) contributions to the dielectric function from infrared active phonons in the 2D layers, (2) screening from homogeneous bulk substrates, and (3) intraband screening from free carriers in doped 2D semiconductor layers. We demonstrate the potential of the extended QEH model by calculating the dispersion of coupled phonons in multilayer stacks of hexagonal boron-nitride (hBN), the strong hybridization of plasmons and optical phonons in graphene/hBN heterostructures, the effect of substrate screening on the exciton series of monolayer MoS2, and the properties of hyperbolic plasmons in a doped phosphorene sheet. The new QEH code is distributed as a Python package with a simple command line interface and a comprehensive library of dielectric building blocks for the most common 2D materials, providing an efficient open platform for dielectric modeling of realistic vdW heterostructures.

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: Electrically gating of trilayer graphene (TLG) opens a band gap offering the possibility to electrically engineer TLG quantum dots. We study the energy levels of such quantum dots and investigate their dependence on a perpendicular magnetic field B and different types of stacking of the graphene layers. The dots are modeled as circular and confined by a truncated parabolic potential which can be realized by nanostructured gates or position-dependent doping. The energy spectra exhibit the intervalley symmetry E-K(e) (m) = -E (h)(K') (m) for the electron (e) and hole (h) states, where m is the angular momentum quantum number and K and K' label the two valleys. The electron and hole spectra for B = 0 are twofold degenerate due to the intervalley symmetry E-K (m) = E-K' [-(m + 1)]. For both ABC [alpha = 1.5 (1.2) for large (small) R] and ABA (alpha = 1) stackings, the lowest-energy levels show approximately a R-alpha dependence on the dot radius R in contrast with the 1/R-3 one for ABC-stacked dots with infinite-mass boundary. As functions of the field B, the oscillator strengths for dipole-allowed transitions differ drastically for the two types of stackings.

Abstract: External fields are a powerful tool to probe optical excitations in a material. The linear energy shift of an excitation in a magnetic field is quantified by its effective g factor. Here we show how exciton g factors and their sign can be determined by converged first-principles calculations. We apply the method to monolayer excitons in semiconducting transition metal dichalcogenides and to interlayer excitons in MoSe2/WSe2 heterobilayers and obtain good agreement with recent experimental data. The precision of our method allows us to assign measured g factors of optical peaks to specific transitions in the band structure and also to specific regions of the samples. This revealed the nature of various, previously measured interlayer exciton peaks. We further show that, due to specific optical selection rules, g factors in van der Waals heterostructures are strongly spin and stacking-dependent. The calculation of orbital angular momenta requires the summation over hundreds of bands, indicating that for the considered two-dimensional materials the basis set size is a critical numerical issue. The presented approach can potentially be applied to a wide variety of semiconductors.

Abstract: The rational design of two-dimensional (2D) piezoelectric materials has recently garnered great interest due to their increasing use in technological applications, including sensor technology, actuating devices, energy harvesting, and medical applications. Several materials possessing high piezoelectric response have been reported so far, but a high-throughput first-principles approach to estimate the piezoelectric potential of layered materials has not been performed yet. In this study, we systematically investigated the piezoelectric (e(11), d(11)) and elastic (C-11 and C-12) properties of 128 thermodynamically stable 2D semiconductor materials by employing first-principle methods. Our high-throughput approach demonstrates that the materials containing Group-V elements produce significantly high piezoelectric strain constants, d(11) > 40 pm V-1, and 49 of the materials considered have the e(11) coefficient higher than MoS2 insomuch as BrSSb has one of the largest d(11) with a value of 373.0 pm V-1. Moreover, we established a simple empirical model in order to estimate the d(11) coefficients by utilizing the relative ionic motion in the unit cell and the polarizability of the individual elements in the compounds.

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: Efficient thermal transport control is a fundamental issue for electronic device applications such as information, communication, and energy storage technologies in modern electronics in order to achieve desired thermal conditions. Structural defects in materials provide a mechanism to adjust the thermal transport properties of these materials on demand. In this context, the effect of structural defects on lattice thermal conductivities of two-dimensional hexagonal binary group-III nitride (XN, X = B, Al, and Ga) semiconductors is systematically investigated by means of classical molecular dynamics simulations performed with recently developed transferable inter-atomic potentials accurately describing defect energies. Here, two different Green-Kubo based approaches and another approach based on non-equilibrium molecular dynamics are compared in order to get an overall understanding. Our investigation clearly shows that defect concentrations of 3% decrease the thermal conductivity of systems containing these nitrites up to 95%. Results hint that structural defects can be used as effective adjustment parameters in controlling thermal transport properties in device applications associated with these materials. Published under an exclusive license by AIP Publishing.

Abstract: We demonstrate the possibility of engineering the optical properties of transition metal dichalcogenide heterobilayers when one of the constitutive layers has a Janus structure. We investigate different MoS2@Janus layer combinations using first-principles methods including excitons and exciton-phonon coupling. The direction of the intrinsic electric field from the Janus layer modifies the electronic band alignments and, consequently, the energy separation between dark interlayer exciton states and bright in-plane excitons. We find that in-plane lattice vibrations strongly couple the two states, so that exciton-phonon scattering may be a viable generation mechanism for interlayer excitons upon light absorption. In particular, in the case of MoS2@WSSe, the energy separation of the low-lying interlayer exciton from the in-plane exciton is resonant with the transverse optical phonon modes (40 meV). We thus identify this heterobilayer as a prime candidate for efficient generation of charge-separated electron-hole pairs.

Abstract: The structural stability of 2D transition metal dichalcogenide (TMD) formations is of particular importance for their reliable device performance in nano-electronics and opto-electronics. Recent observations show that the CVD-grown TMD monolayers are likely to encounter stability problems such as cracking or fracturing when they are kept under ambient conditions. Here, two different growth configurations are investigated and a favorable growth geometry is proposed, which also sheds light onto the growth mechanism and provides a solution for the stability and fracture formation issues for TMDs specifically for MoS2 monolayers. It is shown that 18 months naturally and thermally aged MoS2 monolayer flakes grown using specifically developed conditions, retain their stability. To understand the mechanism of the structural deterioration, two possible effective mechanisms, S vacancy defects and growth-induced tensile stress, are assessed by the first principle calculations where the role of S vacancy defects in obtaining oxidation resistant MoS2 monolayer flakes is revealed to be rather more critical. Hence, these simulations, time-dependent observations and thermal aging experiments show that durability and stability of 2D MoS2 flakes can be controlled by CVD growth configuration.

Abstract: The influence of the nonparabolicity of charge carriers dispersion law (Kane's dispersion) on a magnetoexciton energy spectrum in InSb quantum rings is theoretically investigated The analytical expression for the energy spectrum of exciton in a narrow-gap semiconductor nanoring in a magnetic field is obtained. The Aharonov – Bohm oscillations in the energy of excited states are studied.

Abstract: The structure and thermal boundary conductance of the wurtzite GaN/AlN (0001) interface are investigated using molecular dynamics simulation. Simulation results with three different empirical interatomic potentials have produced similar misfit dislocation networks and dislocation core structures. Specifically, the misfit dislocation network at the GaN/AlN interface is found to consist of pure edge dislocations with a Burgers vector of 1/3(1 (2) over bar 10) and the misfit dislocation core has an eight-atom ring structure. Although different interatomic potentials lead to different dislocation properties and thermal conductance values, all have demonstrated a significant effect of misfit dislocations on the thermal boundary conductance of the GaN/AlN (0001) interface. Published under an exclusive license by AIP Publishing.

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 theoretical investigation on structural, vibrational, and electronic properties of zinc-blende (ZB) AgI were carried out employing first principles density functional theory calculations. Thermoelectric properties then were predicted through semi-classical Boltzmann transport equations within the constant relaxation time approximation. Equilibrium lattice parameter, bulk modulus, elastic constants, and vibrational properties were calculated by using generalized gradient approximation. Calculated properties are in good agreement with available experimental values. Electronic and thermoelectric properties were investigated both with and without considering spin-orbit coupling (SOC) effect which is found to have a strong influence on p-type Seebeck coefficient as well as the power factor of the ZB-AgI. By inclusion of SOC, a reduction of the band-gap and p-type Seebeck coefficients as well as the power factor was found which is the indication of that spin-orbit interaction cannot be ignored for p-type thermoelectric properties of the ZB-AgI. By using deformation potential theory for electronic relaxation time and experimentally predicted lattice thermal conductivity, we obtained aZTvalue 1.69 (0.89) at 400 K for n-type (p-type) carrier concentration of 1.5 x 10(18)(4.6 x10(19)) cm(-3)that makes ZB-AgI as a promising room temperature thermoelectric material.

Abstract: Recent experimental measurements of light absorption in few-layer black phosphorus (BP) revealed a series of high and sharp peaks, interspersed by pairs of lower and broader features. Here, we propose a theoretical model for these excitonic states in few-layer BP within a continuum approach for the in-plane degrees of freedom and a tight-binding approximation that accounts for interlayer couplings. This yields excitonic transitions between different combinations of the subbands created by the coupled BP layers, which leads to a series of high and low oscillator strength excitonic states, consistent with the experimentally observed bright and dark exciton peaks, respectively. The main characteristics of such subband exciton states, as well as the possibility to control their energies and oscillator strengths via applied electric and magnetic fields, are discussed, towards a full understanding of the excitonic spectrum of few-layer BP and its tunability.

Abstract: Two-dimensional systems of inverse patchy colloids modeled as disks with a central charge and having their surface decorated with oppositely pointlike charged patches are investigated using molecular dynamics simulations. The self-assembly of the patchy colloids leads to diverse ground state configurations ranging from crystalline arrangements of monomers to linear clusters, ramified linear clusters and to percolated configurations. Two structural phase diagrams are constructed: (1) as a function of the net charge and area fraction, and (2) as a function of the net charge and the range of the pair interaction potential. An interesting reentrant percolation transition is obtained as a function of the net charge of the colloids. We identify distinct mechanisms that lead to the percolation transition.