pith. sign in

arxiv: 2603.22827 · v4 · submitted 2026-03-24 · ⚛️ physics.comp-ph · physics.app-ph· physics.flu-dyn

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

Kamalendu Ghosh , Bhavesh Shrimali , Subin Jeong This is my paper

Pith reviewed 2026-05-15 00:59 UTC · model grok-4.3

classification ⚛️ physics.comp-ph physics.app-phphysics.flu-dyn
keywords wafer-to-wafer bondingfluid-structure interactionKirchhoff-Love plateReynolds lubricationfinite element simulationbonding dynamicswafer deformation
0
0 comments X

The pith

A coupled plate-Reynolds model reproduces experimental wafer-bonding displacement histories across initial gaps.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives a Kirchhoff-Love plate equation for thin-wafer bending from three-dimensional linear elasticity and couples it to a Reynolds lubrication equation for the thin air film trapped between wafers. The resulting nonlinear system is discretized monolithically with interior-penalty finite elements and solved in FEniCSx using implicit time stepping and Newton iteration. Simulations recover measured probe-displacement curves for several initial gaps and confirm that Reynolds pressure supplies the effective contact reaction at the advancing bond front. Parametric sweeps expose nonlinear and sometimes non-monotonic dependence of front speed on gap size, air viscosity, and interfacial energy.

Core claim

The nonlinear plate-Reynolds system, derived under small-deformation and thin-film assumptions, accurately captures the time-dependent fluid-structure interaction during wafer-to-wafer bonding. When solved monolithically with C0 interior-penalty elements for the fourth-order plate operator and continuous Galerkin elements for pressure, the model reproduces experimentally reported probe-displacement histories for multiple initial gaps and verifies force equilibrium at the bond front, where Reynolds pressure acts as the effective reaction force.

What carries the argument

The nonlinear plate-Reynolds system, which enforces fourth-order bending equilibrium on the wafers while solving second-order lubrication flow for the inter-wafer air film and couples them through pressure and deflection at each time step.

Load-bearing premise

Wafers remain thin and deformations stay small enough for Kirchhoff-Love plate theory, while the air film stays thin with no-slip walls and negligible inertia.

What would settle it

A direct mismatch between simulated and measured probe-displacement histories for an untested combination of initial gap and air viscosity, or a measured force imbalance at the bond front that exceeds the Reynolds pressure prediction.

Figures

Figures reproduced from arXiv: 2603.22827 by Bhavesh Shrimali, Kamalendu Ghosh, Subin Jeong.

