Wafer-to-Wafer Bonding: Part: I -- The Coupled Physics Problem and the 2D Finite Element Implementation

Dit artikel presenteert een wiskundig consistent gereduceerd model voor wafer-tot-wafer-bonding dat de gekoppelde interactie tussen plaatvervorming en luchtfilmdruk beschrijft, en dit numeriek oplost met FEniCSx om de niet-lineaire dynamica van het bondingsproces te analyseren en te optimaliseren.

De kern

De Dans van de Wafers: Hoe Chips Dichtbij elkaar worden Gekleefd

Stel je voor dat je twee enorme, perfect ronde glazen schijven hebt (deze worden "wafers" genoemd in de chipwereld). Je wilt ze zo nauwkeurig tegen elkaar drukken dat ze aan elkaar plakken, zonder dat er ook maar een klein luchtbelletje tussen blijft zitten. Dit proces heet Wafer-to-Wafer Bonding en het is de sleutel tot het maken van superkrachtige, 3D-geprinte computerchips.

Maar hier is het probleem: als je de bovenste schijf neerlaat, zit er een heel dun laagje lucht tussen de twee schijven. Die lucht wil niet weg. Het is alsof je twee natte glazen schijven op elkaar legt; de lucht zit gevangen en duwt de schijven uit elkaar.

Dit artikel beschrijft hoe de auteurs een slim computerprogramma hebben gebouwd om precies te voorspellen wat er gebeurt tijdens dit "plakproces". Ze noemen het een FSI-probleem (Fluïdum-Structuur Interactie), maar laten we het simpel houden: het is een gevecht tussen buigkracht en luchtdruk.

1. Het Gevecht: Buigen vs. Luchtdruk

Stel je de bovenste wafers voor als een heel dun, flexibel trampoline.

• De Buigkracht: Een kleine stamper duwt in het midden van de trampoline. De trampoline buigt naar beneden.
• De Luchtdruk: Als de trampoline naar beneden buigt, wordt de ruimte eronder kleiner. De lucht die daar zit, wordt samengeperst. Die lucht probeert eruit te ontsnappen, maar omdat de opening zo klein is, bouwt er een enorme druk op. Deze druk duwt de trampoline weer omhoog, tegen de stamper in.

Het is een ingewikkeld dansje: de trampoline buigt, de lucht duwt, de trampoline buigt weer, en de lucht moet weer weg. Als je dit niet goed berekent, blijven er luchtbellen achter (wat slecht is voor de chip) of duurt het te lang.

2. De Oplossing: Een Slimme Wiskundige "Korte Weg"

De auteurs hebben een wiskundig model bedacht dat dit hele proces simuleert. In plaats van elke moleculaire beweging van de lucht en het glas te berekenen (wat duizenden jaren zou duren op een computer), hebben ze een reductiemodel gebruikt.

• De Plaat: Ze behandelen de wafers als een klassiek "Kirchhoff-Love plaatje". Denk hierbij niet aan een zware betonnen vloer, maar aan een dun vel papier dat je buigt. Ze hebben bewezen dat je de complexe 3D-wiskunde kunt vervangen door een simpelere 2D-versie die net zo goed werkt voor deze dunne schijven.
• De Lucht: Voor de lucht gebruiken ze een vergelijking die bekendstaat als de Reynolds-vergelijking. Dit is eigenlijk de wiskunde voor smeermiddelen (zoals olie in een motor), maar dan toegepast op lucht. Het beschrijft hoe lucht stroomt door een steeds smaller wordend kanaal.

3. De Simulatie: Een Digitale Danszaal

Ze hebben dit model ingebouwd in een krachtig computerprogramma genaamd FEniCSx.

• Ze hebben de ronde wafers opgedeeld in duizenden kleine driehoekjes (een "mesh"), net als een puzzel.
• Het programma rekent in elke stap van de simulatie uit: "Hoeveel buigt de plaat?" en "Hoeveel druk heeft de lucht?".
• Omdat de lucht de plaat duwt en de plaat de lucht verplaatst, moeten ze dit gelijktijdig oplossen. Het is alsof je twee mensen probeert te laten dansen waarbij ze elkaar voortdurend duwen en trekken; je moet hun bewegingen tegelijk berekenen.

