Atomistic-Continuum Coupling by Homogenization

Technische Samenvatting: Koppeling van atomaire en continuumbenaderingen via Homogenisatie

Probleemstelling
Klassieke atomaire simulaties gebaseerd op interatomaire potentialen bieden een hoge-resolutie weergave van roosterinstabiliteiten, defectnucleatie en microstructuur-evolutie. Hun toegankelijke systeemgroottes zijn echter beperkt tot het nanoschaal, waardoor ze ongeschikt zijn voor de analyse van structurele componenten op micrometerschaal en hun randvoorwaarden. Hoewel gelijktijdige multischaalmethoden (bijvoorbeeld quasicontinuum) atomaire en continuumbereiken ruimtelijk koppelen, kampen ze vaak met algoritmische uitdagingen zoals spookkrachten aan de interface en moeilijkheden bij het handhaven van een globale energieconsistentie. Hiërarchische methoden, die vertrouwen op homogenisatie, vermijden directe ruimtelijke interfaces, maar vereisen traditioneel een continuüm constitutieve wet op microschaal. De uitdaging die in dit werk wordt aangepakt, is het tot stand brengen van een hiërarchische koppeling waarbij de microschaal strikt wordt gedomineerd door discrete atomaire fysica (interatomaire potentialen) in plaats van een fenomenologische continuümwet, terwijl energieconsistentie en robuustheid voor niet-lineaire structurele analyse worden behouden.

Methodologie
De auteurs stellen een tweeschaals atomaire–continuüm raamwerk voor gebaseerd op computationele homogenisatie en het Hill–Mandel energie-equivalentieprincipe. De benadering koppelt een niet-lineair eindige-elementen (EE) randwaardeprobleem op microschaal (structureel niveau) met periodieke moleculaire-statische (MS) celproblemen op de kwadratuurpunten van het EE-net.

• Schaalovergang: De continuüm vervormingsgradiënt ($F$) drijft het atomaire representatieve volume-element (RVE). Periodieke randvoorwaarden (PBC) worden toegepast op de atomaire cel. De atomaire solver voert een energie-minimalisatie bij nul temperatuur uit om de evenwichts-atomaire configuratie te vinden voor de voorgeschreven vervorming.
• Overdracht van Spanning en Raaklijn: In tegenstelling tot standaard EE$^2$-methoden die vertrouwen op een opgeslagen energiefunctie, geeft de atomaire cel volumegemiddelde Cauchy-spanningen terug die zijn afgeleid uit het viriaaltheorema. Deze worden getransformeerd in werkgeconjugeerde eerste Piola–Kirchhoff-spanningen ($P$) en tweede Piola–Kirchhoff-spanningen ($S$) die vereist zijn door de Lagrangiaanse EE-formulering.
• Algoritmische Raaklijn: Om robuuste Newton-iteraties op de grove schaal te waarborgen, wordt een effectieve raaklijnmodulus berekend. Dit wordt bereikt via numerieke differentiatie (voorwaartse differentie) van de spanningsrespons met betrekking tot de vervormingsgradiënt. Dit houdt in dat de vervormingsgradiënt wordt verstoord, onafhankelijke MS-relaxaties worden uitgevoerd voor elke verstoring, en de resulterende spanningsincrementen worden berekend.
• Implementatie: Het raamwerk integreert LAMMPS voor de MS-solver binnen een schaalbaar MPI-gebaseerd workflow. Parallelisme wordt benut over kwadratuurpunten, waarbij elk Gauss-punt een onafhankelijke set atomaire simulaties triggert (één voor de basistoestand en zes voor de raaklijnverstoringen).

Belangrijkste Bijdragen

