Abstract: Electrons in graphene, in addition to their spin, have two pseudospin degrees of freedom: sublattice and valley pseudospin. Valleytronics uses the valley degree of freedom as a carrier of information similarly to the way spintronics uses electron spin. We show how a double-barrier structure consisting of electric and vector potentials can be used to filter massive Dirac electrons based on their valley index. We study the resonant transmission through a finite number of barriers and we obtain the energy spectrum of a superlattice consisting of electric and vector potentials. When a mass term is included, the energy bands and energy gaps at the K and K′ points are different and they can be tuned by changing the potential.

Abstract: Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while nonhomogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced. DOI: 10.1103/PhysRevB.87.134509

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: Bogoliubov-de Gennes theory is used to investigate the effect of the size of a superconducting square on the vortex states in the quantum confinement regime. When the superconducting coherence length is comparable to the Fermi wavelength, the shape resonances of the superconducting order parameter have strong influence on the vortex configuration. Several unconventional vortex states, including asymmetric ones, giant-multivortex combinations, and states comprising giant antivortices, were found as ground states and their stability was found to be very sensitive on the value of k(F)xi(0), the size of the sample W, and the magnetic flux Phi. By increasing the temperature and/or enlarging the size of the sample, quantum confinement is suppressed and the conventional mesoscopic vortex states as predicted by the Ginzburg-Laudau (GL) theory are recovered. However, contrary to the GL results we found that the states containing symmetry-induced vortex-antivortex pairs are stable over the whole temperature range. It turns out that the inhomogeneous order parameter induced by quantum confinement favors vortex-antivortex molecules, as well as giant vortices with a rich structure in the vortex core-unattainable in the GL domain.

Abstract: The scattering of a Gaussian wave packet in armchair and zigzag graphene edges is theoretically investigated by numerically solving the time-dependent Schrodinger equation for the tight-binding model Hamiltonian. Our theory allows us to investigate scattering in reciprocal space, and depending on the type of graphene edge we observe scattering within the same valley, or between different valleys. In the presence of an external magnetic field, the well-known skipping orbits are observed. However, our results demonstrate that in the case of a pseudomagnetic field, induced by nonuniform strain, the scattering by an armchair edge results in a nonpropagating edge state.

Abstract: Experiments have confirmed that double monolayer graphene does not generate finite-temperature electron-hole superfluidity, because of very strong screening of the pairing attraction. The linear dispersing energy bands in monolayer graphene block any attempt to reduce the strength of the screening. We propose a hybrid device with two sheets of monolayer graphene in a modulated periodic perpendicular magnetic field. The field preserves the isotropic Dirac cones of the original monolayers but reduces the slope of the cones, making the monolayer Fermi velocity v(F) smaller. We demonstrate that with current experimental techniques, the reduction in vF can weaken the screening sufficiently to allow electron-hole superfluidity at measurable temperatures.

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: Magnetization reversal is a well-studied problem with obvious applicability in computer hard drives. One can accomplish a magnetization reversal in at least one of two ways: application of a magnetic field or through a spin current. The latter is more amenable to a fully quantum-mechanical analysis. We formulate and solve the problem whereby a spin current interacts with a ferromagnetic Heisenberg spin chain, to eventually reverse the magnetization of the chain. Spin flips are accomplished through both elastic and inelastic scattering. A consequence of the inelastic-scattering channel, when it is no longer energetically possible, is the occurrence of a nonequilibrium bound state, which is an emergent property of the coupled local plus itinerant spin system. For certain definite parameter values the itinerant spin lingers near the local spins for some time, before eventually leaking out as an outwardly diffusing state. This phenomenon results in spin-flip dynamics and filtering properties for this type of system.

Abstract: Correlated, or extended, impurities play an important role in the transport properties of dirty metals. Here, we examine, in the framework of a tight-binding lattice, the transmission of a single electron through an array of correlated impurities. In particular we show that particles transmit through an impurity array in identical fashion, regardless of the direction of traversal. The demonstration of this fact is straightforward in the continuum limit, but requires a detailed proof for the discrete lattice. We also briefly demonstrate and discuss the time evolution of these scattering states, to delineate regions (in time and space) where the aforementioned symmetry is violated.

Abstract: <script type='text/javascript'>document.write(unpmarked('We show that delaminations in bilayer graphene (BLG) with electrostatically induced interlayer symmetry can provide one with ballistic channels for electrons with energies inside the electrostatically induced BLG gap. These channels are formed by a combination of valley-polarized evanescent states propagating along the delamination edges (which persist in the presence of a strong magnetic field) and standing waves bouncing between them inside the delaminated region (in a strong magnetic field, these transform into Landau levels in the monolayers). For inverted stackings in BLGs on the left and right of the delamination (AB-2ML-BA or BA-2ML-AB, where 2ML indicates two decoupled monolayers of graphene), the lowest-energy ballistic channels are gapless, have linear dispersion, and appear to be weakly topologically protected. When BLG stackings on both sides of the delamination are the same (AB-2ML-AB or BA-2ML-BA), the lowest-energy ballistic channels are gapped, with a gap epsilon(g) scaling as epsilon(g) alpha W-1 with delamination width and epsilon(g) alpha delta(-1) with the on-layer energy difference in the delaminated part of the structure. Depending on the width, delaminations may also support several \u0022higher-energy\u0022 waveguide modes. Our results are based on both the analytical study of the wave matching of Dirac states and tight-binding model calculations, and we analyze in detail the dependence of the delamination spectrum on the electrostatic conditions in the structure, such as the vertical displacement field.'));