4. De Verassende Resultaten: Wat Ze Ontdekten

De simulatie leverde enkele verrassende inzichten op die tegen onze intuïtie ingaan:

• Hoe groter de opening, hoe sneller het plakt!
  Je zou denken: "Als de wafers dichter bij elkaar beginnen, plakt het sneller." Maar nee! Als ze heel dicht bij elkaar beginnen (bijvoorbeeld 30 micrometer), zit er al veel lucht opgesloten die moeilijk weg kan. De luchtweerstand is zo groot dat de plaat niet verder kan buigen.
  Als ze verder uit elkaar beginnen (bijvoorbeeld 100 micrometer), is er meer ruimte voor de lucht om te ontsnappen. De plaat kan makkelijker buigen en de "plakrand" (waar de wafers elkaar raken) snelt sneller naar buiten.

• De Lucht is de "Kleefkracht":
  In het gebied waar de wafers al aan elkaar plakken, is er geen lucht meer. Maar de wiskunde toont aan dat de luchtdruk in het niet-geplakte gebied precies de kracht levert die nodig is om de wafers tegen elkaar te drukken. Het is alsof de lucht zelf de "hand" is die de wafers samendrukt.

• Viskeuze Lucht (Dikke Lucht) vertraagt:
  Als je de lucht "dikker" zou maken (hogere viscositeit), duurt het plakproces langer. Het is alsof je probeert een deksel van een pot met honing te draaien in plaats van van een pot met water.

Waarom is dit belangrijk?

Voor chipfabrikanten is dit goud waard. Ze kunnen nu precies voorspellen hoe ze hun machines moeten instellen. Ze weten nu dat ze niet hoeven te proberen de wafers zo dicht mogelijk bij elkaar te starten, maar dat ze de initiële afstand en de snelheid van het proces moeten optimaliseren om luchtbellen te voorkomen.

Kort samengevat:
De auteurs hebben een digitale simulator gebouwd die de dans tussen een buigende glazen plaat en een gevangen luchtlaag perfect nabootst. Ze hebben ontdekt dat "ruimte geven" (een grotere startafstand) paradoxaal genoeg zorgt voor een snellere en betere binding, en dat de luchtweerstand de belangrijkste regelaar is in dit proces. Dit helpt bij het maken van de snellere, kleinere en krachtigere computers van de toekomst.

Technische samenvatting

Technische Samenvatting: Wafer-tot-Wafer Bonding

1. Het Probleem
Wafer-tot-wafer (WxW) bonding is een cruciale technologie voor geavanceerde 3D-integratie, zoals hybride bonding voor Cu-Cu interconnects met een fijne pitch (<10 µm). Het proces wordt echter gehinderd door complexe, sterk gekoppelde fysica. Tijdens het bonden vervormen de wafers en wordt lucht ingesloten tussen de oppervlakken. Dit creëert een Fluid-Structure Interaction (FSI) probleem waarbij:

• De elastische vervorming van de wafer de geometrie van het spleetje (gap) bepaalt.
• De druk van de ingesloten lucht (Reynolds-druk) als een weerstand werkt die de verdere bonding vertraagt.
• Dit leidt tot niet-intuïtieve dynamica, zoals een niet-monotoon gedrag waarbij een grotere initiële spleet soms sneller bondt dan een kleinere spleet, wat tegen de intuïtie ingaat.

Bestaande modellen missen vaak wiskundige consistentie, behandelen de niet-lineariteiten onvoldoende en bieden geen gedetailleerde implementatie voor reproduceerbaarheid.

2. Methodologie
De auteurs ontwikkelen een wiskundig consistent reduced-order model dat de volledige 3D-fysica reduceert tot een 2D-probleem, geïmplementeerd in het high-performance FEniCSx framework.

• Fysische Formulering:

  • Structuur: De 3D lineaire elasticiteit wordt systematisch gereduceerd tot de klassieke Kirchhoff-Love plaatvergelijking door integratie over de dikte. Dit beschrijft de buiging van de wafer ($w$).
  • Vloeistof: De ingesloten lucht wordt gemodelleerd met de Reynolds-smeermiddelvergelijking, die de drukverdeling ($p$) in het spleetje koppelt aan de verandering van de spleethoogte.
  • Koppeling: Het systeem is tweeweg-gekoppeld (non-lineair): de verplaatsing $w$ bepaalt de spleethoogte $h$ in de Reynolds-vergelijking, en de luchtdruk $p$ werkt als een verdeelde transversale last op de plaatvergelijking.

