Evolución de la morfología de sustratos Si(001) con patrones de pozos impulsada por difusión superficial: Un estudio de campo de fase

Resumen

El ordenamiento lateral de islas heteroepitaxiales se logra rutinariamente mediante patrones de pozos en el sustrato antes de la deposición. Controlar la forma, orientación y profundidad de estos pozos no es trivial porque las geometrías metaestables evolucionan durante el recocido o el crecimiento posterior. En este trabajo presentamos un modelo de campo de fase continuo que predice las morfologías típicas de pozos en Si(001) en función de la profundidad y forma inicial. El modelo incorpora difusión superficial y una densidad de energía superficial fuertemente anisotrópica que reproduce la forma de Wulff del silicio. Reproduce con éxito configuraciones metaestables observadas experimentalmente y revela una vía cinética común hacia el aplanamiento global del pozo, independientemente de la profundidad inicial. Al extender el modelo a la heteroepitaxia de dos materiales con facetas de equilibrio distintas, mostramos, por ejemplo, una rotación del contorno del pozo durante el sobrecrecimiento de Ge, en concordancia con observaciones experimentales recientes.

Introducción

La heteroepitaxia de semiconductores con desajuste de red como Ge/Si o InGaAs/GaAs conduce frecuentemente a la formación de islas tridimensionales mediante el modo de crecimiento Stranski-Krastanow (SK). Aunque el autoensamblaje puede producir puntos cuánticos, la nucleación aleatoria inherente a este proceso a menudo produce una amplia distribución de tamaños y formas que limita el rendimiento de dispositivos.

En las últimas décadas, se han desarrollado diversas técnicas para dirigir el crecimiento heteroepitaxial hacia estructuras ordenadas. Entre ellas, los sustratos con patrones de pozos han surgido como uno de los enfoques más versátiles para lograr tanto un alto ordenamiento lateral como un control preciso del tamaño de las islas.

Los patrones de pozos se fabrican usando métodos top-down como litografía por nanoimpresión, litografía de haz de electrones combinada con grabado reactivo-iónico o químico húmedo, e nanoindentación. Estas técnicas permiten el diseño de arrays de pozos ordenados con precisión nanométrica, que bajo condiciones de crecimiento optimizadas dan lugar a un ordenamiento lateral casi perfecto de islas.

Dado que la forma del pozo influye en la energía superficial total y, por consiguiente, en la nucleación de islas, el control preciso de la morfología del pozo es esencial. A temperaturas elevadas, la capilaridad impulsa la evolución de los pozos hacia configuraciones energéticamente favorables, a veces llevando a su curación completa. Por ello, el recocido o la deposición adicional de material se emplea rutinariamente para estabilizar formas metaestables de pozos. Incluso tras la estabilización, puede ocurrir una evolución adicional durante la fase de crecimiento heteroepitaxial.

En este estudio describimos la evolución de sustratos con patrones de pozos impulsada por reducción de energía superficial mediante difusión superficial. Adoptamos una formulación de campo de fase que captura escalas de longitud y tiempo relevantes para experimentos y se ha aplicado previamente a cinética limitada por difusión en sistemas heteroepitaxiales. El modelo es capaz de reproducir formas de equilibrio con energías superficiales anisotrópicas realistas y ha sido validado contra datos experimentales.

Nos centramos en superficies Si(001) con patrones de pozos, un sistema ampliamente investigado en la literatura. El artículo se organiza así: el modelo de campo de fase se introduce en la siguiente sección; el alisado de pozos en Si(001) y la vía cinética hacia el equilibrio se discuten en la Sección 3; la aplicación al sobrecrecimiento de Ge, que induce rotación del pozo, se presenta en la Sección 4; y el artículo concluye con una discusión de las implicaciones de nuestros hallazgos.

Métodos

Modelo de campo de fase

El modelo de campo de fase introduce un parámetro de orden continuo φ que varía suavemente de φ = 1 (sólido) a φ = 0 (vacío). La evolución de φ se rige por una ecuación de Cahn-Hilliard degenerada que captura la difusión superficial respetando la fuerte anisotropía de la energía superficial. El funcional de energía libre total es:

F = ∫Ω γ(n̂) [ ε/2 |∇φ|² + 1/ε B(φ) ] d𝑟 + ∫Ω β/(2ε) [ -ε∇²φ + 1/ε B′(φ) ]² d𝑟,

donde Ω es el dominio de simulación, ε el espesor de interfase, B(φ)=18φ²(1−φ)² el potencial de doble pozo, y β el parámetro de regularización de Willmore que redondea esquinas agudas en el régimen de fuerte anisotropía. El potencial químico μ=δF/δφ impulsa la evolución a través de:

∂φ/∂t = D ∇·[ M(φ) ∇μ ],

con D el coeficiente de difusión y M(φ)=36/ε φ²(1−φ)² la movilidad superficial. La expresión para μ incluye la densidad de energía superficial anisotrópica γ(n̂) y el término de regularización de Willmore, asegurando una descripción físicamente realista de la evolución de facetas.

Energía superficial anisotrópica

