Abstract: We propose that low-strain Si/O superlattices can be constructed by connecting reconstructed Si{001} surfaces by Si-O-Si bridges. Ab initio calculations show that our models are energetically more favorable than all the models proposed so far. The part of our Si/O superlattice model is experimentally accessible just by oxidizing a Si( 001) substrate. To complete our Si/O superlattice model, we propose a three-step method. We also explore the potential of our Si/O superlattice models for new materials used in future Si electronics. We find that the location of the channel where the carriers travel can be controlled between the interfaces and the Si layers by the insertion of O atoms into the Si-Si dimers. By revealing the origins of the interface electron and hole states, we find that similar interface states should be easily achieved for Si slabs and Si substrates. Interestingly, the interface electrons and holes have small effective masses in the direction parallel to the channel and large effective masses in the direction normal to the channel, which makes the Si/O superlattices attractive to be used for channel materials. We also find that the valley splitting of Si is enhanced by the formation of the Si/O/Si interfaces, which is ideal for developing Si-based qubits. Our findings open new perspectives to design and control the electronic properties of Si.

Abstract: When floating on a two-dimensional surface of superfluid He-4, electrons arrange themselves in two-dimensional crystalline structure known as Wigner crystal. In channels, the boundaries interfere the crystalline order and in case of very narrow channels one observes a quasi-one-dimensional (quasi-1D) Wigner crystal formed by just a few rows of electrons and, ultimately, one row in the “quantum wire” regime. Recently, the “quantum wire” regime was accessed experimentally [D.G. Rees, H. Totsuji, K. Kono, Phys. Rev. Lett. 108, 176801 (2012)] resulting in unusual transport phenomena such as, e.g., oscillations in the electron conductance. Using molecular dynamics simulations, we study the nonlinear transport of electrons in channels with various types of constrictions: single and multiple symmetric and asymmetric geometrical constrictions with varying width and length, and saddle-point-type potentials with varying gate voltage. In particular, we analyze the average particle velocity of the particles and the corresponding electron current versus the driving force or the gate voltage. We have revealed a significant difference in the dynamics for long and short constrictions: The oscillations of the average velocity of the particles for the systems with short constrictions exhibit a clear correlation with the transitions between the states with different numbers of rows of particles; on the other hand, for the systems with longer constrictions these oscillations are suppressed. The obtained results qualitatively agree with the experimental observations. Next, we propose a FET-like structure that consists of a channel with asymmetric constrictions. We show that applying a transverse bias results either in increase of the average particle velocity or in its suppression thus allowing a flexible control tool over the electron transport. The advantage of the asymmetric FET is that it does not have a gate and it allows an easy control of relatively large electron flow. Furthermore, the asymmetric device can be used for rectification of an ac-driven electron flow. Our results bring important insights into the dynamics of electrons floating on the surface of superfluid He-4 in channels with constrictions and allow the effective control over the electron transport.

Abstract: The envelope function method traditionally employs a single basis set which, in practice, relates to a single material because the k.p matrix elements are generally only known in a particular basis. In this work, we defined a basis function transformation to alleviate this restriction. The transformation is completely described by the known inter-band momentum matrix elements. The resulting envelope function equation can solve the electronic structure in lattice matched heterostructures without resorting to boundary conditions at the interface between materials, while all unit-cell averaged observables can be calculated as with the standard envelope function formalism. In the case of two coupled bands, this heterostructure formalism is equivalent to the standard formalism while taking position dependent matrix elements. (C) 2015 Elsevier B.V. All rights reserved

Abstract: We numerically investigate the mixed pairing symmetries inmesoscopic superconducting loops in the presence of spin correlations by solving the Bogoliubov-de Gennes equations self-consistently. The spatial variations of the superconducting order parameters and the spontaneous magnetization are determined by the band structure. When the threaded magnetic flux turns on, the charge and spin currents both emerge and depict periodic evolution. In the case of a mesoscopic loop with dominant triplet p(x) +/- ip(y)-wave symmetry, a slight change of the chemical potential may lead to novel flux-dependent evolution patterns of the ground-state energy and the magnetization. The spin-polarized currents show pronounced quantum oscillations with fractional periods due to the appearance of energy jumps in flux, accompanied with a steplike feature of the enhanced spin current. Particularly, at some appropriate flux, the peaks of the zero-energy local density of states clearly indicate the occurrence of the odd-frequency pairing. In the case of a superconducting loop with dominant singlet d(x2-y2)-wave symmetry, the spatial profiles of the zero-energy local density of states and the magnetization show spin-dependent features on different sample diagonals. Moreover, the evolution of the flux-induced spin current always exhibits an hc/e periodicity.

