Abstract: Twisted graphene layers produce a moire pattern (MP) structure with a predetermined wavelength for a given twist angle. However, predicting the membrane corrugation amplitude for any angle other than pure AB-stacked or AA-stacked graphene is impossible using first-principles density functional theory (DFT) due to the large supercell. Here, within elasticity theory, we define the MP structure as the minimum-energy configuration, thereby leaving the height amplitude as the only unknown parameter. The latter is determined from DFT calculations for AB-and AA-stacked bilayer graphene in order to eliminate all fitting parameters. Excellent agreement with scanning tunneling microscopy results across multiple substrates is reported as a function of twist angle.

Abstract: In honeycomb Dirac systems with broken inversion symmetry, orbital magnetic moments coupled to the valley degree of freedom arise due to the topology of the band structure, leading to valley-selective optical dichroism. On the other hand, in Dirac systems with prominent spin-orbit coupling, similar orbital magnetic moments emerge as well. These moments are coupled to spin, but otherwise have the same functional form as the moments stemming from spatial inversion breaking. After reviewing the basic properties of these moments, which are relevant for a whole set of newly discovered materials, such as silicene and germanene, we study the particular impact that these moments have on graphene nanoengineered barriers with artificially enhanced spin-orbit coupling. We examine transmission properties of such barriers in the presence of a magnetic field. The orbital moments are found to manifest in transport characteristics through spin-dependent transmission and conductance, making them directly accessible in experiments. Moreover, the Zeeman-type effects appear without explicitly incorporating the Zeeman term in the models, i.e., by using minimal coupling and Peierls substitution in continuum and the tight-binding methods, respectively. We find that a quasiclassical view is able to explain all the observed phenomena.

Abstract: Quantum thermal transport in armchair and zigzag graphene nanoribbons is investigated in the presence of single atomic vacancies and subject to different boundary conditions. We start with a full comparison of the phonon polarizations and energy dispersions as given by a fifth-nearest-neighbor force-constant model (5NNFCM) and by elasticity theory of continuum membranes (ETCM). For free-edge ribbons, we discuss the behavior of an additional acoustic edge-localized flexural mode, known as fourth acoustic branch (4ZA), which has a small gap when it is obtained by the 5NNFCM. Then, we show that ribbons with supported edges have a sample-size dependent energy gap in the phonon spectrum which is particularly large for in-plane modes. Irrespective to the calculation method and the boundary condition, the dependence of the energy gap for the low-energy optical phonon modes against the ribbon width W is found to be proportional to 1/W for in-plane, and 1/W-2 for out-of-plane phonon modes. Using the 5NNFCM, the ballistic thermal conductance and its contributions from every single phonon mode are then obtained by the nonequilibrium Green's function technique. We found that, while edge and central localized single atomic vacancies do not affect the low-energy transmission function of in-plane phonon modes, they reduce considerably the contributions of the flexural modes. On the other hand, in-plane modes contributions are strongly dependent on the boundary conditions and at low temperatures can be highly reduced in supported-edge samples. These findings could open a route to engineer graphene based devices where it is possible to discriminate the relative contribution of polarized phonons and to tune the thermal transport on the nanoscale.

Abstract: The electronic properties of the two-dimensional material silicene are strongly influenced by the application of a perpendicular electric field E-z and of an exchange field M due to adatoms positioned on the surface or a ferromagnetic substrate. Within the random phase approximation, we investigate how electron-electron interactions are affected by these fields and present analytical and numerical results for the dispersion of plasmons, their lifetime, and their oscillator strength. We find that the combination of the fields E-z and M brings a spin and valley texture to the particle-hole excitation spectrum and allows the formation of spin-and valley-polarized plasmons. When the Fermi level lies in the gap of one spin in one valley, the intraband region of the corresponding spectrum disappears. For zero E-z and finite M the spin symmetry is broken and spin polarization is possible. The lifetime and oscillator strength of the plasmons are shown to depend strongly on the number of spin and valley type electrons that form the electron-hole pairs.

Abstract: Using first principles calculations we showed that the electronic and optical properties of single-layer black phosphorus (BP) depend strongly on the applied strain. Due to the strong anisotropic atomic structure of BP, its electronic conductivity and optical response are sensitive to the magnitude and the orientation of the applied strain. We found that the inclusion of many body effects is essential for the correct description of the electronic properties of monolayer BP; for example, while the electronic gap of strainless BP is found to be 0.90 eV by using semilocal functionals, it becomes 2.31 eV when many-body effects are taken into account within the G(0)W(0) scheme. Applied tensile strain was shown to significantly enhance electron transport along zigzag direction of BP. Furthermore, biaxial strain is able to tune the optical band gap of monolayer BP from 0.38 eV (at -8% strain) to 2.07 eV (at 5.5%). The exciton binding energy is also sensitive to the magnitude of the applied strain. It is found to be 0.40 eV for compressive biaxial strain of -8%, and it becomes 0.83 eV for tensile strain of 4%. Our calculations demonstrate that the optical response of BP can be significantly tuned using strain engineering which appears as a promising way to design novel photovoltaic devices that capture a broad range of solar spectrum.

Abstract: We study the energy and temperature relaxation of electrons in graphene on a piezoelectric substrate. Scattering from the combined potential of extrinsic piezoelectric surface acoustical (PA) phonons of the substrate and intrinsic deformation acoustical phonons of graphene is considered for a (non) degenerate gas of Dirac fermions. It is shown that in the regime of low energies or temperatures the PA phonons dominate the relaxation and change qualitatively its character. This prediction is relevant for quantum metrology and electronic applications using graphene devices and suggests an experimental setup for probing electron-phonon coupling in graphene.

