《Results in Engineering》：Multi-Objective Density Optimization of Cu-Cu Bonded Stacked Dies under Thermal Cycling via Desirability Function Analysis
编辑推荐：
互连致密化对铜-铜堆叠芯片并非单调有益；翘曲抑制仅在中密度区间有所改善；更高密度增加了局部约束和不可逆循环行为；响应从弹性可恢复性向耗散主导转变；通过期望函数法（Desirability Functional Approach），密度64被确定为最稳健的可靠性
互连致密化对铜-铜堆叠芯片并非单调有益；翘曲抑制仅在中密度区间有所改善；更高密度增加了局部约束和不可逆循环行为；响应从弹性可恢复性向耗散主导转变；通过期望函数法（Desirability Functional Approach），密度64被确定为最稳健的可靠性折衷方案。
论文解读：热循环下铜-铜键合堆叠芯片的多目标密度优化
一、研究背景与问题提出
随着人工智能计算、数据中心架构和异构系统集成的快速发展，三维集成电路（3D IC）日益依赖直接的铜-铜（Cu-Cu）键合技术。这种技术能够实现亚微米节距、更低的寄生效应和更紧密的电热集成，成为下一代高带宽存储器、CMOS图像传感器和高性能计算平台的核心。然而，电学与集成优势伴随着严苛的热机械环境。反复的温度波动会在铜、硅及周围介电材料间产生显著的膨胀失配，导致凸起、局部约束、应力重新分布以及键合堆叠的面外变形。现有研究虽已阐明铜微观结构、界面质量、空洞演化、节距缩放和循环损伤等重要方面，但缺乏一个直接关联互连密度与全局翘曲、局部应力及不可逆循环能量分配的框架。因此，本研究旨在填补这一空白，系统考察不同互连密度下Cu-Cu键合堆叠芯片在热循环中的密度依赖性行为。
二、研究方法概述
研究人员采用了一套多尺度数值模拟与实验验证相结合的方法。首先，建立了一个代表性三维有限元模型（Finite Element Model, FEM），包含两个单晶硅芯片和居中对齐的铜互连阵列，并通过改变方形阵列规模（从50×50到70×70）定义五种密度构型。其次，采用顺序耦合瞬态工作流程，先进行瞬态热分析，再将温度历史映射至结构求解器进行变形分析，加载JEDEC标准的三周期热循环。铜的非弹性行为采用多线性随动强化（Multilinear Kinematic Hardening, MLKH）本构模型描述。此外，研究还引入了两个新的循环解析描述符：循环不可逆性分数（Cycle Irreversibility Fraction, CIF）和热塑性适应指数（Thermo-Plastic Accommodation Index, TPAI），用于量化能量分配和塑性适应程度。最后，基于仿真数据构建了监督机器学习代理辅助优化框架，并结合期望函数法（Desirability Function Approach, DFA）进行折衷分析，以识别最优互连密度。实验验证使用了Thermotron S1.2C热循环箱和PicoLog TC-08测温系统，对基准密度试样的温度和面外位移进行了同步监测。
三、研究结果
3.1 数值保真度验证
通过对比加热阶段模拟与实验的平均温度历史和中心点面外变形，验证了数值模型的准确性。模拟结果始终处于热电偶不确定带（±1.25°C）内，平均偏差低于1°C；中心点变形的终末差异为0.20微米，平均绝对偏差为0.14微米。这证实了恒定对流模型足以表征实际热载荷。
3.2 时间分辨的热机械响应
密度扫描表明，热机械优化受控于曲率抑制与约束放大之间的平衡。翘曲从50×50时的0.611微米降至65×65时的0.562微米，但在70×70时回升至0.582微米，表明存在一个临界密度，超过后过度约束反而加剧变形。主应力响应显示，最小主应力在高温保温阶段逐渐变得更压缩，而最大主应力在65×65至70×70间急剧上升，揭示了从有益载荷分担向拉伸应力集中的转变。铜互连的冯·米塞斯应力表现出明显的循环依赖性，首个循环峰值最高，随后迅速稳定，表明主要的塑性重排发生在首次热循环中。
3.3 循环解析的能量响应
CIF从50×70时的0.084单调增至70×70时的0.131，表明密度增加使循环响应从可恢复的硅储能转向不可逆的铜耗散。相反，平均TPAI从0.643降至0.588，意味着密度增大将热失配分散到更广泛的互连网络中，降低了每单位铜的平均塑性适应负担。CIF与TPAI的负相关性进一步支持了它们的配对解释。相关性分析确认，CIF与铜非弹性功和硅应变能摆动强相关，但与峰值铜等效应力呈负相关；TPAI则与铜塑性应变增量和峰值铜等效应力正相关。
四、结论与意义
本研究通过引入CIF和TPAI两个循环解析描述符，揭示了铜互连密度对Cu-Cu键合堆叠芯片热机械可靠性的非单调影响。研究表明，最佳设计并非追求最大密度，而是存在于一个中间密度区间（约64×64附近），在此处翘曲抑制达到最优，同时避免了过高的拉伸应力集中和不可逆能量耗散。期望函数法的分析结果明确指出，密度64是最稳健的可靠性折衷方案。这项工作为高性能三维封装中互连密度的理性设计提供了新的能量视角和多目标优化框架，超越了传统的基于峰值应力或应变的评估方法，具有重要的工程指导意义。
