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: Due to its strong bonds graphene can stretch up to 25% of its original size without breaking. Furthermore, mechanical deformations lead to the generation of pseudo-magnetic fields (PMF) that can exceed 300 T. The generated PMF has opposite direction for electrons originating from different valleys. We show that valley-polarized currents can be generated by local straining of multi-terminal graphene devices. The pseudo-magnetic field created by a Gaussian-like deformation allows electrons from only one valley to transmit and a current of electrons from a single valley is generated at the opposite side of the locally strained region. Furthermore, applying a pressure difference between the two sides of a graphene membrane causes it to bend/bulge resulting in a resistance change. We find that the resistance changes linearly with pressure for bubbles of small radius while the response becomes non-linear for bubbles that stretch almost to the edges of the sample. This is explained as due to the strong interference of propagating electronic modes inside the bubble. Our calculations show that high gauge factors can be obtained in this way which makes graphene a good candidate for pressure sensing.

Abstract: abstract not available

Abstract: We explore the flatness of conduction and valence bands of interlayer excitons in MoS2/WSe2 van der Waals heterobilayers, tuned by interlayer twist angle, pressure, and external electric field. We employ an efficient continuum model where the moire pattern from lattice mismatch and/or twisting is represented by an equivalent mesoscopic periodic potential. We demonstrate that the mismatch moire potential is too weak to produce significant flattening. Moreover, we draw attention to the fact that the quasi-particle effective masses around the Gamma-point and the band flattening are reduced with twisting. As an alternative approach, we show (i) that reducing the interlayer distance by uniform vertical pressure can significantly increase the effective mass of the moire hole, and (ii) that the moire depth and its band flattening effects are strongly enhanced by accessible electric gating fields perpendicular to the heterobilayer, with resulting electron and hole effective masses increased by more than an order of magnitude – leading to record-flat bands. These findings impose boundaries on the commonly generalized benefits of moire twistronics, while also revealing alternative feasible routes to achieve truly flat electron and hole bands to carry us to strongly correlated excitonic phenomena on demand.

Abstract: Polarons are quasiparticles emerging in materials from the interaction of extra charge carriers with the surrounding atomic lattice. They appear in a wide va- riety of compounds and can have a profound impact on their properties, making the concept of a polaron a central and ubiquitous topic in material science. Al- though the concept is known for about 75 years, the origin of polarons is not yet fully elucidated. This thesis focuses on WO 3 as a well-known prototypical system for studying polarons, which inherent polaronic nature is linked to its remark- able electrical and chromic properties. The primary objective of this research is to provide a comprehensive atomistic description and understanding of polaron formation in WO 3 using first-principles density functional theory (DFT) calcula- tions. Additionally, the investigation explores the interactions between polarons and the possibility of bipolaron formation. Following a systematic strategy, we first extensively analyze the dielectric and lattice dynamical properties of WO 3 in both the room-temperature P 2 1 /n and ground-state P 2 1 /c phases. Our specific focus is on characterizing the zone-center phonons, which serve as the founda- tion for identifying the phonon modes involved in the polaron formation and charge localization process. Subsequently, we examine the impact of structural distortions on the electronic structure of WO 3 to elucidate the interplay between structural distortions and electronic properties, thereby laying the groundwork for understanding electron-phonon couplings. By incorporating these critical fac- tors, we address our primary research goals. The most common explanation for the polaron formation is associated with the electrostatic screening of the extra charge by the polarizable lattice. Here, we show that, even in ionic crystals, this is not necessarily the case. We demonstrate that polarons in this compound arise primarily from non-polar atomic distortions. We then unveil that this unexpected behavior originates from the undoing of distortive atomic motions, which lowers the bandgap. As such, we coin the name of anti-distortive polaron and validate its appearance through a simple quantum-dot model, in which charge localization is the result of balancing structural, electronic, and confinement energy costs. Then, we also study the polaron-polaron interaction and present the formation of the antiferromagnetic W 4+ bipolaronic state with relatively large formation energy. Our analysis of the W 4+ bipolaronic distortions on the global structure reveals the same behavior as in experiments where the highly distorted monoclinic phase transforms into a tetragonal phase as a function of doping. Additionally, leveraging our previous findings on asymmetric polaronic distortion and examin- ing different merging orientations, we stabilize the antiferromagnetic W 5+ -W 5+ bipolaronic state with an energy lower than the W 4+ state. This thesis clari- fies the formation of unusual medium-size 2D polarons and bipolarons in WO3,which might be relevant to the whole family of ABO 3 perovskites, to which WO 3 is closely related. The simplicity of the concept provides also obvious guidelines for tracking similar behavior in other families of compounds.

Abstract: In this theoretical study, we assess the impact of a graphene capping layer on the resistivity of defective, extremely scaled interconnects. We investigate the effect of graphene capping on the electronic transport in ultrascaled interconnects, in the presence of grain boundary defects in the metal layer. We compare the results obtained using our quantum mechanical model to a simple parallel-conductor model and find that the parallel-conductor model does not capture the effect of the graphene cap correctly. At 0.5 nm metal thickness, the parallel-conductor model underestimates the conductivity by 3.0% to 4.0% for single-sided and double sided graphene capping, respectively.

Abstract: We calculate the contact resistance for an edge- and top-contacted 2D semiconductor. The contact region consists of a metal contacting a monolayer of MoS2 which is otherwise surrounded by SiO2. We use the quantum transmitting boundary method to compute the contact resistance as a function of the 2D semiconductor doping concentration. An effective mass Hamiltonian is used to describe the properties of the various materials. The electrostatic potentials are obtained by solving the Poisson equation numerically. We incorporate the effects of the image-force barrier lowering on the Schottky barrier and examine the impact on the contact resistance. At low doping concentrations, the contact resistance of the top contact is lower compared to edge contact, while at high doping concentrations, the edge contact exhibits lower resistance.