Simulación ab initio de la respuesta a radiación de baja energía en Si, Ge y superredes Si/Ge

Resumen

Mediante dinámica molecular ab initio, examinamos cómo la irradiación iónica de baja energía afecta al silicio (Si) masivo, germanio (Ge) masivo y una superred (SL) Si/Ge. El estudio revela que los átomos de Ge en la interfaz de la SL mantienen una resistencia al desplazamiento similar al Ge masivo, mientras que los átomos de Si cerca de la interfaz exhiben un umbral de desplazamiento mayor que el Si masivo, indicando una tolerancia a la radiación mejorada. Los mecanismos de generación de defectos difieren entre las estructuras masivas y la SL, produciendo esta última cúmulos de defectos más complejos. Los cálculos de energías de formación y barreras de migración de defectos demuestran que los defectos puntuales en la SL son más difíciles de formar y las vacancias son menos móviles. Estos hallazgos sugieren que las superredes Si/Ge son prometedoras para dispositivos electrónicos y optoelectrónicos que operan en entornos ricos en radiación.

Antecedentes

En años recientes, la superred Si/Ge ha emergido como plataforma versátil para dispositivos electrónicos y optoelectrónicos de próxima generación, impulsada por su estructura de bandas ajustable y movilidad de portadores mejorada. Las aplicaciones abarcan fotodiodos rápidos para comunicaciones ópticas, microelectrónica, células solares y componentes espaciales. Sin embargo, la exposición a iones de alta energía en el espacio puede alterar las propiedades ópticas y electrónicas de estos materiales, llevando a degradación del dispositivo. Por consiguiente, comprender la respuesta a la radiación de las SL Si/Ge bajo condiciones extremas es esencial.

Trabajos experimentales han documentado la superior resistencia a la radiación de las SL Si/Ge comparadas con Si masivo. Por ejemplo, estudios de irradiación electrónica en SL con monocapa de Ge embebida mostraron degradación reducida de fotoluminiscencia relativa al silicio masivo, mientras que la irradiación con protones de 2 MeV en estructuras de puntos cuánticos de Ge en SL reveló tolerancia excepcional. Tendencias similares se observaron en dispositivos de pozo cuántico de Ge sobre multicapas Si/Ge y en mediciones termoeléctricas, donde la fluencia iónica aumentó la figura de mérito al mejorar la dispersión de fonones y transporte electrónico. A pesar de estos conocimientos, falta una comprensión atomística completa de la formación y evolución de defectos en SL Si/Ge.

La dinámica molecular ab initio (AIMD) ofrece una herramienta poderosa para sondear daño por radiación a escala atómica, proporcionando potenciales interatómicos derivados de cálculos de estructura electrónica de primeros principios. A diferencia de MD clásico, AIMD captura efectos de transferencia de carga y apantallamiento electrónico, permitiendo determinación precisa de energías de desplazamiento umbral (E_d) y configuraciones de defectos. En este trabajo, empleamos AIMD para comparar eventos de desplazamiento de baja energía en Si masivo, Ge masivo y una SL Si₂/Ge₂, identificar los mecanismos subyacentes a sus distintos comportamientos ante radiación, y cuantificar la energética de formación y migración de defectos.

Métodos

Las cascadas de desplazamiento se simularon con el código SIESTA. Pseudopotenciales de Troullier-Martins normoconservadores describieron la interacción ión-electrón, mientras que la aproximación de densidad local (LDA) en la parametrización de Ceperley-Alder trató efectos de intercambio-correlación. Los estados de valencia se expandieron en una base de zeta simple más polarización (SZP), con malla de puntos k 1 × 1 × 1 y corte de 60 Ry. La supercelda modelo para la SL comprendió 288 átomos dispuestos en dos capas de Si alternando con dos de Ge (Si₂/Ge₂). La Figura 1 ilustra las geometrías de Si masivo y la SL Si/Ge.

Se seleccionó un átomo primario de retroceso (PKA) y se le dio energía cinética para iniciar un evento de retroceso. Si el PKA retornaba a su sitio de red, la simulación se repetía a mayor energía con incrementos de 5 eV; una vez que el PKA permanecía desplazado, el umbral se refinaba a precisión de 0,5 eV. Para cada elemento, examinamos cuatro a cinco direcciones cristalográficas en el masivo y la SL. Todas las corridas usaron el ensamble NVE y se limitaron a 1,2 ps para evitar inestabilidad numérica.

Resultados y Discusión

Eventos de desplazamiento en silicio y germanio masivos

Nuestros cálculos AIMD arrojaron constantes de red de 5,50 Å para Si masivo (coincidiendo con el experimental 5,43 Å y teórico 5,48 Å) y 5,71 Å para Ge masivo (en línea con experimental 5,77 Å y teórico 5,65 Å). La Tabla 1 resume las energías de desplazamiento umbral (E_d) y las estructuras de defecto resultantes, mientras que las Figuras 2 y 3 muestran las configuraciones de daño para Si y Ge, respectivamente.

