Evolução de Morfologia Dirigida por Difusão Superficial de Substratos Si(001) com Padrão de Cavidades: Estudo por Campo de Fase

Resumo

O ordenamento lateral de ilhas heteroepitaxiais é rotineiramente alcançado por padrão de cavidades no substrato antes da deposição. Controlar a forma, orientação e profundidade dessas cavidades não é trivial porque geometrias metaestáveis evoluem durante recozimento ou crescimento subsequente. Neste trabalho, apresentamos um modelo contínuo de campo de fase que prevê morfologias típicas de cavidades em Si(001) como função da profundidade e forma inicial. O modelo incorpora difusão superficial e densidade de energia superficial fortemente anisotrópica que reproduz a forma de Wulff do silício. Reproduz com sucesso configurações metaestáveis observadas experimentalmente e revela um caminho cinético comum para o achatamento global da cavidade, independentemente da profundidade inicial. Estendendo o modelo para heteroepitaxia de dois materiais com facetas de equilíbrio distintas, mostramos, por exemplo, uma rotação do contorno da cavidade durante sobrecrescimento de Ge, em concordância com observações experimentais recentes.

Introdução

Heteroepitaxia de semicondutores com descompasso de rede como Ge/Si ou InGaAs/GaAs frequentemente leva à formação de ilhas tridimensionais via modo de crescimento Stranski-Krastanow (SK). Embora a auto-organização possa produzir pontos quânticos, a nucleação aleatória inerente a este processo frequentemente resulta em ampla distribuição de tamanho e forma que limita o desempenho do dispositivo.

Nas últimas décadas, várias técnicas foram desenvolvidas para direcionar o crescimento heteroepitaxial para estruturas ordenadas. Entre elas, substratos com padrão de cavidades surgiram como uma das abordagens mais versáteis para alcançar tanto alto ordenamento lateral quanto controle preciso de tamanho de ilhas.

Padrões de cavidades são fabricados usando métodos top-down como nanolitografia por impressão, litografia de feixe de elétrons combinada com gravação iônica reativa ou química úmida, e nanoindentação. Essas técnicas permitem projeto de arrays de cavidades ordenadas com precisão nanométrica, que, sob condições de crescimento otimizadas, dão origem a ordenamento lateral quase perfeito de ilhas.

Como a forma da cavidade influencia a energia superficial total e, consequentemente, a nucleação de ilhas, o controle preciso da morfologia da cavidade é essencial. A temperaturas elevadas, capilaridade impulsiona a evolução das cavidades para configurações energeticamente favoráveis, às vezes levando à cura completa. Portanto, recozimento ou deposição de material adicional é rotineiramente empregado para estabilizar formas de cavidades metaestáveis. Mesmo após estabilização, evolução adicional pode ocorrer durante a fase de crescimento heteroepitaxial.

Neste estudo, descrevemos a evolução de substratos com padrão de cavidades impulsionada por redução de energia superficial via difusão superficial. Adotamos uma formulação de campo de fase que captura escalas de comprimento e tempo relevantes a experimentos e foi previamente aplicada a cinética limitada por difusão em sistemas heteroepitaxiais. O modelo é capaz de reproduzir formas de equilíbrio com energias superficiais anisotrópicas realistas e foi validado contra dados experimentais.

Focamos em superfícies Si(001) com padrão de cavidades, sistema extensivamente investigado na literatura. O artigo está organizado da seguinte forma: o modelo de campo de fase é introduzido na próxima seção; o alisamento de cavidades Si(001) e o caminho cinético para equilíbrio são discutidos na Seção 3; a aplicação ao sobrecrescimento de Ge, que induz rotação da cavidade, é apresentada na Seção 4; e o artigo conclui com discussão das implicações de nossos achados.

Métodos

Modelo de Campo de Fase

O modelo de campo de fase introduz um parâmetro de ordem contínuo φ que varia suavemente de φ = 1 (sólido) a φ = 0 (vácuo). A evolução de φ é governada por uma equação de Cahn-Hilliard degenerada que captura difusão superficial respeitando a forte anisotropia da energia superficial. O funcional de energia livre total lê-se

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

onde Ω é o domínio de simulação, ε a espessura da interface, B(φ)=18φ²(1−φ)² o potencial de poço duplo, e β o parâmetro de regularização Willmore que arredonda cantos agudos no regime de forte anisotropia. O potencial químico μ=δF/δφ impulsiona a evolução através de

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

com D o coeficiente de difusão e M(φ)=36/ε φ²(1−φ)² a mobilidade superficial. A expressão para μ inclui a densidade de energia superficial anisotrópica γ(n̂) e o termo de regularização Willmore, garantindo descrição fisicamente realista da evolução de facetas.

Energia Superficial Anisotrópica

