Visão em Escala Atômica da Deformação por Escorregamento e Maquinabilidade Nanométrica do 6H-SiC

Resumo

Compreender a deformação em escala atômica e a remoção de material do 6H-SiC é essencial para produzir superfícies ultra-lisas, livres de danos, com degraus atômicos. Como a observação direta da zona de usinagem em nanoescala está além do alcance experimental atual, empregamos simulações de dinâmica molecular em larga escala para capturar fenômenos-chave como escorregamento de discordâncias, transições de fase e separação de material. Este estudo explora sistematicamente como a anisotropia cristalográfica influencia a deformação por escorregamento e a maquinabilidade nanométrica, avançando o conhecimento fundamental para processamento de precisão do 6H-SiC.

Introdução

O carbeto de silício (SiC) é pedra angular da tecnologia de semicondutores de terceira geração, valorizado por sua larga banda proibida, alto campo de ruptura, superior tolerância à radiação, alta velocidade de saturação de portadores, condutividade térmica excepcional, baixa constante dielétrica e robusta estabilidade química. Esses atributos tornam o SiC indispensável em aplicações de alta temperatura, alta frequência, alta potência, anti-radiação e optoeletrônica de curto comprimento de onda, incluindo retificadores Schottky, tiratrões e MOSFETs.

Entre os politipos de SiC, 3C, 4H e 6H são os mais usados. Usinagem convencional — retífica, lapidação e polimento — continua sendo a principal rota de fabricação para SiC monocristalino. Contudo, a razão de dureza entre diamante e SiC (~2:1 para profundidades de processamento <50 nm) fica aquém da recomendada 5:1, levando a desgaste acelerado da ferramenta e dano subsuperficial que comprometem a qualidade do wafer. Pesquisa extensiva em 3C-SiC elucidou mecanismos de deformação plástica, desgaste de ferramenta, comportamento de atrito, anisotropia e efeitos de temperatura durante o corte. Em contraste, o 6H-SiC, com seu empilhamento complexo ABCACB e banda proibida de 3 eV, recebeu atenção comparativamente menor, apesar de sua superior maturidade tecnológica e prevalência industrial.

Para melhorar a integridade de superfície e subsuperfície do 6H-SiC, é crucial identificar o plano cristalino ideal (superfície de usinagem) e orientação (direção de usinagem) que minimizem danos enquanto maximizem eficiência de remoção. Experimentos e simulações de risco são ferramentas padrão para sondar remoção de material em nanoescala, oferecendo insights que guiam a usinagem abrasiva. Aqui, empregamos simulações de risco por dinâmica molecular para dissecar a influência da anisotropia cristalográfica na maquinabilidade nanométrica do 6H-SiC.

Metodologia

Todas as simulações de corte foram realizadas com o Large-Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS). Visualização do modelo e identificação de defeitos utilizaram OIVTO e algoritmo de detecção de estrutura diamantina. Peça e ferramenta foram modeladas como corpos totalmente deformáveis, permitindo capturar evolução do desgaste durante o risco. Cada modelo compreendia camada de contorno, camada termostática e região newtoniana. Átomos na camada de contorno na base e lateral direita da peça foram fixados; os demais obedeceram à segunda lei de Newton. Contornos periódicos foram aplicados na direção y. Antes do risco, relaxamento NVE de 50 ps com termostato Berendsen assegurou estado de energia estável.

A ponta abrasiva era uma coroa piramidal triangular esférica com ângulo aresta-aresta de 90°. Selecionamos os três planos mais comuns do 6H-SiC — basal (plano a), prismático (plano m) e plano c — como superfícies de usinagem. Considerando simetria cristalina, examinamos seis combinações superfície/orientação: (0001)/[2-1-1-0], (0001)/[10-1-0], (01-1-0)/[2-1-1-0], (01-1-0)/[0001], (11-2-0)/[1-1-00] e (11-2-0)/[0001]. Parâmetros de simulação estão na Tabela 1. Antes do risco, o abrasivo foi posicionado 50 Å abaixo da superfície superior, 30 Å da peça, bem fora do raio de corte do potencial interatômico. O abrasivo moveu-se ao longo de +x para completar o risco.

