Uma experiência com epitaxia por feixe molecular no Instituto de Física da USP: da organização dos átomos à difração de elétrons e à análise de Arrhenius.
Construindo materiais em escala atômica
Como é possível fabricar um material semicondutor controlando sua composição e espessura em escalas comparáveis às dimensões dos próprios átomos?
Essa pergunta envolve diferentes áreas da Física: estrutura cristalina, mecânica quântica, termodinâmica, difração, física de superfícies e técnicas experimentais de crescimento de materiais.
Em 29 de outubro de 2026, participei de uma atividade experimental no Laboratório de Novos Materiais Semicondutores (LNMS), pertencente ao Departamento de Física dos Materiais e Mecânica (DFMT), do Instituto de Física da Universidade de São Paulo (IFUSP).
O laboratório utiliza uma técnica conhecida como Molecular Beam Epitaxy (MBE), ou epitaxia por feixe molecular, para produzir materiais semicondutores e heteroestruturas com controle preciso de sua composição e espessura. Trata-se de uma infraestrutura utilizada em pesquisas envolvendo poços quânticos, pontos quânticos e dispositivos optoeletrônicos [1].
Durante a atividade, realizamos calibrações de duas fontes de gálio (Ga₁ e Ga₂) e uma de alumínio (Al). Para determinar as taxas de crescimento, utilizamos RHEED (Reflection High-Energy Electron Diffraction), uma técnica baseada na difração de elétrons de alta energia.
Depois, ao analisar os resultados, investiguei a dependência dessas taxas com a temperatura das fontes. Por sugestão do professor Dr. Alain André Quivy, essa investigação foi ampliada utilizando gráficos de Arrhenius, conectando as calibrações experimentais à física da evaporação e dos processos termicamente ativados.
O que começou como uma atividade de calibração revelou uma sequência de perguntas interessantes. O que significa crescer uma monocamada? Como a difração de elétrons permite estudar a superfície de um cristal? Por que a intensidade difratada pode oscilar durante o crescimento? E o que a temperatura de uma fonte de material nos informa sobre a quantidade de átomos que ela fornece?
Para compreender essas questões, precisamos começar pela organização dos próprios cristais.
1. Antes de crescer um cristal, precisamos compreender sua estrutura
Um sólido cristalino não é simplesmente um conjunto de átomos agrupados. Sua característica fundamental é a existência de uma organização periódica no espaço.
Em um cristal ideal, encontramos unidades estruturais que se repetem por meio de translações bem definidas. Essa organização pode ser descrita matematicamente por uma rede de Bravais, à qual associamos uma base atômica.
A rede descreve a geometria da repetição. A base identifica os átomos associados a cada ponto da rede. Juntas, determinam a estrutura cristalina.
Um dos materiais fundamentais na tecnologia de semicondutores é o arseneto de gálio, GaAs. Ele possui estrutura cristalina do tipo zincblenda, constituída por duas sub-redes cúbicas de faces centradas, uma formada por átomos de gálio e outra por átomos de arsênio, deslocadas entre si de um quarto da diagonal espacial da célula cúbica [2].
Essa estrutura é importante porque a disposição dos átomos influencia diretamente as propriedades eletrônicas do material.
A periodicidade do potencial elétrico ao qual os elétrons estão submetidos contribui para a formação das bandas de energia características dos sólidos. Em um semicondutor, a existência de uma região de energias proibidas, conhecida como gap, é essencial para compreender suas propriedades elétricas e ópticas.
O GaAs e o AlGaAs, por exemplo, podem ser combinados para produzir heteroestruturas, isto é, estruturas formadas por regiões de semicondutores com composições diferentes.
Quando essas regiões possuem espessuras suficientemente pequenas, os portadores de carga podem apresentar efeitos de confinamento quântico. Esse é um dos fundamentos físicos de dispositivos baseados em poços quânticos.
Portanto, controlar a deposição de material não é apenas controlar uma quantidade de matéria. É também controlar uma estrutura cuja organização determina suas propriedades.
2. O que diferencia crescimento epitaxial de uma simples deposição?
Imagine que desejamos depositar gálio e arsênio sobre a superfície de um cristal de GaAs.
Não basta que os átomos cheguem ao substrato. Para produzir uma camada epitaxial de boa qualidade, eles precisam incorporar-se à estrutura cristalina de maneira compatível com a organização do material subjacente.
O termo epitaxia descreve justamente o crescimento de uma camada cristalina com uma relação de orientação definida em relação ao substrato.
Durante o crescimento, partículas chegam à superfície, podem ser adsorvidas, deslocam-se por difusão superficial e eventualmente se incorporam à estrutura.
A mobilidade dos átomos, a temperatura do substrato, as condições químicas da superfície e os fluxos incidentes influenciam a maneira como o material cresce.
Em determinadas condições, a deposição ocorre predominantemente camada por camada. Em outras, podem surgir ilhas tridimensionais, mudanças de morfologia ou outros regimes de crescimento.
A diferença é fundamental: queremos não apenas adicionar material, mas produzir uma estrutura cristalina controlada.
É justamente para realizar esse processo com elevado controle experimental que utilizamos uma MBE.
3. Dentro de uma máquina de epitaxia por feixe molecular
Uma MBE opera sob condições de ultra-alto vácuo, permitindo que feixes de partículas provenientes de fontes aquecidas alcancem o substrato com poucas colisões com moléculas do gás residual.
Em sistemas destinados ao crescimento de semicondutores III-V, como GaAs e AlGaAs, essas fontes podem fornecer elementos como Ga, Al e As.
Gálio e alumínio são exemplos de materiais que podem ser evaporados por meio de células de efusão. Nessas células, o aquecimento produz uma pressão de vapor que permite a saída de partículas por uma abertura direcionada à região do substrato.
O processo envolve diferentes grandezas físicas:
- A temperatura das fontes influencia a pressão de vapor e o fluxo de partículas.
- O substrato é aquecido para favorecer os processos necessários ao crescimento cristalino.
- Obturadores, conhecidos como shutters, permitem controlar a exposição do substrato aos diferentes feixes.
- O ambiente de vácuo reduz a contaminação e favorece a propagação dos feixes.
A combinação desses mecanismos permite produzir sequências de camadas com diferentes composições e espessuras.
Entretanto, existe uma dificuldade experimental importante.
Conhecer a temperatura de uma célula de efusão não significa conhecer automaticamente sua taxa de crescimento.
A pressão de vapor depende da temperatura, mas o fluxo efetivo também está relacionado à geometria da fonte, às condições de operação e à incorporação das espécies na superfície.
Além disso, algumas células possuem regiões controladas separadamente. Nas fontes de gálio utilizadas em nossas calibrações, distinguimos as temperaturas de bulk, associadas à região que contém a carga de material, e de tip, associadas à extremidade superior da célula.
Essas temperaturas não devem ser confundidas com a temperatura do substrato, onde o crescimento acontece.
Para utilizar uma fonte com precisão, precisamos estabelecer experimentalmente uma relação entre suas condições de operação e a taxa de crescimento produzida.
Mas como medir o crescimento de camadas tão finas?
4. RHEED: utilizando elétrons para investigar um cristal
O RHEED é uma técnica de difração de elétrons de alta energia em geometria de incidência rasante.
Sua fundamentação está relacionada à natureza ondulatória da matéria.
De acordo com a relação de de Broglie, uma partícula com momento linear p possui comprimento de onda associado:
em que h é a constante de Planck.
Elétrons suficientemente energéticos possuem comprimentos de onda pequenos em comparação com as distâncias características entre átomos de um cristal.
Quando um feixe de elétrons interage com uma estrutura periódica, as amplitudes das ondas espalhadas podem interferir construtiva ou destrutivamente. As condições de interferência dependem da periodicidade cristalina, da geometria do feixe e da transferência de momento envolvida no espalhamento.
A difração, portanto, transporta informações sobre a organização espacial dos átomos.
Na técnica RHEED, os elétrons chegam à superfície com um ângulo pequeno em relação ao plano do cristal. Uma parte dos elétrons espalhados produz um padrão observado em uma tela sensível aos elétrons.
Esse padrão pode conter pontos, linhas alongadas e máximos de diferentes ordens. Suas características dependem da organização da superfície, da sua morfologia e das condições experimentais [3].
Uma ideia essencial é que o RHEED não produz uma fotografia direta dos átomos. O que observamos é um padrão de intensidade resultante do espalhamento eletrônico. Para interpretar esse padrão, precisamos relacioná-lo à estrutura cristalina e à geometria da difração.
Essa distinção entre observação e interpretação é uma parte importante do método científico.
Por que o padrão muda quando giramos o cristal?
Um cristal possui direções e planos identificados por índices cristalográficos, como [100], [110] e (001).
Essas notações descrevem orientações específicas em relação à estrutura periódica.
Em uma superfície cristalina, diferentes direções podem apresentar diferentes periodicidades atômicas. Assim, quando modificamos a orientação azimutal da amostra em relação ao feixe de elétrons, alteramos a geometria de difração observada.
Durante a atividade no LNMS, discutimos direções cristalográficas, rotações e padrões identificados pelas periodicidades 1×, 2× e 4×.
Essas designações estão relacionadas a múltiplos da periodicidade de referência da superfície.
Isso nos conduz a outro fenômeno: as reconstruções superficiais.
A superfície de um cristal não precisa manter exatamente a mesma organização atômica encontrada no seu interior. A interrupção das ligações químicas e as condições de temperatura e composição podem favorecer rearranjos dos átomos superficiais.
No GaAs(001), reconstruções como 2 × 4 e c(4 × 4) são amplamente estudadas. Sua formação depende, entre outros fatores, da temperatura do substrato e das condições de exposição ao arsênio [4].
Durante a atividade, também foram discutidos diagramas de reconstrução com regiões identificadas por letras gregas, como α, β e γ. Existem trabalhos científicos que investigam subfases dessas reconstruções e suas regiões de ocorrência [5].
A identificação precisa dos padrões e diagramas discutidos naquela ocasião, porém, depende de informações experimentais adicionais. Por isso, é importante não associar automaticamente uma determinada posição angular ou periodicidade observada a uma reconstrução específica.
5. Como uma oscilação de intensidade pode medir o crescimento?
Além de fornecer informações estruturais, o RHEED pode ser utilizado para acompanhar a evolução temporal da superfície durante o crescimento.
Considere um regime aproximadamente camada por camada.
Inicialmente, podemos ter uma superfície relativamente plana. Quando o crescimento começa, átomos incorporados à superfície formam ilhas bidimensionais que se expandem progressivamente.
À medida que a camada é preenchida, a distribuição de degraus, ilhas e regiões planas muda.
Essas modificações alteram as condições de espalhamento dos elétrons e podem produzir oscilações periódicas na intensidade de determinadas características do padrão RHEED.
Trabalhos clássicos sobre crescimento por MBE demonstraram a relação entre essas oscilações e a dinâmica de formação das camadas cristalinas [3].
Quando cada período corresponde à deposição de uma unidade de camada definida para a calibração, a taxa pode ser determinada por:
em que:
- R é a taxa de crescimento, em monocamadas por segundo (ML/s);
- τ é o período das oscilações, em segundos.
Por exemplo, se uma oscilação completa corresponde à formação de uma monocamada e dura dois segundos, a taxa será de 0,5 ML/s.
A relação parece simples, mas sua validade depende do regime de crescimento, da característica difratada monitorada e da convenção adotada para a monocamada.
No crescimento de semicondutores compostos, esse último detalhe é particularmente relevante: uma monocamada pode ser definida de maneiras diferentes conforme o contexto cristalográfico e experimental.
É por isso que uma curva de intensidade em função do tempo pode funcionar como instrumento de medição, mas exige interpretação física e condições de aquisição adequadas.
Foi utilizando o RHEED que determinamos as taxas associadas às fontes de Ga₁, Ga₂ e Al durante a atividade de 29 de outubro.
6. Das medições às curvas de calibração
Depois de obter as taxas, o passo seguinte é investigar como elas variam com a temperatura das fontes.
Os registros da atividade contêm 27 medições: dez para Ga₁, dez para Ga₂ e sete para Al.
Nas séries de Ga₁ e Ga₂, dois pontos foram obtidos com temperaturas de tip diferentes das condições principais. Esses pontos foram preservados nos gráficos, mas não incluídos nos ajustes principais, para evitar combinar condições operacionais distintas em um modelo que utiliza apenas a temperatura de bulk como variável independente.
Para cada fonte, foi ajustado um modelo exponencial da forma:
em que Rref é a taxa na temperatura de referência Tref, e b representa a sensibilidade logarítmica da taxa à temperatura.
A utilidade prática é direta: depois de calibrada uma fonte, podemos estimar qual temperatura deverá produzir determinada taxa, desde que respeitemos as condições e os limites de validade da calibração.
Os ajustes foram realizados pelo método dos mínimos quadrados ortogonais ponderados (ODR), que permite considerar incertezas tanto nas temperaturas quanto nas taxas.
Esse cuidado é importante porque nenhuma medição experimental é infinitamente precisa.
Uma curva ajustada pode passar muito próxima dos dados e, ainda assim, não representar exatamente o mecanismo físico responsável pelo comportamento observado.
É necessário distinguir a qualidade matemática do ajuste de sua interpretação física.
7. Da calibração à análise de Arrhenius
Após as calibrações operacionais, o professor Alain André Quivy sugeriu investigar a dependência das taxas com a temperatura utilizando a representação de Arrhenius.
A motivação é física.
Processos como evaporação, difusão e outras transformações termicamente ativadas frequentemente apresentam dependências aproximadamente exponenciais com o inverso da temperatura absoluta.
Um modelo de Arrhenius pode ser escrito como:
em que T deve ser expressa em kelvin, kB é a constante de Boltzmann e Ea,app é uma energia de ativação aparente.
Ao representar o logaritmo da taxa em função do inverso da temperatura absoluta, esse modelo produz uma reta cuja inclinação está relacionada a essa energia.
Utilizando os dados de 29 de outubro, os ajustes de Arrhenius forneceram os seguintes resultados:
| Fonte | Energia de ativação aparente |
|---|---|
| Ga₁ | (2,232 ± 0,033) eV |
| Ga₂ | (2,247 ± 0,048) eV |
| Al | (3,22 ± 0,13) eV |
Esses resultados pertencem à análise v3 realizada posteriormente à atividade experimental. As incertezas indicadas são estimativas estatísticas condicionadas às incertezas atribuídas aos dados e ao modelo utilizado.
Os valores obtidos para as duas fontes de gálio são compatíveis dentro das incertezas reportadas. O alumínio apresentou uma energia aparente maior.
Entretanto, é importante evitar uma interpretação excessiva.
A energia extraída do gráfico não deve ser automaticamente identificada com a entalpia de vaporização ou sublimação do elemento.
A relação entre temperatura, pressão de vapor e fluxo de partículas envolve hipóteses termodinâmicas e cinéticas. Além disso, o que medimos por RHEED é uma taxa de crescimento incorporado, não a pressão de vapor diretamente.
A análise de Arrhenius oferece, portanto, uma maneira de investigar a dependência térmica dos resultados e construir hipóteses físicas. Sua interpretação microscópica exige considerar as características das células, o transporte de partículas e as condições de incorporação no substrato.
Esse é um exemplo de como uma calibração inicialmente operacional pode motivar uma investigação física mais profunda.
8. O que as medições nos dizem — e o que ainda precisamos investigar?
As calibrações demonstram que podemos relacionar quantitativamente as condições de operação das fontes às taxas de crescimento medidas.
Os modelos exponenciais e de Arrhenius descrevem bem os dados nos intervalos investigados. Entretanto, isso não significa que possam ser extrapolados indefinidamente, nem que um bom ajuste estatístico seja suficiente para identificar um mecanismo físico único.
Também é necessário compreender a origem das incertezas, as condições mantidas constantes entre as medições e a reprodutibilidade das taxas previstas.
As próprias observações de RHEED levantam questões adicionais sobre orientação cristalográfica, reconstruções superficiais e dinâmica do crescimento.
Essas questões não diminuem o valor dos resultados. Pelo contrário: ajudam a definir quais conclusões são sustentadas pelas medições e quais precisam ser testadas.
Uma etapa futura especialmente interessante será utilizar as curvas para prever novas condições de operação e verificar experimentalmente se as taxas obtidas correspondem às previsões.
Esse tipo de comparação entre modelo e experimento é essencial para avaliar a utilidade de uma calibração.
9. Uma experiência que conecta diferentes áreas da Física
Uma das características mais interessantes dessa atividade foi a quantidade de conceitos que aparecem conectados em um único procedimento experimental.
A cristalografia descreve a organização dos átomos.
A mecânica quântica permite compreender a natureza ondulatória dos elétrons.
A difração transforma a periodicidade cristalina em padrões observáveis.
A física de superfícies ajuda a interpretar reconstruções, ilhas e degraus atômicos.
A termodinâmica contribui para compreender a dependência da evaporação com a temperatura.
E a análise de dados permite transformar medições em modelos quantitativos, acompanhados de suas incertezas.
Em conjunto, esses conhecimentos tornam possível investigar e controlar o crescimento de materiais semicondutores.
A experiência no LNMS, em 29 de outubro de 2026, tornou-se assim o ponto de partida para um estudo mais amplo sobre a Física da epitaxia por feixe molecular.
Pretendo desenvolver esses temas em publicações posteriores, começando pela geometria dos cristais, avançando para o RHEED e retornando às calibrações e à interpretação de Arrhenius.
Mais do que apresentar o funcionamento de uma máquina ou os resultados de um ajuste, o objetivo é compreender como relacionamos a organização microscópica da matéria a observações experimentais e modelos físicos.
Afinal, medir não significa apenas obter números. Significa estabelecer, com hipóteses e incertezas explícitas, uma conexão verificável entre aquilo que observamos e aquilo que afirmamos sobre a natureza.
Referências
- INSTITUTO DE FÍSICA DA UNIVERSIDADE DE SÃO PAULO. Laboratório de Novos Materiais Semicondutores (LNMS). Informações institucionais e linhas de pesquisa. Disponível em: portal.if.usp.br/mbe/pt-br/node/326.
- ASHCROFT, N. W.; MERMIN, N. D. Solid State Physics. Philadelphia: Holt, Rinehart and Winston, 1976. Especialmente capítulos 4–6, sobre estruturas cristalinas, rede recíproca e difração.
- NEAVE, J. H.; JOYCE, B. A.; DOBSON, P. J.; NORTON, N. Dynamics of film growth of GaAs by MBE from RHEED observations. Applied Physics A, v. 31, p. 1–8, 1983. DOI: 10.1007/BF00617180.
- OHTAKE, A. Surface reconstructions on GaAs(001). Surface Science Reports, v. 63, p. 295–327, 2008. DOI: 10.1016/j.surfrep.2008.03.001.
- LABELLA, V. P.; BULLOCK, D. W.; DING, Z.; EMERY, C.; THIBADO, P. M. Enabling electron diffraction as a tool for determining substrate temperature and surface morphology. Applied Physics Letters, v. 79, p. 3065–3067, 2001. DOI: 10.1063/1.1416477.
- AL-ZEIDAN, Ahmad. Otimização de pontos quânticos de submonocamada de InAs/GaAs crescidos sobre GaAs(001) com uma reconstrução de superfície (2 × 4) para fotodetectores de radiação infravermelha. Tese de Doutorado, Instituto de Física, Universidade de São Paulo, 2023. Especialmente capítulo 3, dedicado ao crescimento por MBE, ao RHEED e às reconstruções superficiais.
- Registros experimentais de calibração MBE — LNMS, 29/10/2026. Dados de temperatura e taxa de crescimento das fontes Ga₁, Ga₂ e Al, reunidos no arquivo
20261029_calibracoes_mbe_dados_com_incertezas_bruto.csv. Tratamento estatístico e calibrações operacionais v2.4; análise de Arrhenius v3, de 31/10/2026. Documentação experimental fornecida pelo autor, ainda não disponibilizada publicamente.