Abstract: Ultrathin superconductors of different materials are becoming a powerful platform to find mechanisms for enhancement of superconductivity, exploiting shape resonances in different superconducting properties. Here, we evaluate the superconducting gap and its spatial profile, the multiple gap components, and the chemical potential, of generic superconducting nanofilms, considering the pairing attraction and its energy scale as tunable parameters, from weak to strong coupling, at fixed electron density. Superconducting properties are evaluated at mean field level as a function of the thickness of the nanofilm, in order to characterize the shape resonances in the superconducting gap. We find that the most pronounced shape resonances are generated for weakly coupled superconductors, while approaching the strong coupling regime the shape resonances are rounded by a mixing of the subbands due to the large energy gaps extending over large energy scales. Finally, we find that the spatial profile, transverse to the nanofilm, of the superconducting gap acquires a flat behavior in the shape resonance region, indicating that a robust and uniform multigap superconducting state can arise at resonance.

Abstract: We show numerically that a harmonically trapped and coherently Rabi-coupled three-component Bose-Einstein condensate can host unconventional vortex lattices in its rotating ground state. The discovered lattices incorporate square and zig-zag patterns, vortex dimers and chains, and doubly quantized vortices, and they can be quantitatively classified in terms of a skyrmionic topological index, which takes into account the multicomponent nature of the system. The exotic ground-state lattices arise due to the intricate interplay of the repulsive density-density interactions and the Rabi couplings as well as the ubiquitous phase frustration between the components. In the frustrated state, domain walls in the relative phases can persist between some components even at strong Rabi coupling, while vanishing between others. Consequently, in this limit the three-component condensate effectively approaches a two-component condensate with only density-density interactions. At intermediate Rabi coupling strengths, however, we face unique vortex physics that occurs neither in the two-component counterpart nor in the purely density-density-coupled three-component system.

Abstract: Electrons on liquid helium can form different phases depending on density, and temperature. Also the electron-ripplon coupling strength influences the phase diagram, through the formation of so-called “ripplonic polarons”, that change how electrons are localized, and that shifts the transition between the Wigner solid and the liquid phase. We use an all-coupling, finite-temperature variational method to study the formation of a ripplopolaron Wigner solid on a liquid helium film for different regimes of the electron-ripplon coupling strength. In addition to the three known phases of the ripplopolaron system (electron Wigner solid, polaron Wigner solid, and electron fluid), we define and identify a fourth distinct phase, the ripplopolaron liquid. We analyse the transitions between these four phases and calculate the corresponding phase diagrams. This reveals a reentrant melting of the electron solid as a function of temperature. The calculated regions of existence of the Wigner solid are in agreement with recent experimental data.

Abstract: van der Waals (vdW) forces play an important role in the adsorption of molecules on the surface of solids. However, the choice of the most suitable vdW functional for different systems is an essential problem which must be addressed for different systems. The lack of a systematic study on the performance of the vdW functionals in the bulk and adsorption properties of metal-oxides motivated us to examine different vdW approaches and compute the bulk and molecular adsorption properties of alpha-Cr2O3, alpha-Fe2O3, and alpha-Al2O3. For the bulk properties, we compared our results for the heat of formation, cohesive energy, lattice parameters and bond distances between the different vdW functionals and available experimental data. Next we studied the adsorption of benzene and CH3 molecules on top of different oxide surfaces. We employed different approximations to exchange and correlation within DFT, namely, the Perdew-Burke-Ernzerhof (PBE) GGA, (PBE)+U, and vdW density functionals [ DFT(vdW-DF/DF2/optPBE/optB86b/optB88)+U] as well as DFT-D2/D3(+U) methods of Grimme for the bulk calculations and optB86b-vdW(+U) and DFT-D2(+U) for the adsorption energy calculations. Our results highlight the importance of vdW interactions not only in the adsorption of molecules, but importantly also for the bulk properties. Although the vdW contribution in the adsorption of CH3 (as a chemisorption interaction) is less important compared to the adsorption of benzene (as a physisorption interaction), this contribution is not negligible. Also adsorption of benzene on ferryl/chromyl terminated surfaces shows an important chemisorption contribution in which the vdW interactions become less significant.