Modelamos la densidad de energía superficial como:

γ(n̂) = γ₀ [ 1 − Σᵢ αᵢ (n̂·mᵢ)^{wᵢ} Θ(n̂·mᵢ) ],

donde los vectores mᵢ corresponden a las direcciones cristalográficas de energía mínima (〈001〉, 〈113〉, 〈110〉, 〈111〉 para silicio). La profundidad de cada mínimo αᵢ se deriva de los valores relativos de energía superficial γᵢ reportados en la literatura, mientras que el parámetro de ancho wᵢ controla la extensión angular de la faceta. Para silicio establecemos γ₀ = 1, α_{〈001〉}=0.15, y valores w de 50 (excepto w_{〈113〉}=100). El parámetro de regularización de Willmore β se elige como 0.005 para resolver esquinas redondeadas manteniendo factibilidad computacional.

Morfología inicial y configuración de simulación

Construimos una geometría inicial de pozo en una superficie (001) definida por una depresión circular de radio L y profundidad H, conectada suavemente al plano circundante mediante una cola gaussiana con ancho σ=L/2. La distancia con signo d(r) desde la interfase se inserta en el perfil tangente hiperbólico φ(r)=½[1−tanh(3 d(r)/ε)] para generar la condición inicial. Se usa un tamaño de dominio típico de L = 1 (en unidades arbitrarias) y espesor de interfase ε = 0.2. Se emplea la caja de herramientas de elementos finitos AMDiS con un esquema de integración temporal semi-implícito y refinamiento de malla en la interfase. Se aplican condiciones de contorno periódicas en direcciones laterales, y condiciones de no flujo en el eje vertical. El coeficiente de difusión D se fija a la unidad, por lo que el tiempo de simulación se reporta en unidades arbitrarias.

Resultados y discusión

Alisado de pozos en Si(001)

La Figura 2 muestra la evolución temprana de pozos con relaciones de aspecto R=0.25 y R=0.5. Los perfiles inicialmente lisos desarrollan espontáneamente facetas correspondientes a los mínimos de la densidad de energía superficial del silicio. Para el pozo somero (R=0.25), una base cuadrada 〈001〉 está rodeada por bordes 〈113〉 y triángulos 〈110〉 estrechos. Para el pozo más profundo (R=0.5), aparecen facetas 〈111〉 adicionales cerca del fondo, y las facetas 〈110〉 se ensanchan. Estas predicciones concuerdan cuantitativamente con observaciones experimentales de superficies Si(001) con pozos.

La Figura 3 traza la evolución a largo plazo del pozo más profundo (R=0.5). La secuencia de vistas en perspectiva y superior ilustra la desaparición progresiva de facetas 〈111〉 empinadas, el engrosamiento de facetas 〈113〉 adyacentes, y la eventual desaparición de facetas 〈110〉, que juntas producen una superficie aplanada. La energía superficial total disminuye monótonamente, como se muestra en la Figura 4, y la cinética sigue una vía común para todo R≤0.5. La escala de tiempo característica para un pozo de decenas de nanómetros a 1100–1200 °C es del orden de horas, consistente con valores de literatura para difusión superficial de silicio.

Imitación del cambio de forma por sobrecrecimiento de Ge

Para evaluar el impacto de depositar un material con un conjunto de facetas de equilibrio diferente, simulamos el sobrecrecimiento de Ge en un sustrato Si(001) con patrón de pozos con R≈0.1. Partiendo de la morfología en t≈5.0 (Fig. 3), reemplazamos la densidad de energía superficial del silicio por una que incluye un mínimo profundo en la dirección 〈105〉, característico de sistemas Ge/Si(001). El parámetro de ancho w_{〈105〉} se aumentó a 500 para capturar las facetas 〈105〉 poco profundas.

La Figura 5 demuestra que la superficie se reorganiza rápidamente: facetas 〈105〉 nuclean entre las facetas 〈113〉 restantes, expandiéndose a expensas de estas. El pozo final está delimitado exclusivamente por facetas 〈105〉, y la vista superior muestra una rotación de 45° del contorno del pozo respecto a la orientación inicial. Esta rotación se ha reportado experimentalmente durante la deposición de Ge y ahora se muestra que surge puramente de la minimización de energía superficial.

Conclusiones

Hemos desarrollado un marco de campo de fase continuo que captura la evolución de sustratos Si(001) con patrones de pozos impulsada por difusión superficial. Al incorporar explícitamente la fuerte anisotropía de la energía superficial del silicio, el modelo reproduce la formación de facetas metaestables, el aplanamiento eventual del pozo, y las vías cinéticas que conducen a estos estados. Extender la descripción de energía superficial para incluir facetas 〈105〉 demuestra que la rotación de contornos de pozos observada durante el sobrecrecimiento de Ge puede explicarse enteramente por reducción de energía superficial. El modelo es general y puede adaptarse fácilmente a otros sustratos re-parameterizando la anisotropía, ofreciendo una herramienta predictiva para diseñar geometrías de pozos que guíen la colocación de islas heteroepitaxiales con precisión nanométrica.