Elastisch gedrag van tantalum einkristallen onder <100> trekkingsbelasting: temperatuur, druk en fasetransitie inzichten uit moleculaire dynamica

Samenvatting

Atomistische simulaties kunnen de mechanismen verhelderen die de mechanische respons van tantalum (Ta) regeren. In deze studie gebruiken we moleculaire dynamica (MD) om te onderzoeken hoe temperatuur en druk de elastische eigenschappen van Ta einkristallen beïnvloeden die onderworpen zijn aan <100> trekkingsbelasting. We vergelijken eerst twee embedded-atom method (EAM) potentiaal – Zhou-EAM en Ravelo-EAM – voor hun vermogen om de elastische constanten van Ta over een brede drukbereik te reproduceren. De Ravelo-EAM potentiaal toont superieure prestaties, met nauwkeurige vastlegging van elastische moduli tot 140 GPa. MD-simulaties gebaseerd op deze potentiaal onthullen een lichaam-gecentreerd-kubisch (BCC) → vlak-gecentreerd-kubisch (FCC) fasetransitie voorafgaand aan breuk bij 1 K, onafhankelijk van modelgrootte of vervormingssnelheid. Over 1–1500 K daalt de <100> elastische modulus E₁₀₀ lineair met temperatuur, terwijl de vloeispanning een kwadratische temperatuurafhankelijkheid volgt. Druk-afhankelijke simulaties (0–140 GPa) tonen een monotone toename in E₁₀₀. Deze resultaten bieden een uitgebreid, atomistisch raamwerk voor het voorspellen van Ta's elastisch gedrag onder extreme omstandigheden.

Achtergrond

Tantal kristalliseert in een lichaam-gecentreerd-kubisch (BCC) rooster bij omstandigheden en is bekend om zijn uitzonderlijke fasestabiliteit, hoge smeltpunt (3269 K) en oxidatiewerstand. Deze eigenschappen maken Ta een ideale kandidaat voor diffusiebarrières in micro-elektronica, slijtvast coatingen en hoge-temperatuur super-alliages. Uitgebreid experimenteel en theoretisch werk – met diamant-stempelcellen, dichtheidsfunctionaaltheorie (DFT) en MD – heeft Ta's gedrag onder hoge druk en temperatuur verkend, consistent bevestigend de persistentie van de BCC-fase tot 135 GPa en onthullend toenemende elastische moduli met druk. Ondanks deze vooruitgang, ontbreekt een systematisch atomistisch onderzoek van Ta's dynamische respons op trekkingsbelasting. Dit werk vult die leemte door grootte, vervormingssnelheid, temperatuur- en drukeffecten op Ta's elastische eigenschappen tijdens <100> trekking te evalueren en de mogelijkheid van een trekking-geïnduceerde fasetransitie te onderzoeken.

Methoden en Experimentele Details

Fysische Modellering

We construeren kubische Ta-cellen door de BCC-eenheidscel te repliceren langs de <001>, <010> en <100> richtingen, met roosterconstanten a = b = c = 3,301 Å. Vier celgroottes (12 a = 3,96 nm, 18 a = 5,94 nm, 24 a = 7,92 nm, 30 a = 9,90 nm) bevatten 3.456 tot 54.000 atomen. Figuur 1 illustreert het 3,96 nm model, de primaire structuur gebruikt in alle simulaties.

MD-simulatieprotocol

Alle simulaties werden uitgevoerd met LAMMPS. Twee EAM-potentialen – Zhou-EAM en Ravelo-EAM – werden geëvalueerd voor hun vermogen om elastische constanten bij 0 K en onder hydrostatische druk te reproduceren. De Ravelo-EAM potentiaal, ontwikkeld door Ravelo et al. (Phys. Rev. B 2013, 88, 134101), presteerde beter dan Zhou-EAM over het drukbereik van belang. We gebruikten de directe methode om C₁₁ en C₁₂ te berekenen en leidden de <100> elastische modulus af via Vergelijking (1). Volgende simulaties gebruikten een 1 fs tijdstap, periodieke randvoorwaarden in alle richtingen, en een 50 ps NVT-ontspanning gevolgd door NPT-vergelijking om hydrostatische drukken (0–140 GPa) op te leggen. Trekkingsbelasting langs x werd toegepast bij vervormingssnelheden van 5 × 10⁸ s⁻¹ tot 7,5 × 10⁹ s⁻¹ tot een 15 % rek bereikt was.

Resultaten en Discussie

Rekproces

Figuur 3 toont de spanning-rek-respons en bijbehorende atomaire configuraties tijdens <100> trekking bij nul druk. Elastische vervorming is lineair tot ~7 % rek, waarna een BCC→FCC-transformatie initieert bij ~7,4 % en voltooid is bij ~9,8 %. De FCC-fase persisteert tot een oppervlaktebreuk optreedt bij ~13 %. De transformatie- en breukmechanismen zijn consistent over alle modelgroottes en vervormingssnelheden, wat verwaarloosbare grootte- of snelheidseffecten op de elastische modulus en vloeispanning aangeeft.

Grootte- en Vervormingssnelheidseffecten

Tabel 2 demonstreert dat modelgrootte (3,96–9,90 nm) de elastische modulus (~139 GPa) of vloeispanning niet beïnvloedt. Op dezelfde manier toont Tabel 3 dat vervormingssnelheden tussen 5 × 10⁸ s⁻¹ en 7,5 × 10⁹ s⁻¹ geen waarneembare veranderingen in deze eigenschappen produceren. Consequent werd een 3,96 nm cel en een vervormingssnelheid van 5 × 10⁸ s⁻¹ gekozen voor alle volgende analyses.

Temperatuurafhankelijkheid

Figuur 4 illustreert de evolutie van de spanning-rek-curves van 1 K tot 1500 K. Zowel de elastische modulus E₁₀₀ als de vloeispanning dalen monotoon met temperatuur. De modulus daalt met ~42 % (136,5 → 76,7 GPa) en de vloeispanning met ~51 % (≈8 → 4 GPa) over dit bereik. Een lineaire fit levert E₁₀₀ = 138,07 – 0,04094 T (GPa) op, geldig tot 1500 K, terwijl een kwadratische fit de vloeispanningsafhankelijkheid vastlegt.

Drukafhankelijkheid

Figuur 6 presenteert E₁₀₀ versus hydrostatische druk (0–140 GPa) bij verschillende temperaturen. Bij 1 K stemmen de dynamische resultaten nauwkeurig overeen met statische berekeningen gebaseerd op Verg. (1). Voor temperaturen ≤ 600 K verschijnt een concave regio tussen 20–60 GPa, consistent met eerdere MD-studies, terwijl hogere-temperatuurcurves monotoon blijven. Over alle temperaturen neemt E₁₀₀ toe met druk boven ~40 GPa, wat bevestigt dat Ta's stijfheid verbetert onder compressie.

Conclusies

We hebben een uitgebreid MD-onderzoek uitgevoerd naar Ta einkristallen onder <100> trekking over een breed spectrum van temperaturen (1–1500 K) en drkken (0–140 GPa). De Ravelo-EAM potentiaal bleek betrouwbaar voor hoge-druk simulaties. Belangrijkste bevindingen omvatten: (1) een BCC→FCC fasetransitie gaat breuk vooraf bij lage temperatuur; (2) elastische modulus en vloeispanning dalen respectievelijk lineair en kwadratisch met temperatuur; (3) druk verhoogt de elastische modulus boven ~40 GPa; en (4) modelgrootte en vervormingssnelheid oefenen verwaarloosbare invloed uit op de gemeten elastische eigenschappen.