Comportement Élastique de Cristaux Uniques de Tantalum sous Chargement Tensile <100> : Aperçus Température, Pression et Transition de Phase par Dynamique Moléculaire

Résumé

Les simulations atomistiques peuvent élucider les mécanismes régissant la réponse mécanique du tantale (Ta). Dans cette étude, nous utilisons la dynamique moléculaire (DM) pour examiner comment la température et la pression influencent les propriétés élastiques de cristaux uniques de Ta soumis à un chargement tensile <100>. Nous comparons d'abord deux potentiels de méthode des atomes embarqués (EAM) — Zhou-EAM et Ravelo-EAM — pour leur capacité à reproduire les constantes élastiques du Ta sur une large gamme de pression. Le potentiel Ravelo-EAM démontre une performance supérieure, capturant précisément les modules élastiques jusqu'à 140 GPa. Les simulations DM basées sur ce potentiel révèlent une transition de phase cubique centré (BCC) → cubique face centré (FCC) précédant la fracture à 1 K, indépendante de la taille du modèle ou du taux de déformation. De 1 à 1500 K, le module élastique <100> E₁₀₀ diminue linéairement avec la température, tandis que la contrainte d'écoulement suit une dépendance quadratique en température. Les simulations dépendantes de la pression (0–140 GPa) montrent une augmentation monotone de E₁₀₀. Ces résultats fournissent un cadre atomistique complet pour prédire le comportement élastique du Ta sous conditions extrêmes.

Contexte

Le tantale cristallise en réseau cubique centré (BCC) aux conditions ambiantes et est renommé pour sa stabilité de phase exceptionnelle, point de fusion élevé (3269 K) et résistance à l'oxydation. Ces propriétés font du Ta un candidat idéal pour barrières de diffusion en microélectronique, revêtements résistants à l'usure et superalliages haute température. De nombreux travaux expérimentaux et théoriques — utilisant cellules à enclume diamant, théorie de la fonctionnelle de densité (DFT), et DM — ont exploré le comportement du Ta sous haute pression et température, confirmant la persistance de la phase BCC jusqu'à 135 GPa et révélant des modules élastiques croissants avec la pression. Malgré ces avancées, une investigation atomistique systématique de la réponse dynamique du Ta au chargement tensile fait défaut. Ce travail comble cette lacune en évaluant les effets de taille, taux de déformation, température et pression sur les propriétés élastiques du Ta pendant tension <100> et en sondant la possibilité d'une transition de phase induite par traction.

Méthodes et Détails Expérimentaux

Modélisation Physique

Nous construisons des cellules cubiques Ta en répliquant la cellule unité BCC le long des directions <001>, <010> et <100>, utilisant constantes de réseau a = b = c = 3,301 Å. Quatre tailles de cellules (12 a = 3,96 nm, 18 a = 5,94 nm, 24 a = 7,92 nm, 30 a = 9,90 nm) contiennent 3 456 à 54 000 atomes. La Figure 1 illustre le modèle 3,96 nm, structure principale utilisée dans toutes simulations.

Protocole Simulation DM

Toutes simulations effectuées avec LAMMPS. Deux potentiels EAM — Zhou-EAM et Ravelo-EAM — évalués pour leur capacité à reproduire constantes élastiques à 0 K et sous pression hydrostatique. Le potentiel Ravelo-EAM, développé par Ravelo et al. (Phys. Rev. B 2013, 88, 134101), a surpassé Zhou-EAM sur la gamme de pression d'intérêt. Nous avons utilisé la méthode directe pour calculer C₁₁ et C₁₂ et dérivé le module élastique <100> via Éq. (1). Simulations ultérieures ont employé pas de temps 1 fs, conditions aux limites périodiques toutes directions, et relaxation NVT 50 ps suivie équilibrage NPT pour imposer pressions hydrostatiques (0–140 GPa). Chargement tensile selon x appliqué à taux déformation 5 × 10⁸ s⁻¹ à 7,5 × 10⁹ s⁻¹ jusqu'à 15 % élongation.

Résultats et Discussion

Processus d'Étirement

Figure 3 montre réponse contrainte-déformation et configurations atomiques correspondantes pendant tension <100> à pression nulle. Déformation élastique linéaire jusqu'à ~7 % déformation, après quoi transformation BCC→FCC s'initie à ~7,4 % et s'achève ~9,8 %. Phase FCC persiste jusqu'à fracture de surface près 13 %. Mécanismes transformation et fracture cohérents toutes tailles modèles et taux déformation, indiquant effets négligeables taille ou taux sur module élastique et contrainte écoulement.

Effets Taille et Taux Déformation

Tableau 2 démontre que taille modèle (3,96–9,90 nm) n'influence pas module élastique (~139 GPa) ni contrainte écoulement. Similairement, Tableau 3 montre taux déformation entre 5 × 10⁸ s⁻¹ et 7,5 × 10⁹ s⁻¹ ne produisent changements discernables ces propriétés. Par conséquent, cellule 3,96 nm et taux déformation 5 × 10⁸ s⁻¹ adoptés pour toutes analyses subséquentes.

Dépendance Température

Figure 4 illustre évolution courbes contrainte-déformation 1 K à 1500 K. Module élastique E₁₀₀ et contrainte écoulement diminuent monotoniquestement avec température. Module décline ~42 % (136,5 → 76,7 GPa) et contrainte écoulement ~51 % (≈8 → 4 GPa) sur cette gamme. Ajustement linéaire donne E₁₀₀ = 138,07 – 0,04094 T (GPa), valide jusqu'à 1500 K, tandis qu'ajustement quadratique capture dépendance contrainte écoulement.

Dépendance Pression

Figure 6 présente E₁₀₀ versus pression hydrostatique (0–140 GPa) à plusieurs températures. À 1 K, résultats dynamiques s'accordent étroitement calculs statiques basés Éq. (1). Pour températures ≤ 600 K, région concave apparaît entre 20–60 GPa, cohérente études DM antérieures, tandis que courbes température supérieure restent monotones. À toutes températures, E₁₀₀ augmente avec pression au-dessus ~40 GPa, confirmant que raideur Ta est renforcée sous compression.

Conclusions

Nous avons effectué une investigation DM complète de cristaux uniques Ta sous tension <100> sur large spectre températures (1–1500 K) et pressions (0–140 GPa). Le potentiel Ravelo-EAM s'est avéré fiable pour simulations haute pression. Résultats clés : (1) transition phase BCC→FCC précède fracture à basse température ; (2) module élastique et contrainte écoulement déclinent linéairement et quadratiquement respectivement avec température ; (3) pression augmente module élastique au-dessus ~40 GPa ; (4) taille modèle et taux déformation influence négligeable sur propriétés élastiques mesurées.