Bond reconstruction and vacancy clustering in monolayer silicon carbide from first principles

Resumen Técnico: Reconstrucción de Enlaces y Agrupamiento de Vacantes en Monocapa de Carburo de Silicio desde Primeros Principios

Planteamiento del Problema
Los materiales bidimensionales (2D), particularmente el carburo de silicio (SiC) monocapa, están emergiendo como plataformas favorables para las tecnologías cuánticas basadas en defectos debido a sus propiedades mecánicas, eléctricas y magnéticas. Aunque los defectos intrínsecos como las vacantes son candidatos potenciales para qubits y emisores cuánticos, su estabilidad y propiedades cuánticas están fuertemente influenciadas por las relajaciones estructurales. Estudios computacionales previos sobre vacantes en SiC monocapa a menudo reportaron estructuras de alta simetría con enlaces colgantes o limitaron su alcance a relajaciones en el plano. Sin embargo, la medida en que las distorsiones de ruptura de simetría, la reconstrucción de enlaces y los desplazamientos fuera del plano afectan las energías de formación, las estructuras electrónicas y la actividad óptica de las vacantes individuales y sus agregados permanece inexplorada. Específicamente, el papel de la reconstrucción de enlaces para estabilizar cúmulos de vacantes compactos y definir sus propiedades cuánticas requiere una investigación exhaustiva desde los primeros principios.

Metodología
Los autores emplearon cálculos de la Teoría de la Funcional de la Densidad (DFT) utilizando el código VASP para investigar monovacantes y agregados de vacantes en SiC monocapa.

• Configuración Computacional: Se utilizó una supercelda periódica de $10 \times 10$ con un acolchado de vacío de 20 Å, lo que permitió un muestreo de un solo punto-k en el punto $\Gamma$. Los electrones del núcleo fueron tratados mediante el formalismo de la Onda Plana de la Función de Onda Aumentada (PAW) con un corte de onda plana de 500 eV.
• Funcionales: Se utilizó el funcional de Perdew-Burke-Ernzerhof (PBE) para las relajaciones de la geometría, los cálculos de fonones y las búsquedas de rutas de reacción de la Banda Más Baja (NEB). El funcional híbrido de Heyd-Scuseria-Ernzerhof (HSE06) se aplicó para calcular las energías de formación, las estructuras de niveles electrónicos y los estados excitados.
• Herramientas de Análisis: Las propiedades de enlace se analizaron mediante el código LOBSTER a través de la Población de Orbitales de Cristal Proyectados Integrados (IpCOHP) y la topología de la Función de Localización Electrónica (ELF). Los momentos dipolares de transición se calcularon con el código VaspBandUnfolding. Las energías de formación y de unión se computaron con correcciones para sistemas cargados (método de Freysoldt) y doble conteo de átomos de la red.
• Estrategia de Búsqueda: Se realizó una búsqueda computacional sistemática para cúmulos de vacantes compactos mediante la eliminación o sustitución de átomos alrededor de semillas de monovacantes, centrándose en estados de carga neutros para modelar condiciones libres de impurezas.

Contribuciones Clave y Resultados

1. Reconstrucción de Enlaces en Monovacantes
El estudio revela que las monovacantes aisladas experimentan distorsiones significativas de ruptura de simetría que involucran tanto la reconstrucción en el plano como el desplazamiento fuera del plano.

• Vacante de Carbono ($V_C$): El estado triplete de alta simetría $D_{3h}$ es inestable. El sistema se relaja a un estado fundamental de singlete de espín con simetría $C_2$ mediante una doble reconstrucción de los enlaces colgantes de silicio, formando un patrón local 5-6-6-5. Esto viene acompañado de una gran distorsión fuera del plano (0.58 Å) y una quiralidad en la disposición de los átomos de Si circundantes. Aunque esta reconstrucción estabiliza la estructura, suprime la actividad óptica; el estado excitado es disociativo, lo que conduce a la relajación no radiativa y a una baja eficiencia cuántica.
• Vacante de Silicio ($V_{Si}$): De manera similar, el defecto $V_{Si}$ se relaja de una simetría $D_{3h}$ a una $C_{1h}$ con una distorsión fuera del plano de 1.05 Å, resultando en un estado fundamental de singlete.
• Energética: Las distorsiones fuera del plano contribuyen significativamente a la estabilización, con energías de relajación total (respecto a $D_{3h}$) que alcanzan 1.03 eV para $V_C$ y 0.31 eV para $V_{Si}$.

2. Formación de Cúmulos de Vacantes Compactos
Los autores investigaron la agregación de vacantes, identificando que la reconstrucción de enlaces actúa como un mecanismo de estabilización clave para los cúmulos compactos.

• Vías de Agregación: Entre varias vías de agregación ($V_{Si} + nV_{Si}$, $V_{Si} + nV_C$, $V_C + nV_C$), la vía $V_{Si} + nV_C$ resultó ser la termodinámicamente más favorable.
• Agregado Estable ($V_{Si}3V_C$): Un agregado específico compuesto por tres vacantes de carbono rodeando una vacante de silio central ($V_{Si}3V_C$) surgió como una estructura altamente estable. Este cúmulo exhibe simetría $D_{3h}$ (coincidiendo con el material prístino), una reconstrucción completa de los enlaces colgantes y una energía de unión excepcional. La reconstrucción del enlace Si-Si en este complejo es casi completa, con fuerzas de enlace comparables al silicio en bulto.

3. Propiedades Cuánticas del Defecto $V_{Si}3V_C$
El defecto $V_{Si}3V_C$ se identifica como un candidato prometedor para un emisor cuántico:

• Estado de Espín: Posee un estado fundamental de triplete ($S=1$), lo que lo convierte en un qubit de espín viable. El parámetro de desdoblamiento de campo cero (ZFS) $D$ se calcula en 214 MHz, consistente con la extensión espacial de la densidad de espín.
• Propiedades Ópticas: El defecto presenta un orbital doblemente degenerado de media ocupación ($E''$) y un orbital $A''_1$ no ocupado. La transición entre ellos está permitida por dipolo eléctrico con un gran momento dipolar de transición (12.5 D). La energía de la Línea de Cero Fonones (ZPL) calculada es de 0.521 eV (infrarrojo cercano).
• Acoplamiento de Fonones: El defecto exhibe dos modos de vibración localizados (LVM) ópticamente activos a 88.66 meV y 101.32 meV. El factor de Huang-Rhys total es 0.59, lo que corresponde a un factor de Debye-Waller de 0.55. Esto indica un acoplamiento débil con los fonones, lo que sugiere transiciones no radiativas suprimidas y una emisión brillante.
• Interacción Hiperfina: Se predicen fuertes interacciones hiperfinas para los átomos de Si del primer vecino más cercano, proporcionando huellas experimentales distintivas para su identificación.

Significancia
Este artículo establece que la reconstrucción de enlaces es un factor decisivo para definir la estabilidad y las propiedades cuánticas de los defectos de vacante en SiC 2D. Al demostrar que las distorsiones fuera del plano y la ruptura de simetría son críticas para la estabilidad de las vacantes individuales, y que agregados específicos de reconstrucción de enlaces pueden producir centros ópticamente activos y estables, este trabajo proporciona una base teórica para la realización de emisores cuánticos en SiC monocapa. La identificación del defecto $V_{Si}3V_C$ como un centro de color robusto en el infrarrojo con un estado fundamental de triplete y características ópticas favorables ofrece un objetivo específico para la futura realización experimental en tecnologías cuánticas.