Comportamento Elástico de Cristais Únicos de Tântalo sob Carga de Tração <100>: Insights de Temperatura, Pressão e Transição de Fase da Dinâmica Molecular

Resumo

Simulações atomísticas podem elucidar os mecanismos que governam a resposta mecânica do tântalo (Ta). Neste estudo, usamos dinâmica molecular (DM) para examinar como temperatura e pressão influenciam as propriedades elásticas de cristais únicos de Ta submetidos a carga de tração <100>. Primeiro comparamos dois potenciais de método de átomo embutido (EAM) — Zhou-EAM e Ravelo-EAM — por sua capacidade de reproduzir constantes elásticas de Ta em ampla faixa de pressão. O potencial Ravelo-EAM demonstra desempenho superior, capturando com precisão módulos elásticos até 140 GPa. Simulações DM baseadas neste potencial revelam transição de fase corpo-centrado-cúbico (BCC) → face-centrado-cúbico (FCC) precedendo fratura a 1 K, independente de tamanho do modelo ou taxa de deformação. Em 1–1500 K, o módulo elástico <100> E₁₀₀ diminui linearmente com a temperatura, enquanto a tensão de escoamento segue dependência quadrática de temperatura. Simulações dependentes de pressão (0–140 GPa) mostram aumento monótono em E₁₀₀. Estes resultados fornecem estrutura atomística abrangente para prever comportamento elástico do Ta sob condições extremas.

Contexto

Tântalo cristaliza em rede corpo-centrado-cúbica (BCC) em condições ambientes e é renomado por sua estabilidade de fase excepcional, alto ponto de fusão (3269 K) e resistência à oxidação. Estas propriedades tornam Ta candidato ideal para barreiras de difusão em microeletrônica, revestimentos resistentes ao desgaste e superligas de alta temperatura. Trabalhos experimentais e teóricos extensivos — usando células de bigorna de diamante, teoria do funcional da densidade (DFT) e DM — exploraram comportamento do Ta sob alta pressão e temperatura, consistentemente confirmando persistência da fase BCC até 135 GPa e revelando aumento de módulos elásticos com pressão. Apesar destes avanços, investigação atomística sistemática da resposta dinâmica do Ta a carga de tração permanece ausente. Este trabalho aborda essa lacuna avaliando efeitos de tamanho, taxa de deformação, temperatura e pressão nas propriedades elásticas do Ta durante tensão <100> e sondando possibilidade de transição de fase induzida por tração.

Métodos e Detalhes Experimentais

Modelagem Física

Construímos células cúbicas de Ta replicando célula unitária BCC ao longo das direções <001>, <010> e <100>, usando constantes de rede a = b = c = 3,301 Å. Quatro tamanhos de célula (12 a = 3,96 nm, 18 a = 5,94 nm, 24 a = 7,92 nm, 30 a = 9,90 nm) contêm 3.456 a 54.000 átomos. Figura 1 ilustra modelo cúbico de 3,96 nm, que é a estrutura primária usada em todas as simulações.

Protocolo de Simulação DM

Todas as simulações foram realizadas com LAMMPS. Dois potenciais EAM — Zhou-EAM e Ravelo-EAM — foram avaliados por sua capacidade de reproduzir constantes elásticas a 0 K e sob pressão hidrostática. O potencial Ravelo-EAM, desenvolvido por Ravelo et al. (Phys. Rev. B 2013, 88, 134101), superou Zhou-EAM na faixa de pressão de interesse. Usamos método direto para calcular C₁₁ e C₁₂ e derivamos módulo elástico <100> via Eq. (1). Simulações subsequentes empregaram passo de tempo de 1 fs, condições de contorno periódicas em todas direções, e relaxamento NVT de 50 ps seguido por equilibração NPT para impor pressões hidrostáticas (0–140 GPa). Carga de tração ao longo de x foi aplicada em taxas de deformação de 5 × 10⁸ s⁻¹ a 7,5 × 10⁹ s⁻¹ até alongamento de 15% ser alcançado.

Resultados e Discussão

Processo de Estiramento

Figura 3 mostra resposta tensão-deformação e configurações atômicas correspondentes durante tensão <100> a pressão zero. Deformação elástica é linear até ~7% deformação, após o que transformação BCC→FCC inicia a ~7,4% e completa por ~9,8%. Fase FCC persiste até fratura de superfície ocorrer próximo a 13%. Mecanismos de transformação e fratura são consistentes em todos tamanhos de modelo e taxas de deformação, indicando efeitos insignificantes de tamanho ou taxa no módulo elástico e tensão de escoamento.

Efeitos de Tamanho e Taxa de Deformação

Tabela 2 demonstra que tamanho do modelo (3,96–9,90 nm) não influencia módulo elástico (~139 GPa) ou tensão de escoamento. Similarmente, Tabela 3 mostra que taxas de deformação entre 5 × 10⁸ s⁻¹ e 7,5 × 10⁹ s⁻¹ não produzem mudanças discerníveis nestas propriedades. Consequentemente, célula de 3,96 nm e taxa de deformação de 5 × 10⁸ s⁻¹ foram adotadas para todas análises subsequentes.

Dependência de Temperatura

Figura 4 ilustra evolução das curvas tensão-deformação de 1 K a 1500 K. Tanto módulo elástico E₁₀₀ quanto tensão de escoamento diminuem monoticamente com temperatura. Módulo declina ~42% (136,5 → 76,7 GPa) e tensão de escoamento ~51% (≈8 → 4 GPa) nesta faixa. Ajuste linear resulta em E₁₀₀ = 138,07 – 0,04094 T (GPa), válido até 1500 K, enquanto ajuste quadrático captura dependência da tensão de escoamento.

Dependência de Pressão

Figura 6 apresenta E₁₀₀ versus pressão hidrostática (0–140 GPa) em várias temperaturas. A 1 K, resultados dinâmicos concordam estreitamente com cálculos estáticos baseados na Eq. (1). Para temperaturas ≤ 600 K, região côncava aparece entre 20–60 GPa, consistente com estudos DM prévios, enquanto curvas de temperatura mais alta permanecem monótonas. Em todas as temperaturas, E₁₀₀ aumenta com pressão acima de ~40 GPa, confirmando que rigidez do Ta é aumentada sob compressão.

Conclusões

Realizamos investigação DM abrangente de cristais únicos de Ta sob tensão <100> em amplo espectro de temperaturas (1–1500 K) e pressões (0–140 GPa). O potencial Ravelo-EAM provou confiável para simulações de alta pressão. Descobertas-chave incluem: (1) transição de fase BCC→FCC precede fratura a baixa temperatura; (2) módulo elástico e tensão de escoamento declinam linear e quadraticamente, respectivamente, com temperatura; (3) pressão aumenta módulo elástico acima de ~40 GPa; e (4) tamanho do modelo e taxa de deformação exercem influência insignificante nas propriedades elásticas medidas.