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 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 report the formation of nanobubbles on graphene with a radius of the order of 1 nm, using ultralow energy implantation of noble gas ions (He, Ne, Ar) into graphene grown on a Pt(111) surface. We show that the universal scaling of the aspect ratio, which has previously been established for larger bubbles, breaks down when the bubble radius approaches 1 nm, resulting in much larger aspect ratios. Moreover, we observe that the bubble stability and aspect ratio depend on the substrate onto which the graphene is grown (bubbles are stable for Pt but not for Cu) and trapped element. We interpret these dependencies in terms of the atomic compressibility of the noble gas as well as of the adhesion energies between graphene, the substrate, and trapped atoms.

Abstract: A twin boundary in one of the layers of a twisted van der Waals heterostructure separates regions with near opposite inter-layer twist angles. In a MoS₂/WSe₂bilayer, the regions with<inline-formula><tex-math><?CDATA $Rh^h$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>h</mi></msubsup></math><inline-graphic href=“tdmac529dieqn1.gif” type=“simple” /></inline-formula>and<inline-formula><tex-math><?CDATA $Rh^X$?></tex-math><math overflow=“scroll”><msubsup><mi>R</mi><mi>h</mi><mi>X</mi></msubsup></math><inline-graphic href=“tdmac529dieqn2.gif” type=“simple” /></inline-formula>stacking registry that defined the sub-lattices of the moiré honeycomb pattern would be mirror-reflected across such a twist boundary. In that case, we demonstrate that topologically protected chiral moiré exciton states are confined at the twist boundary. These are one-dimensional and uni-directional excitons with opposite velocities for excitons composed by electronic states with opposite valley/spin character, enabling intrinsic, guided, and far reaching valley-polarized exciton currents.

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: 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 examine from first principles the flexoelectric properties of phosphorene nanoribbons under mechanical bending along armchair and zigzag directions. In both cases we find that the radial polarization depends linearly on the strain gradient. The flexoelectricity along the armchair direction is over 40% larger than along the zigzag direction. The obtained flexoelectric coefficients of phosphorene are four orders of magnitude larger than those of graphene and comparable to transition metal dichalcogenides. Analysis of charge density shows that the flexoelectricity mainly arises from the pz orbitals of phosphorus atoms. The electron mobilities in bent phosphorene can be enhanced by over 60% along the armchair direction, which is significantly higher than previous reports of mobility tuned by uniaxial strain. Our results indicate phosphorene is a candidate for a two-dimensional material applicable in flexible-electronic devices.

Abstract: We numerically investigate the mixed pairing symmetries inmesoscopic superconducting loops in the presence of spin correlations by solving the Bogoliubov-de Gennes equations self-consistently. The spatial variations of the superconducting order parameters and the spontaneous magnetization are determined by the band structure. When the threaded magnetic flux turns on, the charge and spin currents both emerge and depict periodic evolution. In the case of a mesoscopic loop with dominant triplet p(x) +/- ip(y)-wave symmetry, a slight change of the chemical potential may lead to novel flux-dependent evolution patterns of the ground-state energy and the magnetization. The spin-polarized currents show pronounced quantum oscillations with fractional periods due to the appearance of energy jumps in flux, accompanied with a steplike feature of the enhanced spin current. Particularly, at some appropriate flux, the peaks of the zero-energy local density of states clearly indicate the occurrence of the odd-frequency pairing. In the case of a superconducting loop with dominant singlet d(x2-y2)-wave symmetry, the spatial profiles of the zero-energy local density of states and the magnetization show spin-dependent features on different sample diagonals. Moreover, the evolution of the flux-induced spin current always exhibits an hc/e periodicity.

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: 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: 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: 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: 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: 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: 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.