Abstract: The split-operator technique for wave packet propagation in quantum systems is expanded here to the case of propagatingwave functions describing Schrodinger particles, namely, charge carriers in semiconductor nanostructures within the effective mass approximation, in the presence of Zeeman effect, as well as of Rashba and Dresselhaus spin-orbit interactions. We also demonstrate that simple modifications to the expanded technique allow us to calculate the time evolution of wave packets describing Dirac particles, which are relevant for the study of transport properties in graphene.

Abstract: The Tomasch effect (TE) is due to quasiparticle interference (QPI) as induced by a nonuniform superconducting order parameter, which results in oscillations in the density of states (DOS) at energies above the superconducting gap. Quantum confinement in nanoscale superconductors leads to an inhomogenerous distribution of the Cooperpair condensate, which, as we found, triggers the manifestation of a new TE. We investigate the electronic structure of nanoscale superconductors by solving the Bogoliubov-de Gennes (BdG) equations self-consistently and describe the TE determined by two types of processes, involving two-or three-subband QPIs. Both types of QPIs result in additional BCS-like Bogoliubov-quasiparticles and BCS-like energy gaps leading to oscillations in the DOS and modulated wave patterns in the local density of states. These effects are strongly related to the symmetries of the system. A reduced 4 x 4 inter-subband BdG Hamiltonian is established in order to describe analytically the TE of two-subband QPIs. Our study is relevant to nanoscale superconductors, either nanowires or thin films, Bose-Einsten condensates, and confined systems such as two-dimensional electron gas interface superconductivity.

Abstract: To determine and control the spatial confinement of charge carriers is of importance for nanoscale optoelectronic device applications. Using first-principles calculations, we investigate the tunability of band alignment and Charge localization in lateral and combined lateral vertical heterostructures of MoS2 and WS2. First, we Show that a type-II to type-I band alignment transition takes place when tensile strain is applied on the WS2 region. This band alignment transition is a result of the different response of the band edge states with strain and is caused by their different wave function characters. Then we show that the presence of the grain boundary introduces localized in-gap states. The boundary at the armchair interface significantly modifies the charge distribution of the valence band maximum (VBM) state, whereas in a heterostructure with tilt grain domains both conducation band maximum (CBM) and VBM are found to be localized around the grain boundary. We also found that the thickness of the constituents in a lateral heterostructure also determines how the electrons and holes are confined. Creating combined lateral vertical heterostructures of MOS2/WS2 provides another way cif tuning the charge confinement. These results provide possible ways to tune the carrier confinement in MoS2/WS2 heterostructures, which are interesting for its practical: applications in the future.

Abstract: We propose systems that allow a tuning of the phonon transmission function T(omega) in graphene nanoribbons by using C-13 isotope barriers, antidot structures, and distinct boundary conditions. Phonon modes are obtained by an interatomic fifth-nearest neighbor force-constant model (5NNFCM) and T(omega) is calculated using the non-equilibrium Green's function formalism. We show that by imposing partial fixed boundary conditions it is possible to restrict contributions of the in-plane phonon modes to T(omega) at low energy. On the contrary, the transmission functions of out-of-plane phonon modes can be diminished by proper antidot or isotope arrangements. In particular, we show that a periodic array of them leads to sharp dips in the transmission function at certain frequencies omega(nu) which can be pre-defined as desired by controlling their relative distance and size. With this, we demonstrated that by adequate engineering it is possible to govern the magnitude of the ballistic transmission functions T(omega) in graphene nanoribbons. We discuss the implications of these results in the design of controlled thermal transport at the nanoscale as well as in the enhancement of thermo-electric features of graphene-based materials.

Abstract: Propagation of an electron wave packet through a quantum point contact (QPC) defined by electrostatic gates in bilayer graphene is investigated. The gates provide a bias between the layers, in order to produce an energy gap. If the gates on both sides of the contact produce the same bias, steps in the electron transmission probability are observed, as in the usual QPC. However, if the bias is inverted on one of the sides of the QPC, only electrons belonging to one of the Dirac valleys are allowed to pass, which provides a very efficient valley filtering.

Abstract: Different allotropes of phosphorene are possible of which black and blue phosphorus are the most stable. While blue phosphorus has isotropic properties, black phosphorus is strongly anisotropic in its electronic and optical properties due to its anisotropic crystal structure. In this work, we systematically investigated the lattice thermal properties of black and blue phosphorene by using first-principles calculations based on the quasiharmonic approximation approach. Similar to the optoelectronic and electronic properties, we predict that black phosphorene has highly anisotropic thermal properties, in contrast to the blue phase. The linear thermal expansion coefficients along the zigzag and armchair direction differ up to 20% in black phosphorene. The armchair direction of black phosphorene is more expandable as compared to the zigzag direction and the biaxial expansion of blue phosphorene under finite temperature. Our comparative analysis reveals that the inclusion of finite-temperature effects makes the blue phase thermodynamically more stable over the black phase above 135 K.

