Abstract: The alpha-T3 lattice model interpolates a honeycomb (graphene-like) lattice and a T3 (also known as dice) lattice via the parameter alpha. These lattices are made up of three atoms per unit cell. This gives rise to an additional dispersionless flat band touching the conduction and valence bands. Electrons in this model are analogous to Dirac fermions with an enlarged pseudospin, which provides unusual tunneling features like omnidirectional Klein tunneling, also called super-Klein tunneling (SKT). However, it is unknown how small deviations in the equivalence between the atomic sites, i.e., variations in the alpha parameter, and the number of tunnel barriers changes the transmission properties. Moreover, it is interesting to learn how tunneling occurs through regions where the energy spectrum changes from linear with a middle flat band to a hyperbolic dispersion. In this paper we investigate these properties, its dependence on the number of square barriers and the alpha parameter for either gapped and gapless cases. Furthermore, we compare these results to the case where electrons tunnel from a region with linear dispersion to a region with a bandgap. In the latter case, contrary to tunneling through a potential barrier, the SKT is no longer observed. Finally, we find specific cases where transmission is allowed due to a symmetry breaking of sublattice equivalence.

Abstract: The moire pattern observed in stacked noncommensurate crystal lattices, such as heterobilayers of transition metal dichalcogenides, produces a periodic modulation of their band gap. Excitons subjected to this potential landscape exhibit a band structure that gives rise to a quasiparticle dubbed the moire exciton. In the case of MoS2/WSe2 heterobilayers, the moire trapping potential has honeycomb symmetry and, consequently, the moire exciton band structure is the same as that of a Dirac-Weyl fermion, whose mass can be further tuned down to zero with a perpendicularly applied field. Here we show that, analogously to other Dirac-like particles, the moire exciton exhibits a trembling motion, also known as Zitterbewegung, whose long timescales are compatible with current experimental techniques for exciton dynamics. This promotes the study of the dynamics of moire excitons in van der Waals heterostructures as an advantageous solid-state platform to probe Zitterbewegung, broadly tunable by gating and interlayer twist angle.

Abstract: The time evolution of a low-energy two-dimensional Gaussian wave packet in ABC-stacked n-layer graphene (ABC-NLG) is investigated. Expectation values of the position (x, y) of center-of-mass and the total probability densities of the wave packet are calculated analytically using the Green's function method. These results are confirmed using an alternative numerical method based on the split-operator technique within the Dirac approach for ABC-NLG, which additionally allows to include external fields and potentials. The main features of the zitterbewegung (trembling motion) of wave packets in graphene are demonstrated and are found to depend not only on the wave packet width and initial pseudospin polarization, but also on the number of layers. Moreover, the analytical and numerical methods proposed here allow to investigate wave packet dynamics in graphene systems with an arbitrary number of layers and arbitrary potential landscapes.

Abstract: The energy spectrum and local current patterns in graphene quantum dots (QD) are investigated for different geometries in the presence of an external perpendicular magnetic field. Our results demonstrate that, for specific geometries and edge configurations, the QD exhibits vortex and anti-vortex patterns in the local current density, in close analogy to the vortex patterns observed in the probability density current of semiconductor QD, as well as in the order parameter of mesoscopic superconductors.

Abstract: Using the tight-binding approach coupled with mean-field Hubbard model, we theoretically study the effect of mechanical deformations on the magnetic properties of bilayer graphene (BLG) quantum dots (QDs). Results are obtained for AA-and AB(Bernal)-stacked BLG QDs, considering different geometries (hexagonal, triangular and square shapes) and edge types (armchair and zigzag edges). In the absence of strain, our results show that (i) the magnetization is affected by taking different dot sizes only for hexagonal BLG QDs with zigzag edges, exhibiting different critical Hubbard interactions, and (ii) the magnetization does not depend on the interlayer hopping energies, except for the geometries with zigzag edges and AA stacking. In the presence of in-plane and uniaxial strain, for all geometries we obtain two different magnetization regimes depending on the applied strain amplitude. The appearance of such different regimes is due to the breaking of layer and sublattice symmetries in BLG QDs.

