Publications

Thin adhesive films can be removed from substrates, torn, and folded in distinct geometries under external driving forces. In two-dimensional materials, however, these processes can be self-driven as shown in previous studies on folded twisted bilayer graphene nanoribbons produced by spontaneous tearing and peeling from a substrate. Here, we use atomic force microscopy techniques to generate and characterize the geometrical structure of naturally self-grown folded nanoribbon structures. Measurements of nanoribbon width and interlayer separation reveal similar twist-angle dependences possibly caused by the anisotropy in the bilayer potential. In addition, analysis of the data shows an unexpected correlation between the height of the folded arc edge—parameterized by a radius R—, and the ribbon width, suggestive of a self-growth process driven by a variable cross-sectional shape. These observations are well described by an energy minimization model that includes the bilayer adhesion energy density as represented by a distance dependent Morse potential. We obtain an analytical expression for the radius R versus the ribbon width that predicts a renormalized bending rigidity and stands in good agreement with experimental observations. The newly found relation between these geometrical parameters suggests a mechanism for tailored growth of folded twisted bilayer graphene- a platform for many intriguing physics phenomena.

One of the intriguing characteristics of honeycomb lattices is the appearance of a pseudomagnetic field as a result of mechanical deformation. In the case of graphene, the Landau quantization resulting from this pseudomagnetic field has been measured using scanning tunneling microscopy. Here we show that a signature of the pseudomagnetic field is a local sublattice symmetry breaking observable as a redistribution of the local density of states. This can be interpreted as a polarization of graphene’s pseudospin due to a strain induced pseudomagnetic field, in analogy to the alignment of a real spin in a magnetic field. We reveal this sublattice symmetry breaking by tunably straining graphene using the tip of a scanning tunneling microscope. The tip locally lifts the graphene membrane from a SiO ₂  support, as visible by an increased slope of the  I ( z ) curves. The amount of lifting is consistent with molecular dynamics calculations, which reveal a deformed graphene area under the tip in the shape of a Gaussian. The pseudomagnetic field induced by the deformation becomes visible as a sublattice symmetry breaking which scales with the lifting height of the strained deformation and therefore with the pseudomagnetic field strength. Its magnitude is quantitatively reproduced by analytic and tight-binding models, revealing fields of 1000 T. These results might be the starting point for an effective THz valley filter, as a basic element of valleytronics.

An Anderson model for a magnetic impurity in a two-dimensional electron gas with bulk Rashba spin-orbit interaction is solved using the numerical renormalization group under two different experimental scenarios. For a fixed Fermi energy, the Kondo temperature 𝑇𝐾 varies weakly with Rashba coupling 𝜆𝑅, as reported previously. If instead the band filling is low and held constant, increasing 𝜆𝑅 can drive the system into a helical regime with exponential enhancement of 𝑇𝐾. Under either scenario, thermodynamic properties at low temperatures 𝑇 exhibit the same dependencies on 𝑇/𝑇𝐾 as are found for 𝜆𝑅=0. Unlike the conventional Kondo effect, however, the impurity exhibits static spin correlations with conduction electrons of nonzero orbital angular momentum about the impurity site. We also consider a magnetic field that Zeeman splits the conduction band but not the impurity level, an effective picture that arises under a proposed route to access the helical regime in a driven system. The impurity contribution to the system's ground-state angular momentum is found to be a universal function of the ratio of the Zeeman energy to a temperature scale that is not 𝑇𝐾 (as would be the case in a magnetic field that couples directly to the impurity spin), but rather is proportional to 𝑇𝐾 divided by the impurity hybridization width. This universal scaling is explained via a perturbative treatment of field-induced changes in the electronic density of states.

We study the Ruderman-Kittel-Kasuya-Yosida interaction in three-dimensional Dirac semimetals. Using retarded Green's functions in real space, we obtain and analyze asymptotic expressions for the interaction, with magnetic impurities at different distances and relative angle with respect to high symmetry directions on the lattice. We show that the Fermi velocity anisotropy in these materials produces a strong renormalization of the magnitude of the interaction, as well as a correction to the frequency of oscillation in real space. Hybridization of the impurities to different conduction electron orbitals are shown to result in interesting anisotropic spin-spin interactions which can generate spiral spin structures in doped samples.

