Investigación desde primeros principios de la estabilidad estructural e imagen STM de borofeno en Ag(111)

Resumen

Trabajos experimentales recientes han demostrado el crecimiento exitoso de láminas atómicamente delgadas de boro —borofeno— en Ag(111) por epitaxia de haz molecular. Se observaron dos fases de superficie distintas, pero los arreglos atómicos precisos y su estabilidad termodinámica siguen debatidos. En este estudio empleamos teoría de funcional de densidad para examinar los polimorfos de borofeno más plausibles —triangular abovedado, β12 y χ3— cuando están soportados en Ag(111). Nuestros cálculos revelan que las tres láminas independientes son termodinámicamente inestables y metálicas, mientras que el sustrato Ag(111) las estabiliza. Las imágenes de microscopía de túnel de barrido (STM) simuladas reproducen las fases rayada y homogénea observadas experimentalmente e identifican inequívocamente las estructuras de red subyacentes.

Antecedentes

Desde el descubrimiento del grafeno, los materiales bidimensionales (2D) han atraído intenso interés investigador debido a sus propiedades electrónicas, mecánicas y catalíticas únicas. 1–7 Durante la última década, se han reportado varias propuestas teóricas para láminas 2D de boro —colectivamente denominadas borofeno— con diversos arreglos de vacancias hexagonales y patrones de abovedamiento. 8–21 Aunque extensos estudios desde primeros principios han sugerido que redes triangulares planares de boro con huecos hexagonales periódicos son energéticamente favorecidas, la forma 2D sigue siendo más alta en energía que el boro masivo, implicando un impulso termodinámico hacia la cristalización 3D. 22–24 Consecuentemente, se requiere un sustrato suficientemente "pegajoso" para suprimir la nucleación 3D y favorecer la formación de capas atómicamente delgadas. 25–27

Avances experimentales recientes han realizado borofeno en superficies metálicas. Mannix y colaboradores depositaron boro exitosamente en Ag(111) y observaron dos fases de superficie distintas —una fase rayada y una fase homogénea— usando STM de alta resolución. 28 Feng et al. reportaron fases similares y asignaron la estructura homogénea a una red χ3, mientras que la fase rayada se interpretó como una red triangular abovedada o rectangular β12. 29–31 Estas interpretaciones divergentes subrayan la necesidad de una investigación teórica sistemática de las configuraciones atómicas posibles y su estabilidad en Ag(111).

Más allá de la identificación estructural, los polimorfos de borofeno exhiben propiedades notables. La lámina triangular abovedada se predice como un metal altamente anisotrópico con módulo de Young superior al grafeno en dirección zigzag, y muestra conductividad térmica de red fuertemente anisotrópica. 32,33 La superconductividad también se ha predicho para borofeno β12 y χ3, potencialmente rivalizando con películas delgadas de MgB2. 34–36 Sin embargo, la estabilidad termodinámica de β12 y χ3 sigue siendo controvertida, con algunos estudios reportando frecuencias fonónicas imaginarias cerca del punto Γ. 37,38

Métodos computacionales

Todos los cálculos se realizaron con el paquete Vienna Ab-Initio Simulation Package (VASP) usando teoría de funcional de densidad (DFT). 38 El método de onda aumentada proyectada describió interacciones electrón-ión, 39 y la aproximación de gradiente generalizado (GGA) con el funcional Perdew-Burke-Ernzerhof (PBE) se empleó para intercambio-correlación. 40 Los conjuntos de base de onda plana se truncaron a 500 eV, y la zona de Brillouin se muestreó con mallas k-punto 25 × 15 × 1, 15 × 9 × 1 y 11 × 11 × 1 para las fases triangular abovedado, β12 y χ3, respectivamente. Un espaciado de vacío de al menos 20 Å en dirección z minimizó interacciones entre imágenes periódicas. Las optimizaciones geométricas convergieron a fuerzas bajo 0,02 eV Å⁻¹, mientras que las dos capas inferiores de Ag se fijaron. Las dispersiones fonónicas se obtuvieron usando el método de desplazamiento finito implementado en PHONOPY. 41