Abstract: The interaction of monolayer graphene with specific substrates may break its sublattice symmetry and results in unidirectional chiral states with opposite group velocities in the different Dirac cones (Zarenia et al 2012 Phys. Rev. B 86 085451). Taking advantage of this feature, we propose a valley filter based on a transversal mass kink for low energy electrons in graphene, which is obtained by assuming a defect region in the substrate that provides a change in the sign of the substrate-induced mass and thus creates a non-biased channel, perpendicular to the kink, for electron motion. By solving the time-dependent Schrodinger equation for the tight-binding Hamiltonian, we investigate the time evolution of a Gaussian wave packet propagating through such a system and obtain the transport properties of this graphene-based substrate-induced quantum point contact. Our results demonstrate that efficient valley filtering can be obtained, provided: (i) the electron energy is sufficiently low, i.e. with electrons belonging mostly to the lowest sub-band of the channel, and (ii) the channel length (width) is sufficiently long (narrow). Moreover, even though the transmission probabilities for each valley are significantly affected by impurities and defects in the channel region, the valley polarization in this system is shown to be robust against their presence.

Abstract: A pseudomagnetic field kink can be realized along a graphene nanoribbon using strain engineering. Electron transport along this kink is governed by snake states that are characterized by a single propagation direction. Those pseudomagnetic fields point towards opposite directions in the K and K' valleys, leading to valley polarized snake states. In a graphene nanoribbon with armchair edges this effect results in a valley filter that is based only on strain engineering. We discuss how to maximize this valley filtering by adjusting the parameters that define the stress distribution along the graphene ribbon.

Abstract: Using the tight-binding approach, we study the influence of a perpendicular magnetic field on the energy levels of hexagonal, triangular, and circular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We obtain the energy levels for AB (Bernal)-stacked BLG QDs in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). We find different regions in the spectrum of biased QDs with respect to the crossing point between the lowest-electron and -hole Landau levels of a biased BLG sheet. Those different regions correspond to electron states that are localized at the center, edge, or corner of the BLG QD. Quantum Hall corner states are found to be absent in circular BLG QDs. The spatial symmetry of the carrier density distribution is related to the symmetry of the confinement potential, the position of zigzag edges, and the presence or absence of interlayer inversion symmetry.

Abstract: Propagation of an electron wave packet through a quantum point contact (QPC) defined by electrostatic gates in bilayer graphene is investigated. The gates provide a bias between the layers, in order to produce an energy gap. If the gates on both sides of the contact produce the same bias, steps in the electron transmission probability are observed, as in the usual QPC. However, if the bias is inverted on one of the sides of the QPC, only electrons belonging to one of the Dirac valleys are allowed to pass, which provides a very efficient valley filtering.

Abstract: Using the four-band continuum model we derive a general expression for the infinite-mass boundary condition in bilayer graphene. Applying this new boundary condition we analytically calculate the confined states and the corresponding wave functions in a bilayer graphene quantum dot in the absence and presence of a perpendicular magnetic field. Our results for the energy spectrum show an energy gap between the electron and hole states at small magnetic fields. Furthermore the electron (e) and hole (h) energy levels corresponding to the K and K' valleys exhibit the E-K(e(h)) (m) = E-K'(e(h)) (m) symmetry, where m is the angular momentum quantum number. (C) 2014 Elsevier Ltd. All rights reserved.

Abstract: The scattering of a Gaussian wave packet in armchair and zigzag graphene edges is theoretically investigated by numerically solving the time-dependent Schrodinger equation for the tight-binding model Hamiltonian. Our theory allows us to investigate scattering in reciprocal space, and depending on the type of graphene edge we observe scattering within the same valley, or between different valleys. In the presence of an external magnetic field, the well-known skipping orbits are observed. However, our results demonstrate that in the case of a pseudomagnetic field, induced by nonuniform strain, the scattering by an armchair edge results in a nonpropagating edge state.