Para Si masivo, los valores E_d calculados están cerca de mediciones experimentales: 21 eV para [001], ~47,6 eV para [110], y ~12,9 eV para [111]. El tipo de defecto primario en todos los casos es un par Frenkel (vacancia más intersticial). Nuestro E_d para [110] (47 eV) supera el valor de MD clásico (24 eV), reflejando la inclusión de efectos de transferencia de carga en AIMD. Los mecanismos de desplazamiento varían: a lo largo de [111], el PKA se dispersa a un sitio intersticial mientras el átomo desplazado retorna a su sitio de red; a lo largo de [001], una secuencia más compleja de intercambios produce dos vacancias y dos intersticiales; y a lo largo de [110], múltiples colisiones llevan a un cúmulo de pares Frenkel. Estos resultados confirman que la dirección [110] es la más resistente al desplazamiento en Si masivo.

El Ge masivo muestra valores E_d que se alinean bien con experimento y teoría: ~18 eV para [001] y 9,5 eV para [111] y [1 1 1]. El desplazamiento a lo largo de [111] produce una simple formación intersticial, mientras que a lo largo de [110] el PKA intercambia posiciones con un átomo de Ge vecino, creando una vacancia y un intersticial. Como en Si, la dirección [110] es la más difícil de desplazar. El panorama general de defectos en Si y Ge masivos está dominado por pares Frenkel.

Eventos de desplazamiento en superred Si/Ge

Investigamos la SL Si₂/Ge₂ seleccionando átomos PKA adyacentes a la interfaz Si/Ge. La Tabla 2 lista los valores E_d y defectos asociados, mientras que las Figuras 4 y 5 ilustran las configuraciones resultantes.

En la SL, los átomos de Si son altamente resistentes a lo largo de la dirección [1 1 1] (E_d = 10 eV), mientras que Si[001] y Si[00 1̅] requieren 46,5 eV y 42,5 eV respectivamente, y generan cúmulos complejos de antisitio. Los átomos de Ge en la SL retienen bajos valores E_d a lo largo de [111] y [1 1 1̅] (≈16 eV) similares al Ge masivo, pero muestran umbrales mayores a lo largo de [001] (≈17 eV) y [00 1̅] (≈18 eV). Importante, Ge[110] en la SL tiene un E_d reducido (8,5 eV) y produce un cúmulo de defecto mixto (vacancia, antisitio, intersticial). Estos hallazgos confirman que la SL mejora la resistencia a la radiación para Si y, en menor medida, para Ge, de acuerdo con observaciones experimentales previas.

Energía de formación de defectos y barrera de migración en Si masivo, Ge masivo y superred Si/Ge

Calculamos energías de formación de defectos para vacancias, intersticiales y antisitios en Si masivo, Ge masivo y la SL usando una supercelda de 64 átomos, muestreo de puntos k 6 × 6 × 6 y corte de 500 eV. Los resultados, mostrados en la Tabla 3, indican que V_Si (3,60 eV) y V_Ge (2,23 eV) son los defectos más favorables en sus respectivos materiales masivos, mientras que la SL muestra energías de formación ligeramente mayores para la mayoría de defectos, con V_Ge (2,73 eV) y V_Si (2,85 eV) siendo aún los más bajos.

Las barreras de migración se evaluaron para V_Ge y V_Si adyacentes a la interfaz SL usando el método de banda elástica nudged (Tabla 4). Las vacancias migran más fácilmente a lo largo de la dirección [111] tanto en masivo como en SL, pero las barreras son consistentemente menores en el masivo (p. ej., V_Si a lo largo de [111]: 0,11 eV en masivo vs 0,12 eV en SL). En contraste, la migración a lo largo de [100] y [110] muestra diferencias modestas. Estos resultados sugieren que la movilidad de vacancias se suprime en la SL, lo que podría reducir la agrupación de defectos e hinchamiento bajo irradiación.

Conclusiones

Nuestra investigación AIMD demuestra que las energías de desplazamiento umbral en Si y Ge masivos son fuertemente anisotrópicas, siendo la dirección [110] la más resistente. En la SL Si/Ge, los átomos de Si exhiben resistencia mejorada a lo largo de direcciones [111] y [001], mientras que los átomos de Ge mantienen umbrales tipo masivo excepto por un valor reducido a lo largo de [110]. Las energías de formación de defectos confirman que los defectos puntuales son más difíciles de crear en la SL, y las barreras de migración indican movilidad de vacancias reducida, especialmente a lo largo de la dirección [111] más favorable. Estos efectos combinados sustentan la superior tolerancia a la radiación de las superredes Si/Ge, haciéndolas atractivas para dispositivos que deben soportar entornos de radiación severos.