电镀模拟
理论与物理
电镀仿真是什么
电镀仿真是预测什么的？
简单来说，是预测镀层厚度均匀性的技术。我们希望在工件（被镀件）整个表面镀上相同厚度的镀层，但实际上会出现有的地方厚、有的地方薄的情况。电镀仿真就是事先通过计算来识别这种情况。
为什么不同位置的厚度会不同呢？
因为电流倾向于走“电阻更小的路径”。例如，在汽车保险杠上镀铬时，靠近阳极的凸起部位电流集中，镀层变厚；而凹陷部位或孔洞深处电流密度低，镀层变薄。这个“分散能力（均一电着性）”问题，正是催生电镀仿真的最大动机。
具体在哪些场景下使用呢？
列举一些典型应用——
- PCB的过孔内电镀：在直径0.1〜0.3 mm的孔洞深处均匀镀铜。通孔的高宽比越大越困难
- 半导体的铜互连电镀：用铜完全填充宽度10 nm以下的沟槽（超填充）
- 汽车的装饰性电镀：在复杂3D形状的保险杠或格栅上均匀镀铬
- 飞机零件的硬铬电镀：对要求耐磨性的零件进行膜厚管理
- 电镀治具/挂具设计：优化辅助阴极或屏蔽板的布局以提高均匀性
原来在很多常见产品中都有应用！那么，理论上是如何计算的呢？
电镀仿真根据精度分为三个层次的模型：一次电流分布、二次电流分布、三次电流分布。层次越高，物理过程考虑得越精细，但计算成本也越高。我们按顺序来看。
一次电流分布（拉普拉斯方程）
首先，一次电流分布是什么？
一次电流分布是完全忽略电极反应电阻，仅根据电解液中的电位分布来求电流的模型。假设电解液电导率 $\kappa$ 均匀且无化学反应影响，则电位 $\phi$ 遵循拉普拉斯方程：
电流密度由电位梯度得到。这是欧姆定律的微分形式：
诶，就这些？相当简单啊。
是的，方程本身很简单。但一次分布给出的是最不均匀的分布，是“最坏情况”。它将电极面视为等电位面，因此电流会一味地试图走几何上的最短路径。凸起边缘的电流会极端集中。实际电镀中，由于电极反应存在电阻，会朝着比这更均匀的方向作用。所以一次分布可以作为“均匀性的下限值”来使用。
一次电流分布的边界条件
- 阴极面：$\phi = \phi_{\text{cathode}}$（等电位）
- 阳极面：$\phi = \phi_{\text{anode}}$（等电位）
- 绝缘壁面：$\frac{\partial \phi}{\partial n} = 0$（电流法向分量为零）
忽略电极/电解液界面的过电位 $\eta$ 是一次分布的本质假设。
二次电流分布（巴特勒-沃尔默方程）
二次分布有什么不同？
二次电流分布考虑了电极表面的电荷转移反应（电极反应动力学）。也就是说，金属离子获得电子变成金属原子的反应存在一定的“阻力”。描述这种关系的是巴特勒-沃尔默方程（Butler-Volmer equation）：
出现了好多符号…请一个一个解释。
| 符号 | 含义 | 典型值·备注 |
|---|---|---|
| $i$ | 电极面的电流密度 | [A/m²] |
| $i_0$ | 交换电流密度 | 取决于反应体系（Cu: 约1〜10 A/m²、Ni: 约0.1〜1 A/m²） |
| $\alpha_a, \alpha_c$ | 阳极侧·阴极侧的传递系数 | 通常 $\alpha_a + \alpha_c = 1$（单电子反应时） |
| $F$ | 法拉第常数 | 96,485 C/mol |
| $\eta$ | 过电位（活化过电位） | $\eta = \phi_m - \phi_s - E^{\text{eq}}$（金属电位 - 溶液电位 - 平衡电位） |
| $R$ | 气体常数 | 8.314 J/(mol·K) |
| $T$ | 绝对温度 | [K] |
过电位 $\eta$ 是关键参数呢。这个值越大反应越快吗？
没错。过电位 $\eta$ 是偏离平衡状态的驱动力，$|\eta|$ 越大电流越大。但重要的是，这个方程作为非线性边界条件起作用。电解液内的拉普拉斯方程与一次分布相同，但在电极面，不是将 $\phi$ 设为等电位，而是用巴特勒-沃尔默方程将 $i$ 和 $\eta$（即 $\phi$）联系起来。由于这种非线性，试图集中在凸起部位的电流会撞上“反应阻力的墙”，结果比一次分布更均匀。
原来如此！反应慢（$i_0$ 小）的电镀液反而镀得更均匀，是吗？
很敏锐。没错，$i_0$ 越小，电极反应成为控制步骤，电流分布就越趋向均匀。氰化物系镀铜液均匀性优越，很大程度上就是因为这种反应动力学效应。
塔菲尔近似（高过电位区域）
当 $|\eta| \gg RT/F$（25°C时约25 mV）时，巴特勒-沃尔默方程的阴极侧可简化为一个指数函数：
$$ i \approx -i_0 \exp\left(-\frac{\alpha_c F \eta}{RT}\right) $$取对数可得线性关系（塔菲尔图），这是从实验数据求取 $i_0$ 和 $\alpha_c$ 的经典方法：
$$ \eta = a + b \log|i| $$其中 $a = -\frac{RT}{\alpha_c F}\ln i_0$、$b = -\frac{2.303RT}{\alpha_c F}$（塔菲尔斜率）。
三次电流分布（能斯特-普朗克方程）
三次电流分布还要考虑什么？
三次电流分布考虑了物质传输。电镀进行时，电极附近的金属离子浓度会下降（浓差极化），供应跟不上。支配离子迁移的是能斯特-普朗克方程：
右边的三项分别是：
- $-D_i \nabla c_i$：扩散（由浓度梯度引起的菲克定律）
- $-\frac{z_i F D_i c_i}{RT} \nabla \phi$：电迁移（电场引起的带电离子移动）
- $c_i \mathbf{v}$：对流（液体流动引起的传输）
其中 $D_i$ 是扩散系数，$c_i$ 是物种 $i$ 的浓度，$z_i$ 是电荷数，$\mathbf{v}$ 是流速矢量。
对流也考虑进来，那是不是也需要计算流体？
是的。要严格进行三次分布计算，就需要用纳维-斯托克斯方程求解流场，并将其代入能斯特-普朗克方程的对流项。计算量会相当大。不过在实际工作中，常用能斯特扩散层模型来近似电极附近的扩散边界层，只将层厚 $\delta_N$（典型值10〜500 $\mu$m）作为参数给定的简化方法。
比如PCB的过孔，孔深处液体不流动所以浓度容易下降，是这样吧？
正是如此。过孔深处几乎处于停滞状态，没有对流。离子仅靠扩散补充，因此高宽比（孔深/直径）越大，浓差极化越严重，镀层越薄。精确预测这种效应正是三次电流分布发挥作用的地方。
极限电流密度
当电极表面的离子浓度降为零时，电流无法再增大。这个上限就是极限电流密度$i_L$：
$$ i_L = \frac{nFD c_{\infty}}{\delta_N} $$$c_{\infty}$ 是本体浓度，$\delta_N$ 是能斯特扩散层厚度。当 $i > i_L$ 时，沉积质量会急剧恶化，产生烧焦、树枝状沉积等现象。
电中性条件
假设电解液的本体区域局部电荷呈中性：
$$ \sum_i z_i c_i = 0 $$将这个条件与各种能斯特-普朗克方程、电极面的巴特勒-沃尔默方程联立，就是三次电流分布的完整公式化。
Wagner数与均一电着性
一次、二次、三次，有选择的标准吗？全部用三次是不是最保险？
计算成本完全不同，所以选择必要且精度足够的模型很重要。用于判断的是Wagner数（瓦格纳数）Wa：
$\kappa$ 是电解液电导率，$L$ 是工件的特征尺寸，$\partial\eta/\partial i$ 是极化曲线的斜率（电极反应的微分电阻）。物理意义上是“电极反应电阻 / 电解液电阻”之比。
- Wa ≪ 1：电解液电阻占主导 → 接近一次分布（不均匀）
- Wa ≫ 1：电极反应电阻占主导 → 电流分布均匀
- Wa ≈ 1：两者影响程度相当 → 必须进行二次分布计算
也就是说，看Wa数就能判断是用一次分布还是需要二次分布，对吧？实际应用中大概是什么值？
| 电镀液 | 典型的Wa | 均匀性倾向 |
|---|---|---|
| 硫酸铜电镀（酸性） | 0.1〜1 | 低〜中 |
| 氰化铜电镀 | 5〜50 | 高（优异的均匀性） |
| 瓦特镍电镀 | 0.5〜5 | 中 |
| 硬铬电镀 | 0.01〜0.1 | 非常低 |
| 锡电镀 | 1〜10 | 中〜高 |
硬铬电镀的Wa是0.01…看起来会非常不均匀啊。
没错。所以铬电镀中会使用辅助阴极（罗伯）或屏蔽板来控制电流分布。优化这些治具的布局，正是仿真最擅长的地方。
膜厚生长模型
知道了电流分布，那膜厚怎么计算呢？
根据法拉第电解定律，由局部电流密度求膜的生长速度：
- $\delta$：膜厚 [m]
- $M$：沉积金属的摩尔质量 [kg/mol]（Cu: 0.0636, Ni: 0.0587, Cr: 0.0520）
- $i_n$：电极面的法向电流密度 [A/m²]
- $n$：电子转移数（Cu²⁺→Cu: $n=2$, Cr⁶⁺→Cr: $n=6$）
- $\rho$：沉积金属的密度 [kg/m³]
- $\eta_{\text{CE}}$：电流效率（扣除副反应部分。铬电镀约10〜25%，极低）
例如，硫酸铜电镀中 $i = 300$ A/m²、电流效率100%时，生长速度约为3.7 $\mu$m/min。
智能手机中的电镀仿真
您智能手机的处理器中，布满了数十亿条宽度10 nm以下的铜互连线。这些全部都是通过电镀（大马士革工艺）形成的。尤其困难的是，从高宽比（又细又深）的沟槽底部开始生长铜的“超填充”。PEG（聚乙二醇）等抑制剂、SPS（双(3-磺丙基)二硫化物）等加速剂、JGB（杰纳斯绿B）等整平剂——通过物质传输模型计算这三种有机添加剂的吸附、扩散、消耗的竞争，探索实现自底向上填充的工艺条件。据说，制造一片半导体芯片所需的电镀仿真计算量，是汽车保险杠一个零件的数千倍。
数值解法与实现
FEM与BEM的选用
电镀仿真使用什么数值解法？是有限元法吗？
主要有两种方法。
- 有限元法（FEM）：将整个电解液区域划分为单元。擅长处理二次、三次电流分布及非线性材料（如电导率的浓度依赖性等）。COMSOL、Ansys等使用此法
- 边界元法（BEM）：只对电极和绝缘壁的表面进行网格划分。对于一次分布（拉普拉斯方程），不需要内部的体网格，网格生成容易。对3D复杂形状的电镀治具设计有利。ElSyCA PlatingMaster是代表
なった
詳しく
報告