Efecto Hall Anómalo Cuántico en Monocapas de Miel Ti/Zr/Hf-Bi/Sb: Predicción por Primeros Principios

Resumen

Los avances en investigación de materiales bidimensionales revelan continuamente oportunidades para nuevas fases topológicas. Aunque el efecto Hall anómalo cuántico (QAHE) se ha demostrado recientemente en aislantes topológicos dopados magnéticamente, permanece confinado a temperaturas criogénicas. Usando teoría de funcional de densidad, predecimos que monocapas de miel compuestas por metales de transición Ti, Zr u Hf combinados con pnictógenos Bi o Sb (MPn, M = Ti, Zr, Hf; Pn = Bi, Sb) albergan QAHE intrínseco. Entre los sistemas estudiados, TiBi y HfBi planares emergen como aislantes QAH con brechas de banda de 15 meV y 7 meV, respectivamente, bajo tensión tracción moderada. Revelamos cómo variar constantes de red y distancias de abultamiento impulsa transiciones de fase topológica, observamos inversiones de banda inducidas por acoplamiento espín-órbita, y confirmamos modos de borde quirales en nanocintas zigzag mediante modelos de enlace fuerte derivados de Wannier. Estos hallazgos expanden el catálogo de materiales QAH realistas y sugieren plataformas viables para dispositivos espintrónicos de baja energía.

Antecedentes

Los aislantes topológicos bidimensionales (2D TIs) exhiben estados de borde sin brecha polarizados en espín mientras mantienen un bulk aislante, sello del efecto Hall cuántico de espín (QSH) [1]. Cuando la simetría de inversión temporal (TRS) se rompe —ya sea por magnetismo intrínseco o perturbaciones externas— los canales de borde helicoidales se convierten en quirales, dando origen al efecto Hall anómalo cuántico (QAHE) [3]. El transporte de carga disipación-libre y sin campo magnético del QAHE lo hace atractivo para espintrónica de ultra-bajo consumo [4], sin embargo, realizaciones experimentales se han limitado a aislantes topológicos magnéticos diluidos [8]. El trabajo teórico indica que introducir ferromagnetismo y fuerte acoplamiento espín-órbita (SOC) en un material QSH puede abrir una brecha topológica de banda y generar número de Chern no nulo [9,10]. Mientras muchos sistemas QSH —como monocapas de grupo IV (Sn) y grupo V (Bi, Sb)— se han estudiado extensamente, el potencial de mieles basadas en metales de transición para QAHE permanece largamente inexplorado. Aquí llenamos este vacío explorando mieles MPn (M = Ti, Zr, Hf; Pn = Bi, Sb), basándonos en predicciones recientes de que mieles III-V pueden albergar fases QSH y que dopaje con metales de transición puede inducir magnetismo [17,25-27]. Nuestro trabajo emplea cálculos de primeros principios para evaluar si la sustitución de Ti, Zr u Hf en mieles de Bi/Sb puede estabilizar una fase QAH intrínseca.

Resultados y Discusiones

Reemplazar la mitad de átomos de Bi en una red de miel prístina con un metal de transición (Ti, Zr u Hf) produce monocapas MPn que pueden adoptar geometrías de bajo abultamiento o planares. La Figura 1a muestra vista superior de la red M-Pn, mientras Figs. 1b y 1c ilustran vistas laterales de configuraciones abultadas y planares, respectivamente. La correspondiente primera zona de Brillouin (BZ) con puntos de alta simetría se muestra en Figura 1d.

a Estructura cristalina de miel M-Sb/Bi. b, c Vistas laterales de estructuras abultada y planar, respectivamente. d Primera zona de Brillouin (BZ) con puntos de alta simetría

Evaluamos estabilidad estructural y efectos de tensión variando constante de red y relajando posiciones atómicas para fases abultadas y planares. Curvas de energía total para TiBi, ZrBi y HfBi se grafican en Figura 2a-c. Estructuras MPn de bajo abultamiento consistentemente poseen energías menores que sus contrapartes planares. No obstante, el QAHE aparece predominantemente en régimen planar tensionado, con solo una estrecha ventana de constante de red soportando fase QAH en HfBi abultado (ver Fig. 2c).

Diagrama de fase de a TiBi, b ZrBi, y c HfBi mostrando energía total a diferentes constantes de red. El diagrama se divide en regiones etiquetadas QAH (fase Hall anómalo cuántico), I (aislante), y SM (semi-metal). Círculos azules y triángulos rojos representan casos abultados y planares, respectivamente

La Tabla 1 lista constantes de red de equilibrio, momentos magnéticos, brechas de banda y números de Chern calculados para todos los sistemas MPn. Números de Chern enteros no nulos identifican aislantes QAH: TiBi planar (C = +1, Δ = 15 meV) y HfBi planar (C = -1, Δ = 7 meV) son los candidatos más prometedores. La transición de fase en TiBi puede impulsarse reduciendo distancia de abultamiento (Fig. 3), mientras que tensión tracción induce transición en HfBi (Fig. 4).

