Verkenning van tip-gebaseerde nanomachining met moleculaire dynamica: een uitgebreide review

Samenvatting

Tip-gebaseerde nanomachining (TBN) heeft zich beweten als een krachtige en haalbare techniek voor de fabricage van microstructuren. Moleculaire dynamica (MD)-simulaties worden breed toegepast in TBN-onderzoek om mechanismen te ontrafelen die experimenteel niet volledig te achterhalen zijn. Deze review geeft een overzicht van de recente wetenschappelijke vooruitgang in MD-simulaties van TBN. Allereerst worden de opzetmethoden van simulatiemodellen voor diverse materialen besproken. Vervolgens wordt de analyse van het machinisme behandeld, waaronder snijkraftanalyse, materiaalverwijdering en defectanalyse in de onderlaag. Tot slot worden huidige tekortkomingen en toekomstperspectieven van de TBN-methode in MD-simulaties gegeven. Deze review biedt een waardevolle referentie voor vervolgonderzoek.

Inleiding

Micro-/nanofabricagetechnologie vindt brede toepassing in sectoren als milieu, energie, biologie, geneeskunde en landverdediging, en speelt een steeds belangrijkere rol in nationale ontwikkeling en maatschappelijke vooruitgang [1-4]. Voor hoogprecize micro-/nanofabricage is de op atoomkrachtmicroscopie (AFM) gebaseerde mechanische benadering, bekend als TBN, aantrekkelijk dankzij eenvoudige apparatuur, nanometerprecisie en lage omgevingsvereisten. Tot nu toe zijn nanodotten, lijnen/groeven, 2D/3D-structuren en zelfs nanostructuren op gekromde oppervlakken met succes gerealiseerd via TBN [5-7].

Om nanostructuren met nog hogere precisie te fabriceren, is een diepgaand begrip van het TBN-proces essentieel. Veel onderzoekers hebben het mechanisme experimenteel bestudeerd in metalen [8], halfgeleiders [9,10] en polymeren [11]. Echter, experimenten bieden beperkt inzicht in de interne mechanismen. MD-simulaties zijn een krachtig hulpmiddel om het nanomachiningsproces te verduidelijken [12,13]. In vergelijking met experimenten kan MD het micromechanisme van materiaalverwijdering en oppervlaktevorming beter verklaren via analyse van snijkrachten, spanningsstaat, energie-dissipatie en oppervlaktetopografie [12].

In het verleden hebben diverse auteurs MD-simulaties gebruikt voor TBN-processen. Fang et al. onderzochten de invloed van voedersnelheid en kegelhoek op nanolithografie [14,15]. Isono en Tanaka analyseerden temperatuur, bewerkbaarheid en interatoomkrachten van nikkel [16,17]. Yan et al. bestudeerden tipgeometrie in AFM-gebaseerde lithografie [18]. Recent zijn nieuwe TBN-processen met MD onderzocht. Xiao et al. vergelijken statische en dynamische ploeging en tonen aan dat dynamische ploeging kleinere nanostructuren oplevert [19]. Geng et al. simuleerden last-gecontrolleerd nanokrabben door een constante normalkracht op de tip toe te passen [20]. Voor hogere dichtheid nanogroeven is de minimale voeder (MF) cruciaal. Ren et al. presenteerden een grof-naar-fijn criterium om MF te bepalen met MD [21]. Eerdere reviews [12,13,22-24] omvatten deze nieuwe vooruitgang niet. Deze review richt zich daarom op de recente state-of-the-art van MD-simulaties voor TBN, inclusief modelleringstechnieken, nieuwe machiningsmethodes, mechanisme-analyse (snijkraft, materiaalverwijdering, onderlaagdefecten) en resterende uitdagingen.

Simulatiemethoden

Opzet van het MD-model

