Investigação por Dinâmica Molecular do Corte Diamantado de Precisão Ultra do Cério

Resumo

O acoplamento entre transformações de fase estruturais e dinâmica de discordâncias complica o entendimento da deformação nanoescalas do cério. Aqui empregamos dinâmica molecular (DM) para desvendar o mecanismo do cério durante corte diamantado de precisão ultra. Um modelo DM incorporando potenciais Ce–Ce e Ce–C derivados empiricamente é construído, e a resposta elástica de duas fases fcc de cério é avaliada. Simulações mostram que escorregamento de discordâncias domina deformação plástica, enquanto análise de função de distribuição radial revela apenas mudanças menores de fase γ→δ na superfície usinada e cavaco. Estudos de parâmetros identificam condições ideais de usinagem para alcançar superfície de alta qualidade.

Introdução

Cério (Ce, Z=58) é um metal lantânido valorizado por suas propriedades mecânicas, físicas e químicas únicas. Morfologia de superfície influencia criticamente desempenho, resistência à corrosão e comportamento de armazenamento de hidrogênio de peças de cério [1–3]. Corte diamantado de precisão ultra pode produzir superfícies ultra-lisas com dano subsuperficial mínimo [6,7], mas pouco trabalho experimental ou teórico abordou corte de cério. Como o raio de aresta da ferramenta em corte de precisão ultra é comparável à profundidade de corte, as propriedades do material dominam o processo, tornando o mecanismo de corte de cério particularmente desafiador.

Ce exibe diagrama de fase pressão-temperatura rico impulsionado por deslocalização de elétrons 4f. A pressão ambiente e <110 K a fase α-Ce (fcc) é estável, transformando a β-Ce (dhcp) entre 45–275 K, γ-Ce (fcc) entre 270–999 K, e δ-Ce (bcc) entre 999–1071 K [8–11]. Notavelmente, a transformação γ→α a 295 K e 8 kbar colapsa volume em ~20% [12–14], enquanto a transição γ→δ pode ser desencadeada por alta temperatura e pressão. Essas transformações alteram ligações e propriedades elásticas, afetando deformação. Ce metálico também é dúctil, governado por atividade de discordâncias [15]. No entanto, a interação entre mudanças de fase e discordâncias durante corte diamantado permanece largamente inexplorada.

Simulação DM oferece rota poderosa para sondar mecanismos fundamentais de usinagem em escala atômica. Estudos prévios examinaram nanorisco e corte ortogonal de Cu, Al, Fe e ligas de alta entropia, revelando insights sobre nucleação de discordâncias, acúmulo superficial e comportamento de atrito [16–27]. Tal trabalho destaca importância de orientação cristalina, ângulo de inclinação da ferramenta e potenciais de interação. Ainda nenhum estudo DM investigou usinagem mecânica de cério, motivando o presente trabalho.

Métodos

Modelo DM de Corte Diamantado

A Figura 1 ilustra a configuração da simulação: peça de trabalho de Ce monocristalino (41 × 25 × 31 nm) contendo ~1 × 10⁶ átomos de γ-Ce, e ferramenta de corte diamantada (0,1 × 10⁶ átomos C) com ângulo de alívio de 9°. A peça é dividida em camada fixa inferior de 2 nm e camada móvel que segue dinâmica newtoniana com passo de tempo de 1 fs. Contornos periódicos são aplicados apenas na direção longitudinal. Três orientações cristalográficas — (010), (110) e (111) — são consideradas. Sete ângulos de inclinação (−30°, −20°, −10°, 0°, 10°, 20°, 30°) são testados. A ferramenta diamantada é tratada como corpo rígido; assim interações C–C são omitidas.

Interações atômicas são modeladas da seguinte forma: forças Ce–Ce usam parametrização do método do átomo embutido (EAM) de Sheng et al. [29], reproduzindo com precisão constantes elásticas de Ce fcc. Interações Ce–C adotam potencial de Morse: E_total = Σ D₀[exp(-2α(r-r₀)) - 2exp(-α(r-r₀))] com D₀=0,087 eV, α=5,14, r₀=2,93 Å; corte de 1,0 nm aplicado [30].

O sistema é primeiro equilibrado a 30 K, 0 bar no ensemble NPT, então submetido a corte diamantado a velocidade constante de 100 m/s e profundidade de corte 4 nm no ensemble NVT. Forças de corte são resolvidas na direção de corte, com direção indicada por setas vermelhas nas figuras. Embora dimensões simuladas e velocidades de corte sejam ordens de magnitude menores que condições experimentais de precisão ultra, o modelo DM captura a física essencial do processo. Análise de vizinhos comuns (CNA) distingue átomos fcc, hcp, bcc e defeito [31]. Simulações realizadas com LAMMPS (passo 1 fs) e visualizadas usando OVITO [32,33].

Caracterização das Fases do Cério

Cinco fases de cério — γ, α, β, ε e δ — são modeladas usando parâmetros de rede da literatura [8–11]. Configurações em massa passam por tração uniaxial, cisalhamento e compressão para extrair constantes elásticas e propriedades mecânicas (Tabela 2). O módulo de Young calculado de γ-Ce (24,17 GPa) alinha-se com dados experimentais de nanoindentação (36,7 GPa) [10], e constantes elásticas concordam com resultados de espalhamento de nêutrons inelástico [34]. α-Ce mais denso exibe rigidez significativamente maior que γ-Ce.

