Deformationsbeteende och skademekanism för 3C-SiC under nanoindentation +

PAM-XIAMEN levererar 3C-SiC-skivor för forskningsändamål. För mer information, vänligen kontakta:[email protected]

Kubisk kiselkarbid (3C-SiC) uppvisar stor potential inom banbrytande områden som mikroelektromekaniska system, kvantberäkningar och optoelektroniska enheter, på grund av dess breda bandgap, höga värmeledningsförmåga, höga kritiska nedbrytningsfältstyrka och utmärkta kemiska tröghet. Den ultrahöga hårdheten och den inneboende sprödheten hos 3C-SiC gör dock dess ultraprecisionsbearbetning extremt utmanande, och bearbetningsinducerade skador på ytan/underytan har blivit en viktig flaskhals som begränsar dess avancerade applikationer. Därför är en djup förståelse av deformationsbeteendet och skademekanismerna hos 3C-SiC under kontaktbelastning av stor teoretisk och teknisk betydelse för att optimera dess bearbetningsteknik.

Under de senaste åren har forskare i stor utsträckning utforskat det mekaniska svaret av 3C-SiC genom att kombinera experiment med nanoindentation med simuleringar av molekylär dynamik. Zhu et al. undersökte systematiskt nanoindentationsbeteendet hos 3C-SiC på (010), (110) och (111) kristallplanen med hjälp av simuleringar av molekylär dynamik (MD), avslöjade de anisotropa egenskaperna hos elastisk deformation och dislokationsevolution, och introducerade Thompsons tetrahedrons tetraedermekanism för att förklara den plastiska deformationsmekanismen. Samtidigt påverkar temperaturen, som en nyckelfaktor i faktiska servicemiljöer, avsevärt de mekaniska egenskaperna hos 3C-SiC. Kourbani et al. funnit genom MD-simuleringar att ökande temperatur minskar hårdheten och Youngs modul för materialet och främjar amorfisering. Å andra sidan har traditionella kontinuummekaniska metoder begränsningar när det gäller att hantera diskontinuerliga problem såsom sprickinitiering och fortplantning. För att ta itu med detta använde Yang och Zhu peridynamik för att simulera skadeutvecklingen av 3C-SiC under indrag och repning, och framgångsrikt reproducera pop-in-fenomenet och bildandet av radiella sprickor. Dessa studier ger referenser i flera skalor för att förstå det mekaniska beteendet och bearbetningsteorin för 3C-SiC.

1. Anisotropisk elastisk deformation och tidigt plastiskt beteende hos 3C-SiC

3C-SiC har en zinkblandningsstruktur och dess mekaniska respons är starkt beroende av kristallorienteringen. Genom MD-simuleringar, Zhu et al. fann att under nanointryckning på (010), (110) och (111) planen, uppvisar de elastiska deformationsområdena distinkt olika symmetrier: (010) planet uppvisar fyrfaldig symmetri, med en pyramidformad deformationszon som expanderar längs <110>-riktningen; planet (110) uppvisar central symmetri, med deformationszonen delad i två delar; planet (111) visar trefaldig symmetri, med tre elastiska distorsionszoner som också expanderar längs <110>-riktningen. Detta resultat indikerar att <110>-riktningen är den optimala vägen för elastisk spänningsöverföring.

Fig. 1 Morfologi för den elastiska distorsionszonen av 3C-SiC vid ett intryckningsdjup av 1,1 nm: (a,d) (010) plan; (b,e) (110) plan; (c,f) (111) plan; röda pilar indikerar riktningen för den elastiska deformationens utbredning

Från belastningsdjupkurvorna varierar längden på det elastiska steget avsevärt mellan olika kristallplan. Ph-kurvan för (010)-planet matchar Hertz-kontaktmodellen mycket väl för h<2.0nm, medan de elastiska stegen för (110) och (111) planen är betydligt kortare, vilket indikerar tidigare plastisk eftergivenhet. Noterbart visas pop-in-händelser i kurvorna för (010)- och (110)-planen, men observeras inte för (111)-planet. Genom atomförskjutningsanalys, Zhu et al. bekräftade att pop-in-händelsen motsvarar en riktningsrörelse av ett stort antal atomer längs <110>-riktningen, vilket leder till en plötslig frisättning av lagrad elastisk energi. Till exempel, när fördjupningsdjupet på (110)-planet ökar från 2,2 nm till 2,3 nm, ökar antalet atomer med förskjutning större än 0,1 nm från 9 558 till 13 767; denna kollektiva atommigrering är den direkta mikroskopiska mekanismen för pop-in.