Voor nauwkeurige voorspellingen is optimalisatie van het simulatiemodel en de verwerkingsconfiguratie essentieel. Het model omvat de atoomconfiguratie en potentiaalfuncties. Atoomconfiguraties zijn kristallijn (bijv. koper, ijzer, silicium) of amorfs (polymers). Enkeltkristallen worden opgebouwd door herhaling van de eenheidscel; de oppervlakteoriëntatie beïnvloedt de bewerkingseigenschappen [29-31]. Poly kristallijne structuren met verschillende korrelgroottes worden samengesteld via de Voronoi-tessellatiemethode [32]. Nanotwined (NT) materialen, een speciale vorm van polykristallen, trekken aandacht door hun uitstekende mechanische eigenschappen [33-36]. Als voorbeeld wordt de constructie van NT-polykristallijn Cu beschreven [25]: eerst wordt een multicouche van gelijke enkeltkristallijnen Cu-lagen gebouwd met twin boundaries (TB) tussen opeenvolgende lagen. Vervolgens worden korrelhoeken berekend en wordt de NT-polykristallijne structuur met periodieke randvoorwaarden gegenereerd via Voronoi-constructie. Figuur 1a toont de atoomconfiguraties van enkeltkristallijn Cu, polykristallijn Cu en NT-polykristallijn Cu.

Polymers worden gemodelleerd met coarse-grained modellen. Voor polyethyleen (PE) wordt een initiëel systeem van 10 moleculaire ketens gegenereerd via een Monte Carlo self-avoiding random walk algoritme [39]. De simulatiebox is ongeveer 80,06 × 80,06 × 80,06 ų. Het MD-model van PE is weergegeven in Fig. 1d.

Tipgeometrieën variëren (kegel, driehoekige piramide, bolvormig) en beïnvloeden het TBN-proces sterk. De piramidale tip komt overeen met AFM-machining, de bolvormige met afschuiving. Meestal wordt de tip als oneindig stijf beschouwd, door de relatieve atoomposities constant te houden en met constante snelheid te verplaatsen [41].

De keuze van potentiaalfuncties bepaalt de betrouwbaarheid. De Embedded Atom Method (EAM) is geschikt voor metalen (Cu [42], Fe [43]). Tersoff [44] en Stillinger-Weber (SW) [45] potentiealen zijn geschikt voor diamantkubische structuren. Voor silicium in TBN is de SW-potentieel superieur in het beschrijven van dislocatie-nucleatie en fasetransformaties [46,47]. Analytical Bond Order Potential (ABOP), Reactive Empirical Bond Order (REBO) en Adaptive Intermolecular Reactive Empirical Bond Order (AIREBO) zijn uitbreidingen van Tersoff. ABOP is ideaal voor Si-C interacties (SiC). REBO is populair voor koolstof en koolstofnanobuizen. AIREBO verbetert REBO en modelleert complexere interacties [41]. Polystyreen wordt beschreven met AIREBO [49].

De meeste simulaties draaien in LAMMPS [50]. Complexe modellen (polykristallen, NT) worden gegenereerd met Atomsk [51]. Polymermodellen komen uit Material Studio [52]. Visualisatie geschiedt met OVITO [53] en VMD [54].

Ondanks de effectiviteit blijven meeste werkstukmodellen kleiner dan 50 nm³, wat afwijkingen van de werkelijkheid kan geven [40,55,56]. Bovendien ontbreken potentiealen voor materialen als GaAs en Lu₂O₃. Verdere optimalisatie van MD-modellen en potentiealen is noodzakelijk.

Verwerkingsconfiguratie

In TBN-experimenten is de last-gecontrolleerde modus (constante normalkracht) gebruikelijk en biedt hogere precisie, vooral op hellende of gekromde oppervlakken. Veel MD-simulaties gebruiken echter verplaatsingsgestuurde modi, wat verschillen met experimenten kan opleveren [18,57]. Recent worden last-gecontrolleerde nanokrab-simulaties uitgevoerd door direct een constante normalkracht toe te passen [20,58]. Het proces heeft drie fasen: relaxatie, penetratie en krabben. De tip is niet perfect loodrecht door balkbuiging van de AFM, substraatafwijking en oppervlakteruwheid; tip-kanteling mag niet genegeerd worden. Liu et al. onderzochten de invloed van tip-inclinatiehoeken [59]. Tip-slijtage is cruciaal voor machiningskwaliteit. Waterlaag-smering kan slijtage verminderen, maar het mechanisme is nog onvoldoende begrepen. Ren et al. simuleerden het effect van waterlaag-smering op monokristallijn koper [60].