Abstract: We investigate mesoscopic Josephson-junction arrays created by patterning superconducting disks on monolayer graphene, concentrating on the high-T/T-c regime of these devices and the phenomena which contribute to the superconducting glass state in diffusive arrays. We observe features in the magnetoconductance at rational fractions of flux quanta per array unit cell, which we attribute to the formation of flux-quantized vortices. The applied fields at which the features occur are well described by Ginzburg-Landau simulations that take into account the number of unit cells in the array. We find that the mean conductance and universal conductance fluctuations are both enhanced below the critical temperature and field of the superconductor, with greater enhancement away from the graphene Dirac point.

Abstract: The angular momentum quantization of chiral gapless modes confined to a circularly shaped interface between two different topological phases is investigated. By examining several different setups, we show analytically that the angular momentum of the topological modes exhibits a highly chiral behavior, and can be coupled to spin and/or valley degrees of freedom, reflecting the nature of the interface states. A simple general one-dimensional model, valid for arbitrarily shaped loops, is shown to predict the corresponding energies and the magnetic moments. These loops can be viewed as building blocks for artificial magnets with tunable and highly diverse properties.

Abstract: A Hg(Cd)Te insulator heterostructure hosts a two-dimensional electron system that can simulate the physics of Dirac fermions with only a single valley. We investigate the magnetotransport properties of this structure and show that, unlike most two-dimensional crystals with spin and valley coupled levels, the Shubnikov-de Haas oscillations exhibit a high spin polarization in the absence of any valley degree of freedom. This effect can be observed using magnetospectroscopy measurements for quantum well thicknesses corresponding to either the topologically trivial or quantum spin Hall phases. The pseudospin texture of the electrons near the Fermi level is also studied and we show that a tunable pseudospin-polarized quantum Hall liquid can only be observed for thicknesses corresponding to the inverted regime.

Abstract: The electronic properties of monolayer phosphorene are exotic due to its puckered structure and large intrinsic direct band gap. We derive and discuss its band structure in the presence of a perpendicular magnetic field. Further, we evaluate the magneto-optical Hall and longitudinal optical conductivities as functions of temperature, magnetic field, and Fermi energy, and show that they are strongly influenced by the magnetic field. The imaginary part of the former and the real part of the latter exhibit regular interband oscillations as functions of the frequency omega in the range (h) over bar omega similar to 1.5-2 eV. Strong intraband responses in the latter and weak ones in the former occur at much lower frequencies. The magneto-optical response can be tuned in the microwave-to-terahertz and visible frequency ranges in contrast with a conventional two-dimensional electron gas or graphene in which the response is limited to the terahertz regime. This ability to isolate carriers in an anisotropic structure may make phosphorene a promising candidate for new optical devices.

Abstract: The structural transitions of the ground state of a system of repulsively interacting particles confined in a quasi-one-dimensional channel, and the effect of the interparticle interaction as well as the functional form of the confinement potential on those transitions are investigated. Although the nonsequential ordering of transitions (non-SOT), i.e., the 1 – 2 – 4 – 3 – 4 – 5 – 6 – ... sequence of chain configurations with increasing density, is widely robust as predicted in a number of theoretical studies, the sequential ordering of transitions (SOT), i.e., the 1 – 2 – 3 – 4 – 5 – 6 – ... chain, is found as the ground state for long-ranged interparticle interaction and hard-wall-like confinement potentials. We found an energy barrier between every two different phases around its transition point, which plays an important role in the preference of the system to follow either a SOT or a non-SOT. However, that preferential transition requires also the stability of the phases during the transition. Additionally, we analyze the effect of a small structural disorder on the transition between the two phases around its transition point. Our results show that a small deformation of the triangular structure changes dramatically the picture of the transition between two phases, removing in a considerable region the non-SOT in the system. This feature could explain the fact that the non-SOT is, up to now, not observed in experimental systems, and suggests a more advanced experimental setup to detect the non-SOT.

Abstract: Recent experimental studies showed that formation of a two-dimensional crystal structure of nitrogenated holey graphene (NHG) is possible. Similar to graphene, NHGs have an atomically thin and strong crystal structure. Using first-principles calculations, we investigate the structural, phononic, and thermal properties of monolayer NHG crystal. Our charge analysis reveals that the charged holey sites of NHG provide a reactive ground for further functionalization by adatoms or molecules. We also found that similar to graphene, the NHG structure has quite high-frequency phonon modes and the presence of nitrogen atoms leads to the emergence of additional vibrational modes. Our phonon analysis reveals the presence of three characteristic Raman-active modes of NHG. Furthermore, the analysis of constant-volume heat capacity showed that the NHG structure has a linear temperature dependence in the low-temperature region. The strong lattice structure and unique thermal properties of the NHG crystal structure are desirable in nanoscale device applications.

