Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “twisted bilayer graphene”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Spin magnetometry as a probe of stripe superconductivity in twisted bilayer graphene

The discovery of alternating superconducting and insulating ground states in magic angle graphene has suggested an intriguing analogy with cuprate high-T c materials. In this paper we argue that the network states of small angle twisted bilayer graphene (TBG) afford a further perspective on the cuprates by emulating their stripe-ordered phases, as in La 1.875 Ba 0.125 CuO 4 . We show that the spin and valley quantum numbers of stripes in TBG graphene fractionalize, developing characteristic signatures in the tunneling density of states and the magnetic noise spectrum of impurity spins. By examining the coupling between the charge rivers we determine the superconducting transition temperature. Our study suggests that twisted bilayer graphene can be used for a controlled emulation of stripe superconductivity and quantum sensing experiments of emergent anyonic excitations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Moiré commensurability and the quantum anomalous Hall effect in twisted bilayer graphene on hexagonal boron nitride

The quantum anomalous Hall (QAH) effect is sometimes observed in twisted bilayer graphene (tBG) when it is nearly aligned with an encapsulating hexagonal boron nitride (hBN) layer. We propose that the appearance or absence of the QAH effect in individual devices could be related to commensurability between the graphene/graphene and graphene/hBN moiré patterns. We identify a series of points in the (θ GG , θ GBN ) twist-angle space at which the two moiré patterns are commensurate, allowing moiré band theory to be applied, and we show that the band Chern numbers are in this case sensitive to a rigid in-plane hBN displacement. Given this property, we argue that the QAH effect is likely only when (i) the (θ GG , θ GBN ) twist-angle-pair is close enough to a commensurate point that the two moiré patterns yield a supermoiré pattern with a sufficiently long length scale, and (ii) the supermoiré has a percolating topologically nontrivial QAH phase. For twist angles far from commensurability, the hBN layer acts as a source of disorder that can destroy the QAH effect. Our proposal can explain a number of current experimental observations. Finally, further experimental studies that can test this proposal more directly are suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Top-Down Fabrication of Atomic Patterns in Twisted Bilayer Graphene

Atomic-scale engineering typically involves bottom-up approaches, leveraging parameters such as temperature, partial pressures, and chemical affinity to promote spontaneous arrangement of atoms. These parameters are applied globally, resulting in atomic-scale features scattered probabilistically throughout the material. In a top-down approach, different regions of the material are exposed to different parameters, resulting in structural changes varying on the scale of the resolution. In this work, the application of global and local parameters is combined in an aberration-corrected scanning transmission electron microscope (STEM) to demonstrate atomic-scale precision patterning of atoms in twisted bilayer graphene. The focused electron beam is used to define attachment points for foreign atoms through the controlled ejection of carbon atoms from the graphene lattice. The sample environment is staged with nearby source materials such that the sample temperature can induce migration of the source atoms across the sample surface. Under these conditions, the electron-beam (top-down) enables carbon atoms in the graphene to be replaced spontaneously by diffusing adatoms (bottom-up). Using image-based feedback control, arbitrary patterns of atoms and atom clusters are attached to the twisted bilayer graphene with limited human interaction. Finally, the role of substrate temperature on adatom and vacancy diffusion is explored by first-principles simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spin-orbit–driven ferromagnetism at half moiré filling in magic-angle twisted bilayer graphene

Strong electron correlation and spin-orbit coupling (SOC) can have a profound influence on the electronic properties of materials. We examined their combined influence on a two-dimensional electronic system at the atomic interface between magic-angle twisted bilayer graphene and a tungsten diselenide crystal. We found that strong electron correlation within the moiré flatband stabilizes correlated insulating states at both quarter and half filling, and that SOC transforms these Mott-like insulators into ferromagnets, as evidenced by a robust anomalous Hall effect with hysteretic switching behavior. The coupling between spin and valley degrees of freedom could be demonstrated through control of the magnetic order with an in-plane magnetic field or a perpendicular electric field. Our findings establish an experimental knob to engineer topological properties of moiré bands in twisted bilayer graphene and related systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate tight-binding model for twisted bilayer graphene describes topological flat bands without geometric relaxation

