Elastisches Verhalten von Tantal-Einkristallen unter <100>-Zugbelastung: Temperatur-, Druck- und Phasenübergangs-Einblicke aus Molekulardynamik

Zusammenfassung

Atomistische Simulationen können die Mechanismen aufklären, die das mechanische Verhalten von Tantal (Ta) bestimmen. In dieser Studie verwenden wir Molekulardynamik (MD), um zu untersuchen, wie Temperatur und Druck die elastischen Eigenschaften von Ta-Einkristallen unter <100>-Zugbelastung beeinflussen. Zunächst vergleichen wir zwei eingebettete-Atom-Methode (EAM) Potenziale – Zhou-EAM und Ravelo-EAM – hinsichtlich ihrer Fähigkeit, die elastischen Konstanten von Ta über einen weiten Druckbereich wiederzugeben. Das Ravelo-EAM-Potenzial zeigt überlegene Leistung und erfasst die Elastizitätsmoduln präzise bis 140 GPa. MD-Simulationen basierend auf diesem Potenzial offenbaren einen körperzentriert-kubischen (BCC) → flächenzentriert-kubischen (FCC) Phasenübergang vor dem Bruch bei 1 K, unabhängig von Modellgröße oder Dehnrate. Im Bereich 1–1500 K nimmt der <100> Elastizitätsmodul E₁₀₀ linear mit der Temperatur ab, während die Fließspannung einer quadratischen Temperaturabhängigkeit folgt. Druckabhängige Simulationen (0–140 GPa) zeigen einen monotonen Anstieg von E₁₀₀. Diese Ergebnisse liefern einen umfassenden, atomistischen Rahmen zur Vorhersage von Tas elastischem Verhalten unter extremen Bedingungen.

Hintergrund

Tantal kristallisiert unter Umgebungsbedingungen in einem körperzentriert-kubischen (BCC) Gitter und ist bekannt für seine außergewöhnliche Phasenstabilität, hohen Schmelzpunkt (3269 K) und Oxidationsbeständigkeit. Diese Eigenschaften machen Ta zu einem idealen Kandidaten für Diffusionsbarrieren in der Mikroelektronik, verschleißfeste Beschichtungen und Hochtemperatur-Superlegierungen. Umfassende experimentelle und theoretische Arbeiten – mit Diamantstempelzellen, Dichtefunktionaltheorie (DFT) und MD – haben Tas Verhalten unter hohem Druck und Temperatur erforscht und konsistent die Beständigkeit der BCC-Phase bis 135 GPa bestätigt sowie steigende Elastizitätsmoduln mit Druck aufgezeigt. Trotz dieser Fortschritte fehlt eine systematische atomistische Untersuchung von Tas dynamischer Antwort auf Zugbelastung. Diese Arbeit schließt diese Lücke, indem sie Größen-, Dehnrate-, Temperatur- und Druckeffekte auf Tas elastische Eigenschaften während <100>-Zug bewertet und die Möglichkeit eines zuginduzierten Phasenübergangs untersucht.

Methoden und Experimentelle Details

Physikalisches Modell

Wir konstruieren kubische Ta-Zellen durch Replikation der BCC-Einheitszelle entlang der <001>, <010> und <100> Richtungen, mit Gitterkonstanten a = b = c = 3,301 Å. Vier Zellengrößen (12 a = 3,96 nm, 18 a = 5,94 nm, 24 a = 7,92 nm, 30 a = 9,90 nm) enthalten 3456 bis 54.000 Atome. Abbildung 1 illustriert das 3,96 nm Modell, das die primäre Struktur in allen Simulationen ist.

MD-Simulationsprotokoll