Modelamos a densidade de energia superficial como

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

onde os vetores mᵢ correspondem às direções cristalográficas de energia mínima (〈001〉, 〈113〉, 〈110〉, 〈111〉 para silício). A profundidade de cada mínimo αᵢ é derivada dos valores relativos de energia superficial γᵢ reportados na literatura, enquanto o parâmetro de largura wᵢ controla a extensão angular da faceta. Para silício definimos γ₀ = 1, α_{〈001〉}=0,15, e valores w de 50 (exceto w_{〈113〉}=100). O parâmetro de regularização Willmore β é escolhido como 0,005 para resolver cantos arredondados mantendo viabilidade computacional.

Morfologia Inicial e Configuração de Simulação

Construímos geometria inicial de cavidade em superfície (001) definida por depressão circular de raio L e profundidade H, conectada suavemente ao plano circundante usando cauda gaussiana com largura σ=L/2. A distância assinada d(r) da interface é inserida no perfil tangente hiperbólico φ(r)=½[1−tanh(3 d(r)/ε)] para gerar condição inicial. Tamanho de domínio típico L = 1 (unidades arbitrárias) e espessura de interface ε = 0,2 são usados. Toolbox de elementos finitos AMDiS empregado com esquema de integração temporal semi-implícito e refinamento de malha na interface. Condições de contorno periódicas nas direções laterais e sem fluxo no eixo vertical. Coeficiente de difusão D definido como unidade, então tempo de simulação reportado em unidades arbitrárias.

Resultados e Discussão

Alisamento de Cavidades Si(001)

Figura 2 mostra evolução em tempo inicial de cavidades com razões de aspecto R=0,25 e R=0,5. Perfis inicialmente sem feições desenvolvem espontaneamente facetas correspondentes aos mínimos da densidade de energia superficial do silício. Para cavidade rasa (R=0,25), base quadrada 〈001〉 é cercada por arestas 〈113〉 e triângulos 〈110〉 estreitos. Para cavidade mais profunda (R=0,5), facetas 〈111〉 adicionais aparecem próximo ao fundo, e facetas 〈110〉 se alargam. Essas previsões concordam quantitativamente com observações experimentais de superfícies Si(001) com cavidades.

Figura 3 traça evolução em longo tempo da cavidade mais profunda (R=0,5). Sequência de vistas em perspectiva e topo ilustra desaparecimento progressivo de facetas 〈111〉 íngremes, coalescência de facetas 〈113〉 adjacentes, e eventual desaparecimento de facetas 〈110〉, que juntas produzem superfície achatada. A energia superficial total diminui monotonicamente, como mostrado na Figura 4, e a cinética segue caminho comum para todos R≤0,5. A escala de tempo característica para cavidade de dezenas de nanômetros a 1100–1200 °C é da ordem de horas, consistente com valores da literatura para difusão superficial de silício.

Mimetizando Mudança de Forma devida a Sobrecrescimento de Ge

Para avaliar impacto de depositar material com conjunto de facetas de equilíbrio diferente, simulamos sobrecrescimento de Ge em substrato Si(001) com padrão de cavidade R≈0,1. Partindo da morfologia em t≈5,0 (Fig. 3), substituímos a densidade de energia superficial do silício por outra que inclui mínimo profundo na direção 〈105〉, característico de sistemas Ge/Si(001). O parâmetro de largura w_{〈105〉} foi aumentado para 500 para capturar facetas 〈105〉 rasas.

Figura 5 demonstra que a superfície se reorganiza rapidamente: facetas 〈105〉 nucleam entre as facetas 〈113〉 remanescentes, expandindo-se em detrimento destas. A cavidade final é limitada exclusivamente por facetas 〈105〉, e a vista superior mostra rotação de 45° do contorno da cavidade relativa à orientação inicial. Essa rotação foi reportada experimentalmente durante deposição de Ge e agora mostrada como resultante puramente de minimização de energia superficial.

Conclusões

Desenvolvemos uma estrutura contínua de campo de fase que captura a evolução dirigida por difusão superficial de substratos Si(001) com padrão de cavidades. Ao incorporar explicitamente a forte anisotropia da energia superficial do silício, o modelo reproduz a formação de facetas metaestáveis, o eventual achatamento da cavidade e os caminhos cinéticos que levam a esses estados. Estender a descrição de energia superficial para incluir facetas 〈105〉 demonstra que a rotação de contornos de cavidades observada durante sobrecrescimento de Ge pode ser explicada inteiramente por redução de energia superficial. O modelo é geral e pode ser prontamente adaptado a outros substratos reparametrizando a anisotropia, oferecendo ferramenta preditiva para projetar geometrias de cavidades que guiem posicionamento de ilhas heteroepitaxiais com precisão nanométrica.