Deformations in graphene systems are central elements in the novel field of  straintronics . Various strain geometries have been proposed to produce specific properties, but their experimental realization has been limited. Because strained folds can be engineered on graphene samples on appropriate substrates, we study their effects on graphene transport properties. We show the existence of an enhanced local density of states (LDOS) along the direction of the strained fold that originates from localization of higher energy states and provides extra conductance channels at lower energies. In addition to exhibiting sublattice symmetry breaking, these states are valley polarized, with quasiballistic properties in smooth disorder potentials. We confirmed that these results persist in the presence of strong edge disorder, making these geometries viable electronic waveguides. These findings could be tested in properly engineered experimental settings.

We study transport properties of hexagonal zigzag graphene quantum rings connected to semi-infinite nanoribbons. Open two-fold symmetric structures support localized states that can be traced back to those existing in the isolated six-fold symmetric rings. Using a tight-binding Hamiltonian within the Green's function formalism, we show that an external magnetic field promotes these localized states to Fano resonances with robust signatures in transport. Local density of states and current distributions of the resonant states are calculated as a function of the magnetic flux intensity. For structures on corrugated substrates we analyze the effect of strain by including an out-of-plane centro-symmetric deformation in the model. We show that small strains shift the resonance positions without further modifications, while high strains introduce new ones.

We calculate the local density of states (LDOS) for an infinite graphene sheet with a single centrosymmetric out-of-plane deformation, in order to investigate measurable strain signatures on graphene. We focus on the regime of small deformations and show that the strain-induced pseudomagnetic field induces an imbalance of the LDOS between the two triangular graphene sublattices in the region of the deformation. Real-space imaging reveals a characteristic sixfold symmetry pattern where the sublattice symmetry is broken within each fold, consistent with experimental and tight-binding observations. The open geometry we study allows us to make use of the usual continuum model of graphene and to obtain results independent of boundary conditions. We provide an analytic perturbative expression for the contrast between the LDOS of each sublattice, showing a scaling law as a function of the amplitude and width of the deformation. We confirm our results by a numerically exact iterative scattering matrix method.

We propose a quantum electronic device based on a strained graphene nanoribbon. Mechanical strain, internal exchange field, and spin-orbit couplings (SOCs) have been exploited as principle parameters to tune physical properties of the device. We predict a remarkable zero-field topological quantum phase transition between the time-reversal-symmetry-broken quantum spin Hall (QSH) and quantum anomalous Hall states, which was previously thought to take place only in the presence of finite magnetic field. We illustrate as intrinsic SOC is tuned, how two different helicity edge states located in the opposite edges of the nanoribbon exchange their locations. Our results indicate that the pseudomagnetic field induced by the strain could be coupled to the spin degrees of freedom through the SOC responsible for the stability of a QSH state. The controllability of this zero-field phase transition with strength and direction of the strain is also demonstrated. Our prediction offers a tempting prospect of strain, electric, and magnetic manipulation of the QSH effect.

We study the Kondo screening of a magnetic impurity adsorbed in graphene in the presence of Rashba spin-orbit interaction. The system is described by an effective single-channel Anderson impurity model, which we analyze using the numerical renormalization group. The nontrivial energy dependence of the host density of states gives rise to interesting behaviors under variation of the chemical potential or the spin-orbit coupling. Varying the Rashba coupling produces strong changes in the Kondo temperature characterizing the many-body screening of the impurity spin, and at half filling allows an approach to a quantum phase transition separating the strong-coupling Kondo phase from a free-moment phase. Tuning the chemical potential close to sharp features of the hybridization function results in striking features in the temperature dependencies of thermodynamic quantities and in the frequency dependence of the impurity spectral function.

We study the Ruderman-Kittel-Kasuya-Yosida interaction in 𝑝-doped transition-metal dichalcogenides such as MoS2and WS2. We consider magnetic impurities hybridized to the Mo 𝑑 orbitals characteristic of the valence bands. Using the Matsubara Green's-function formalism, we obtain the two-impurity interaction vs their separation and chemical potential of the system, accounting for the important angular dependence which reflects the underlying triangular lattice symmetry. The inclusion of the valence-band valley at the 𝛤 point results in a strong enhancement of the interaction. Electron-scattering processes transferring momentum between valleys at different symmetry points give rise to complex spatial oscillation patterns. Variable doping would allow the exploration of rather interesting behavior in the interaction of magnetic impurities on the surfaces of these materials, including the control of the interaction symmetry, which can be directly probed in scanning tunneling microscopy experiments.