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: The ground-state properties of a paramagnetic Mott insulator at half-filling are investigated in the presence of an external electric field using the inhomogeneous Gutzwiller approximation for a single-band Hubbard model in a slab geometry. We find that the metal-insulator transition is shifted toward higher Hubbard repulsions by applying an electric field perpendicular to the slab. The main reason is the accumulation of charges near the surface. The spatial distribution of site-dependent quasiparticle weight shows that it is maximal in a few layers beneath the surface, while the central sites where the field is screened have a very low quasiparticle weight. Our results show that above a critical-field value, states near the surface will be metallic, while the bulk quasiparticle weight is extremely suppressed but never vanishing, even for large Hubbard repulsions above the bulk zero-field critical value. Below the critical-field value, our results hint toward an insulating state in which the electric field is totally screened and the slab is again at half-filling.

Abstract: Surface effects of a doped thin film made of a strongly correlated material are investigated both in the absence and presence of a perpendicular electric field. We use an inhomogeneous Gutzwiller approximation for a single-band Hubbard model in order to describe correlation effects. For low doping, the bulk value of the quasiparticle weight is recovered exponentially deep into the slab, but with increasing doping, additional Friedel oscillations appear near the surface. We show that the inverse correlation length has a power-law dependence on the doping level. In the presence of an electrical field, considerable changes in the quasiparticle weight can be realized throughout the system. We observe a large difference (as large as five orders of magnitude) in the quasiparticle weight near the opposite sides of the slab. This effect can be significant in switching devices that use the surface states for transport. DOI: 10.1103/PhysRevB.87.035131

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.

Abstract: Recent experiments showed that nonuniform strain can be produced by depositing graphene over pillars. We employed atomistic calculations to study the nonuniform strain and the induced pseudomagnetic field in graphene on top of nanopillars. By decreasing the distance between the nanopillars a complex distribution for the pseudomagnetic field can be generated. Furthermore, we performed tight-binding calculations of the local density of states (LDOS) by using the relaxed graphene configuration obtained from atomistic calculations. We find that the quasiparticle LDOS are strongly modified near the pillars, both at low energies showing sublattice polarization and at high energies showing shifts of the van Hove singularity. Our study shows that changing the specific pattern of the nanopillars allows us to create a desired shape of the pseudomagnetic field profile while the LDOS maps provide an input for experimental verification by scanning tunneling microscopy.

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: 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 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: Using highly efficient simulations of the tight-binding Bogoliubov-de-Gennes model, we solved self-consistently for the pair correlation and the Josephson current in a superconducting-bilayer graphene-superconducting Josephson junction. Different doping levels for the non-superconducting link are considered in the short- and long-junction regimes. Self-consistent results for the pair correlation and superconducting current resemble those reported previously for single-layer graphene except at the Dirac point, where remarkable differences in the proximity effect are found, as well as a suppression of the superconducting current in the long-junction regime. Inversion symmetry is broken by considering a potential difference between the layers and we found that the supercurrent can be switched if the junction length is larger than the Fermi length.

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 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: 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 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: Two-dimensional atomic crystals can radically change their properties in response to external influences, such as substrate orientation or strain, forming materials with novel electronic structure(1-5). An example is the creation of weakly dispersive, 'flat' bands in bilayer graphene for certain 'magic' angles of twist between the orientations of the two layers(6). The quenched kinetic energy in these flat bands promotes electron-electron interactions and facilitates the emergence of strongly correlated phases, such as superconductivity and correlated insulators. However, the very accurate fine-tuning required to obtain the magic angle in twisted-bilayer graphene poses challenges to fabrication and scalability. Here we present an alternative route to creating flat bands that does not involve fine-tuning. Using scanning tunnelling microscopy and spectroscopy, together with numerical simulations, we demonstrate that graphene monolayers placed on an atomically flat substrate can be forced to undergo a buckling transition(7-9), resulting in a periodically modulated pseudo-magnetic field(10-14), which in turn creates a 'post-graphene' material with flat electronic bands. When we introduce the Fermi level into these flat bands using electrostatic doping, we observe a pseudogap-like depletion in the density of states, which signals the emergence of a correlated state(15-17). This buckling of two-dimensional crystals offers a strategy for creating other superlattice systems and, in particular, for exploring interaction phenomena characteristic of flat bands. Buckled monolayer graphene superlattices are found to provide an alternative to twisted bilayer graphene for the study of flat bands and correlated states in a carbon-based material.

Abstract: We perform the wavepacket dynamics simulation of a graphene-based device where propagating electron trajectories are tamed by an applied magnetic field toward a normal/superconductor interface. The magnetic field controls the incidence angle of the incoming electronic wavepacket at the interface, which results in the tunable electron-hole ratio in the reflected wave function due to the angular dependence of the Andreev reflection. Here, mapped control of the quasiparticle trajectories by the external magnetic field not only defines an experimental probe for fundamental studies of the Andreev reflection in graphene but also lays the foundation for further development of magnetic focusing devices based on nanoengineered superconducting two-dimensional materials.