Transición de fase tras variar distancia abultada. a Diagrama de fase de TiBi a a=4.6 Å. La flecha muestra el camino de la transición. b–f Transición de estructura de banda al reducir distancia de abultamiento (δ) de 0.44 a 0.4 Å. La transición ocurre en δ=0.41 Å

Transición de fase tras variar constante de red. a Diagrama de fase de HfBi abultado. La flecha muestra el camino de la transición. b–h Transición de estructura de banda al aumentar constante de red de 4.7 a 5.1 Å

Estructuras de banda electrónicas para geometrías de equilibrio se presentan en Figura 5. En geometría planar, TiBi exhibe clara inversión de banda y brecha de 15 meV, mientras HfBi planar muestra brecha de 7 meV con C = -1. HfBi abultado, sin embargo, es semi-metal con alto número de Chern (C = -3). Otras configuraciones, como ZrBi abultado y TiSb, permanecen aislantes triviales.

Estructuras de banda electrónicas de M-Pn (M=Ti, Zr, Hf; Pn=Sb, Bi) en sus constantes de red de equilibrio para a planar y b abultado. Círculos rojos y azules indican contribuciones +s_z y -s_z, respectivamente

Para diseccionar el papel de SOC y magnetismo, nos enfocamos en TiBi planar (a = 4.76 Å). Figura 6 muestra que cálculo no magnético produce estado metálico, mientras ordenamiento ferromagnético (momento magnético ≈ 1.05 µB por celda, dominado por Ti) produce espectro semi-metálico. Inclusión de SOC abre brecha de 15 meV, confirmando que inversión de banda es impulsada por SOC y responsable del QAHE.

Estructuras de banda electrónicas de película TiBi planar a a=4.76 Å para cálculos no magnéticos (a) sin SOC y (c) con SOC así como ferromagnéticos (b) sin SOC y (d) con SOC. Círculos rojos y azules indican contribuciones +s_z y -s_z, respectivamente, para (c) no magnético (d) cálculos ferromagnéticos con SOC

Estados de borde son huellas digitales esenciales de una fase topológica. Usando Hamiltonianos de enlace fuerte derivados de Wannier, construimos nanocinta zigzag de HfBi de 127 Å de ancho (Fig. 7). La estructura de banda revela modos de borde quirales cruzando nivel de Fermi un número impar de veces, igual a |C| = 1, confirmando la naturaleza QAH de HfBi planar.

Estructura de banda a lo largo del borde de nanocinta zigzag de HfBi abultado con a=4.9 Å y ancho 127 Å. Círculos azules (rojos) indican contribución de bordes izquierdo (derecho). Bandas de bulk denotadas por región rellena naranja

Cálculos fonónicos muestran frecuencias imaginarias pequeñas para monocapas libres, sugiriendo que un sustrato puede requerirse para realización experimental. Adicionalmente, dentro de supercelda 2×2, la configuración ferromagnética sigue siendo energéticamente favorable comparada a orden antiferromagnético en casos abultados, mientras sistemas planares exhiben cuasi-degeneración entre estados FM y AFM.

Conclusiones

Nuestro estudio de primeros principios demuestra que sustituir Ti, Zr u Hf en redes de miel de Bi o Sb puede inducir QAHE intrínseco. Mientras geometrías abultadas son energéticamente preferidas, tensión tracción moderada o reducción de distancia de abultamiento impulsa transición a configuración planar que alberga número de Chern no nulo. TiBi y HfBi planares emergen como aislantes QAH robustos con brechas de banda de 15 meV y 7 meV, respectivamente. Estos materiales se suman al creciente catálogo de candidatos QAH bidimensionales y son plataformas prometedoras para dispositivos espintrónicos de baja potencia.

Métodos/Experimental

Todos los cálculos se realizaron dentro del marco de teoría de funcional de densidad usando aproximación de gradiente generalizado (GGA) y potenciales de onda aumentada proyectada (PAW) como implementados en VASP 5.3. El corte de onda plana se fijó a 350 eV, y estructuras se relajaron hasta fuerzas residuales por debajo de 5×10⁻³ eV/Å. Convergencia del bucle auto-consistente electrónico se alcanzó a 10⁻⁶ eV. Espaciado de vacío de 20 Å a lo largo del eje z previno interacciones entre capas, y malla Monkhorst-Pack 24×24×1 centrada en Gamma muestreó la zona de Brillouin 2D. Funciones de Wannier maximamente localizadas se generaron con WANNIER90, permitiendo construcción de modelos de enlace fuerte y cálculo de estados de borde. Invariantes topológicos se obtuvieron via paquete Z2Pack, que rastrea centros de carga híbridos de Wannier para obtener número de Chern.