Enkelvoudig krabben produceert nanogroeven, maar beperkt de afmetingen. Meervoudig krabben vergroot de groef [62]. Geng et al. vergelijken enkelvoudig en meervoudig krabben met MD en experimenten [20]. Beide zijn statische ploeglithografie (contactmodus AFM). Lange statische ploeging veroorzaakt significante tipslijtage. Tapping-modus machining (dynamische ploeglithografie) kan slijtage verminderen. Figuur 2 toont schematisch statische en dynamische ploeging op enkeltkristallijn koper. Bij statische ploeging daalt de diamanttip, ploegt langs de negatieve x-as en stijgt weer op. Bij dynamische ploeging beweegt de tip langs een sinusoïdaal pad [19]. Op basis van tapping-modus ontwikkelde Sundaram et al. Vibration-Assisted Nano Impact machining by Loose Abrasives (VANILA): de AFM-tip trilt, geeft kinetische energie aan nanoafschuivende deeltjes in een slurry, wat nanometer-schaal materiaalverwijdering veroorzaakt [61]. Figuur 2b,c toont het VANILA-proces.

Analyse van simulateresultaten

Snijkraftanalyse

De snijkraft weerspiegelt het materiaalverwijderingsproces en defectevolutie, een cruciale fysieke parameter voor het begrijpen van het snijfenomeen [19,63,64]. Deze kracht hangt samen met snijvermogen en toollevensduur, en biedt begeleiding voor TBN [65]. In TBN wordt de kracht gemeten via krachtkalibratie, maar laterale krachten blijven onzichtbaar [66,67]. MD maakt realtime waarneming van zowel laterale als tangentiële krachten mogelijk gedurende het hele proces [68], evenals gemiddelde en resultante krachten.

Uit snijkraftvariatie worden verschillen in materiaalstructuur en invloed van bewerkingsparameters zichtbaar. Li et al. vonden hogere krachten bij polykristallijn Cu ten opzichte van enkeltkristallijn Cu, toegeschreven aan anisotropie van korreloriëntaties, terwijl enkeltkristallijn Cu stabiele plastische vloeistof toont [25]. Snijkracht neemt toe met snij snelheid door meer chipsvorming [69]. In SiC gedraagt de kracht zich tegenovergesteld: hogere snelheid genereert meer amorfe structuur, waardoor het materiaal duktiler en makkelijker te verwijderen is [70]. Yan et al. vonden bij Cu/Ni-bilayers dat Ni-Cu de hoogste kracht vertoont, doordat de interface dislocatiepropagatie belemmert [71]. Bij amorfe polymers beïnvloedt krab snelheid de machiningskracht op drie manieren [72]: hogere stapelhoogte voor de tip verhoogt de tangentiële kracht; hogere vervormingssnelheid veroorzaakt strain hardening, wat krachten verhoogt; thermische verzachting bij hogere snelheid verlaagt krachten. In polymers domineren stapelvorming en strain hardening. Tipgeometrie speelt ook een grote rol [18,73,74]. Ren et al. toonden aan dat krachten stijgen met toenemende half-tophoek van kegelvormige tips door grotere contactoppervlakte [21]. Wrijvingscoëfficiënt daalt sterk met half-tophoek, hardheid neemt toe [75].

Specifieke machiningscondities zijn ook onderzocht. Ren et al. analyseerden de correlatie tussen waterlaagdikte en krabbekrachten (Fig. 3). In tegenstelling tot macrokrabben, waar water smeert, domineert de weerstand van de waterlaag en leidt een dikker laag tot grotere totale krabbekrachten [60]. Tip-inclinatie effecten op hardheid (normaalkracht per atoom) en wrijvingscoëfficiënt tonen aan dat normaalhardheid gevoeliger is voor voorwaarts/achterwaarts richting, laterale inclinatie effecten zijn verwaarloosbaar. De invloed op de normaaalkracht is de hoofdoorzaak van wrijvingscoëfficiëntverandering [59].