Abstract: The intrinsic field effect, the change in surface conductance with an applied transverse electric field, of prototypal strongly correlated VO(2) has remained elusive. Here we report its measurement enabled by epitaxial VO(2) and atomic layer deposited high-kappa dielectrics. Oxygen migration, joule heating, and the linked field-induced phase transition are precluded. The field effect can be understood in terms of field-induced carriers with densities up to approximately 5x10(13) cm(-2) which are trongly localized, as shown by their low, thermally activated mobility ( approximately 1x10(-3) cm(2)/V s at 300 K). These carriers show behavior consistent with that of Holstein polarons and strongly impact the (opto)electronics of VO(2).

Abstract: The electron states in axially symmetric quantum wires are computed by means of the effective-mass Schrodinger equation, which is written in cylindrical coordinates phi, rho, and z. We show that a direct discretization of the Schrodinger equation by central finite differences leads to a non-symmetric Hamiltonian matrix. Because diagonalization of such matrices is more complex it is advantageous to transform it in a symmetric form. This can be done by the Liouville-like transformation proposed by Rizea et al. (2008), which replaces the wave function psi(rho) with the function F(rho) = psi(rho)root rho and transforms the Hamiltonian accordingly. Even though a symmetric Hamiltonian matrix is produced by this procedure, the computed wave functions are found to be inaccurate near the origin, and the accuracy of the energy levels is not very high. In order to improve on this, we devised a finite-difference scheme which discretizes the Schrodinger equation in the first step, and then applies the Liouville-like transformation to the difference equation. Such a procedure gives a symmetric Hamiltonian matrix, resulting in an accuracy comparable to the one obtained with the finite element method. The superior efficiency of the new finite-difference scheme (FDM) is demonstrated for a few p-dependent one-dimensional potentials which are usually employed to model the electron states in free-standing and core shell quantum wires. The new scheme is compared with the other FDM schemes for solving the effective-mass Schrodinger equation, and is found to deliver energy levels with much smaller numerical error for all the analyzed potentials. It also gives more accurate results than the scheme of Rizea et al., except for the ground state of an infinite rectangular potential in freestanding quantum wires. Moreover, the PT symmetry is invoked to explain similarities and differences between the considered FDM schemes. (C) 2015 Elsevier B.V. All rights reserved.

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: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

Abstract: Motivated by the recent synthesis of single layer TiSe2, we used state-of-the-art density functional theory calculations, to investigate the structural and electronic properties of zigzag and armchair-edged nanoribbons (NRs) of this material. Our analysis reveals that, differing from ribbons of other ultra-thin materials such as graphene, TiSe2 NRs have some distinctive properties. The electronic band gap of the NRs decreases exponentially with the width and vanishes for ribbons wider than 20 angstrom. For ultranarrow zigzag-edged NRs we find odd-even oscillations in the band gap width, although their band structures show similar features. Moreover, our detailed magnetic-ground-state analysis reveals that zigzag and armchair edged ribbons have non-magnetic ground states. Passivating the dangling bonds with hydrogen at the edges of the structures influences the band dispersion. Our results shed light on the characteristic properties of T phase NRs of similar crystal structures.

Abstract: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

Abstract: Fluorinated graphene is one of the most important derivatives of graphene and has been found to have great potential in optoelectronic and photonic nanodevices. However, the stability of F atoms on fluorinated graphene under different conditions, which is essential to maintain the desired properties of fluorinated graphene, is still unclear. In this work, we investigate the diffusion of F atoms on pristine graphene, graphene with defects, and at graphene/fluorographene interfaces by using density functional theory calculations. We find that an isolated F atom diffuses easily on graphene, but those F atoms can be localized by inducing vacancies or absorbates in graphene and by creating graphene/fluorographene interfaces, which would strengthen the binding energy of F atoms on graphene and increase the diffusion energy barrier of F atoms remarkably.

Abstract: Using density functional theory, we study the electronic properties of a graphene-silicene bilayer (GSB). A single layer of silicene binds to the graphene layer with adhesion energy of about 25 meV/atom. This adhesion energy between the two layers follows accurately the well-known -1/z(2) dispersion energy as found between two infinite parallel plates. In small flakes of GSB with hydrogenated edges, negative charge is transferred from the graphene layer to the silicene layer, producing a permanent and a switchable polar bilayer, while in an infinite GSB, the negative charge is transferred from the silicene layer to the graphene layer. The graphene-silicene bilayer is a good candidate for a nanocapacitor with piezoelectric capabilities. We found that the permanent dipole of the bilayer can be tuned by an external perpendicular electric field.

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.