Las imágenes STM se simularon dentro del marco Tersoff-Hamann, que relaciona la corriente de túnel con la densidad local de estados (LDOS) de la muestra. 42 Se usó el modo de corriente constante para generar imágenes a varios voltajes de sesgo, proporcionando comparación directa con datos STM experimentales.

Resultados y discusión

La Figura 1 muestra las geometrías optimizadas de los tres polimorfos de borofeno considerados: triangular abovedado (a), β12 (b) y χ3 (c). La lámina triangular abovedada exhibe un único abovedamiento fuera del plano a lo largo del eje b, mientras que β12 y χ3 son planares con arreglos distintos de vacancias hexagonales. Las constantes de red —1,613 Å × 2,866 Å para triangular abovedado, 2,916 Å × 5,075 Å para β12 y 4,448 Å para χ3— concuerdan bien con reportes teóricos y experimentales previos. 23,25,27

La concentración de vacancias η es 1/6 para β12 y 1/5 para χ3, indicando un mayor grado de ordenamiento de defectos en esta última. La estabilidad relativa se evaluó calculando la energía de enlace promedio por átomo de boro, E_FB, definida como E_boron/n. Los resultados (Tabla 2) muestran que la lámina β12 independiente es la más estable, mientras que χ3 es la menos estable, diferenciándose por 0,08 eV por átomo. Sin embargo, los cálculos de dispersión fonónica (Figura 2) revelan frecuencias imaginarias cerca de Γ para las tres estructuras, confirmando que son termodinámicamente inestables en aislamiento.

Las estructuras de banda electrónicas (Figura 3) confirman que los tres polimorfos son metálicos. La lámina triangular abovedada exhibe fuerte anisotropía, con cruces de bandas a lo largo de S–Y y G–X, y huecos considerables a lo largo de X–S (9,63 eV) e Y–G (4,32 eV). Estas características sugieren conductividad direccional confinada al eje a no abovedado.

Al colocarse en Ag(111), el desajuste de red es mínimo (~1 %) para β12 y moderado (~3 %) para triangular abovedado. La lámina χ3 adopta dos registros distintos (χ3 y χ3′) en el sustrato, con constantes de red 8,67 Å y 2,89 Å × 25,02 Å, respectivamente. La separación vertical entre la capa de borofeno y la superficie de Ag varía de 2,4 a 3,3 Å, indicando fisisorción débil que preserva el carácter planar de β12, χ3 y χ3′. Estas distancias intercapas son consistentes con los espesores experimentales reportados para borofeno en Ag(111). 23,28

Las imágenes STM simuladas (Figura 5) reproducen las fases rayada (Mannix et al.) y homogénea (Feng et al.) observadas experimentalmente. El patrón rayado para triangular abovedado aparece como crestas brillantes en forma de huso, mientras que β12 produce protuberancias ovaladas rodeadas por esquinas brillantes. La lámina χ3 muestra un patrón de mancuerna romboédrica coincidente con la fase S2 reportada por Feng et al. Simulaciones dependientes del sesgo (Figura 6) distinguen aún más las dos configuraciones rayadas: la raya triangular abovedada invierte contraste entre sesgo positivo y negativo, mientras que la raya β12 retiene su forma ovalada independientemente del sesgo, apoyando la asignación de la raya de Mannix a triangular abovedado y la de Feng a β12.

En resumen, cálculos desde primeros principios demuestran que borofeno triangular abovedado, β12 y χ3 independientes son termodinámicamente inestables y metálicos. El sustrato Ag(111) estabiliza estas láminas, reduciendo sus energías en 0,1–0,2 eV por átomo de boro. Las imágenes STM simuladas proporcionan un mapeo definitivo entre fases de superficie observadas y estructuras de red subyacentes, resolviendo el debate prolongado sobre la configuración atómica de borofeno en plata. Estos insights sientan las bases para el diseño racional de dispositivos basados en borofeno y para guiar síntesis experimental futura en sustratos metálicos.

Abreviaturas

2D

Bidimensional

3D

Tridimensional

STM

Microscopía de túnel de barrido