Vergelijking van statische en dynamische ploeging onthult mechanistische verschillen. De snijkracht in dynamische ploeging oscilleert sterk met de periodieke tapping, terwijl statische ploeging lichte fluctuaties toont. De gemiddelde machiningskracht in dynamische ploeging is minder dan de helft van die in statische ploeging, wat minder tipslijtage impliceert [19]. In dynamische ploeging beïnvloedt tiporiëntatie de snijkracht sterk. Yan et al. combineerden groefdiepte en ploegkracht voor efficiëntievergelijking: sideface-forward heeft de beste efficiëntie (vergelijkbare kracht als face-forward, maar meer dan tweemaal de groefdiepte) [63]. Snijkraftanalyse biedt inzicht in TBN-mechanismen en richtlijnen voor procesoptimalisatie.

Analyse van materiaalverwijderingstoestand

Materiaalverwijderingstoestand is essentieel voor mechanismebegrip. Experimenteel (SEM, AFM) is alleen na-bewerking morfologie zichtbaar; het dynamische verwijderingsproces blijft verborgen. MD-simulatie maakt realtime observatie mogelijk van materiaalverwijdering en oppervlaktegeneratie op nanoschaal.

Uit MD-resultaten blijkt dat vervormingstoestanden vallen in ploeging en snijden. Tijdens snijden worden atomen voor de tip samengeperst tot continue chips; tegelijkertijd stapelen atomen zich aan de groefzijden op [57]. De verhouding snijden/niet-snijden bij variërende tipradius en krabdiepte bepaalt de dominante toestand [77]. Enkeltkristallijne materialen vertonen anisotropie in TBN, wat materiaalverwijdering sterk beïnvloedt [20]. Bij polykristallen is korreloriëntatie dominanter dan korrelgrootte voor stapelprofiel [78]. Tipgeometrie heeft sterke invloed. Kegelvormige tips tonen duidelijke afhankelijkheid van half-tophoek: grotere hoek geeft meer chipvolume, gladder oppervlak, maar hogere kracht, wrijving, temperatuur en onderlaagschade [40]. Alhafez toonde aan dat bij kleine half-tophoek stapeling voornamelijk voorwaarts is, bij grote hoek lateraal [75]. Driehoekige piramide tips (consistent met AFM) worden veel bestudeerd [19,20,63,79]. Drie machiningsrichtingen: edge-forward, face-forward, sideface-forward. Bij SiC kan materiaalverwijdering worden gecontroleerd via krabrichting; edge-forward geeft stabielere processen en betere grootteconsistentie [79]. Polymers zijn gevoeliger voor temperatuur dan metalen/halfgeleiders. Lokale temperatuur in krabzone overschrijdt glasovergangstemperatuur, wat duktiele verwijdering mogelijk maakt [72]. Zhan et al. vonden dat moleculaire ketens bij hogere temperatuur flexibeler worden, wat makkelijkere vervorming en herstel toelaat [56]. Du et al. toonden aan dat machiningsnelheid polystyrene resultaten sterk beïnvloedt: lage snelheid geeft intra-keten veranderingen (permanente vervorming), hoge snelheid geeft inter-keten glijden [80].

Onder droge condities geconcentreerd, tonen studies met waterlaag aan dat atoomopstapeling voor de tool afneemt, graten minder duidelijk zijn, en oppervlakteruwheid afneemt met toenemende waterlaagdikte, wat leidt tot gladde oppervlaktes [60]. De waterlaag fungeert als smeermiddel, vermindert het plakgebied tussen tool en verse chip, en verlaagt de wrijving aan de chip-tool grens. In nieuwe benaderingen onderzocht Shockly et al. parameterinvloeden (impactsnelheid, hoek, temperatuur) in vibratie-ondersteunde nano-impact machining op nanocaviteitsvorming [61] (Fig. 4). Xiao et al. vonden dat dynamische ploeging kleinere groefdiepte en -breedte oplevert dan statische ploeging, wat fabricage van kleinere kenmerken mogelijk maakt; groefafmetingen zijn controleerbaar via de dryf-amplituderatio [19].