A major hurdle in understanding the phase diagram of twisted bilayer graphene is the roles of lattice relaxation and electronic structure on isolated band flattening near magic twist angles. Here in this work, the authors develop an accurate local environment tight-binding model fit to tight-binding parameters computed from ab initio density-functional theory calculations across many atomic configurations. With the accurate parametrization, it is found that the magic angle shifts to slightly lower angles than often quoted, from around 1.05° to around 0.99°, and that isolated flat bands appear for rigidly rotated graphene layers, with enhancement of the flat bands when the layers are allowed to distort. Study of the orbital localization supports the emergence of fragile topology in the isolated flat bands without the need for lattice relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cascades between Light and Heavy Fermions in the Normal State of Magic-Angle Twisted Bilayer Graphene

We present a framework for understanding the cascade transitions and the Landau level degeneracies of twisted bilayer graphene. The Coulomb interaction projected onto narrow bands causes the charged excitations at an integer filling to disperse, forming new bands. If the excitation moves the filling away from the charge neutrality point, then it has a band minimum at the moiré Brillouin zone center with a small mass that compares well with the experiment; if towards the charge neutrality point, then it has a much larger mass and a higher degeneracy. At a nonzero density away from an integer filling the excitations interact. The system on the small mass side has a large bandwidth and forms a Fermi liquid. On the large mass side the bandwidth is narrow, the compressibility is negative and the Fermi liquid is likely unstable. This explains the observed sawtooth features in compressibility, the Landau fans pointing away from charge neutrality and their degeneracies. The framework sets the stage for superconductivity at lower temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topological and nematic superconductivity mediated by ferro-SU(4) fluctuations in twisted bilayer graphene

We propose an SU(4) spin-valley-fermion model to investigate the superconducting instabilities of twisted bilayer graphene (TBG). In this approach, bosonic fluctuations associated with an emergent SU(4) symmetry, corresponding to combined rotations in valley and spin spaces, couple to the low-energy fermions that comprise the flat bands. These fluctuations are peaked at zero wave vector, reflecting the “ferromagnetic-like” SU(4) ground state recently found in strong-coupling solutions of microscopic models for TBG. Focusing on electronic states related to symmetry-imposed points of the Fermi surface, dubbed here “valley hot-spots” and “van Hove hot-spots,” we find that the coupling to the itinerant electrons partially lifts the huge degeneracy of the ferro-SU(4) ground state manifold, favoring intervalley order, spin-valley coupled order, ferromagnetic order, spin-current order, and valley-polarized order, depending on details of the band structure. Here, these fluctuations, in turn, promote attractive pairing interactions in a variety of closely competing channels, including a nodeless f-wave state, a nodal i-wave state, and topological d+id and p+ip states with unusual Chern numbers 2 and 4, respectively. Nematic superconductivity, although not realized as a primary instability of the system, appears as a consequence of the near-degeneracy of superconducting order parameters that transform as one-dimensional and two-dimensional irreducible representations of the point group D 6 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Electron– K -phonon interaction in twisted bilayer graphene

We develop an analytic theory to describe the interaction between electrons and K phonons and study its influence on superconductivity in the bare bands of twisted bilayer graphene (TBG). We find that, due to symmetry and the two-center approximation, only one optical K phonon (160meV) of graphene is responsible for the intervalley electron-phonon interaction. Here, this phonon has recently been found in angular-resolved photoemission spectroscopy to be responsible for replicas of the TBG flat bands. By projecting the interaction to the TBG flat bands, we perform the full symmetry analysis of the phonon-mediated attractive interaction and pairing channels in the Chern basis, and show that several channels are guaranteed to have gapless order parameters. From the linearized gap equations, we find that the highest T c pairing induced by this phonon is a singlet gapped s-wave inter-Chern-band order parameter, followed closely by a gapless nematic d-wave intra-Chern-band order parameter. We justify these results analytically, using the topological heavy-fermion mapping of TBG which has allowed us to obtain an analytic form of a phonon-mediated attractive interaction and to analytically solve the linearized and T = 0 gap equations. For the intra-Chern-band channel, the nematic state with nodes is shown to be stabilized in the chiral flat-band limit. While the flat-band Coulomb interaction can be screened sufficiently enough—around the Van Hove singularities—to allow for electron-phonon based superconductivity, it is unlikely that this effect can be maintained in the lower density of states excitation bands around the correlated insulator states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Symmetric Kondo Lattice States in Doped Strained Twisted Bilayer Graphene

