Disorder-induced symmetry breaking in moiré bands of marginally twisted bilayer MoS$_2$

Planteamiento del Problema
Los dicalcogenuros de metales de transición retorcidos (tTMDs) albergan potenciales de moiré altamente ajustables y bandas planas que sustentan fases electrónicas fuertmente correlacionadas. Mientras que los modelos teóricos predicen ricos fenómenos cuánticos emergentes en estos sistemas, el papel específico del desorden intrínseco en la configuración de estas propiedades permanece en gran medida sin resolver. A diferencia del grafeno bicapa retorcido, que requiere un "ángulo mágico" específico para las bandas planas, los tTMDs de ángulo pequeño exhiben anchuras de banda que disminuyen continuamente cerca del apilamiento paralelo (0°) y antiparalelo (60°), ofreciendo una amplia plataforma para investigar fases cuánticas. Sin embargo, el desorden intrínseco —como los defectos puntuales y los charcos de carga— es una característica inevitable de los materiales bidimensionales. La naturaleza microscópica del desorden en los tTMDs y su impacto en la estructura electrónica de las superredes de moiré, particularmente en lo que respecta a la ruptura de simetría, ha recibido una atención limitada.

Metodología
Los autores investigaron el disulfuro de molibdeno bicapa ligeramente retorcido (tb-MoS$_2$) con un ángulo de giro $\theta \approx 0.95^\circ$ utilizando microscopía de efecto túnel de barrido (STM) y espectroscopía (STS).

• Experimental: Las muestras se fabricaron utilizando una técnica de "desgarro y apilamiento" (tear-and-stack) para crear bicapas sobre un sustrato de grafito/SiO$_2$. Se utilizó la topografía de STM para caracterizar la superred de moiré, mientras que se empleó STS para mapear la densidad local de estados (LDOS). Los investigadores midieron espectros de conductancia diferencial ($dI/dV$) en regiones de apilamiento de alta simetría (MM, MX, XM y paredes de dominio) y generaron mapas espaciales de "inicio" (onset maps) para visualizar la energía a la cual la LDOS se eleva por encima de un umbral de ruido tanto para las bandas de valencia como de conducción.
• Teórica: El estudio utilizó un modelo de continuo para el MoS$_2$ retorcido que incorpora efectos de relajación atómica derivados de cálculos de la teoría del funcional de la densidad (DFT) de primera escala. Para modelar las observaciones experimentales, los autores desarrollaron un modelo de desorden electrostático tratando el sistema como capas dieléctricas alternas sobre un sustrato metálico de grafito. Este modelo contabilizó defectos cargados distribuidos aleatoriamente (específicamente vacantes de azufre) para calcular el potencial electrostático resultante y su efecto en las bandas de moiré.

Resultados Clave

1. Desplazamientos de inicio por ruptura de simetría: Las mediciones de STS revelaron un desdoblamiento inesperado en las energías de inicio de las bandas de valencia y de conducción entre las regiones de apilamiento simétricas de espejo MX y XM. Específicamente, se observó una diferencia de $15 \pm 4$ meV para la banda de valencia y de $15 \pm 6$ meV para la banda de conducción. En un sistema prístino y simétrico, estas regiones deberían exhibir energías de inicio idénticas.
2. Desorden espacialmente correlacionado: Los mapas de inicio mostraron variaciones de energía espacialmente correlacionadas en escalas de longitud superiores al periodo de moiré. El análisis de autocorrelación de estos mapas mostró una dependencia radial consistente con un potencial de $1/r^3$, característico de los dipolos formados por defectos cargados y sus cargas imagen en el sustrato.
3. Estimación de la densidad de defectos: Al ajustar los datos de autocorrelación al modelo de desorden electrostático, los autores infirieron una densidad de carga aleatoria de fondo de aproximadamente $1 \times 10^{12}$ cm$^{-2}$ a $3 \times 10^{11}$ cm$^{-2}$. Este rango concuerda con las densidades de vacantes de azufre reportadas en monocapas de MoS$_2$ exfoliadas mecánicamente.
4. Acuerdo Teórico: Los cálculos del modelo de continuo incorporando un campo de desplazamiento eléctrico entre capas (simulando el potencial proveniente del desorden inferido) reprodujeron los resultados experimentales. Los cálculos demostraron que el campo eléctrico induce una polarización de capa en las bandas de moiré, causando un desdoblamiento de las energías de inicio entre los sitios MX y XM. La densidad local de estados calculada, al ser ensanchada para considerar el desorden y los efectos térmicos, coincidió cuantitativamente con los mapas de inicio y los desplazamientos espectrales experimentales.

Significancia y Reivindicaciones
El artículo sostiene que incluso niveles modestos de desorden electrostático, inherentes a los TMDs exfoliados, desempeñan un papel crítico en la determinación de la estructura electrónica de los materiales de moiré a escala nanométrica. La principal significancia de estos hallazgos es la demostración de que el desorden puede inducir la ruptura de simetría en las bandas de moiré, transformando el carácter de los estados electrónicos. Específicamente, las bandas de valencia superiores, que teóricamente se predice que sean no polarizadas y de tipo grafeno en un sistema prístino, demuestran estar completamente polarizadas en capa y ampliamente separadas en energía debido a la presencia de defectos cargados.

Los autores concluyen que estos efectos inducidos por el desorden ocurren en una escala de energía de decenas de meV, la cual es comparable a las anchuras de banda de las propias bandas planas. Por consiguiente, el artículo afirma que se debe tener cuidado al interpretar los resultados experimentales en los tTMDs, ya que el desorden nativo puede alterar fundamentalmente la realización de topologías de banda plana específicas y fases correlacionadas. El trabajo destaca que la relajación estructural y el desorden electrostático deben contabilizarse simultáneamente para lograr un acuerdo cuantitativo entre la teoría y el experimento en estos sistemas.