Tijdens materiaalverwijdering treedt tipslijtage op. Meestal wordt de tip als rigide lichamen gemodelleerd, waardoor slijtage niet direct observeerbaar is [13]. Slijtage wordt indirect bestudeerd via materiaalverwijdering, spanningen en temperatuurverdeling [70,77]. Meng et al. modeleerden de tip als deformeerbaar en observeerden slijtagegedrag [81]: initiële adhesieve slijtage, gevolgd door atoom-voor-atoom slijtage in stabiele fase; atoom-voor-atoom bijdrage wordt nauwelijks beïnvloed door snijsnelheid. Er zijn nog weinig rapporten over tipslijtage; deformeerbare tipmodellen kunnen dit onderzoek bevorderen.

Materiaalverwijderingstoestandanalyse ondersteunt realtime procesmonitoring. Door lengteschaalbeperkingen blijft analyse kwalitatief; kwantitatieve voorspelling van bewerkingsresultaten is moeilijk.

Defectevolutieprocesanalyse

Defectinformatie wordt experimenteel verkregen via TEM na FIB-monsterbereiding, wat complex, duur en materiaalafhankelijk is. MD-simulaties maken defecten onder het oppervlak en hun evolutie tijdens krabben direct zichtbaar. Beschikbare algoritmes voor defectextractie: Common Neighbor Analysis (CNA) [82], Centro-Symmetry Parameter (CSP) [83], slip vector analyse [84], Ackland-Jones analyse [85]. Deze zijn ongeschikt voor dislocatiepropagatie tracing bij grote aantallen dislocaties [86]. Chen et al. introduceerden slip-analyse om atomaire glijding te visualiseren, zelfs bij grote atomaantallen [86]. Xiao et al. pasten dit toe op dynamische en statische ploeging [19] (Fig. 5). Dynamische ploeging toont naar beneden en zijwaartse materiaalstroom; statische ploeging toont geen duidelijke neerwaartse stroom. Dislocatiepropagatie hangt af van tiporiëntatie, wat leidt tot verschillende groefmorfologieën. Dislocation Extraction Algorithm (DXA) identificeert rooster-dislocaties en hun Burgers-vector [87]. Gao et al. pasten DXA toe op nanokrabben van ijzer [88]: dislocatienetwerk herorganiseert zich; plastische zone groeit lineair met krabafstand, daarna neemt dislocatiedichtheid snel af door reacties; plastische activiteit concentreert zich aan de krabfront; weinig dislocaties blijven in het midden; vacancies ontstaan door dislocatiereacties; deformatie-twinning treedt op.

Bij polykristallen wordt vervorming gemedieerd door dislocatienucleatie in korrelinterieurs en korrelgrensdislocaties [43]. Bij grotere korrels leidt dislocatiepropagatie tot nanovoids en interstitiële clusters; bij kleinere korrels domineert twin-vorming aan korrelgrenzen. Dit wordt toegeschreven aan overgang van gladde (grote korrel) naar ruwe (kleine korrel) dislocatiebeweging. Li vergelijkde vervormingsmechanismen van enkeltkristallijn Cu, polykristallijn Cu en NT-polykristallijn Cu [25]: enkeltkristallijn Cu wordt gedomineerd door dislocatie-interacties; polykristallijn Cu door dislocaties en korrelgrenzen (GB); NT-polykristallijn Cu door interacties van dislocaties, GB en twin boundaries (TB) met twinning/detwinning.

Veel studies richten zich op halfgeleiders (Si, SiC, GaAs, AlN). Naast dislocatieglijding is fasetransformatie een belangrijke vervormingsmodus [89]. Interactie tussen dislocatie en fasetransformatie varieert met kristallografische oriëntatie. Si (010) ondergaat inelastische vervorming met fasetransformatie Si-I naar Si-III/Si-XII vóór "Pop-In"; dit gebeurt niet in Si (110) en Si (111). Fasetransformatie Si-I naar bct-5 is dominant voor incipient plasticity in alle oriëntaties; dislocatienucleatie is ook actief in Si (010). Dai et al. analyseerden onderlaagschade in enkeltkristallijn Si: kristallijne fase-evolutie is consistent met hydrostatische spanning en temperatuurdistributie [40]. SiC verwijdering gebeurt via fasetransformatie van zincblende naar amorfe structuur met weinig hexagonale diamantstructuren. Hogere krab snelheid genereert meer amorfe atomen, minder hexagonale diamant atomen en minder dislocatie-atomen [90]. Meng et al. toonden aan dat wanneer fasetransformatie niet dominant is, de Schmidt-coëfficiëntmethode effectief glijbeweging van 3C-SiC voorspelt [79]. Verder onderzoek van Meng et al. toont aan dat strain rate en thermische verzachting direct materiaalverwijdering en onderlaagschade (SSD) beïnvloeden. Thermische verzachting onder hoge strain rate overtreft de remming van dislocatiegeneratie. Polykristallijn SiC wordt gedomineerd door amorfe fasetransformatie; hexagonale diamant atomen en dislocaties in GB's. Hogere snelheid en diepere snede bevorderen amorfe structuur. In vergelijking met monokristallijn SiC is polykristallijn SiC zachter (minder normaaalkracht, dunner SSD) [91]. Dislocatiepropagatie en fasetransformatieanalyse verklaren machinismechanismen, maar complexe materialen zijn nog onderbelicht.

