Abstract: Static and dynamic properties of superconducting vortices in a superconducting stripe with a periodic array of weakly-superconducting (or normal metal) regions are studied in the presence of external magnetic and electric fields. The time-dependent Ginzburg-Landau theory is used to describe the electronic transport, where the anisotropy is included through the spatially-dependent critical temperature T-c. Superconducting vortices penetrating into the weak-superconducting region with smaller T-c are more mobile than the ones in the strong superconducting regions. We observe periodic entrance and exit of vortices which reside in the weak link for some short interval. The mobility of the weakly-pinned vortices can be reduced by increasing the uniform applied magnetic field leading to distinct features in the voltage vs. magnetic field response of the system.

Abstract: We investigate the effects of coupling to multiple-phonon modes on the properties of a Holstein polaron. To this end, we generalize the Momentum Average approximations MA((0)) and MA((1)) to deal with multiple-phonon modes. As for a single-phonon mode, these approximations are found to be numerically very efficient. They become exact for very weak or very strong couplings, and are highly accurate in the intermediate regimes, e.g. the spectral weights obey exactly the first six, respectively eight, sum rules. Our results show that the effect on ground-state properties is cumulative in nature. As a result, if the effective coupling to one mode is much larger than to all the others, this mode effectively determines the ground-state properties. However, even very weak coupling to a second phonon mode has important non-perturbational effects on the higher-energy spectrum, in particular on the dispersion and the phonon statistics of the polaron band. This has important consequences on the analysis and interpretation of data for real materials.

Abstract: A periodic spatial modulation, as created by a moire pattern, has been extensively studied with the view to engineer and tune the properties of graphene. Graphene encapsulated by hexagonal boron nitride (hBN) when slightly misaligned with the top and bottom hBN layers experiences two interfering moire patterns, resulting in a so-called supermoire (SM). This leads to a lattice and electronic spectrum reconstruction. A geometrical construction of the nonrelaxed SM patterns allows us to indicate qualitatively the induced changes in the electronic properties and to locate the SM features in the density of states and in the conductivity. To emphasize the effect of lattice relaxation, we report band gaps at all Dirac-like points in the hole doped part of the reconstructed spectrum, which are expected to be enhanced when including interaction effects. Our result is able to distinguish effects due to lattice relaxation and due to the interfering SM and provides a clear picture on the origin of recently experimentally observed effects in such trilayer heterostuctures.

Abstract: The mechanical and electronic properties of a graphene membrane placed on top of a triangular superlattice of nanopillars are investigated. We use molecular dynamics simulations to access the deformation fields and the tight-binding approaches to calculate the electronic properties. Ripples form in the graphene layer that span across the unit cell, connecting neighboring pillars, in agreement with recent experiments. We find that the resulting pseudo-magnetic field (PMF) varies strongly across the unit cell. We investigate the dependence of PMF on unit cell boundary conditions, height of the pillars, and the strength of the van der Waals interaction between graphene and the substrate. We find direct correspondence with typical experiments on pillars, showing intrinsic “slack” in the graphene membrane. PMF values are confirmed by the local density of states calculations performed at different positions of the unit cell showing pseudo-Landau levels with varying spacings. Our findings regarding the relaxed membrane configuration and the induced strains are transferable to other flexible 2D membranes.

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.

Abstract: The conductance through graphene nanoribbons (GNR) connected to a partially unzipped carbon nanotube (CNT) is studied in the presence of an external magnetic field applied parallel to the long axis of the tube by means of non-equilibrium Green's function technique. We consider CNTs that are partially unzipped to form armchair-GNR/zigzag-CNT/armchair-GNR or zigzag-GNR/armchair-CNT/zigzag-GNR junctions. We find that the inclusion of a longitudinal magnetic field affects the electronic states only in the CNT region, leading to the suppression of the conductance at low energies. We demonstrate that both types of junctions can be used as magnetic field sensors. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4726039]

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: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasiparticles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low-energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition-metal dichalcogenides. However, the method is easily adapted for other cases such as electrons in few-layer black phosphorus or any Schrodinger quasiparticles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in mono-, bi-, and trilayer graphene, where specular-and retro-reflection cases are observed for electrons scattered by a steplike superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.

