Comportamento di deformazione e meccanismo di danno del 3C-SiC sotto nanoindentazione +
PAM-XIAMEN fornisce wafer 3C-SiC per scopi di ricerca. Per i dettagli, si prega di contattare:[email protected]
Il carburo di silicio cubico (3C-SiC) presenta un grande potenziale in campi all'avanguardia come i sistemi microelettromeccanici, l'informatica quantistica e i dispositivi optoelettronici, grazie all'ampio gap di banda, all'elevata conduttività termica, all'elevata intensità del campo di rottura critica e all'eccellente inerzia chimica. Tuttavia, la durezza ultraelevata e la fragilità intrinseca del 3C-SiC rendono la sua lavorazione ultraprecisa estremamente impegnativa e i danni superficiali/sottosuperficiali indotti dalla lavorazione sono diventati un collo di bottiglia chiave che ne limita le applicazioni di fascia alta. Pertanto, una profonda comprensione del comportamento di deformazione e dei meccanismi di danneggiamento del 3C-SiC sotto carico di contatto è di grande importanza teorica e ingegneristica per ottimizzarne le tecniche di lavorazione.
Negli ultimi anni, i ricercatori hanno esplorato ampiamente la risposta meccanica del 3C-SiC combinando esperimenti di nanoindentazione con simulazioni di dinamica molecolare. Zhu et al. hanno studiato sistematicamente il comportamento di nanoindentazione di 3C-SiC sui piani cristallini (010), (110) e (111) utilizzando simulazioni di dinamica molecolare (MD), rivelando le caratteristiche anisotrope della deformazione elastica e dell'evoluzione delle dislocazioni, e hanno introdotto il tetraedro di Thompson per spiegare il meccanismo intrinseco della deformazione plastica. Nel frattempo, la temperatura, in quanto fattore chiave negli ambienti di servizio reali, influisce in modo significativo sulle proprietà meccaniche del 3C-SiC. Kourbani et al. scoperto attraverso simulazioni MD che l'aumento della temperatura riduce la durezza e il modulo di Young del materiale e promuove l'amorfizzazione. D'altra parte, i metodi tradizionali della meccanica del continuo presentano limitazioni nell'affrontare problemi discontinui come l'innesco e la propagazione delle cricche. Per risolvere questo problema, Yang e Zhu hanno utilizzato la peridinamica per simulare l'evoluzione del danno del 3C-SiC sotto indentazione e graffi, riproducendo con successo il fenomeno del pop-in e la formazione di crepe radiali. Questi studi forniscono riferimenti multiscala per comprendere il comportamento meccanico e la teoria della lavorazione del 3C-SiC.
1. Deformazione elastica anisotropa e comportamento plastico iniziale del 3C-SiC
Il 3C-SiC ha una struttura di blenda di zinco e la sua risposta meccanica dipende fortemente dall'orientamento del cristallo. Attraverso simulazioni MD, Zhu et al. ha scoperto che durante la nanoindentazione sui piani (010), (110) e (111), le regioni di deformazione elastica mostrano simmetrie nettamente diverse: il piano (010) mostra una simmetria quadrupla, con una zona di deformazione piramidale che si espande lungo la direzione <110>; il piano (110) presenta una simmetria centrale, con la zona di deformazione divisa in due parti; il piano (111) mostra una tripla simmetria, con tre zone di distorsione elastica che si espandono anche lungo la direzione <110>. Questo risultato indica che la direzione <110> è il percorso ottimale per il trasferimento elastico della sollecitazione.
Fig. 1 Morfologia della zona di distorsione elastica di 3C-SiC ad una profondità di indentazione di 1,1 nm: piano (a,d) (010); (b,e) (110) aereo; (c,f) (111) piano; le frecce rosse indicano la direzione di propagazione della deformazione elastica
Dalle curve carico-profondità, la lunghezza dello stadio elastico varia significativamente tra i diversi piani cristallini. La curva Ph per il piano (010) corrisponde molto bene al modello di contatto Hertz per h<2,0 nm, mentre gli stadi elastici per i piani (110) e (111) sono significativamente più brevi, indicando uno snervamento plastico precedente. In particolare, gli eventi pop-in appaiono nelle curve per i piani (010) e (110), ma non sono osservati per il piano (111). Attraverso l'analisi dello spostamento atomico, Zhu et al. ha confermato che l'evento pop-in corrisponde ad un movimento direzionale di un gran numero di atomi lungo la direzione <110>, portando ad un improvviso rilascio di energia elastica immagazzinata. Ad esempio, quando la profondità di indentazione sul piano (110) aumenta da 2,2 nm a 2,3 nm, il numero di atomi con spostamento maggiore di 0,1 nm aumenta da 9.558 a 13.767; questa migrazione atomica collettiva è il meccanismo microscopico diretto del pop-in.
Fig. 2 (ac) Curve di profondità di carico della nanoindentazione su superfici 3C-SiC (010), (110) e (111); (d) Confronto tra tre curve Ph e modello di contatto elastico Hertz
Per descrivere quantitativamente l'influenza della distribuzione dello stress sull'iniziazione plastica, Zhu et al. calcolato la distribuzione dello sforzo di taglio di von Mises. Ad una profondità di indentazione di 2,0 nm, il gradiente di sollecitazione diminuisce all'interno della zona di distorsione elastica; quando la profondità raggiunge i 2,1 nm e le dislocazioni iniziano a nuclearsi, lo stress nella regione della dislocazione aumenta a 8–10 GPa, mentre lo stress tra le regioni distorte e non distorte scende a 5–8 GPa. Ciò indica che l’espansione della dislocazione dissipa efficacemente l’elevato stress locale, consentendo a parte della struttura distorta di ritornare ad una perfetta struttura a diamante.
2. Nucleazione, propagazione delle dislocazioni e analisi del tetraedro di Thompson
La deformazione plastica del 3C-SiC è principalmente dominata dallo scorrimento delle dislocazioni. Utilizzando un algoritmo di estrazione della dislocazione, Zhu et al. identificato che la dislocazione iniziale è un anello di dislocazione perfetto con il vettore di Burgers b = 1/2[1-10], nucleante dalla regione ad alto stress di taglio sotto il penetratore. Lo scorrimento delle dislocazioni trasforma la sequenza di impilamento atomico dall'originale ABCABC ad ABABC, dove gli atomi di silicio deviano dalle loro posizioni originali e formano nuovi legami con gli atomi di carbonio adiacenti lungo la direzione <110>. Questo processo segna un passaggio irreversibile dall’elasticità alla plasticità.
Fig. 3 Cross-sectional views of the dislocation zone during nanoindentation: (a) local view; (b) cross-section; red dashed arrows indicate the slip of silicon atoms from their original bonding positions to new positions
As the indentation depth increases, embryonic dislocations gradually evolve into stacking fault ribbons and prismatic dislocation loops. Taking the (010) plane as an example: at a depth of 2.1nm, two stacking faults (SF1 and SF2) located on the (111) planes intersect, forming a 1/2[1-10] stair-rod dislocation (SRD), which then connects into a ribbon dislocation loop; when the depth increases to 2.6nm, the dislocation “punches off”, forming isolated prismatic dislocation loops. Similarly, on the (110) and (111) planes, the evolution sequence from embryonic dislocation → ribbon dislocation → prismatic dislocation is also observed, but the evolution speed and geometric morphology differ due to anisotropy.
Fig. 4 Evolution of dislocations with indentation depth on the (010) plane: (a) h = 2.1nm; (b) h = 2.2nm; (c) h = 2.5nm; (d) h = 2.6nm
Per spiegare sistematicamente questa anisotropia, Zhu et al. ha introdotto il modello del tetraedro di Thompson. Per l'indentazione sul piano {100}, i difetti di impilamento sui quattro piani di scorrimento {111} (ACE, AEG, ACG, CEG) si intersecano, ciascuno producendo dislocazioni delle aste lungo la direzione <110>, con la reazione rappresentata come δβ = δE + Eβ, ecc. Per il piano {110}, alcuni piani di scorrimento sono perpendicolari alla superficie di indentazione, causando l'espansione degli anelli di dislocazione prismatica perpendicolare alla superficie; per il piano {111}, tre dislocazioni del nastro mostrano una simmetria tripla, formando infine tre anelli di dislocazioni prismatiche. Il tetraedro di Thompson fornisce una struttura geometrica concisa ma potente per comprendere l'anisotropia plastica delle ceramiche a struttura cubica diamantata.
3. Effetto della temperatura sulle proprietà meccaniche e sul comportamento deformativo
Nelle applicazioni pratiche, i dispositivi 3C-SiC possono funzionare in un ampio intervallo di temperature, quindi gli effetti della temperatura non possono essere ignorati. Kourbani et al. hanno studiato la risposta alla nanoindentazione di 3C-SiC sul piano (010) in un intervallo di temperature compreso tra 50K e 1000K utilizzando simulazioni MD. Hanno scoperto che all’aumentare della temperatura, il carico massimo di indentazione diminuisce da 2555,4 nN (50 K) a 2197,4 n N (1000 K). Utilizzando il metodo Oliver-Pharr, è stato riscontrato che il modulo di Young diminuisce linearmente da 328 GPa (50 K) a 293,5 GPa (1000 K), con un valore della temperatura ambiente (300 K) di circa 318 GPa, vicino al valore sperimentale di 314 GPa; allo stesso modo la durezza diminuisce da 45,3 GPa a 37,9 GPa. Questo effetto ammorbidente deriva dall'aumento della spaziatura interatomica e dalla riduzione dell'energia del legame covalente Si–C con l'aumento della temperatura.
Anche la temperatura influenza in modo significativo il comportamento della dislocazione. Kourbani et al. analizzato le morfologie delle lussazioni a diverse temperature utilizzando DXA. Dislocazioni perfette con b = 1/2<110> e dislocazioni parziali di Shockley con b = 1/6<112> sono state osservate a tutte le temperature, con dominanza di dislocazioni perfette. All’aumentare della temperatura, la lunghezza totale della dislocazione aumenta notevolmente e la formazione di anse di dislocazione prismatica di tipo “lazo” diventa più attiva. Ciò indica che l’attivazione termica promuove la nucleazione e lo scorrimento delle lussazioni, migliorando la deformabilità plastica.
Fig. 5 Analisi DXA e IDS: strutture di dislocazione di 3C-SiC dopo indentazione alla massima profondità (h = 2,5 nm) a diverse temperature (50K, 100K, 300K, 500K, 800K, 1000K)
Di particolare interesse è il ruolo della temperatura nella regolazione della trasformazione di fase e dell'amorfizzazione. Utilizzando il metodo Identifica struttura del diamante, Kourbani et al. hanno scoperto che gli atomi con una struttura a diamante esagonale (cioè la struttura della wurtzite B4) compaiono nella zona di rientranza, come mostrato in Fig. 6. All'aumentare della temperatura, il numero di "altri" atomi (che rappresentano lo stato amorfo) diminuisce, mentre anche il numero di atomi HD (struttura della wurtzite) diminuisce a 1000K. L'analisi della funzione di distribuzione radiale mostra che dopo l'indentazione, i picchi dell'ordine a lungo raggio scompaiono, confermando che l'amorfizzazione avviene effettivamente e diventa più pronunciata a temperature più elevate. L'analisi dei numeri di coordinazione rivela inoltre che all'aumentare della temperatura, la proporzione di atomi con CN = 4 diminuisce, mentre il numero di atomi a bassa coordinazione (CN = 1, 2, 3) aumenta, indicando che la temperatura promuove la transizione degli atomi dal reticolo perfetto a una struttura disordinata. Questi risultati rivelano il duplice effetto della temperatura sulla stabilità strutturale del 3C-SiC: da un lato, l'attivazione termica promuove lo scorrimento della dislocazione; d'altro canto, esacerba il danno da amorfizzazione.
Fig. 6 Distribuzione degli atomi deformati nel 3C-SiC dopo indentazione a diverse temperature (esclusi gli atomi con struttura cubica del diamante): HD, HD (1° vicino), HD (2° vicino) rappresentano atomi legati alla struttura della wurtzite; “altro” rappresenta gli atomi amorfi
4. Simulazione di danni e fratture basata sulla peridinamica
Sebbene le simulazioni MD possano risolvere la deformazione plastica su scala atomica, la loro scala computazionale è tipicamente limitata alla gamma dei nanometri (decine di nanometri), rendendo difficile il confronto diretto con gli esperimenti su scala micrometrica. Inoltre, l’MD è inefficiente nel gestire la propagazione delle cricche durante la frattura fragile. Per affrontare questi problemi, Yang e Zhu hanno simulato i processi di indentazione e graffiatura del 3C-SiC utilizzando un approccio peridinamico (PD). La peridinamica è una teoria non locale che sostituisce le equazioni differenziali parziali con equazioni integrali, rendendola naturalmente adatta a gestire problemi discontinui come l'innesco e la propagazione delle cricche.
Nelle simulazioni di indentazione, Yang e Zhu hanno osservato un chiaro fenomeno di pop-in nella curva carico-profondità, che è altamente coerente con il comportamento di scoppio della plastica osservato nelle simulazioni MD. Le curve di carico a diverse velocità di rientro (1–20 m/s) si sovrappongono ampiamente, ma durante lo scarico, velocità più elevate portano a cadute di carico più rapide. Ancora più importante, le simulazioni PD hanno riprodotto con successo la morfologia della fessura indotta dall'indentazione: quando la profondità dell'indentazione raggiunge i 15μm a una velocità di 5 m/s, si possono osservare chiaramente crepe radiali, laterali e mediane, coerenti con la morfologia sperimentale della fessura dell'indentazione. Ciò dimostra che la PD può catturare le caratteristiche di frattura dei materiali fragili su scala continua.
Fig. 7 Morfologia del danno da indentazione di 3C-SiC a una velocità di indentazione di 5 m/s: (a–c) viste dall'alto; (d–i) viste in sezione trasversale
Nelle simulazioni di graffio, Yang e Zhu hanno esaminato le variazioni della forza di attrito e della forza normale a diverse profondità di graffio (1–4μm). A profondità di 1μm e 2μm non compaiono crepe sui lati della scanalatura antigraffio e la rimozione del materiale è prevalentemente plastica; quando la profondità aumenta a 3μm e 4μm, compaiono crepe radiali sui bordi della scanalatura e sia l'attrito che le forze normali mostrano fluttuazioni periodiche. L'energia di taglio specifica (energia consumata per unità di volume di materiale rimosso) aumenta in modo non lineare al diminuire della profondità di graffiatura, raggiungendo un massimo di 22GPa a 1μm, dimostrando un effetto dimensionale significativo.
Yang e Zhu hanno simulato ulteriormente il processo di graffiatura a doppio penetratore per imitare l'interazione di più abrasivi nella rettifica effettiva. Quando la spaziatura tra i due penetratori è di 50μm, la larghezza del graffio è di circa 92μm, più del doppio della larghezza di un graffio a penetratore singolo e il volume di rimozione del materiale è il più grande; ad una distanza di 60μm, i due graffi sono indipendenti e il volume di rimozione è solo il doppio di quello di un singolo penetratore; con una distanza di 40μm le zone danneggiate si sovrappongono, ma il volume totale di rimozione diminuisce. Questo fenomeno è attribuito all'effetto di sovrapposizione del campo di spostamento: ad una distanza adeguata, le fessure si collegano tra i due penetratori, portando alla rimozione cooperativa del materiale; una spaziatura troppo piccola provoca la soppressione reciproca dei campi di spostamento, mentre una spaziatura troppo grande non produce alcuna interazione. Questa scoperta fornisce una base teorica per ottimizzare i parametri del processo di macinazione.
Fig. 8 Morfologia del danno da graffio a doppio penetratore (profondità di indentazione 3μm): (a,d) spaziatura 40μm, graffi irregolari; (b,e) spaziatura 50μm, larghezza del graffio notevolmente aumentata, volume massimo di rimozione del materiale; (c,f) spaziatura 60μm, due graffi indipendenti
Fig. 9 Viste laterali in sezione trasversale del graffio con doppio penetratore: (a) spaziatura 40μm; (b) spaziatura 50μm; (c) spaziatura 60μm
Attraverso una combinazione di simulazioni MD e PD, i ricercatori hanno rivelato la deformazione anisotropa, l'evoluzione delle dislocazioni, gli effetti della temperatura e i meccanismi di danno/frattura del 3C-SiC sottoposto a nanoindentazione. Lo studio rileva che i comportamenti elastici e plastici dipendono fortemente dall'orientamento dei cristalli, con dislocazioni che si nucleano lungo il sistema di scorrimento {111<110> e si evolvono in anelli prismatici; l'aumento della temperatura riduce la durezza e il modulo promuovendo l'amorfizzazione. Il metodo PD riproduce con successo il fenomeno del pop-in nell'indentazione e la propagazione della cricca nel graffio e identifica una spaziatura ottimale tra il doppio penetratore per migliorare l'efficienza di rimozione del materiale. Il lavoro futuro potrebbe sviluppare ulteriormente modelli di accoppiamento multiscala e potenziali funzioni apprese dalle macchine per ottenere previsioni quantitative dai meccanismi atomici ai processi di lavorazione.
Se hai bisogno di wafer 3C-SiC per la ricerca o per applicazioni industriali, contattaci via e-mail all'indirizzo[email protected] e [email protected].
Riferimenti:
- Zhu, B., Zhao, D., Zhang, Z., Niu, Y., Zhang, Z., Zhao, J., … & Zhao, H. (2024). Studio atomico sul comportamento di deformazione e sull'effetto anisotropico del 3C-SiC sotto nanoindentazione. Giornale di ricerca e tecnologia dei materiali, 28, 2636-2647.
- Kourbani, H., Bogtob, S., Samiri, A., Hassani, A., & Hasnaoui, A. (2025). Nanoindentazione termica del carburo di silicio (3C) 3C-SiC utilizzando la simulazione della dinamica molecolare. Comunicazioni a stato solido, 403, 116013.
- Xu, Y. e Zhu, P. (2022). Simulazioni peridinamiche di danno nell'indentazione e nel graffio del 3C-SiC. Giornale di ricerca sui materiali, 37(23), 4381-4391.