Abstract: We utilize an exact variational numerical procedure to calculate the ground-state properties of a polaron in the presence of Rashba and linear Dresselhaus spin-orbit coupling. We find that when the linear Dresselhaus spin-orbit coupling approaches the Rashba spin-orbit coupling, the Van Hove singularity in the density of states will be shifted away from the bottom of the band and finally disappear when the two spin-orbit couplings are tuned to be equal. The effective mass will be suppressed; the trend will become more significant for low phonon frequency. The presence of two dominant spin-orbit couplings will make it possible to tune the effective mass with more varied observables.

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 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: 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: 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: We provide a comprehensive analysis of the prominent tight-binding (TB) models for transition metal dichalcogenides (TMDs) available in the literature. We inspect the construction of these TB models, discuss their parameterization used and conduct a thorough comparison of their effectiveness in capturing important electronic properties. Based on these insights, we propose a novel TB model for TMDs designed for enhanced computational efficiency. Utilizing MoS2 as a representative case, we explain why specific models offer a more accurate description. Our primary aim is to assist researchers in choosing the most appropriate TB model for their calculations on TMDs.

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: When using hexagonal boron-nitride (hBN) as a substrate for graphene, the resulting moire pattern creates secondary Dirac points. By encapsulating a multilayer graphene within aligned hBN sheets the controlled moire stacking may offer even richer benefits. Using advanced tight-binding simulations on atomistically-relaxed heterostructures, here we show that the gap at the secondary Dirac point can be opened in selected moire-stacking configurations, and is independent of any additional vertical gating of the heterostructure. On the other hand, gating can broadly tune the gap at the principal Dirac point, and may thereby strongly compress the first moire mini-band in width against the moire-induced gap at the secondary Dirac point. We reveal that in hBN-encapsulated bilayer graphene this novel mechanism can lead to isolated bands flatter than 10 meV under moderate gating, hence presenting a convenient pathway towards electronically-controlled strongly-correlated states on demand.

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: 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: The flexoelectric and electronic properties of zig-zag graphene nanoribbons are explored under mechanical bending using state of the art first principles calculations. A linear dependence of the bending induced out of plane polarization on the applied strain gradient is found. The inferior flexoelectric properties of graphene nanoribbons can be improved by more than two orders of magnitude by hydrogen and fluorine functionalization (CH and CF nanoribbons). A large out of plane flexoelectric effect is predicted for CF nanoribbons. The origin of this enhancement lies in the electro-negativity difference between carbon and fluorine atoms, which breaks the out of plane charge symmetry even for a small strain gradient. The flexoelectric effect can be further improved by co-functionalization with hydrogen and fluorine (CHF Janus-type nanoribbon), where a spontaneous out of plane dipole moment is formed even for flat nanoribbons. We also find that bending can control the charge localization of valence band maxima and therefore enables the tuning of the hole effective masses and band gaps. These results present an important advance towards the understanding of flexoelectric and electronic properties of hydrogen and fluorine functionalized graphene nanoribbons, which can have important implications for flexible electronic applications. (C) 2020 Elsevier Ltd. All rights reserved.

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