Análise de função de distribuição radial (RDF) caracteriza diferenças estruturais entre fases. Figura 2a mostra picos RDF distintos para cada fase, correspondendo a distâncias de vizinhos mais próximos (ex: γ-Ce: 3,64, 5,13, 6,30 Å; α-Ce: 3,41, 4,85, 5,92 Å). Compressão uniforme de γ-Ce para colapso de volume de 20% resulta em RDF correspondente ao de α-Ce (Figura 2b), confirmando transição γ→α [12–14].

Resultados e Discussão

Mecanismos de Usinagem do Cério

Simulações DM de corte diamantado em Ce(010) com ângulo de inclinação 0° revelam três regimes de força (Figura 3). Inicialmente, forças adesivas produzem forças de corte e normais negativas; deformação elástica segue, levando a rápido aumento nas forças. Em ~2,3 nm de comprimento de corte, queda aguda indica início de plasticidade, dominada por discordâncias parciais de Shockley 1/6<112> nucleadas próximo à superfície livre e deslizando em planos {111}. Forças de corte e normal estabilizam entre 10–35 nm, enquanto discordâncias se acumulam à frente e sob a ferramenta. Além de 35 nm, ferramenta separa da peça, causando redução de força. Ao longo do processo, força normal permanece menor que força de corte.

Figura 4 retrata evolução de estruturas de defeito e morfologia de superfície em comprimentos de corte chave. Parciais de Shockley geram falhas de empilhamento limitadas por núcleos de discordância. À medida que corte progride, nuvem densa de parciais emana da superfície frontal, formando cavacos ao longo da face de inclinação (Figura 4f). Discordâncias atrás da ferramenta ascendem à superfície superior, promovendo recuperação superficial. Quando ferramenta aproxima-se do limite esquerdo, propagação de discordâncias é dificultada, aumentando volume de cavaco (Figura 4g). Após separação do cavaco, densidade de discordâncias diminui devido a aniquilação na superfície livre (Figura 4h).

Análise RDF de zonas de defeito confirma transformação γ→δ limitada. RDF da região de defeito sob superfície alinha-se com picos de δ-Ce, indicando transição γ→δ menor impulsionada por aquecimento local (Figura 5a). Assinaturas similares de δ-Ce aparecem no cavaco (Figura 5b), mas quantidade geral permanece trivial, sugerindo que mudanças de fase não são dominantes durante corte.

Influência do Ângulo de Inclinação

Figura 6 mostra que tanto forças de corte quanto normais diminuem com aumento do ângulo de inclinação, de acordo com teoria de Merchant de energetica do plano de cisalhamento [37]. Ângulos de inclinação maiores também reduzem diferencial de força e melhoram estabilidade de usinagem.

Análise de defeitos (Figura 7a,b) revela que ângulo de inclinação negativo (−30°) resulta em maior densidade de discordâncias que ângulo positivo (30°). Número de segmentos de discordância (Figura 7c) declina monoticamente com ângulo de inclinação, refletindo deformação plástica menos complexa e menor endurecimento por trabalho.

Morfologia de superfície (Figura 8) corrobora achados: ângulos de inclinação negativos produzem acúmulo superficial assimétrico e pronunciado, enquanto inclinação de 30° gera menor volume de acúmulo e perfil simétrico. Consequentemente, ângulo de inclinação de 30° é ótimo para alcançar baixa força de corte, densidade mínima de discordâncias e acabamento superficial superior.

Influência da Orientação Cristalográfica

Sob inclinação ótima de 30°, comparamos orientações Ce(010), Ce(110) e Ce(111) (Figura 9). Ce(010) exibe menor força de corte e contagem de discordâncias. Embora Ce(111) ofereça menos discordâncias, seu plano {111} paralelo à superfície promove acúmulo superficial, aumentando resistência de corte. Morfologia de superfície usinada (Figura 10) mostra que Ce(010) resulta no menor volume de acúmulo e perfil mais simétrico, enquanto Ce(110) sofre maior acúmulo.

Conclusões

Desenvolvemos estrutura DM incorporando potenciais EAM e Morse para modelar corte diamantado de precisão ultra de cério. Constantes elásticas, propriedades mecânicas e comportamento de transição de fase de fases de Ce são validadas, confirmando capacidade preditiva do modelo. Simulações revelam que nucleação e deslizamento de discordâncias dominam plasticidade, com apenas mudanças menores de fase γ→δ observadas. Condições ideais de usinagem — ângulo de inclinação 30° e orientação cristalográfica (010) — produzem menores forças de corte, densidade mínima de discordâncias e maior qualidade de superfície.

Referências

[1]–[38] (ver artigo original para lista completa de citações)

Figuras

Figura 1. Modelo DM de corte diamantado de Ce. (a) Vista frontal, (b) vista superior. Vermelho = átomos Ce inferiores, azul = átomos Ce móveis, cinza = átomos C.

Figura 2. (a) RDF de fases do cério. (b) RDF antes e depois de compressão de γ-Ce para colapso de volume 20%.

Figura 3. Variação das forças de corte e normal com comprimento de corte. Subfiguras exibem configurações de defeito por zona.

Figura 4. Estruturas de defeito e morfologia de superfície usinada em vários comprimentos de corte.

Figura 5. Análise RDF de zona de defeito (a) e cavaco (b) após corte.

Figura 6. Dependência do ângulo de inclinação das forças médias de corte e normal.

Figura 7. Estruturas de defeito para ângulos de inclinação −30° (a) e 30° (b); contagem de discordâncias (c).

Figura 8. Morfologias de superfície usinada para ângulos de inclinação −10°, −20°, −30°, 10°, 20°, 30°.

Figura 9. Força de corte e número de discordâncias versus orientação cristalográfica.

Figura 10. Morfologias de superfície usinada para Ce(010), Ce(110), Ce(111).