Comportamiento Elástico de Cristales Simples Tantalio Bajo Carga Tracción <100>: Perspectivas de Temperatura, Presión y Transición de Fase desde Dinámica Molecular

Resumen

Las simulaciones atomísticas pueden elucidar mecanismos que gobiernan respuesta mecánica del tantalio (Ta). En este estudio, usamos dinámica molecular (MD) para examinar cómo temperatura y presión influyen propiedades elásticas de cristales simples Ta sometidos a carga tracción <100>. Primero comparamos dos potenciales método átomo embebido (EAM) —Zhou-EAM y Ravelo-EAM— por su capacidad reproducir constantes elásticas Ta en amplio rango presión. El potencial Ravelo-EAM demuestra rendimiento superior, capturando accuradamente módulos elásticos hasta 140 GPa. Simulaciones MD basadas en este potencial revelan transición fase cuerpo-céntrico-cúbico (BCC) → cara-céntrico-cúbico (FCC) precediendo fractura a 1 K, independiente de tamaño modelo o tasa deformación. En 1–1500 K, módulo elástico <100> E₁₀₀ disminuye linealmente con temperatura, mientras tensión cedencia sigue dependencia cuadrática temperatura. Simulaciones dependientes presión (0–140 GPa) muestran aumento monótono E₁₀₀. Estos resultados proveen marco atomístico integral para predecir comportamiento elástico Ta bajo condiciones extremas.

Antecedentes

El tantalio cristaliza en red cuerpo-céntrico-cúbico (BCC) a condiciones ambiente y es renombrado por su excepcional estabilidad fase, alto punto fusión (3269 K) y resistencia oxidación. Estas propiedades hacen Ta candidato ideal para barreras difusión en microelectrónica, recubrimientos resistentes desgaste y superaleaciones alta temperatura. Extenso trabajo experimental y teórico —usando celdas yunque diamante, teoría funcional densidad (DFT) y MD— ha explorado comportamiento Ta bajo alta presión y temperatura, consistentemente confirmando persistencia fase BCC hasta 135 GPa y revelando aumento módulos elásticos con presión. A pesar estos avances, investigación atomística sistemática de respuesta dinámica Ta a carga tracción permanece ausente. Este trabajo aborda esa brecha evaluando efectos tamaño, tasa deformación, temperatura y presión sobre propiedades elásticas Ta durante tensión <100> y sondeando posibilidad transición fase inducida tracción.

Métodos y Detalles Experimentales

Modelado Físico

Construimos celdas cúbicas Ta replicando celda unitaria BCC entlang direcciones <001>, <010> y <100>, usando constantes red a = b = c = 3.301 Å. Cuatro tamaños celda (12 a = 3.96 nm, 18 a = 5.94 nm, 24 a = 7.92 nm, 30 a = 9.90 nm) contienen 3456 a 54000 átomos. Figura 1 ilustra modelo 3.96 nm, estructura primaria usada en todas simulaciones.

Protocolo Simulación MD

Todas simulaciones se realizaron con LAMMPS. Dos potenciales EAM —Zhou-EAM y Ravelo-EAM— se evaluaron por su capacidad reproducir constantes elásticas a 0 K y bajo presión hidrostática. El potencial Ravelo-EAM, desarrollado por Ravelo et al. (Phys. Rev. B 2013, 88, 134101), superó Zhou-EAM en rango presión interés. Usamos método directo para computar C₁₁ y C₁₂ y derivamos módulo elástico <100> vía Ec. (1). Simulaciones subsiguientes emplearon paso temporal 1 fs, condiciones frontera periódicas todas direcciones, y relajación NVT 50 ps seguida equilibración NPT para imponer presiones hidrostáticas (0–140 GPa). Carga tracción lungo x se aplicó a tasas deformación 5 × 10⁸ s⁻¹ a 7.5 × 10⁹ s⁻¹ hasta 15% elongación alcanzada.

Resultados y Discusión

Proceso Estiramiento

Figura 3 muestra respuesta tensión-deformación y configuraciones atómicas correspondientes durante tensión <100> a presión cero. Deformación elástica es lineal hasta ~7% deformación, tras cual transformación BCC→FCC inicia ~7.4% y completa ~9.8%. Fase FCC persiste hasta fractura superficial cerca 13%. Mecanismos transformación y fractura son consistentes en todos tamaños modelo y tasas deformación, indicando efectos despreciables tamaño o tasa sobre módulo elástico y tensión cedencia.

Efectos Tamaño y Tasa Deformación

Tabla 2 demuestra que tamaño modelo (3.96–9.90 nm) no influye módulo elástico (~139 GPa) o tensión cedencia. Similarmente, Tabla 3 muestra tasas deformación entre 5 × 10⁸ s⁻¹ y 7.5 × 10⁹ s⁻¹ no producen cambios discernibles en estas propiedades. Consecuentemente, celda 3.96 nm y tasa deformación 5 × 10⁸ s⁻¹ se adoptaron para todos análisis subsiguientes.

Dependencia Temperatura

Figura 4 ilustra evolución curvas tensión-deformación 1 K a 1500 K. Tanto módulo elástico E₁₀₀ como tensión cedencia disminuyen monótonamente con temperatura. El módulo declina ~42% (136.5 → 76.7 GPa) y tensión cedencia ~51% (≈8 → 4 GPa) en este rango. Ajuste lineal produce E₁₀₀ = 138.07 – 0.04094 T (GPa), válido hasta 1500 K, mientras ajuste cuadrático captura dependencia tensión cedencia temperatura.

Dependencia Presión

Figura 6 presenta E₁₀₀ versus presión hidrostática (0–140 GPa) a varias temperaturas. A 1 K, resultados dinámicos concuerdan estrechamente con cálculos estáticos basados Ec. (1). Para temperaturas ≤600 K, región cóncava aparece entre 20–60 GPa, consistente con estudios MD previos, mientras curvas mayor temperatura permanecen monótonas. En todas temperaturas, E₁₀₀ aumenta con presión sobre ~40 GPa, confirmando rigidez Ta se realza bajo compresión.

Conclusiones

Realizamos investigación MD integral de cristales simples Ta bajo tensión <100> en amplio espectro temperaturas (1–1500 K) y presiones (0–140 GPa). El potencial Ravelo-EAM probó confiable para simulaciones alta presión. Hallazgos clave incluyen: (1) transición fase BCC→FCC precede fractura a baja temperatura; (2) módulo elástico y tensión cedencia declinan lineal y cuadráticamente, respectivamente, con temperatura; (3) presión realza módulo elástico sobre ~40 GPa; y (4) tamaño modelo y tasa deformación ejercen influencia despreciable en propiedades elásticas medidas.