Abstract: Chiral p-wave superconductors in applied magnetic field can exhibit more complex topological defects than just conventional superconducting vortices, due to the two-component order parameter (OP) and the broken time-reversal symmetry. We investigate the electronic properties of those exotic states, some of which contain clusters of one-component vortices in chiral components of the OP and/or exhibit skyrmionic character in the relative OP space, all obtained as a self-consistent solution of the microscopic Bogoliubov-de Gennes equations. We reveal the link between the local density of states (LDOS) of the novel topological states and the behavior of the chiral domain wall between the OP components, enabling direct identification of those states in scanning tunneling microscopy. For example, a skyrmion always contains a closed chiral domain wall, which is found to be mapped exactly by zero-bias peaks in LDOS. Moreover, the LDOS exhibits electron-hole asymmetry, which is different from the LDOS of conventional vortex states with same vorticity. Finally, we present the magnetic field and temperature dependence of the properties of a skyrmion, indicating that this topological defect can be surprisingly large in size, and can be pinned by an artificially indented nonsuperconducting closed path in the sample. These features are expected to facilitate the experimental observation of skyrmionic states, thereby enabling experimental verification of chirality in emerging superconducting materials.

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: Recent measurements of the Josephson critical current through LSCO/LCO/LSCO thin films showed an unusually large proximity effect. Using the Bogoliubov-de Gennes equations for a tight-binding Hamiltonian we describe the proximity effect in weak links between a superconductor with critical temperature T-c and one with critical temperature T-c('), where T-c > T-c('). The weak link (N-') is therefore a superconductor above its own critical temperature and the superconducting regions are considered to have either s-wave or d-wave symmetry. We note that the proximity effect is enhanced due to the presence of superconducting correlations in the weak link. The dc Josephson current is calculated, and we obtain a nonzero value for temperatures greater than T-c(') for sizes of the weak links that can be almost an order of magnitude greater than the conventional coherence length. Considering pockets of superconductivity in the N-' layer, we show that this can lead to an even larger effect on the Josephson critical current by effectively shortening the weak link.

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: Recent scanning tunneling microscopy (STM) measurements of the proximity effect in Au/La2−xSrxCuO4 and La1.55Sr0.45CuO4/La2−xSrxCuO4 bilayers showed a proximity-induced pseudogap [O. Yuli, I. Asulin, Y. Kalcheim, G. Koren, and O. Millo, Phys. Rev. Lett. 103, 197003 (2009)]. We describe the proximity effect in mesoscopic superconductor/normal-metal bilayers by using the Bogoliubov-de Gennes equations for a tight-binding Hamiltonian with competing antiferromagnetic and d-wave superconductivity orders. The temperature-dependent local density of states is calculated as a function of the distance from the interface. Bound state due to both d-wave and spin-density wave gaps are formed in the normal metal for energies less than the respective gaps. If there is a mismatch between the Fermi velocities in the two layers we observe that these states will shift in energy when spin-density wave order is present, thus inducing a minigap at finite energy. We conclude that the STM measurement in the proximity structures is able to distinguish between the two scenarios proposed for the pseudogap (competing or precursor to superconductivity).

Abstract: We utilize an exact variational numerical procedure to calculate the ground state properties of a polaron in the presence of a Rashba-like spin-orbit interaction. Our results corroborate previous work performed with the momentum average approximation and with weak-coupling perturbation theory. We find that spin-orbit coupling increases the effective mass in the regime with weak electron-phonon coupling, and decreases the effective mass in the regimes of intermediate and strong electron-phonon coupling. Analytical strong-coupling perturbation theory results confirm our numerical results in the small-polaron regime. A large amount of spin-orbit coupling can lead to a significant lowering of the polaron effective mass.