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.

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: We describe an efficient numerical approach to calculate the longitudinal and transverse Kubo conductivities of large systems using Bastin's formulation. We expand the Green's functions in terms of Chebyshev polynomials and compute the conductivity tensor for any temperature and chemical potential in a single step. To illustrate the power and generality of the approach, we calculate the conductivity tensor for the quantum Hall effect in disordered graphene and analyze the effect of the disorder in a Chern insulator in Haldane's model on a honeycomb lattice.

Abstract: Vortex matter in mesoscopic superconductors is known to be strongly affected by the geometry of the sample. Here we show that in nanoscale superconductors with coherence length comparable to the Fermi wavelength the shape resonances of the order parameter results in an additional contribution to the quantum topological confinement-leading to unconventional vortex configurations. Our Bogoliubov-de Gennes calculations in a square geometry reveal a plethora of asymmetric, giant multivortex, and vortex-antivortex structures, stable over a wide range of parameters and which are very different from those predicted by the Ginzburg-Landau theory. These unconventional states are relevant for high-T-c nanograins, confined Bose-Einstein condensates, and graphene flakes with proximity-induced superconductivity.

Abstract: A Josephson phase shift can be induced in a Josephson junction by a strategically nearby pinned Abrikosov vortex (AV). For an asymmetric distribution of an imprinted phase along the junction (controlled by the position of the AV) such a simple system is capable of rectification of ac current in a broad and tunable frequency range. The resulting rectified voltage is a consequence of the directed motion of a Josephson antivortex which forms a pair with the AV when at local equilibrium. The proposed realization of the ratchet potential by an imprinted phase is more efficient than the asymmetric geometry of the junction itself, is easily realizable experimentally, and provides rectification even in the absence of an applied magnetic field.

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: We study the effects of a nearby surface on the spectral weight of a Holstein polaron, using the inhomogeneous momentum average approximation which is accurate over the entire range of electron-phonon (e-ph) coupling strengths. The broken translational symmetry is taken into account exactly. We find that the e-ph coupling gives rise to a large additional surface potential, with strong retardation effects, which may bind surface states even when they are not normally expected. The surface, therefore, has a significant effect and bulk properties are recovered only very far away from it. These results demonstrate that interpretation in terms of bulk quantities of spectroscopic data sensitive only to a few surface layers is not always appropriate.

Abstract: We study the effects of the Rashba spin-orbit coupling on polaron formation, using a suitable generalization of the momentum average approximation. While previous work on a parabolic band model found that spin-orbit coupling increases the effective mass, we show that the opposite holds for a tight-binding model, unless both the spin-orbit and the electron-phonon couplings are weak. It is thus possible to lower the effective mass of the polaron by increasing the spin-orbit coupling. We also show that when the spin-orbit coupling is large as compared to the phonon energy, the polaron retains only one of the spin-polarized bands in its coherent spectrum. This has major implications for the propagation of spin-polarized currents in such materials, and thus for spintronic applications.

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: We study the effects of the rippling of a graphene sheet on quasiparticle dispersion. This is achieved using a generalization to the honeycomb lattice of the momentum average approximation, which is accurate for all coupling strengths and at all energies. We show that even though the position of the Dirac points may move and the Fermi speed can be renormalized significantly, quasiparticles with very long lifetimes survive near the Dirac points even for very strong couplings.

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: <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.'));

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.