Toekomstige onderzoeksrichtingen en uitdagingen

1. Nieuwe TBN-technologieën als AFM-tip nanomilling (uitstekende prestaties, groottecontrole) zijn veelbelovend, maar het mechanisme is onvolledig begrepen door detectielimities. MD kan snijkraftvariatie, defectdynamica en materiaalverwijdering tijdens tiprotatie onthullen. Ook sample-vibratie-ondersteund nanokrabben is nog niet gerapporteerd; MD kan vroege richting geven.

2. Door lengte- en tijdschaallimitaties kan MD experimentele processen nog niet kwantitatief volledig beschrijven. In gecombineerde studies blijft MD vaak kwalitatief [19,20,94]. Voor tijdsafhankelijke materialen (amorfe polymers) beïnvloedt tipsnelheid de verwijderingstoestand sterk. Algoritme- en rekenkrachtverbetering is onmisbaar voor kwantitatieve voorspelling.

3. Veel nanostructuren in polymers zijn gerealiseerd met TBN [95-97]. Mechanische machining met constante normaaalkracht garandeert precisie [98]. MD-simulaties van polymers in TBN zijn echter beperkt. Thermische scanning probelithografie in polymers kan TBN-precisie verbeteren. Het onthullen van het materiaalverwijderingsmechanisme onder mechanisch-thermische effecten is essentieel.

Conclusie

MD-simulaties spelen een onmisbare rol in nanomachining om tot nu toe onbekende fenomeenen te onthullen [99]. Deze review vat de recente vooruitgang samen:

1. Opzet van MD-modellen voor diverse materialen en gerelateerde potentiaalfuncties, met nadruk op NT-polykristallijnen en amorfe polymers. Modelnauwkeurigheid is cruciaal voor simulateresultaten.
2. Nieuwe TBN-technologieën: meervoudig krabben, dynamische ploeging, VANILA. MD onthult essentiële verschillen tussen statische en dynamische ploeging qua interne defecten, morfologie en snijkrachten.
3. Analyse van MD-simulaties: snijkraft, materiaalverwijderingstoestand, defectanalyse. Nieuwe analysemethoden (slip-analyse) zijn geïntroduceerd. Deze methoden tonen de materiaalafhankelijkheid van TBN-machining.

Afkortingen

TBN:

Tip-gebaseerde nanomachining

AFM:

Atoomkrachtmicroscopie

MD:

Moleculaire dynamica

NT:

Nanotwined

PE:

Polyethyleen

EAM:

Embedded atom method

SW:

Stillinger-Weber

ABOP:

Analytical bond order potential

REBO:

Reactive Empirical Bond Order

AIREBO:

Adaptive Intermolecular Reactive Empirical Bond Order

LAMMPS:

Large-scale Atomic/Molecular Massively Parallel Simulator

MS:

Material Studio

OVITO:

Open Visualization Tool

VMD:

Visual molecular dynamics

VANILA:

Vibration-Assisted Nano Impact machining by Loose Abrasives

SEM:

Scanning elektronenmicroscopie

CNA:

Common neighbor analysis

CSP:

Centro-symmetry parameter

DXA:

Dislocation extraction algorithm

GB:

Korrelgrens

TB:

Twin boundary

SSD:

Onderlaagschade