We use the topological heavy fermion (THF) model and its Kondo lattice (KL) formulation to study the possibility of a symmetric Kondo (SK) state in twisted bilayer graphene. Via a large-N approximation, we find a SK state in the KL model at fillings v = 0,±1,±2 where a KL model can be constructed. In the SK state, all symmetries are preserved and the local moments are Kondo screened by the conduction electrons. At the mean-field level of the THF model at v = 0,±1,±2,±3 we also find a similar symmetric state that is adiabatically connected to the symmetric Kondo state. We study the stability of the symmetric state by comparing its energy with the ordered (symmetry-breaking) states found in [H. Hu et al., Phys. Rev. Lett. 131, 026502 (2023)., Z.-D. Song and B. A. Bernevig, Phys. Rev. Lett. 129, 047601 (2022).] and find the ordered states to have lower energy at v = 0,±1,±2. However, moving away from integer fillings by doping the light bands, our mean-field calculations find the energy difference between the ordered state and the symmetric state to be reduced, which suggests the loss of ordering and a tendency toward Kondo screening. In order to include many-body effects beyond the mean-field approximation, we also performed dynamical mean-field theory calculations on the THF model in the nonordered phase. The spin susceptibility follows a Curie behavior at v = 0,±1,±2 down to ~2 K where the onset of screening of the local moment becomes visible. This hints to very low Kondo temperatures at these fillings, in agreement with the outcome of our mean-field calculations. At noninteger filling v = ±0.5,±0.8,±1.2 dynamical mean-field theory shows deviations from a 1/T susceptibility at much higher temperatures, suggesting a more effective screening of local moments with doping. Lastly, we study the effect of a C 3⁢z -rotational-symmetry-breaking strain via mean-field approaches and find that a symmetric phase (that only breaks C 3⁢z symmetry) can be stabilized at sufficiently large strain at v = 0,±1,±2. Our results suggest that a symmetric Kondo phase is strongly suppressed at integer fillings, but could be stabilized either at noninteger fillings or by applying strain.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Critical magnetic fields and electron pairing in magic-angle twisted bilayer graphene

The velocities of the quasiparticles that form Cooper pairs in a superconductor are revealed by the upper critical magnetic field. Here we use this property to assess superconductivity in magic-angle twisted bilayer graphene (MATBG), which has been observed over a range of moiré band-filling, twist angle, and screening environment conditions. An average Fermi velocity can be defined as $v$$^{*}_{F}$ ≡ k B T c ℓ c /ℏ, where T c and ℓ c are the critical temperature and magnetic length, respectively. An advantage of this definition is that $v$$^{*}_{F}$ can be directly extracted from the existing experimental data. Mean-field theory calculations of upper critical fields in model superconductors are consistent with the expectation that Fermi velocities defined in this way are nearly independent of the strength of pairing interaction. Moreover, for fixed strength pairing interaction, minima in $v$$^{*}_{F}$ as a function of band filling are coincident with maxima in T c , as expected from the McMillan formula. Since no association between T c maxima and $v$$^{*}_{F}$ minima is present in MATBG experimental data, we argue that the pairing interaction in MATBG is strongly filling-factor dependent. Furthermore, any theory of MATBG superconductivity must explain this dependence, which is apparently primarily responsible for the observed superconducting domes.

2-dimensional systems↗
DOE OSTI2023

Insulators at fractional fillings in twisted bilayer graphene partially aligned to hexagonal boron nitride