Fig. 2 (ac) Belastningsdjupkurvor för nanoindentation på 3C-SiC (010), (110) och (111) ytor; (d) Jämförelse av tre Ph-kurvor och Hertz elastisk kontaktmodell

För att kvantitativt beskriva inverkan av spänningsfördelning på plastisk initiering, Zhu et al. beräknade von Mises skjuvspänningsfördelning. Vid ett fördjupningsdjup av 2,0 nm minskar spänningsgradienten inom den elastiska distorsionszonen; när djupet når 2,1 nm och dislokationer börjar bilda kärnor stiger spänningen i dislokationsregionen till 8–10GPa, medan spänningen mellan de förvrängda och oförvrängda områdena sjunker till 5–8GPa. Detta indikerar att dislokationsexpansion effektivt avleder lokal hög spänning, vilket tillåter en del av den förvrängda strukturen att återgå till en perfekt diamantstruktur.

2. Dislokationskärnbildning, -förökning och Thompson-tetraederanalys

Den plastiska deformationen av 3C-SiC domineras huvudsakligen av dislokationsglidning. Genom att använda en dislokationsextraktionsalgoritm, Zhu et al. identifierade att den initiala dislokationen är en perfekt dislokationsslinga med Burgers vektor b = 1/2[1-10], kärnor från området med hög skjuvspänning under indentern. Dislokationsglidning omvandlar den atomära staplingssekvensen från den ursprungliga ABCABC till ABABC, där kiselatomer avviker från sina ursprungliga positioner och bildar nya bindningar med intilliggande kolatomer längs <110>-riktningen. Denna process markerar en irreversibel övergång från elasticitet till plasticitet.

Fig. 3 Tvärsnittsvyer av dislokationszonen under nanoindentation: (a) lokal vy; (b) tvärsnitt. röda streckade pilar indikerar glidningen av kiselatomer från deras ursprungliga bindningspositioner till nya positioner

När fördjupningsdjupet ökar, utvecklas embryonala dislokationer gradvis till staplande förkastningsband och prismatiska dislokationsslingor. Om man tar (010)-planet som ett exempel: vid ett djup av 2,1nm skärs två staplingsfel (SF1 och SF2) på (111) planen och bildar en 1/2[1-10] trappstångsdislokation (SRD), som sedan ansluts till en banddislokationsslinga; när djupet ökar till 2,6 nm "slår dislokationen av" och bildar isolerade prismatiska dislokationsslingor. På liknande sätt, på (110) och (111) plan, observeras även evolutionssekvensen från embryonal dislokation → banddislokation → prismatisk dislokation, men utvecklingshastigheten och den geometriska morfologin skiljer sig på grund av anisotropi.

Fig. 4 Utveckling av dislokationer med intryckningsdjup på (010)-planet: (a) h = 2,1 nm; (b) h = 2,2 nm; (c) h = 2,5 nm; (d) h = 2,6 nm

För att systematiskt förklara denna anisotropi, Zhu et al. introducerade Thompson tetrahedronmodellen. För indragning på {100}-planet skär staplingsfel på de fyra {111}-glidplanen (ACE, AEG, ACG, CEG), var och en producerar trappstavsförskjutningar längs <110>-riktningen, med reaktionen representerad som δβ = δE + Eβ, etc. För {110}-planen är indragen, en del av ytan, indragen orsakar prismatiska dislokationsslingor att expandera vinkelrätt mot ytan; för {111}-planet uppvisar tre banddislokationer trefaldig symmetri och bildar så småningom tre prismatiska dislokationsslingor. Thompson-tetraedern ger en kortfattad men kraftfull geometrisk ram för att förstå den plastiska anisotropin hos kubisk diamantstrukturkeramik.

3. Temperaturens inverkan på mekaniska egenskaper och deformationsbeteende