Abstract: When the charge Z of an atom exceeds the critical value of 170, it will undergo a process called atomic collapse which triggers the spontaneous creation of electron-positron pairs. The high charge requirements have prevented the observation of this phenomenon with real atomic nuclei. However, thanks to the relativistic nature of the carriers in graphene, the same physics is accessible at a much lower scale. The atomic collapse analogue in graphene is realized using artificial nuclei which can be created via the deposition of impurities on the surface of graphene or using charged vacancies. These supercritically charged artificial nuclei trap electrons in a sequence of quasi-bound states which can be observed experimentally as resonances in the local density of states.

Abstract: Fermi's golden rule underpins the investigation of mobile carriers propagating through various solids, being a standard tool to calculate their scattering rates. As such, it provides a perturbative estimate under the implicit assumption that the effect of the interaction Hamiltonian which causes the scattering events is sufficiently small. To check the validity of this assumption, we present a general framework to derive simple validity criteria in order to assess whether the scattering rates can be trusted for the system under consideration, given its statistical properties such as average size, electron density, impurity density et cetera. We derive concrete validity criteria for metallic nanowires with conduction electrons populating a single parabolic band subjected to different elastic scattering mechanisms: impurities, grain boundaries and surface roughness.

Abstract: Chalcopyrite semiconductors are of considerable interest for application as absorber layers in thin-film photovoltaic cells. When growing films of these compounds, however, they are often found to contain CuAu-like domains, a metastable phase of chalcopyrite. It has been reported that for CuInS2, the presence of the CuAu-like phase improves the short circuit current of the chalcopyrite-based photovoltaic cell. We investigate the thermodynamic stability of both phases for a selected list of I-III-VI2 materials using a first-principles density functional theory approach. For the CuIn-VI2 compounds, the difference in formation energy between the chalcopyrite and CuAu-like phase is found to be close to 2 meV per atom, indicating a high likelihood of the presence of CuAu-like domains. Next, we calculate the spectroscopic limited maximum efficiency (SLME) of the CuAu-like phase and compare the results with those of the corresponding chalcopyrite phase. We identify several candidates with a high efficiency, such as CuAu-like CuInS2, for which we obtain an SLME of 29% at a thickness of 500 nm. We observe that the SLME can have values above the Shockley-Queisser (SQ) limit, and show that this can occur because the SQ limit assumes the absorptivity to be a step function, thus overestimating the radiative recombination in the detailed balance approach. This means that it is possible to find higher theoretical efficiencies within this framework simply by calculating the J-V characteristic with an absorption spectrum. Finally, we expand our SLME analysis to indirect band gap absorbers by studying silicon, and find that the SLME quickly overestimates the reverse saturation current of indirect band gap materials, drastically lowering their calculated efficiency.

Abstract: In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the ul-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Abstract: Ginzburg-Landau theory is one of the most powerful phenomenological theories in physics, with particular predictive value in superconductivity. The formalism solves coupled nonlinear differential equations for both the electronic and magnetic responsiveness of a given superconductor to external electromagnetic excitations. With order parameter varying on the short scale of the coherence length, and the magnetic field being long-range, the numerical handling of 3D simulations becomes extremely challenging and time-consuming for realistic samples. Here we show precisely how one can employ graphics-processing units (GPUs) for this type of calculations, and obtain physics answers of interest in a reasonable time-frame – with speedup of over 100x compared to best available CPU implementations of the theory on a 2563grid. (C) 2016 Elsevier Inc. All rights reserved.

Abstract: Recently penta-graphene was proposed as a stable two-dimensional carbon allotrope consisting of a single layer of interconnected carbon pentagons [Zhang et al., PNAS, 2015, 112, 2372]. Its silicon counterpart, penta-silicene, however, is not stable. In this work, we show that multilayers of penta-silicene form stable materials with semiconducting or metallic properties, depending on the stacking mode. We demonstrate their dynamic stability through their phonon spectrum and using molecular dynamics. A particular type of bilayer penta-silicene is found to have lower energy than all of the known hexagonal silicene bilayers and forms therefore the most stable bilayer silicon material predicted so far. The electronic and mechanical properties of these new silicon allotropes are studied in detail and their behavior under strain is investigated. We demonstrate that strain can be used to tune its band gap.