Abstract: The properties of a new class of phosphorene nanotubes (PNT) are investigated by performing first-principles calculations. We demonstrate that it is advantageous to use blue phosphorene in order to make small nanotubes and propose a way to create low-energy PNTs by the inclusion of defect lines. Five different types of defect lines are investigated and incorporated in various combinations. The resulting defect-induced faceted PNTs have negligible bending stresses which leads to a reduction in the formation energy with respect to round PNTs. Our armchair faceted PNTs have similar formation energies than the recently proposed multiphase faceted PNTs, but they have a larger variety of possible structures. Our zigzag faceted PNTs have lower formation energies than round tubes and multiphase faceted nanotubes. The electronic properties of the defect-induced faceted PNTs are determined by the defect lines which control the band gap and the shape of the electronic states at the band edges. These band gaps increase with the radius of the nanotubes and converge to those of isolated defect lines.

Abstract: The Majorana zero modes and persistent spin current in mesoscopic d-wave-superconducting loops with spin-orbit (SO) interaction are investigated by numerically solving the spin-generalized Bogoliubov-de Gennes equations self-consistently. For some appropriate strength of the SO coupling, Majorana zero-energy states and sharp jumps of the spin-polarized currents can be observed when the highest energy levels cross the Fermi energy in the spectrum, leading to spin currents with opposite chirality flowing near the inner and outer edges of the sample. When the threaded magnetic flux turns on, four flux-dependent patterns of the persistent spin current with step-like features show up, accompanied by Majorana edge modes at flux values where the energy gap closes. Moreover, the Majorana zero mode is highly influenced by the direction of the Zeeman field. A finite in-plane field can lead to the gap opening since the inversion symmetry is broken. Remarkably, multiple Majorana zero-energy states occur in the presence of an out-of-plane field h(z), and the number of steps in the spin current evolution can be effectively tuned by the field strength due to the shift of Majorana zero modes. Finally, when the loop sample contains surface indentation defects, zero-energy modes can always show up in the presence of an appropriate h(z). Interestingly, multiple Majorana states may be present in the system with a corner defect even if h(z) = 0.

Abstract: The effect of the number of stacking layers and the type of stacking on the electronic and optical properties of bilayer and trilayer black phosphorus are investigated by using first-principles calculations within the framework of density functional theory. We find that inclusion of many-body effects (i.e., electron-electron and electron-hole interactions) modifies strongly both the electronic and optical properties of black phosphorus. While trilayer black phosphorus with a particular stacking type is found to be a metal by using semilocal functionals, it is predicted to have an electronic band gap of 0.82 eV when many-body effects are taken into account within the G(0)W(0) scheme. Though different stacking types result in similar energetics, the size of the band gap and the optical response of bilayer and trilayer phosphorene are very sensitive to the number of layers and the stacking type. Regardless of the number of layers and the type of stacking, bilayer and trilayer black phosphorus are direct band gap semiconductors whose band gaps vary within a range of 0.3 eV. Stacking arrangements that are different from the ground state structure in both bilayer and trilayer black phosphorus exhibit significant modified valence bands along the zigzag direction and result in larger hole effective masses. The optical gap of bilayer (trilayer) black phosphorus varies by 0.4 (0.6) eV when changing the stacking type. The calculated binding energy of the bound exciton hardly changes with the type of stacking and is found to be 0.44 (0.30) eV for bilayer (trilayer) phosphorous.

Abstract: Hydrogen adsorption on boron and nitrogen doped graphene is investigated in detail by means of first-principles calculations. A comprehensive study is performed of the structural, electronic, and magnetic properties of chemisorbed hydrogen atoms and atom pairs near the dopant sites. The main effect of the substitutional atoms is charge doping which is found to greatly affect the adsorption process by increasing the binding energy at the sites closest to the substitutional species. It is also found that doping does not induce magnetism despite the odd number of electrons per atom introduced by the foreign species, and that it quenches the paramagnetic response of chemisorbed H atoms on graphene. Overall, the effects are similar for B and N doping, with only minor differences in the adsorption energetics due to different sizes of the dopant atoms and the accompanying lattice distortions.

Abstract: The structural properties of partially hydrogenated and fluorinated graphene with different percentages of H/F atoms are investigated using molecular dynamics simulations based on reactive force field (ReaxFF) potentials. We found that the roughness of graphene varies with the percentage (p) of H or F and in both cases is maximal around p = 50%. Similar results were obtained for partially oxidized graphene. The two-dimensional area size of partially fluorinated and hydrogenated graphene exhibits a local minimum around p = 35% coverage. The lattice thermal contraction in partially functionalized graphene is found to be one order of magnitude larger than that of fully covered graphene. We also show that the armchair structure for graphene oxide (similar to the structure of fully hydrogenated and fluorinated graphene) is unstable. Our results show that the structure of partially functionalized graphene changes nontrivially with the C : H and C : F ratio as well as with temperature.