[논문 리뷰] Kekulé spiral order in magic-angle graphene: a density matrix renormalization group study

Tianle Wang, Daniel E. Parker외 5명|arXiv (Cornell University)|2022. 11. 04.

Graphene research and applications참고 문헌 75인용 수 6

한 줄 요약

이 연구는 빈도 ν=−3에서 마법의 각도로 얽힌 이중층 그래핀(TBG)에서 이방성 잔류변형의 영향을 조사하기 위해 고 bond 차원 χ=24576를 사용한 대규모 비편향 밀도 행렬 분리군(DMRG) 계산을 수행한다. 0.05% 이방성 잔류변형에서 양자화된 비정상 홀 절연체에서 이방성 키클레 스피RAL(IKS) 상으로의 전이가 발생하며, 이는 밸리 및 전이 대칭성을 깨뜨리지만 수정된 전이 대칭성 T′을 유지한다. 이후 더 높은 변형에서 정상 금속으로의 전이가 발생하여, 정렬되지 않은 샘플에서 이론과 실험 간의 모순을 해결한다.

ABSTRACT

When the two layers of a twisted moiré system are subject to different degrees of strain, the effect is amplified by the inverse twist angle, e.g., by a factor of 50 in magic angle twisted bilayer graphene (TBG). Samples of TBG typically have heterostrains of 0.1-0.7%, increasing the bandwidth of the "flat'' bands by as much as tenfold, placing TBG in an intermediate coupling regime. Here we study the phase diagram of TBG in the presence of heterostrain with unbiased, large-scale density matrix renormalization group calculations (bond dimension $χ=24576$), including all spin and valley degrees of freedom. Working at filling $ν= -3$, we find a strain of $0.05\%$ drives a transition from a quantized anomalous Hall insulator into an incommensurate-Kekulé spiral (IKS) phase. This peculiar order, proposed and studied at mean-field level by Kwan et al (PRX 11, 041063), breaks both valley conservation and translation symmetry $\hat{T}$, but preserves a modified translation symmetry $\hat{T}'$ with moiré-incommensurate phase modulation. Even higher strains drive the system to a fully symmetric metal.

연구 동기 및 목표

• 강상호작용 이론이 ν=−3에서 양자화된 비정상 홀(QAH) 절연체를 예측하는 데 반해, 정렬되지 않은 TBG 샘플에서 실험적으로 관측된 C=0 절연체 또는 금속과의 괴리를 해결하기 위함.
• 작은 그래핀 체적 기준의 현실적인 이방성 잔류변형이 역-twist 각도에 의해 증폭됨에 따라 TBG의 저에너지 물리학이 어떻게 변화하는지 조사하기 위함.
• 모든 스핀 및 밸리 자유도를 포함한 비편향 대규모 DMRG 시뮬레이션을 통해 ν=−3에서 TBG의 진동 상태의 성격을 규명하기 위함.
• 강상호작용 근사법으로는 다루기 어려운 중간 상호작용 영역에서 TBG의 위상도를 명확히 하기 위함.

제안 방법

• ν=−3에서 휘어진 이중층 그래핀의 허버드 해밀토니안을 시뮬레이션하기 위해 고 bond 차원 χ=24576를 사용한 대규모 밀도 행렬 분리군(DMRG) 계산을 수행하였다.
• 다체 상관관계를 비편향적으로 처리하기 위해 DMRG 시뮬레이션에 모든 스핀 및 밸리 자유도를 통합하였다.
• 두 그래핀 층에 이방성 잔류변형 εGr=0.1–0.7%를 적용하였으며, 이는 역-twist 각도에 의해 증폭되어 모리 수반격자와 평탄한 밴드의 대역폭을 크게 변화시켰다.
• DMRG 기저 상태 에너지 및 순서 매개변수(예: 밸리 및 스핀 밀도파, 얽힘 엔트로피)를 사용하여 다양한 변형 값에서의 상을 식별하였다.
• 상의 분류 및 대칭성 깨짐 패tern(예: 이방성 키클레 스피RAL(IKS) 순서)을 규명하기 위해 양자수(τz, sz)를 분석하였다.
• 전역 기저 상태 및 상 경계를 결정하기 위해 다양한 양자수 영역 간의 기저 상태 에너지를 비교하였다.

실험 결과

연구 질문

• RQ1 그래프레인 체적 기준으로 작은 크기임에도 불구하고, 이방성 잔류변형은 마법의 각도로 얽힌 이중층 그래핀의 다체 기저 상태에 어떻게 영향을 미치는가?
• RQ2 정렬되지 않은 TBG 샘플에서 ν=−3에서 관측된 절연체 상태의 성격은 무엇이며, 현실적인 변형 효과로 설명될 수 있는가?
• RQ3 기존에 평균장 수준에서 제안된 이방성 키클레 스피RAL(IKS) 순서가 비편향 DMRG 시뮬레이션에서 안정된 상으로 나타나는가?
• RQ4 변형이 증가함에 따라 TBG의 위상도는 강상호작용 영역에서 중간 상호작용 영역으로 어떻게 변화하는가?
• RQ5 IKS 상에서 대칭성 깨짐의 역할은 무엇이며, 기존의 양자 홀 페로자성체와 어떻게 다를까?

주요 결과

• 이방성 잔류변형 0.05%에서 양자화된 비정상 홀(QAH) 절연체에서 이방성 키클레 스피RAL(IKS) 상으로의 전이가 발생한다.
• IKS 상은 밸리 보존 대칭성과 전이 대칭성(T)을 깨뜨리지만, 모리-비공진 단계 조절을 갖는 수정된 전이 대칭성 T′을 유지한다.
• 더 높은 변형(εGr ≥ 0.2%)에서는 시스템이 완전히 대칭적인 정상 금속으로 전이되며, (0,0) 양자수 영역이 전역 기저 상태가 된다.
• IKS 상은 스핀 비대칭의 간섭상관 및 일시적인 스핀 밀도파(SDW) 순서를 보이지만, 순수한 스핀 극화는 관측되지 않는다.
• 상 간의 에너지 차이는 작지만 수치 정밀도를 초월하여 안정적인 세 상의 구조(저변형에서 QAH, 중간 변형에서 IKS, 고변형에서 정상 금속)를 나타낸다.
• 상 경계는 시스템 크기 및 유전율에 민감하지만, 전체적인 세 상도 그림은 강건할 것으로 예상된다.