I praktiska tillämpningar kan 3C-SiC-enheter fungera över ett brett temperaturområde, så temperatureffekter kan inte ignoreras. Kourbani et al. undersökte nanoindentationsresponsen för 3C-SiC på (010)-planet över ett temperaturområde på 50K till 1000K med hjälp av MD-simuleringar. De fann att när temperaturen ökar, minskar den maximala intryckningsbelastningen från 2555,4nN (50K) till 2197,4nN (1000K). Med Oliver-Pharr-metoden visade sig Youngs modul minska linjärt från 328GPa (50K) till 293,5GPa (1000K), med ett rumstemperaturvärde (300K) på cirka 318GPa, nära experimentvärdet på 314GPa; hårdheten minskar på liknande sätt från 45,3 GPa till 37,9 GPa. Denna mjukgörande effekt uppstår från ökningen av interatomärt avstånd och minskningen av Si-C kovalent bindningsenergi med stigande temperatur.

Temperaturen påverkar också avsevärt dislokationsbeteendet. Kourbani et al. analyserade dislokationsmorfologier vid olika temperaturer med användning av DXA. Perfekta dislokationer med b = 1/2<110> och Shockley partiella dislokationer med b = 1/6<112> observerades vid alla temperaturer, med perfekta dislokationer som dominerade. När temperaturen ökar ökar den totala dislokationslängden markant och bildandet av prismatiska dislokationsslingor av "lasso"-typ blir mer aktiva. Detta indikerar att termisk aktivering främjar dislokationskärnbildning och glidning, vilket förbättrar plastisk deformerbarhet.

Fig. 5 DXA- och IDS-analys: dislokationsstrukturer av 3C-SiC efter intryckning till maximalt djup (h = 2,5 nm) vid olika temperaturer (50K, 100K, 300K, 500K, 800K, 1000K)

Av särskilt intresse är temperaturens roll vid reglering av fastransformation och amorfisering. Med hjälp av metoden Identify Diamond Structure, Kourbani et al. fann att atomer med en hexagonal diamantstruktur (dvs. wurtzite B4-struktur) uppträder i fördjupningszonen, som visas i fig. 6. När temperaturen ökar minskar antalet "andra" atomer (som representerar det amorfa tillståndet), medan antalet HD-atomer (wurtzite-strukturen) också minskar vid 1000K. Radiell distributionsfunktionsanalys visar att efter indragning försvinner ordningstopparna på lång räckvidd, vilket bekräftar att amorfisering verkligen sker och blir mer uttalad vid högre temperaturer. Koordinationstalsanalys avslöjar vidare att när temperaturen ökar, minskar andelen atomer med CN = 4, medan antalet lågkoordinerade atomer (CN = 1, 2, 3) ökar, vilket indikerar att temperaturen främjar övergången av atomer från det perfekta gittret till en oordnad struktur. Dessa fynd avslöjar den dubbla effekten av temperatur på den strukturella stabiliteten av 3C-SiC: å ena sidan främjar termisk aktivering dislokationsglidning; å andra sidan förvärrar det amorfiseringsskador.

Fig. 6 Fördelning av deformerade atomer i 3C-SiC efter intryckning vid olika temperaturer (kubiska diamantstrukturatomer exkluderade): HD, HD (första granne), HD (andra granne) representerar atomer relaterade till wurtzitstrukturen; "annat" representerar amorfa atomer

4. Skador och fraktursimulering baserad på peridynamik

Även om MD-simuleringar kan lösa plastisk deformation på atomär skala, är deras beräkningsskala vanligtvis begränsad till nanometerintervallet (tiotals nanometer), vilket gör direkt jämförelse med experiment i mikrometerskala svår. Dessutom är MD ineffektivt när det gäller att hantera sprickutbredning under sprödbrott. För att lösa dessa problem simulerade Yang och Zhu indragnings- och repningsprocesserna för 3C-SiC med hjälp av en peridynamik (PD)-metod. Peridynamics är en icke-lokal teori som ersätter partiella differentialekvationer med integralekvationer, vilket gör den naturligt lämplig för att hantera diskontinuerliga problem som sprickinitiering och fortplantning.

I indragningssimuleringar observerade Yang och Zhu ett tydligt pop-in-fenomen i belastningsdjupkurvan, vilket är mycket överensstämmande med det plastiska burstbeteende som observerats i MD-simuleringar. Lastkurvorna vid olika intryckningshastigheter (1–20m/s) överlappar till stor del, men vid lossning leder högre hastigheter till snabbare lastfall. Ännu viktigare, PD-simuleringarna reproducerade framgångsrikt den indragningsinducerade sprickmorfologin: när intryckningsdjupet når 15μm med en hastighet av 5m/s, kan radiella sprickor, laterala sprickor och mediansprickor tydligt observeras, i överensstämmelse med den experimentella indragningssprickmorfologin. Detta visar att PD kan fånga sprickegenskaperna hos spröda material på kontinuumskalan.

