Molecular Dynamics Study of Irradiation-Induced Defect and Dislocation Evolution in Strained Nickel

Resumen Técnico: Estudio de Dinámica Molecular de la Evolución de Defectos e Inducción por Irradiación y Deformación en Níquel Tensionado

Planteamiento del Problema
Los materiales en entornos nucleares enfrentan desafíos simultáneos de irradiación de neutrones y tensiones termomecánicas, lo que conduce a la degradación microestructural como el endurecimiento, la fragilización y la hinchazón por huecos (void swelling). Si bien los exámenes post-irradiación y los experimentos de irradiación iónica proporcionan datos macroscópicos, estos son costosos y lentos. Además, la comprensión fundamental de la formación y evolución de los defectos primarios depende en gran medida del modelado computacional. Existe una brecha crítica en la cuantificación de cómo la deformación de tracción elástica altera específicamente el daño por irradiación en metales de estructura cúbica centrada en las caras (fcc). Estudios previos se han centrado mayoritariamente en sistemas sin deformación, dejando insuficientemente caracterizada la interacción entre el estrés externo y la dinámica de defectos —particularmente durante las fases transitorias de las cascadas de colisión—.

Metodología
Este estudio emplea simulaciones de Dinámica Molecular (MD) para investigar el acoplamiento entre la deformación mecánica y la evolución de defectos inducidos por irradiación en monocristales de níquel (Ni) puro.

• Configuración del Sistema: Se simuló un cristal de Ni fcc orientado en [100] que contiene más de 250,000 átomos utilizando el Simulador Atómico de Gran Escala Masivamente Paralelo (LAMMPS). Las interacciones interatómicas se modelaron mediante potenciales de Método de Átomo Embebido Modificado (MEAM), unidos suavemente al potencial repulsivo universal de Ziegler–Biersack–Littmark (ZBL) para describir con precisión las colisiones de alta energía.
• Protocolo de Irradiación: El daño por irradiación se simuló iniciando cascadas de desplazamiento acumulativas y superpuestas utilizando Átomos Primarios de Choque (PKA) de 5 keV a 300 K. Cada cascada se permitió evolucionar durante 20 ps, una duración validada para capturar completamente la fase balística, la formación del pico de calor (heat-spike) y la recuperación de defectos.
• Aplicación de la Deformación: Se aplicó deformación de tracción a lo largo de la dirección [001] a una tasa de $10^8 s^{-1}$ para alcanzar deformaciones objetivo dentro del régimen elástico (por ejemplo, 2%). La deformación se aplicó para establecer un campo de deformación mecánicamente equilibrado antes de cada evento de irradiación, asegurando que la cascada evolucionara dentro de una red pre-tensionada estática.
• Rigor Estadístico: Para abordar la estocasticidad, se ejecutó una secuencia de diez cascadas superpuestas. Esta secuencia se repitió 10 veces con proyectiles PKA seleccionados al azar, resultando en 100 cascadas por condición. Las poblaciones de defectos y las densidades de dislocación se reportaron como promedios de conjunto.
• Herramientas de Análisis: La identificación de defectos utilizó el análisis de Wigner–Seitz para pares de Frenkel y el Algoritmo de Extracción de Dislocaciones (DXA) mediante OVITO para la caracterización de dislocaciones. La evolución de la densidad de dislocaciones se analizó utilizando el modelo de Kocks–Mecking.

Resultos Clave

1. Producción y Recuperación de Defectos: Cada cascada de 5 keV produjo un promedio de 54 pares de Frenkel, consistente con la estimación NRT. En el Ni sin deformación, los procesos de recuperación de la fase balística y del pico de calor se completaron en $\lesssim 10$ ps. La deformación de tracción aumentó sistemáticamente la población de defectos retenidos en comparación con el caso sin deformación.
2. Dinámica del Pico de Calor (Heat-Spike): La duración y el desorden máximo del pico de calor disminuyeron monótnicamente con el aumento de la deformación. Esto indica una disipación de energía local más rápida bajo tensión. Sin embargo, a pesar de la disipación más rápida, la población de defectos supervivientes fue mayor en los sistemas deformados. Esto se atribuye a la movilidad de defectos asistida por estrés, donde la tensión de tracción sesga la migración de intersticiales y estabiliza configuraciones ricas en intersticiales, reduciendo la probabilidad de recombinación vacío-intersticial.
3. Evolución de Dislocaciones: El análisis de dislocaciones reveló un fuerte sesgo hacia bucles de dislocación parcial de Shockley. Los componentes de dislocación de Hirth, stair-rod, Frank y perfectos permanecieron como secundarios.
   • En el Ni sin deformación, las parciales de Shockley dominaron, pero su densidad era menor.
   • En el Ni deformado, la formación de dislocaciones de Shockley se incrementó notablemente, mientras que las dislocaciones de tipo Hirth disminuyeron aproximadamente un orden de magnitud.
   • Tanto la densidad total como la de dislocaciones de Shockley convergieron a una saturación de estado estacionario de aproximadamente $1 \times 10^{16} m^{-2}$ (normalizado a $\sim 0.1$ en unidades específicas) conforme aumentó la deformación.
4. Ajuste del Modelo: La evolución de la densidad de dislocaciones totales dependiente de la deformación fue capturada exitosamente por el modelo de Kocks–Mecking, que describe el equilibrio entre el almacenamiento de dislocaciones y la recuperación dinámica. Los parámetros ajustados confirmaron un comportamiento de saturación en estado estacionario consistente con los datos de la simulación.

Significancia y Reivindicaciones
Los autores afirman que este trabajo proporciona el primer vínculo cuantitativo entre los mecanismos de irradiación a escala atómica y la estabilización de la microestructura dependiente de la deformación en sistemas metálicos relevantes para aplicaciones nucleares. El estudio demuestra un robusto acoplamiento de irradiación-deformación en metales fcc a dosis bajas, revelando que la carga mecánica puede sesgar significativamente el estado de daño primario incluso a dosis extremadamente bajas (aprox. 0.002 dpa).

Las contribuciones clave incluyen:

• Validar que la deformación de tracción altera la morfología y la vida útil de los picos de calor, promoviendo una movilidad de defectos asistida por estrés que favorece la retención de intersticiales sobre la recombinación.
• Establecer que la respuesta microestructural está dominada por la nucleación y estabilización de bucles de dislocación parcial de Shockley, una tendencia que se intensifica bajo deformación de tracción.
• Demostrar que el modelo de Kocks–Mecking puede describir eficazmente la evolución de las dislocaciones bajo irradiación y deformación simultáneas, proporcionando una base predictiva para comprender el acoplamiento mecánico-radiación.

El artículo concluye que estos mecanismos —el acoplamiento elástico de los defectos puntuales a los campos de estrés y la estabilización preferencial de cúmulos ricos en intersticiales— son características fundamentales de la cristalografía fcc. Por consiguiente, aunque la retención cuantitativa puede variar, las tendencias cualitativas de mayor supervivencia de intersticiales y la prevalencia de bucles de Shockley se espera que persistan en aleaciones basadas en Ni y otros sistemas fcc expuestos a carga mecánica e irradiación concurrentes.