Figure 1
Figure 1. Figure 1: (a) Cross-sectional view of the initial domain configuration showing the top wafer [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Quarter-domain ωq exploiting four-fold symmetry of the circular wafer mid-surface ω. (a) Red boundaries: dashed lines denote symmetry edges Γx, Γy (zero normal slope); solid arc denotes the outer boundary ΓR (Dirichlet conditions). The striker region ωS is shown near the origin. (b) The corresponding finite-element mesh (42,449 triangular elements) used in the DOLFINx computation. 3.2. Time discretization … view at source ↗
Figure 3
Figure 3. Figure 3: Top wafer displacement w(r, t) versus time at a fixed radial location r = mid-region for h0 = 30, 70, and 100 µm. Solid lines (orange) show the finite-element solution; dashed lines (blue) show the experimental measurement. 14 [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Snapshots of the displacement field w at nine time instances during bonding (h = 70 µm gap). The out-of￾plane displacement is warped by a factor of 500 to enhance the visualisation. The bottom wafer shape is represented by the gray shape, and the colormap indicates the out-of-plane displacement of the top wafer. 15 [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Displacement w(r) (left) and pressure p(r) (right) at the final time step (t = 50 s) for a 70 µm gap. The average pressure p¯ computed by integrating the Reynolds pressure over the bonded area (r ≤ Rpin) equals the applied pressure |pe|, verifying force equilibrium across the bonding front. This confirms that the Reynolds pressure field correctly captures the contact mechanics at the bond front. 4.3. Effec… view at source ↗
Figure 6
Figure 6. Figure 6: Bond speed diagnostics for three bond gaps: (a) bonded radius versus time, (b) bond speed versus time, [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Effect of air viscosity on: (a) bond speed over time, (b) bond speed versus radius, and (c) temporal growth [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Interfacial energy sensitivity diagnostics: (a) bond speed over time, (b) bond speed versus radius, and (c) [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Spatio-temporal evolution of the radial displacement field [PITH_FULL_IMAGE:figures/full_fig_p019_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Air-pressure diagnostics for h0 = 30, 70, and 100 µm: (a) radial profiles pair(r) in the unbonded region, r > rb(t), evaluated at the common time t = 9 s; and (b) radial profiles pair(r, tc) evaluated at the time tc when the bond front reaches rb(tc) = 65 mm. the monolithic integration of Si-CMOS and III–V-on-Si wafers. Journal of Semiconductors, 42 (2):023106, 2021. doi: 10.1088/1674-4926/42/2/023106. [2… view at source ↗
read the original abstract

Wafer-to-wafer (WxW) bonding is a key enabler for three-dimensional integration, including hybrid bonding for fine-pitch Cu-Cu interconnects. During bonding, wafer deformation and the air entrapped between the wafers interact through a strongly coupled, time-dependent fluid-structure interaction (FSI) that can produce non-intuitive bonding dynamics and process sensitivities. This paper develops a mathematically consistent reduced-order model for WxW bonding by deriving a Kirchhoff-Love plate equation for wafer bending from three-dimensional linear elasticity and coupling it to a Reynolds lubrication equation for the inter-wafer air film. The resulting nonlinear plate-Reynolds system is discretized and solved monolithically in the high-performance FEniCSx framework using a $C^0$ interior-penalty formulation for the fourth-order plate operator, standard continuous Galerkin discretization for the pressure field, implicit time integration, and a Newton solver with automatic differentiation. Simulations reproduce experimentally reported probe-displacement histories for multiple initial gaps and verify force equilibrium at the bond front, where the Reynolds pressure acts as an effective contact reaction. Parametric studies reveal nonlinear, and in some cases non-monotonic, sensitivities of bonding-front kinetics to the initial gap, air viscosity, and interfacial energy, providing actionable trends for process optimization.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The paper develops a reduced-order model for wafer-to-wafer bonding by deriving the Kirchhoff-Love plate equation from 3D linear elasticity and coupling it to the Reynolds lubrication equation for the entrapped air film. The resulting nonlinear FSI system is discretized monolithically in FEniCSx with a C0 interior-penalty scheme for the fourth-order plate operator, continuous Galerkin elements for pressure, implicit time stepping, and a Newton solver with automatic differentiation. Simulations are shown to reproduce experimental probe-displacement histories across multiple initial gaps and to satisfy force equilibrium at the bond front, where Reynolds pressure acts as an effective contact reaction. Parametric studies explore sensitivities to initial gap, viscosity, and interfacial energy.

Significance. If the reproduction of experimental displacement histories holds under the stated assumptions, the work provides a computationally tractable, parameter-free predictive tool for bonding-front kinetics that can inform process optimization in 3D integration. The monolithic discretization and verification of internal force balance are strengths that support the claim of mathematical consistency without ad-hoc contact terms.

minor comments (3)
  1. §2.1: clarify the precise range of wafer thickness-to-radius ratios for which the Kirchhoff-Love reduction remains accurate, and state the corresponding error estimate relative to 3D elasticity.
  2. Figure 4: the legend for the multiple initial-gap curves is difficult to read at the printed size; increase font size or add a table of gap values.
  3. §4.3: the statement that Reynolds pressure 'acts as an effective contact reaction' would benefit from an explicit plot or table showing the integrated pressure force equaling the applied load at the bond front for at least two time instants.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary and recommendation of minor revision. The assessment accurately reflects the derivation of the coupled Kirchhoff-Love/Reynolds model, its monolithic discretization in FEniCSx, and the verification against experimental probe-displacement data. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation starts from the standard reduction of 3D linear elasticity to the Kirchhoff-Love plate equation under the explicit thin-wafer and small-deformation assumptions, then couples it to the Reynolds lubrication equation under the standard thin-film, no-slip, and negligible-inertia limits. These are textbook continuum reductions that do not reference the target results or any fitted quantities from the present simulations. The monolithic C0-IP + CG discretization, implicit time stepping, and Newton solver are conventional numerical methods for the resulting system. The reported reproduction of experimental probe-displacement histories across multiple initial gaps and the internal force-equilibrium check at the bond front follow directly from solving the coupled PDE system; they are not obtained by fitting parameters to the same data or by any self-citation chain that collapses the claim to its inputs. The paper is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard assumptions from linear elasticity and lubrication theory with no new free parameters or invented entities introduced beyond physical inputs such as viscosity and interfacial energy.

axioms (2)
  • domain assumption Kirchhoff-Love plate theory derived from 3D linear elasticity
    Invoked to obtain the fourth-order bending equation for thin wafers under small deflections.
  • domain assumption Reynolds lubrication equation for thin air film
    Used to model viscous flow in the inter-wafer gap with no-slip boundaries.

pith-pipeline@v0.9.0 · 5547 in / 1358 out tokens · 46145 ms · 2026-05-15T00:59:31.135661+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages

  1. [1]

    Shuyu Bao, Yue Wang, Khaw Lina, Li Zhang, Bing Wang, Wardhana Aji Sasangka, Ken- neth Eng Kian Lee, Soo Jin Chua, Jurgen Michel, Eugene Fitzgerald, Chuan Seng Tan, and Kwang Hong Lee. A review of silicon-based wafer bonding processes, an approach to realize 19 0 20 40 60 80 100 120 140 r (mm) 0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 pair(r) (...

    work page doi:10.1088/1674-4926/42/2/023106 2021

  2. [2]

    Moriceau, F

    H. Moriceau, F. Rieutord, F. Fournel, Y. Le Tiec, L. Di Cioccio, C. Morales, A. M. Charvet, and C. Deguet. Physics of direct bonding: Applications to 3D heterogeneous or homogeneous integration.Microelectronics Reliability, 50(8):1154–1162, 2010. doi: 10.1016/j.microrel.2010. 01.011

    work page doi:10.1016/j.microrel.2010 2010

  3. [3]

    Rieutord, B

    F. Rieutord, B. Bataillou, and H. Moriceau. Dynamics of a bonding front.Physical Review Letters, 94:236101, 2005. doi: 10.1103/PhysRevLett.94.236101

    work page doi:10.1103/physrevlett.94.236101 2005

  4. [4]

    Navarro et al

    E. Navarro et al. Direct silicon bonding dynamics: a coupled fluid/structure analysis.Applied Physics Letters, 103:034104, 2013. doi: 10.1063/1.4813843

    work page doi:10.1063/1.4813843 2013

  5. [5]

    Radisson, F

    D. Radisson, F. Fournel, and E. Charlaix. Modelling of the direct bonding wave.Microsystem Technologies, 21:969–971, 2015. doi: 10.1007/s00542-014-2101-y

    work page doi:10.1007/s00542-014-2101-y 2015

  6. [6]

    Multi- physics simulation of wafer-to-wafer bonding dynamics

    Nathan Ip, Nima Nejadsadeghi, Carlos Fonseca, Norifumi Kohama, and Kimio Motoda. Multi- physics simulation of wafer-to-wafer bonding dynamics. In2022 IEEE 72nd Electronic Com- ponents and Technology Conference (ECTC), pages 502–506. IEEE, 2022

    work page 2022

  7. [7]

    Designandsimulationofsymmetricwafer-to-waferbondingcompensating agravityeffect

    KyeongbinLimetal. Designandsimulationofsymmetricwafer-to-waferbondingcompensating agravityeffect. In2020 IEEE 70th Electronic Components and Technology Conference (ECTC), pages 1480–1485. IEEE, 2020. doi: 10.1109/ECTC32862.2020.00233

    work page doi:10.1109/ectc32862.2020.00233 2020

  8. [8]

    Numerical and mechanical analysis of direct wafer bonding considering non-uniform impurity particle distri- butions.Microsystems & Nanoengineering, 11:151, 2025

    Feixiang Tang, Siyu He, Yuhan Li, Wenjin Liu, Fang Dong, and Sheng Liu. Numerical and mechanical analysis of direct wafer bonding considering non-uniform impurity particle distri- butions.Microsystems & Nanoengineering, 11:151, 2025. doi: 10.1038/s41378-025-00994-4

    work page doi:10.1038/s41378-025-00994-4 2025

  9. [9]

    An experimental and numerical study on glass frit wafer-to-wafer bonding.Micromachines, 14(1): 165, 2023

    Seyed Amir Fouad Farshchi Yazdi, Matteo Garavaglia, Aldo Ghisi, and Alberto Corigliano. An experimental and numerical study on glass frit wafer-to-wafer bonding.Micromachines, 14(1): 165, 2023

    work page 2023

  10. [10]

    Simulation of Cu pad 20 expansion in wafer-to-wafer Cu/SiCN hybrid bonding.Microelectronics Reliability, 138:114716, 2022

    Yan Wen Tsau, Joke De Messemaeker, Abdellah Salahouelhadj, Mario Gonzalez, Liesbeth Witters, Boyao Zhang, Marc Seefeldt, Eric Beyne, and Ingrid De Wolf. Simulation of Cu pad 20 expansion in wafer-to-wafer Cu/SiCN hybrid bonding.Microelectronics Reliability, 138:114716, 2022

    work page 2022

  11. [11]

    and Dean, Joseph P

    Igor A Baratta, Joseph P Dean, Jørgen S Dokken, Michal Habera, Jack Hale, Chris N Richardson, Marie E Rognes, Matthew W Scroggs, Nathan Sime, and Garth N Wells. DOLFINx: the next generation FEniCS problem solving environment.Zenodo, 2023. doi: 10.5281/zenodo.10447666. URLhttps://doi.org/10.5281/zenodo.10447666. Preprint

    work page doi:10.5281/zenodo.10447666 2023

  12. [12]

    A hybrid finite element implementation of two- potential constitutive model of dielectric elastomers.Finite Elements in Analysis and Design, 247:104348, 2025

    Kamalendu Ghosh and Bhavesh Shrimali. A hybrid finite element implementation of two- potential constitutive model of dielectric elastomers.Finite Elements in Analysis and Design, 247:104348, 2025

    work page 2025

  13. [13]

    The remarkable bending properties of perforated plates.Journal of the Mechanics and Physics of Solids, 154:104514, 2021

    Bhavesh Shrimali, Matteo Pezzulla, Samuel Poincloux, Pedro M Reis, and Oscar Lopez-Pamies. The remarkable bending properties of perforated plates.Journal of the Mechanics and Physics of Solids, 154:104514, 2021

    work page 2021

  14. [14]

    Constructionof arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes

    MatthewWScroggs, JørgenSDokken, ChrisNRichardson, andGarthNWells. Constructionof arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes. ACM Transactions on Mathematical Software, 48(2):1–23, 2022

    work page 2022

  15. [15]

    Basix: a runtime finite element basis evaluation library.Journal of Open Source Software, 7(73):3982, 2022

    Matthew W Scroggs, Igor A Baratta, Chris N Richardson, and Garth N Wells. Basix: a runtime finite element basis evaluation library.Journal of Open Source Software, 7(73):3982, 2022

    work page 2022

  16. [16]

    Unified form language: A domain-specific language for weak formulations of partial differential equations.ACM Transactions on Mathematical Software, 40(2):1–37, 2014

    Martin S Alnæs, Anders Logg, Kristian B Ølgaard, Marie E Rognes, and Garth N Wells. Unified form language: A domain-specific language for weak formulations of partial differential equations.ACM Transactions on Mathematical Software, 40(2):1–37, 2014

    work page 2014

  17. [17]

    Kamalendu Ghosh, Bhavesh Shrimali, Aditya Kumar, and Oscar Lopez-Pamies. The nonlinear viscoelastic response of suspensions of rigid inclusions in rubber: I—Gaussian rubber with constant viscosity.Journal of the Mechanics and Physics of Solids, 154:104544, 2021

    work page 2021

  18. [18]

    The nonlinear viscoelastic re- sponseofsuspensionsofvacuousbubblesinrubber: I—Gaussianrubberwithconstantviscosity

    Bhavesh Shrimali, Kamalendu Ghosh, and Oscar Lopez-Pamies. The nonlinear viscoelastic re- sponseofsuspensionsofvacuousbubblesinrubber: I—Gaussianrubberwithconstantviscosity. Journal of Elasticity, 153(4):479–508, 2023

    work page 2023

  19. [19]

    On the two-potential constitutive modeling of dielectric elastomers.Meccanica, 56(6):1505–1521, 2021

    Kamalendu Ghosh and Oscar Lopez-Pamies. On the two-potential constitutive modeling of dielectric elastomers.Meccanica, 56(6):1505–1521, 2021

    work page 2021

  20. [20]

    On Nitsche’s method for elastic contact problems.SIAM Journal on Scientific Computing, 42(2):B425–B446, 2020

    Tom Gustafsson, Rolf Stenberg, and Juha Videman. On Nitsche’s method for elastic contact problems.SIAM Journal on Scientific Computing, 42(2):B425–B446, 2020. 21

    work page 2020