Fig. 7 Intryckningsskadas morfologi för 3C-SiC vid en intryckningshastighet på 5m/s: (a–c) vyer ovanifrån; (d–i) tvärsnittsvyer

I skrapsimuleringar undersökte Yang och Zhu variationerna i friktionskraft och normalkraft vid olika repningsdjup (1–4μm). På djupen 1μm och 2μm uppstår inga sprickor på sidorna av repspåret, och borttagning av material är huvudsakligen plast; när djupet ökar till 3μm och 4μm uppstår radiella sprickor i spårets kanter, och både friktion och normalkrafter uppvisar periodiska fluktuationer. Den specifika skärenergin (energi som förbrukas per volymenhet borttaget material) ökar olinjärt när repningsdjupet minskar och når så högt som 22GPa vid 1μm, vilket visar en signifikant storlekseffekt.

Yang och Zhu simulerade ytterligare repningsprocessen med dubbla indrag för att efterlikna interaktionen mellan flera slipmedel i verklig slipning. När avståndet mellan de två indragen är 50 μm är repbredden ungefär 92 μm, mer än dubbelt så bred som en repa med enkel indragning, och materialavlägsningsvolymen är den största; vid ett avstånd på 60 μm är de två reporna oberoende av varandra, och borttagningsvolymen är bara dubbelt så stor som för en enda indragare; vid ett avstånd på 40 μm överlappar skadezonerna, men den totala borttagningsvolymen minskar. Detta fenomen tillskrivs överlagringseffekten av förskjutningsfältet: på lämpligt avstånd ansluter sprickor mellan de två indragen, vilket leder till samverkande materialavlägsnande; för litet avstånd resulterar i ömsesidigt undertryckande av förskjutningsfälten, medan ett för stort avstånd inte ger någon interaktion. Detta fynd ger en teoretisk grund för att optimera malningsprocessparametrar.

Fig. 8 Skademorfologi av repor med dubbel indrag (indragningsdjup 3μm): (a,d) avstånd 40μm, oregelbundna repor; (b,e) avstånd 50μm, avsevärt ökad repbredd, maximal materialavlägsningsvolym; (c,f) avstånd 60μm, två oberoende repor

Fig. 9 Sidovyer i tvärsnitt av repor med dubbel indrag: (a) avstånd 40 μm; (b) avstånd 50 μm; (c) avstånd 60μm

Genom en kombination av MD- och PD-simuleringar har forskarna avslöjat den anisotropa deformationen, dislokationsutvecklingen, temperatureffekterna och mekanismerna för skada/fraktur av 3C-SiC under nanoindentation. Studien finner att elastiska och plastiska beteenden är starkt beroende av kristallorientering, med dislokationer som kärnor längs {111}<110> glidsystemet och utvecklas till prismatiska slingor; ökad temperatur minskar hårdhet och modul samtidigt som amorfisering främjas. PD-metoden reproducerar framgångsrikt pop-in-fenomenet vid indragning och sprickutbredning vid repning, och identifierar ett optimalt dubbelindragsavstånd för att förbättra materialavlägsningseffektiviteten. Framtida arbete kan ytterligare utveckla flerskaliga kopplingsmodeller och maskininlärda potentiella funktioner för att uppnå kvantitativa förutsägelser från atomära mekanismer till bearbetningsprocesser.

Oavsett om du behöver 3C-SiC wafer för forskning eller för industriella tillämpningar, vänligen kontakta oss via e-post på[email protected] och [email protected].

 

Referenser:

1. Zhu, B., Zhao, D., Zhang, Z., Niu, Y., Zhang, Z., Zhao, J., … & Zhao, H. (2024). Atomstudie om deformationsbeteende och anisotropieffekt av 3C-SiC under nanoindentation. Journal of Materials Research and Technology, 28, 2636-2647.
2. Kourbani, H., Bogtob, S., Samiri, A., Hassani, A., & Hasnaoui, A. (2025). Termisk nanoindentation av (3C) kiselkarbid 3C-SiC med simulering av molekylär dynamik. Solid State Communications, 403, 116013.
3. Xu, Y., & Zhu, P. (2022). Peridynamiska simuleringar av skador i indrag och repning av 3C-SiC. Journal of Materials Research, 37(23), 4381-4391.