Alle Simulationen wurden mit LAMMPS durchgeführt. Zwei EAM-Potenziale – Zhou-EAM und Ravelo-EAM – wurden auf ihre Fähigkeit bewertet, elastische Konstanten bei 0 K und unter hydrostatischem Druck wiederzugeben. Das Ravelo-EAM-Potenzial, entwickelt von Ravelo et al. (Phys. Rev. B 2013, 88, 134101), übertraf Zhou-EAM über den interessierenden Druckbereich. Wir nutzten die direkte Methode zur Berechnung von C₁₁ und C₁₂ und leiteten den <100> Elastizitätsmodul über Gl. (1) ab. Nachfolgende Simulationen verwendeten 1 fs Zeitschritt, periodische Randbedingungen in alle Richtungen und 50 ps NVT-Relaxation gefolgt von NPT-Gleichgewicht, um hydrostatische Drücke (0–140 GPa) aufzuerlegen. Zugbelastung entlang x wurde bei Dehnraten von 5 × 10⁸ s⁻¹ bis 7,5 × 10⁹ s⁻¹ bis zu 15 % Dehnung appliziert.

Ergebnisse und Diskussion

Dehnungsprozess

Abbildung 3 zeigt die Spannungs-Dehnungs-Antwort und entsprechende Atomkonfigurationen während <100>-Zug bei null Druck. Elastische Deformation ist linear bis ~7 % Dehnung, wonach eine BCC→FCC-Transformation bei ~7,4 % einsetzt und bis ~9,8 % abgeschlossen ist. Die FCC-Phase persistiert bis zum Oberflächenbruch nahe 13 %. Transformations- und Bruchmechanismen sind über alle Modellgrößen und Dehnraten konsistent, was vernachlässigbare Größen- oder Rateeffekte auf Elastizitätsmodul und Fließspannung anzeigt.

Größen- und Dehnrate-Effekte

Tabelle 2 zeigt, dass Modellgröße (3,96–9,90 nm) weder Elastizitätsmodul (~139 GPa) noch Fließspannung beeinflusst. Ebenso zeigt Tabelle 3, dass Dehnraten zwischen 5 × 10⁸ s⁻¹ und 7,5 × 10⁹ s⁻¹ keine erkennbaren Änderungen in diesen Eigenschaften bewirken. Folglich wurden eine 3,96 nm Zelle und eine Dehnrate von 5 × 10⁸ s⁻¹ für alle nachfolgenden Analysen übernommen.

Temperaturabhängigkeit

Abbildung 4 illustriert die Entwicklung der Spannungs-Dehnungs-Kurven von 1 K bis 1500 K. Sowohl Elastizitätsmodul E₁₀₀ als auch Fließspannung nehmen monoton mit Temperatur ab. Der Modul sinkt um ~42 % (136,5 → 76,7 GPa) und die Fließspannung um ~51 % (≈8 → 4 GPa) über diesen Bereich. Ein linearer Fit ergibt E₁₀₀ = 138,07 – 0,04094 T (GPa), gültig bis 1500 K, während ein quadratischer Fit die Fließspannungsabhängigkeit erfasst.

Druckabhängigkeit

Abbildung 6 stellt E₁₀₀ versus hydrostatischem Druck (0–140 GPa) bei verschiedenen Temperaturen dar. Bei 1 K stimmen die dynamischen Ergebnisse eng mit statischen Berechnungen basierend auf Gl. (1) überein. Für Temperaturen ≤ 600 K erscheint ein konkaver Bereich zwischen 20–60 GPa, konsistent mit früheren MD-Studien, während Kurven bei höheren Temperaturen monoton bleiben. Über alle Temperaturen hinweg steigt E₁₀₀ mit Druck über ~40 GPa, was bestätigt, dass Tas Steifigkeit unter Kompression zunimmt.

Schlussfolgerungen

Wir führten eine umfassende MD-Untersuchung von Ta-Einkristallen unter <100>-Zug über ein breites Spektrum von Temperaturen (1–1500 K) und Drücken (0–140 GPa) durch. Das Ravelo-EAM-Potenzial erwies sich als zuverlässig für Hochdrucksimulationen. Wesentliche Erkenntnisse umfassen: (1) ein BCC→FCC-Phasenübergang geht dem Bruch bei niedriger Temperatur voraus; (2) Elastizitätsmodul und Fließspannung nehmen linear bzw. quadratisch mit Temperatur ab; (3) Druck erhöht den Elastizitätsmodul über ~40 GPa; und (4) Modellgröße und Dehnrate üben vernachlässigbaren Einfluss auf die gemessenen elastischen Eigenschaften aus.