Abstract: We explore the flatness of conduction and valence bands of interlayer excitons in MoS2/WSe2 van der Waals heterobilayers, tuned by interlayer twist angle, pressure, and external electric field. We employ an efficient continuum model where the moire pattern from lattice mismatch and/or twisting is represented by an equivalent mesoscopic periodic potential. We demonstrate that the mismatch moire potential is too weak to produce significant flattening. Moreover, we draw attention to the fact that the quasi-particle effective masses around the Gamma-point and the band flattening are reduced with twisting. As an alternative approach, we show (i) that reducing the interlayer distance by uniform vertical pressure can significantly increase the effective mass of the moire hole, and (ii) that the moire depth and its band flattening effects are strongly enhanced by accessible electric gating fields perpendicular to the heterobilayer, with resulting electron and hole effective masses increased by more than an order of magnitude – leading to record-flat bands. These findings impose boundaries on the commonly generalized benefits of moire twistronics, while also revealing alternative feasible routes to achieve truly flat electron and hole bands to carry us to strongly correlated excitonic phenomena on demand.

Abstract: Despite its apparent simplicity, the idealized model of a particle constrained to move on a circle has intriguing dynamic properties and immediate experimental relevance. While a rotor is rather easy to set up classically, the quantum regime is harder to realize and investigate. Here we demonstrate that the quantum dynamics of quasiparticles in certain classes of nanostructured superconductors can be mapped onto a quantum rotor. Furthermore, we provide a straightforward experimental procedure to convert this nanoscale superconducting rotor into a regular or inverted quantum pendulum with tunable gravitational field, inertia, and drive. We detail how these novel states can be detected via scanning tunneling spectroscopy. The proposed experiments will provide insights into quantum dynamics and quantum chaos.

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 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 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: 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: 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: 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: 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: The time evolution of a wave packet in strained graphene is studied within the tight-binding model and continuum model. The effect of an external magnetic field, as well as a strain-induced pseudomagnetic field, on the wave-packet trajectories and zitterbewegung are analyzed. Combining the effects of strain with those of an external magnetic field produces an effective magnetic field which is large in one of the Dirac cones, but can be practically zero in the other. We construct an efficient valley filter, where for a propagating incoming wave packet consisting of momenta around the K and K' Dirac points, the outgoing wave packet exhibits momenta in only one of these Dirac points while the components of the packet that belong to the other Dirac point are reflected due to the Lorentz force. We also found that the zitterbewegung is permanent in time in the presence of either external or strain-induced magnetic fields, but when both the external and strain-induced magnetic fields are present, the zitterbewegung is transient in one of the Dirac cones, whereas in the other cone the wave packet exhibits permanent spatial oscillations.

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