At partial fillings of its flat electronic bands, magic-angle twisted bilayer graphene (MATBG) hosts a rich variety of competing correlated phases that show sample-to-sample variations. Divergent phase diagrams in MATBG are often attributed to the sublattice polarization energy scale, tuned by the degree of alignment of the hexagonal boron nitride (hBN) substrates typically used in van der Waals devices. Unaligned MATBG exhibits unconventional superconductor and correlated insulator phases, while nearly perfectly aligned MATBG/hBN exhibits zero-field Chern insulating phases and lacks superconductivity. Here we use scanning tunneling microscopy and spectroscopy (STM/STS) to observe gapped phases at partial fillings of the flat bands of MATBG in a new intermediate regime of sublattice polarization, observed when MATBG is only partially aligned (θ Gr-hBN ≈ 1.65°) to the underlying hBN substrate. Under this condition, MATBG hosts not only phenomena that naturally interpolate between the two sublattice potential limits, but also unexpected gapped phases absent in either of these limits. At charge neutrality, we observe an insulating phase with a small energy gap (Δ < 5 meV) likely related to weak sublattice symmetry breaking from the hBN substrate. In addition, we observe new gapped phases near fractional fillings ν = ±1/3 and ν = ±1/6, which have not been previously observed in MATBG. Importantly, energy-resolved STS unambiguously identifies these fractional filling states to be of single-particle origin, possibly a result of the super-superlattice formed by two moiré superlattices. Finally, our observations emphasize the power of STS in distinguishing single-particle gapped phases from many-body gapped phases in situations that could be easily confused in electrical transport measurements, and demonstrate the use of substrate engineering for modifying the electronic structure of a moiré flat-band material.

Physics↗
DOE OSTI2021

Twisted bilayer graphene. I. Matrix elements, approximations, perturbation theory, and a k · p two-band model

We investigate the twisted bilayer graphene (TBG) model of Bistritzer and MacDonald (BM) [Bistritzer and MacDonald, Proc. Natl. Acad. Sci. 108, 12233 (2011)] to obtain an analytic understanding of its energetics and wave functions needed for many-body calculations. We provide an approximation scheme for the wave functions of the BM model, which first elucidates why the BM K M -point centered original calculation containing only four plane waves provides a good analytical value for the first magic angle (θ M ≈ 1°). The approximation scheme also elucidates why most of the many-body matrix elements in the Coulomb Hamiltonian projected to the active bands can be neglected. By applying our approximation scheme at the first magic angle to a Γ M -point centered model of six plane waves, we analytically understand the reason for the small Γ M -point gap between the active and passive bands in the isotropic limit w 0 = w 1 . Furthermore, we analytically calculate the group velocities of the passive bands in the isotropic limit, and show that they are almost doubly degenerate, even away from the Γ M point, where no symmetry forces them to be. Furthermore, moving away from the Γ M and K M points, we provide an explicit analytical perturbative understanding as to why the TBG bands are flat at the first magic angle, despite the first magic angle is defined by only requiring a vanishing K M -point Dirac velocity. We derive analytically a connected “magic manifold” w 1 = $2\sqrt{1 + w^{2}_{0}}$ $-\sqrt{2 + 3w^2_0}$, on which the bands remain extremely flat as w 0 is tuned between the isotropic (w 0 = w 1 ) and chiral (w 0 = 0) limits. We analytically show why going away from the isotropic limit by making w 0 less (but not larger) than w 1 increases the Γ M -point gap between the active and the passive bands. Finally, by perturbation theory, we provide an analytic Γ M point k ∙ p two-band model that reproduces the TBG band structure and eigenstates within a certain w 0 , w 1 parameter range. Further refinement of this model are discussed, which suggest a possible faithful representation of the TBG bands by a two-band Γ M point k ∙ p model in the full w 0 , w 1 parameter range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Junctions and Superconducting Symmetry in Twisted Bilayer Graphene

Junctions provide a wealth of information on the symmetry of the order parameter of superconductors. Here, we analyze junctions between a scanning tunneling microscope (STM) tip and superconducting twisted bilayer graphene (TBG) and TBG Josephson junctions (JJs). We compare superconducting phases that are even or odd under valley exchange (s- or ƒ-wave). The critical current in mixed (s and ƒ) JJs strongly depends on the angle between the junction and the lattice. In STM-TBG junctions, due to Andreev reflection, the ƒ-wave leads to a prominent peak in subgap conductance, as seen in experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Topological Mixed Valence Model for Twisted Bilayer Graphene