O potencial analítico de ordem de ligação (ABOP) desenvolvido por Erhart e Albe foi escolhido para descrever interações Si-Si, C-C e Si-C devido à sua precisão comprovada para propriedades mecânicas e mecanismos de remoção do SiC. Parâmetros na Tabela 2.

Resultados e Discussão

Análise da Maquinabilidade Nanométrica

As figuras demonstram que a anisotropia do 6H-SiC afeta fortemente a maquinabilidade nanométrica — profundidade usinada, modo de remoção, quantidade de material removido e profundidade de dano subsuperficial (SSD). Relacionando superfície e direção de usinagem ao plano basal/eixo c, classificamos os modos de processamento em três grupos: (i) plano basal como superfície de usinagem; (ii) plano basal perpendicular à superfície com eixo c perpendicular à direção de usinagem; (iii) eixo c paralelo à direção de usinagem.

Observações-chave:

Modos (0001)/[2-1-1-0] e (0001)/[10-1-0] produziram maquinabilidade similar, mas o último favoreceu modo de remoção frágil.
Modos (01-1-0)/[2-1-1-0] e (11-2-0)/[1-1-00] atingiram maior profundidade usinada e taxas de remoção. Para profundidade teórica de 5 nm, (11-2-0)/[1-1-00] removeu 3,4× mais material que (01-1-0)/[0001], com razão remoção-desgaste de 10,1, embora SSD fosse 2,3× maior que (0001)/[10-1-0]. Este modo adequa-se a processamento de nanocanais de alta taxa onde dano subsuperficial é aceitável. Em contraste, (01-1-0)/[0001] oferece taxa de remoção comparável mas apenas 50% do SSD, resultando em canal mais liso e mecanicamente robusto — ideal para usinagem de precisão.
Quando a direção de usinagem alinha-se com o eixo c, a ponta abrasiva desgasta-se fortemente no início, levando a menores taxas de remoção (≈1,0) e maior SSD. Embora exibam excelente resistência ao desgaste, não são recomendados para fabricação de micro/nano-canais, mas podem ser valiosos para projeto de faces de avanço de ferramentas de corte.

Análise do Movimento de Escorregamento e Distribuição de Dano Subsuperficial

Distribuição do Fator de Schmid no Processo de Risco

O SiC hexagonal apresenta sistemas de escorregamento basal, prismático e piramidal. A propensão ao escorregamento depende da energia generalizada de falha de empilhamento e do fator de Schmid. Escorregamento basal tipicamente ocorre no conjunto shuffle <11-20>, enquanto escorregamento prismático ativa-se ao longo das direções e . As figuras ilustram os sistemas de escorregamento e fatores de Schmid para cada direção de carregamento.

A tensão de cisalhamento resolvida τ_ss em um plano de deslizamento é τ_ss = σ_cont · m, onde σ_cont é a tensão de contato e m o fator de Schmid. O sistema de coordenadas global é fixo enquanto o local segue a orientação cristalina.

Durante o risco, a face de avanço domina o contato até que desgaste severo da ponta desloque o contato primário para um composto de face de avanço e topo da ponta, alterando a direção de carregamento. As direções de carregamento são f₁=(1,0,-√2) para face de avanço e f₂=(0,0,-1) para topo da ponta.

Usando parâmetros cristalográficos (c/a≈4,901), escorregamento basal requer menor tensão de cisalhamento crítica que escorregamento prismático. Assim, ao usinar o plano basal (0001)/[2-1-1-0] ou (0001)/[10-1-0], escorregamento basal inicia primeiro. Para orientações (01-1-0)/[2-1-1-0] e (11-2-0)/[1-1-00], escorregamento basal é suprimido e escorregamento prismático domina. Em modos paralelos ao eixo c, desgaste severo da ponta eleva influência de f₂, levando a escorregamento prismático simétrico no plano YZ.