1. Energie-consistent Raamwerk: De auteurs vestigen een overgang van nano- naar microschaal die strikt gebaseerd is op het Hill–Mandel-postulaat, waardoor energie-equivalentie wordt gewaarborgd tussen het macroscopische virtuele vermogen en het volumegemiddelde microscopische virtuele vermogen.
2. Interface van Discreet naar Continuüm: Een specifieke interface wordt gedefinieerd waarbij de continuümvervorming een periodiek atomaire RVE aandrijft, en de atomaire solver volumegemiddelde spanningen in werkgeconjugeerde vorm teruggeeft, samen met effectieve raaklijnmoduli.
3. Robuustheid in Niet-Gladde Regimes: Het werk toont aan dat een statisch-impliciete Newton-solver robuust kan blijven, zelfs wanneer de atomaire respons niet-glad is als gevolg van defectnucleatie en instabiliteit. De methode herstelt succesvol een bijna-kwadratische convergentie in de laatste iteraties van kritieke laststappen.
4. Praktische Implementatie: De studie schetst een strategie voor het integreren van gevestigde simulatieblokken (LAMMPS) in een schaalbaar workflow, waarbij de rekenkosten van meerdere gelijktijdige atomaire oplossingen per macroscopische iteratie worden beheerd.

Numerieke Resultaten
Het raamwerk werd getest op enkelkristallijn koper (EAM-potentiaal) in verschillende scenario's:

• Uniaxiale Trek en Druk: Het model ving een uitgesproken asymmetrie tussen trek en druk op. Bij trek leidde instabiliteit tot de vorming van stapelfouten en door tweeling gerelateerde lamellen. Bij druk evolueerde de microstructuur naar brede, gelaagde bundels van stapelfouten en hexagonaal dichtgestapelde (hcp) gebieden. De resultaten toonden aan dat het begin van instabiliteit afhankelijk is van de grootte van de simulatiebox, terwijl de elastische respons voorafgaand aan instabiliteit grootte-onafhankelijk is.
• Cyclische Belasting: Onder omgekeerde cyclische belasting vertoonde het systeem een overgang van een initiële elastische regime met hoge kracht naar een toestand die wordt gedomineerd door een "getrainde" metastabiele defectarchitectuur. Dit resulteerde in brede hysterese-lussen en snelle cyclische stabilisatie na het eerste instabiliteitsgebeuren.
• Kragbalkbuiging: Een balk op micrometerschaal (10 µm lengte) werd gesimuleerd. Het raamwerk loste een gelokaliseerde elastisch-plastische overgang op waarbij defectnucleatie abrupt plaatsvond in nanodomains nabij het ingeklemde uiteinde (onderste vezels), terwijl de rest van de balk elastisch bleef. De Newton-solver verwerkte de scherpe stijfheidsveranderingen die gepaard gaan met deze gelokaliseerde instabiliteiten effectief.

Betekenis en Beweringen
Het artikel beweert de eerste atomaire-naar-continuüm opschaling te bieden die tegelijkertijd: (i) gebaseerd is op het strikte energieprincipe van Hill–Mandel, (ii) compatibel is met een statisch-impliciete eindige-elementenformulering, en (iii) is uitgerust met een effectieve raaklijn die geschikt is voor Newton-solvers in sterk niet-lineaire regimes.

De auteurs benadrukken dat hun benadering potential-gebaseerde atomaire modellering uitbreidt naar structurele lengteschalen (micrometers) die ontoegankelijk zijn voor directe atomaire simulatie in een quasi-statische, athermale setting. Een kritieke bevinding is dat de Newton-solver op grove schaal robuust blijft en een bijna-kwadratische convergentie herstelt, zelfs wanneer de respons op fijne schaal abrupte stijfheidsveranderingen ondergaat als gevolg van defectnucleatie. Dit suggereert dat energie-consistente spannings- en raaklijninformatie uit atomaire cellen nuttige hulpmiddelen blijven voor structurele analyse, zelfs in sterk niet-lineaire, instabiele scenario's.

Het werk is expliciet beperkt tot quasi-statische, athermale responsen bij nul temperatuur. De auteurs merken op dat toekomstige uitbreidingen modellen bij eindige temperaturen, alternatieve randvoorwaarden en de integratie van interatomaire potentialen op basis van machine learning kunnen omvatten.