Abstract: Extended defects are known to strongly affect nanoscale superconductors. Here, we report the properties of superconducting nanoribbons with a constriction formed between two adjacent step edges by solving the Bogoliubov-de Gennes equations self-consistently in the regime where quantum confinement is important. Since the quantum resonances of the superconducting gap in the constricted area are different from the rest of the nanoribbon, such constriction forms a quantum-confined S-S'-S Josephson junction, with a broadly tunable performance depending on the length and width of the constriction with respect to the nanoribbon, and possible gating. These findings provide an intriguing approach to further tailor superconducting quantum devices where Josephson effect is of use.

Abstract: The strain fields of periodically buckled graphene induce a periodic pseudomagnetic field (PMF) that modifies the electronic band structure. From the geometry, amplitude, and period of the periodic pseudomagnetic field, we determine the necessary conditions to access the regime of correlated phases by examining the band flattening. As compared to twisted bilayer graphene the proposed system has the advantages that (1) only a single layer of graphene is needed, (2) one is not limited to hexagonal superlattices, and (3) narrower flat bandwidth and larger separation between flat bands can be induced. We, therefore, propose that periodically strained graphene single layers can become a platform for the exploration of exotic many-body phases.

Abstract: We present KITE, a general purpose open-source tight-binding software for accurate real-space simulations of electronic structure and quantum transport properties of large-scale molecular and condensed systems with tens of billions of atomic orbitals (N similar to 10(10)). KITE's core is written in C++, with a versatile Python-based interface, and is fully optimized for shared memory multi-node CPU architectures, thus scalable, efficient and fast. At the core of KITE is a seamless spectral expansion of lattice Green's functions, which enables large-scale calculations of generic target functions with uniform convergence and fine control over energy resolution. Several functionalities are demonstrated, ranging from simulations of local density of states and photo-emission spectroscopy of disordered materials to large-scale computations of optical conductivity tensors and real-space wave-packet propagation in the presence of magneto-static fields and spin-orbit coupling. On-the-fly calculations of real-space Green's functions are carried out with an efficient domain decomposition technique, allowing KITE to achieve nearly ideal linear scaling in its multi-threading performance. Crystalline defects and disorder, including vacancies, adsorbates and charged impurity centres, can be easily set up with KITE's intuitive interface, paving the way to user-friendly large-scale quantum simulations of equilibrium and non-equilibrium properties of molecules, disordered crystals and heterostructures subject to a variety of perturbations and external conditions.

Abstract: We have grown epitaxial Cr-doped V2O3 thin films with Cr concentrations between 0% and 20% on (0001)-Al2O3 by oxygen-assisted molecular beam epitaxy. For the highly doped samples (>3%), a regular and monotonous increase of the resistance with decreasing temperature is measured. Strikingly, in the low doping samples (between 1% and 3%), a collapse of the insulating state is observed with a reduction of the low temperature resistivity by up to 5 orders of magnitude. A vacuum annealing at high temperature of the films recovers the low temperature insulating state for doping levels below 3% and increases the room temperature resistivity towards the values of Cr-doped V2O3 single crystals. It is well-know that oxygen excess stabilizes a metallic state in V2O3 single crystals. Hence, we propose that Cr doping promotes oxygen excess in our films during deposition, leading to the collapse of the low temperature insulating state at low Cr concentrations. These results suggest that slightly Cr-doped V2O3 films can be interesting candidates for field effect devices. (C) 2015 AIP Publishing LLC.

Abstract: Using an extension of the Gutzwiller approximation for an inhomogeneous system, we study the two-band Hubbard model with unequal band widths for a slab geometry. The aim is to investigate the mutual effect of individual bands on the spatial distribution of quasi-particle weight and charge density, especially near the surface of the slab. The main effect of the difference in band width is the presence of two different length scales corresponding to the quasi-particle profile of each band. This is enhanced in the vicinity of the critical interaction of the narrow band where an orbitally selective Mott transition occurs and a surface dead layer forms for the narrow band. For the doped case, two different regimes of charge transfer between the surface and the bulk of the slab are revealed. The charge transfer from surface/ center to center/ surface depends on both the doping level and the average relative charge accumulated in each band. Such effects could also be of importance when describing the accumulation of charges at the interface between structures made of multi-band strongly correlated materials.