Abstract: We present a systematic study of the energy spectra of graphene quantum rings having different geometries and edge types in the presence of a perpendicular magnetic field. Results are obtained within the tight-binding (TB) and Dirac models and we discuss which features of the former can be recovered by using the approximations imposed by the latter. Energy levels of graphene quantum rings obtained by diagonalizing the TB Hamiltonian are demonstrated to be strongly dependent on the rings geometry and the microscopical structure of the edges. This makes it difficult to recover those spectra by the existing theories that are based on the continuum (Dirac) model. Nevertheless, our results show that both approaches (i.e., TB and Dirac model) may provide similar results, but only for very specific combinations of ring geometry and edge types. The results obtained by a simplified model describing an infinitely thin circular Dirac ring show good agreement with those obtained for hexagonal and rhombus armchair graphene rings within the TB model. Moreover, we show that the energy levels of a circular quantum ring with an infinite mass boundary condition obtained within the Dirac model agree with those for a ring defined by a ring-shaped staggered potential obtained within the TB model.

Abstract: We present a theoretical study of the plasmon and coupled plasmon-phonon modes induced by intraband electron-electron interaction in graphene in the presence of driving dc electric field. We find that the electric field dependence of these collective excitation modes in graphene differs significantly from that in a conventional two-dimensional electron gas with a parabolic energy spectrum. This is due mainly to the fact that graphene has a linear energy spectrum and the Fermi velocity of electrons in graphene is much larger than the drift velocity of electrons. The obtained results demonstrate that the plasmon and coupled plasmon-phonon modes in graphene can be tuned by applying not only the gate voltage but also the source-to-drain field. The manipulation of plasmon and coupled plasmon-phonon modes by source-to-drain voltage can let graphene be more conveniently applied as an advanced plasmonic material.

Abstract: Diffusive properties of interacting magnetic dipoles confined in a parabolic narrow channel and in the presence of a periodic modulated (corrugated) potential along the unconfined direction are studied using Brownian dynamics simulations. We compare our simulation results with the analytical result for the effective diffusion coefficient of a single particle by Festa and d'Agliano [Physica A 90, 229 (1978)] and show the importance of interparticle interaction on the diffusion process. We present results for the diffusion of magnetic dipoles as a function of linear density, strength of the periodic modulation and commensurability factor.

Abstract: We calculate the spectrum of states that are localized at the interface between two regions of opposite bias in trilayer graphene (TLG). These potential profiles, also known as potential kinks, have been predicted to support two different branches of localized states for the case of bilayer graphene, and show similarities to the surface states of topological insulators. On the other hand, we found that ABC stacked TLG exhibits three different unidimensional branches of states in each valley that are confined to the kink interface. They have the property E(k(y)) = -E(-k(y)) when belonging to the same valley and E-K(k(y)) = -E-K' (-k(y)). A kink-antikink potential profile opens a gap in the spectrum of these one-dimensional states.

Abstract: We modeled the electron and hole states in Si/SiO2 quantum wells within a basis of standing waves using the 30-band k . p theory. The hard-wall confinement potential is assumed, and the influence of the peculiar band structure of bulk silicon on the quantum-well sub-bands is explored. Numerous spurious solutions in the conduction-band and valence-band energy spectra are found and are identified to be of two types: (1) spurious states which have large contributions of the bulk solutions with large wave vectors (the high-k spurious solutions) and (2) states which originate mainly from the spurious valley outside the Brillouin zone (the extravalley spurious solutions). An algorithm to remove all those nonphysical solutions from the electron and hole energy spectra is proposed. Furthermore, slow and oscillatory convergence of the hole energy levels with the number of basis functions is found and is explained by the peculiar band mixing and the confinement in the considered quantum well. We discovered that assuming the hard-wall potential leads to numerical instability of the hole states computation. Nonetheless, allowing the envelope functions to exponentially decay in a barrier of finite height is found to improve the accuracy of the computed hole states.

Abstract: We theoretically demonstrate that the transport inefficiency recently found experimentally for branched-out mesoscopic networks can also be observed in a quantum ring of finite width with an attached central horizontal branch. This is done by investigating the time evolution of an electron wave packet in such a system. Our numerical results show that the conductivity of the ring does not necessary improve if one adds an extra channel. This ensures that there exists a quantum analog of the Braess paradox, originating from quantum scattering and interference.

Abstract: We investigate the effect of Rashba and intrinsic spin-orbit couplings on the electronic properties and spin configurations of Dirac fermions confined in: (i) a flat graphene sheet, (ii) a graphene wire with p-n-p structure, and (iii) a superlattice of graphene wires. The interplay between the spin-orbit interaction mechanisms breaks the electron-hole symmetry and the spin configuration induced by Rashba spin-orbit coupling lacks inversion symmetry in k space. We show that the Rashba spin-orbit interaction doubles the Fabry-Perot resonant modes in the transmission spectrum of a graphene wire and opens new channels for the electron transmission. Moreover, it leads to the appearance of spin split extra Dirac cones in the energy spectrum of a graphene superlattice. It is shown that the spin of the electrons and holes confined in a flat graphene sheet is always perpendicular to their motion while this is not the case for the other nanostructures.