Song and Bernevig (SB) have recently proposed a topological heavy-fermion description of the physics of magic angle twisted bilayer graphene (MATBG), involving the hybridization of flat-band electrons with a relativistic conduction sea. Here, we explore the consequences of this model, seeking a synthesis of understanding drawn from heavy-fermion physics and MATBG experiments. Our work identifies a key discrepancy between measured and calculated on-site Coulomb interactions, implicating renormalization effects that are not contained in the current model. With these considerations in mind, we consider a SB model with a single, renormalized on-site interaction between the f electrons, containing a phenomenological heavy-fermion binding potential on the moiré A A sites. This feature allows the simplified model to capture the periodic reset of the chemical potential with filling and the observed stability of local moment behavior. We argue that a two-stage Kondo effect will develop in MATBG as a consequence of the relativistic conduction band: Kondo I occurs at high temperatures, establishing a coherent hybridization at the Γ points and a non-Fermi liquid of incoherent fermions at the moiré K points; at much lower temperatures, Kondo II leads to a Fermi liquid in the flat band. Utilizing an auxiliary-rotor approach, we formulate a mean-field treatment of MATBG that captures this physics, describing the evolution of the normal state across a full range of filling factors. By contrasting the relative timescales of phonons and valence fluctuations in bulk heavy-fermion materials with that of MATBG, we are led to propose a valley-polaron origin to the Coulomb renormalization and the heavy-fermion binding potential identified from experiment. We also discuss the possibility that the two-fluid, non-Fermi liquid physics of the relativistic Kondo lattice is responsible for the strange-metal physics observed in MATBG. Published by the American Physical Society 2025

Lau, Liam L. H. (ORCID:0000000166039088)↗
DOE OSTI2020

Current-Induced Reversal of Anomalous Hall Conductance in Twisted Bilayer Graphene

It is observed experimentally that the sign of the Hall resistance can be flipped by a dc electric current in the twisted bilayer graphene (TBG) at 3/4 filling of the fourfold degenerate conduction flat bands. The experiment implies a switching of the valley polarization (VP) and topology in TBG. In this work, we present a theory on the current-induced switching of VP and topology. The presence of current in the bulk causes the redistribution of electron occupation in bands near the Fermi energy, which then deforms and shifts the band dispersion due to the Coulomb interaction. Above a critical current, the original occupied and empty bands can be swapped, resulting in the switching of VP and topology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Twisted bilayer graphene as topological heavy fermion: II. Analytical approximations of the model parameters

The recently-introduced topological heavy fermion (THF) model [1] of twisted bilayer graphene (TBG) aims to reconcile the quantum-dot-like electronic structure of the latter observed by scanning tunneling microscopy, with its electron delocalization seen in transport measurements. The THF model achieves this by coupling localized (heavy) fermions with anomalous conduction electrons. Originally, the parameters of the THF model were obtained numerically from the Bistritzer–Macdonald (BM) model of TBG [1]. In this work, we derive analytical expressions for the THF model parameters as a function of the twist angle, the ratio between the tunneling amplitudes at the AA and AB regions (w 0 /w 1 ), and the screening length of the interaction potential. By numerically computing the THF model parameters across an extensive experimentally-relevant parameter space, we show that the resulting approximations are remarkably good, i.e., within the 30% relative error for almost the entire parameter space. At the single-particle level, the THF model accurately captures the energy spectrum of the BM model over a large phase space of angles and tunneling amplitude ratios. When interactions are included, we also show that the THF description of TBG is good around the magic angle for realistic values of the tunneling amplitude ratios (0.6 ≤ w 0 /w 1 ≤ 1.0), for which the hybridization between the localized and conduction fermions γ is smaller than the onsite repulsion of the heavy fermions U 1 (i.e., |γ| < U 1 ).

Approximation methods↗
DOE OSTI2024

Twisted bilayer graphene’s gallery of phases

The simultaneous occurrence of exotic phases, and the ability to easily tune them, has positioned magic-angle twisted bilayer graphene as one of the richest materials platforms in condensed-matter physics.

Physics↗
DOE OSTI2023

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund’s Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

For this work, we apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)], to obtain its Kondo lattice limit. In this limit localized $\mathcal{f}$ electrons on a triangular lattice interact with topological conduction $\mathcal{c}$ electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0, ±1, ±2 are controlled by the heavy $\mathcal{f}$ electrons, while ν=±3 is at the border of a phase transition between two $\mathcal{f}$-electron fillings. For ν=0, ±1, ±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the $\mathcal{f}$ moments in the full model and analytically prove the SU(4) Hund’s rule for the ground state which maintains that two $\mathcal{f}$ electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite “bond-operators” method. We then compute the excitation spectrum of the $\mathcal{f}$ moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the $\mathcal{f}$ moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