Abstract: Apart from the most studied tin-oxide compounds, SnO and SnO2, intermediate states have been claimed to exist for more than a hundred years. In addition to the known homologous series (Seko et al., Phys. Rev. Lett. 100, 045702 (2008)), we here predict the existence of several new compounds with an O concentration between 50 % (SnO) and 67 % (SnO2). All these intermediate compounds are constructed from removing one or more (101) oxygen layers of SnO2. Since the van der Waals (vdW) interaction is known to be important for the Sn-Sn interlayer distances, we use a vdW-corrected functional, and compare these results with results obtained with PBE and hybrid functionals. We present the electronic properties of the intermediate structures and we observe a decrease of the band gap when (i) the O concentration increases and (ii) more SnO-like units are present for a given concentration. The contribution of the different atoms to the valence and conduction band is also investigated.

Abstract: A weakly coupled system of two crossed graphene nanoribbons exhibits direct tunneling due to the overlap of the wavefunctions of both ribbons. We apply the Bardeen transfer Hamiltonian formalism, using atomistic band structure calculations to account for the effect of the atomic structure on the tunneling process. The strong quantum-size confinement of the nanoribbons is mirrored by the one-dimensional character of the electronic structure, resulting in properties that differ significantly from the case of inter-layer tunneling, where tunneling occurs between bulk two-dimensional graphene sheets. The current-voltage characteristics of the inter-ribbon tunneling structures exhibit resonance, as well as stepwise increases in current. Both features are caused by the energetic alignment of one-dimensional peaks in the density-of-states of the ribbons. Resonant tunneling occurs if the sign of the curvature of the coupled energy bands is equal, whereas a step-like increase in the current occurs if the signs are opposite. Changing the doping modulates the onset-voltage of the effects as well as their magnitude. Doping through electrostatic gating makes these structures promising for application towards steep slope switching devices. Using the atomistic empirical pseudopotentials based Bardeen transfer Hamiltonian method, inter-ribbon tunneling can be studied for the whole range of two-dimensional materials, such as transition metal dichalcogenides. The effects of resonance and of step-like increases in the current we observe in graphene ribbons are also expected in ribbons made from these alternative two-dimensional materials, because these effects are manifestations of the one-dimensional character of the density-of-states. Published by AIP Publishing.

Abstract: Density-functional theory in combination with the nonequilibrium Green's function formalism is used to study the effect of substitutional doping on the electronic transport properties of hydrogen passivated zig-zag graphene nanoribbon devices. B, N and Si atoms are used to substitute carbon atoms located at the center or at the edge of the sample. We found that Si -doping results in better electronic transport as compared to the other substitutions. The transmission spectrum also depends on the location of the substitutional dopants: for single atom doping the largest transmission is obtained for edge substitutions, whereas substitutions in the middle of the sample give larger transmission for double carbon substitutions. The obtained results are explained in terms of electron localization in the system due to the presence of impurities. (C) 2016 Elsevier B.V. All rights reserved.

Abstract: The interaction between the surface of a 3D topological insulator and a skyrmion/anti-skyrmion structure is studied in order to investigate the possibility of electron confinement due to the skyrmion presence. Both hedgehog (Neel) and vortex (Bloch) skyrmions are considered. For the hedgehog skyrmion, the in-plane components cannot be disregarded and their interaction with the surface state of the topological insulator (TI) has to be taken into account. A semi-classical description of the skyrmion chiral angle is obtained using the variational principle. It is shown that both the hedgehog and the vortex skyrmion can induce bound states on the surface of the TI. However, the number and the properties of these states depend strongly on the skyrmion type and the skyrmion topological number N-Sk. The probability densities of the bound electrons are also derived where it is shown that they are localized within the skyrmion region. Published by AIP Publishing.