• Numerieke Implementatie:

  • Discretisatie: De plaatvergelijking (vierde orde PDE) wordt opgelost met een $C_0$-interne-straf (C0IP) formulering, wat het gebruik van standaard continue Galerkin-elementen mogelijk maakt zonder de complexiteit van $C_1$-continuïteit. De drukveld wordt gemodelleerd met standaard continue Galerkin-elementen.
  • Oplossing: Het gekoppelde systeem wordt monolithisch (gelijktijdig) opgelost met een Newton-oplosser en automatische differentiatie (via UFL in FEniCSx).
  • Randvoorwaarden: Er worden contactkrachten gemodelleerd via een strafmethode (penalty method) voor de bovenste chuck, terwijl de Reynolds-druk in het gebonden gebied fungeert als een effectieve contactreactie om doordringing te voorkomen.

3. Belangrijkste Bijdragen

• Wiskundige Rigor: Voor het eerst wordt een WxW-bondmodel strikt afgeleid van de 3D lineaire elasticiteit naar een 2D plaat-Reynolds-systeem, waarbij alle aannames expliciet worden gemaakt.
• Monolithische FSI-oplossing: Een robuuste implementatie van het sterk gekoppelde, niet-lineaire FSI-probleem in FEniCSx, wat reproduceerbaarheid garandeert.
• Ontmaskering van Niet-Lineariteiten: Het model onthult en analyseert de oorzaken van tegen-intuïtief gedrag, zoals de niet-monotone relatie tussen de initiële spleet en de bonding-snelheid.
• Validatie: Het model is gevalideerd tegen experimentele data (verplaatsingsgeschiedenis) en toont een uitstekende overeenkomst.

4. Resultaten
De simulaties leveren de volgende inzichten op:

• Validatie: De gesimuleerde verplaatsing van de wafer komt goed overeen met experimentele metingen voor verschillende initiële spleten (30, 70 en 100 µm).
• Krachtenbalans: De Reynolds-druk fungeert effectief als een contactkracht in het reeds gebonden gebied; de gemiddelde druk in het gebonden gebied balanceert de aangebrachte kracht van de striker.
• Invloed van de Initiële Spleet: In tegenstelling tot een zuiver buigingsprobleem (zonder lucht), verloopt het bonden sneller bij een grotere initiële spleet dan bij een kleinere. Een kleinere spleet genereert hogere luchtdruk, wat de voortplanting van de bonding-front meer weerstand biedt.
• Viscositeit en Interfaciale Energie:
  • Lagere luchtviscositeit leidt tot een snellere bonding-front.
  • Hogere interfaciale energie (oppervlakte-energie) versnelt de bonding-front en vergroot het gebonden oppervlak.
• Niet-monotoon Gedrag: De snelheid van de bonding-front is niet-lineair afhankelijk van de radius en de tijd, gedreven door de dynamiek van de Reynolds-vergelijking.

5. Betekenis en Toekomstperspectief
Dit werk biedt een voorspellend en computerefficiënt instrument voor het optimaliseren van WxW-bondprocessen. Het helpt bij het begrijpen van faalmechanismen (zoals holtes door deeltjes) en biedt richtlijnen voor procesoptimalisatie.

Toekomstige uitbreidingen omvatten:

• Integratie van anisotrope buigstijfheid voor kristallijne wafers en gestapelde lagen.
• Uitbreiding naar een shell-model om ook in-vlakke vervormingen (rek, schuif) en residualle vervormingen (bow/warp) na het bonden te voorspellen.
• Modelkalibratie en onzekerheidskwantificatie met behulp van uitgebreide metrologie-data.

Kortom, dit artikel legt de basis voor een nieuwe generatie van fysica-gedreven simulaties die essentieel zijn voor de ontwikkeling van geavanceerde halfgeleiderverpakkingen.