Abstract: Simulations are carried out based on the dynamical mean-field theory (DMFT) in order to investigate the properties of correlated thin films for various values of the chemical potential, temperature, interaction strength, and applied transverse electric field. Application of a sufficiently strong field to a thin film at half filling leads to the appearance of conducting regions near the surfaces of the film, whereas in doped slabs the application of a field leads to a conductivity enhancement on one side of the film and a gradual transition to the insulating state on the opposite side. In addition to the inhomogeneous DMFT, a local density approximation (LDA) is considered in which the particle density n, quasiparticle residue Z, and spectral weight at the Fermi level A(ω=0) of each layer are approximated by a homogeneous bulk environment. A systematic comparison between the two approaches reveals that the less expensive LDA results are in good agreement with the DMFT approach, except close to the metal-to-insulator transition points and in the layers immediately at the film surfaces. LDA values for n are overall more reliable than those for Z and A(ω=0). The hysteretic behavior (memory effect) characteristic of the bulk doping driven Mott transition persists in the slab.

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: A tight-binding approach based on the Chebyshev-Bogoliubov-de Gennes method is used to describe disordered single-layer graphene Josephson junctions. Scattering by vacancies, ripples, or charged impurities is included. We compute the Josephson current and investigate the nature of multiple Andreev reflections, which induce bound states appearing as peaks in the density of states for energies below the superconducting gap. In the presence of single-atom vacancies, we observe a strong suppression of the supercurrent, which is a consequence of strong intervalley scattering. Although lattice deformations should not induce intervalley scattering, we find that the supercurrent is still suppressed, which is due to the presence of pseudomagnetic barriers. For charged impurities, we consider two cases depending on whether the average doping is zero, i.e., existence of electron-hole puddles, or finite. In both cases, short-range impurities strongly affect the supercurrent, similar to the vacancies scenario.

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: When a tunneling barrier between two superconductors is formed by a normal material that would be a superconductor in the absence of phase fluctuations, the resulting Josephson effect can undergo an enormous enhancement. We establish this novel proximity effect by a general argument as well as a numerical simulation and argue that it may underlie recent experimental observations of the giant proximity effect between two cuprate superconductors separated by a barrier made of the same material rendered normal by severe underdoping.

Abstract: Using a numerical tight-binding approach based on the Chebyshev–Bogoliubov–de Gennes method we describe Josephson junctions made of multilayer graphene contacted by top superconducting gates. Both Bernal (ABA) and rhombohedral (ABC) stacking are considered and we find that the type of stacking has a strong effect on the proximity effect and the supercurrent flow. For both cases the pair amplitude shows a polarization between dimer and nondimer atoms, being more pronounced for rhombohedral stacking. Even though the proximity effect in nondimer sites is enhanced when compared to single-layer graphene, we find that the supercurrent is suppressed. The spatial distribution of the supercurrent shows that for Bernal stacking the current flows only in the topmost layers while for rhombohedral stacking the current flows throughout the whole structure.

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: 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: The manner in which spin-polarized electrons interact with a magnetized thin film is currently described by a semi-classical approach. This in turn provides our present understanding of the spin transfer, or spin torque phenomenon. However, spin is an intrinsically quantum-mechanical quantity. Here, we make the first strides towards a fully quantum-mechanical description of spin transfer through spin currents interacting with a Heisenberg-coupled spin chain. Because of quantum entanglement, this requires a formalism based on the density matrix approach. Our description illustrates how individual spins in the chain time-evolve as a result of spin transfer.

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: <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: The in-plane dc conductivity of twisted bilayer graphene is calculated using an expansion of the real-space Kubo-Bastin conductivity in terms of Chebyshev polynomials. We investigate within a tight-binding approach the transport properties as a function of rotation angle, applied perpendicular electric field, and vacancy disorder. We find that for high-angle twists, the two layers are effectively decoupled, and the minimum conductivity at the Dirac point corresponds to double the value observed in monolayer graphene. This remains valid even in the presence of vacancies, hinting that chiral symmetry is still preserved. On the contrary, for low twist angles, the conductivity at the Dirac point depends on the twist angle and is not protected in the presence of disorder. Furthermore, for low angles and in the presence of an applied electric field, we find that the chiral boundary states emerging between AB and BA regions contribute to the dc conductivity, despite the appearance of localized states in the AA regions. The results agree qualitatively with recent transport experiments in low-angle twisted bilayer graphene.

Abstract: The interplay between quantum Hall states and Cooper pairs is usually hindered by the suppression of the superconducting state due to the strong magnetic fields needed to observe the quantum Hall effect. From this point of view, graphene is special since it allows the creation of strong pseudomagnetic fields due to strain. We show that in a Josephson junction made of strained graphene, Cooper pairs will diffuse into the strained region. The pair correlation function will be sublattice polarized due to the polarization of the local density of states in the zero pseudo-Landau level. We uncover two regimes: (1) one in which the cyclotron radius is larger than the junction length, in which case the supercurrent will be enhanced, and (2) the long junction regime where the supercurrent is strongly suppressed because the junction becomes an insulator. In the latter case quantized Hall states form and Andreev scattering at the normal/superconducting interface will induce edge states. Our numerical calculation has become possible due to an extension of the Chebyshev-Bogoliubovde Gennes method to computations on video cards (GPUs).

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: <script type='text/javascript'>document.write(unpmarked('Spin-triplet chiral p-wave superconductivity is typically described by a two-component order parameter, and as such is prone to unique emergent effects when compared to the standard single-component superconductors. Here we present the equilibrium phase diagram for small mesoscopic chiral p-wave superconducting disks in the presence of magnetic field, obtained by solving the microscopic Bogoliubov-de Gennes equations self-consistently. In the ultrasmall limit, the cylindrically symmetric giant-vortex states form the ground state of the system. However, with increasing sample size, the cylindrical symmetry is broken as the two components of the order parameter segregate into domains, and the number of fragmented domain walls between them characterizes the resulting states. Such domain walls are topological defects unique for the p-wave order, and constitute a dominant phase in the mesoscopic regime. Moreover, we find two possible types of domain walls, identified by their chirality-dependent interaction with the edge states.'));