結論
ボンディングワイヤの回路シミュレータ用モデル
ワイヤ長が$$$\lambda/4$$$以下の範囲では,
ボンディングワイヤは単なる$$$\pi$$$型LPFとみなすことができます.
LPFの対地キャパシタンスは
パッド面積$$$S$$$がパッドからグランドまでの距離$$$h$$$より
十分大きい場合,すなわち$$$\sqrt{S} \gg h$$$な領域では
$$ C = \varepsilon_r \varepsilon_0 \frac{S}{h} $$
そうでない領域では概ね
$$ C = \varepsilon_r \varepsilon_0 ( \frac{S}{h} + 2 \sqrt{S} ) $$
で計算できます.
直列インダクタンス$$$L$$$は
ボンディングワイヤの直径$$$d$$$とループ高さ$$$h_L$$$から概ね
$$L \approx h_L(7.04\times10^{-4} \cdot \ln{\frac{h}{d}}+3.8\times10^{-4})(\mathrm{nH})$$
で計算できます.
複数のボンディングワイヤが間隔$$$D$$$を隔てて並行している場合の ボンディングワイヤ間の結合係数$$$k$$$は概ね $$k \approx e^{-1.29\cdot\frac{D}{h_L}-1}$$ です.
ボンディングワイヤのSonnet向けモデル
ワイヤ長が$$$\lambda/2$$$程度までは, ボンディングワイヤを他の回路パターンと同時に Sonnetで解析することができます. その場合,図のように ワイヤ径$$$d$$$,ループ高さ$$$h_L$$$のボンディングワイヤを 幅$$$1.75 d$$$,高さ$$$0.85 h_L$$$のストリップ線路で 近似した構造に置き換えれば,計算時間を大幅に節約できます.
考察
ボンディングワイヤが図のように円弧を描くとすると, ボンディングワイヤを流れる電流と, そのリターン電流の経路が電流ループを構成する. この電流ループの大きさにより三つの部分に分けて考える必要があります.[1]
波長より十分小さい領域
ボンディングワイヤ自体は単なるインダクタとみなすことができ,
そのインダクタンス$$$L$$$は概ね
$$L \propto l \cdot \log{\frac{1}{d}}$$
です.
ボンディングパッドは平行平板キャパシタンスとして計算できますが,
パッドが小さい場合は周囲の漏れ電界を考慮した補正が必要になるでしょう.
電流ループが波長と同程度の領域
ワイヤ長が波長に対して無視できない周波数領域では, ボンディングワイヤのインダクタンスがワイヤ長に比例しません. ワイヤ長$$$l=\lambda /2$$$となる周波数で ボンディングワイヤのインピーダンスは極大になります. おそらくその付近でボンディングワイヤからの放射も極大になるでしょう.
Sonnetによる電磁界解析で ボンディングワイヤの特性を知ることは可能ですが その伝送特性は周波数に強く依存し, 広帯域の伝送路としては使用できないでしょう.
電流ループが波長より大きな領域
電磁界解析では,ボンディングワイヤの形状をきちんと再現する必要があり Sonnetでの解析は難しくなります.
たとえ解析ができたとしても, 伝送特性と放射は製造上のばらつきにも強く依存し, 狭帯域の伝送路としても使用できなくなるでしょう. バンプか,あるいは別の特別なフィードスルー構造を 検討しなければならないでしょう.
Sonnetでの解析範囲の検討
以上の考察をSonnetの解析結果で確認します.
右図のボンディングワイヤモデルは ループ高さ$$$h_L=1\mathrm{mm}$$$,基板厚さ$$$h=1\mathrm{mm}$$$,ループ長$$$l=約3\mathrm{mm}$$$で
- 25GHzで
- ワイヤ長$$$l=\lambda/4$$$
- 50GHzで
- ワイヤ長$$$l=\lambda/2$$$になります.
このボンディングワイヤの$$$\mathrm{Im}(Z)$$$は 概ね25GHzまでは周波数に対して直線的に変化しており, その領域でインダクタンスとみなすことができます. そして 60GHzに並列共振点があり,この周波数ではボンディングワイヤは もはや伝送路として使うことができません. 以上の解析結果から
- ワイヤ長が$$$\lambda/4$$$以下の領域では
- 等価インダクタンスを
- ワイヤ長が$$$\lambda/4$$$から$$$\lambda/2$$$程度の領域では
- 共振周波数を
調べれば ボンディングワイヤの電気的特性を十分把握できるはずです.
λ/4以下の領域での等価インダクタンスの導出
単独のボンディングワイヤのインダクタンス
このモデルで
ワイヤ径をd=10um,31um,100um
ループ高さをhl=100um,316um,1000umと変化させて
それぞれ30GHzまで解析しました
いずれの場合も20GHz程度まで
インピーダンスは周波数に比例していることが確認できます.
右図は
ループ高さ$$$h_L$$$とワイヤ径$$$d$$$に対する
15.9GHzでの$$$\mathrm{Im}(Z_{11})$$$をプロットしたグラフです.
15.9GHzでは$$$\mathrm{Im}(Z_{11})\approx 0.01\cdot L(\mathrm{nH})$$$
なのでグラフからインダクタンスを容易に読み取ることができます.
このデータから近似式
$$L \approx h_L(7.04\times10^{-4} \cdot \ln{\frac{h}{d}}+3.8\times10^{-4})(\mathrm{nH})$$
が得られます
結合したボンディングワイヤの結合係数
グラフは二つのワイヤを距離$D$まで接近させた
モデルbw_pair.zonで
二つのワイヤの距離D をループ高さ$$$h_L$$$の0.2倍から5倍まで
ループ高さを100$$$\mu$$$mと316$$$\mu$$$mに変化させた結果です.
結合係数$$$k$$$は $$$\mathrm{Im}(Z_{21})/\mathrm{Im}(Z_{11})$$$ なので,このデータから近似式 $$k \approx e^{-1.29\cdot\frac{D}{h_L}-1}$$ $$M=k \cdot L$$ が得られます
$$$\lambda/4$$$から$$$\lambda/2$$$の領域でのSonnet向け簡易モデル
階段近似モデル
上記のSonnetの解析モデルでは 円弧を描いたボンディングワイヤを, 右図のように 階段状のモデルで近似しました. 円弧を$$$\frac{\pi}{2N}$$$毎の角度で分割し, 円弧に内接する階段と外接する階段が考えられますが, 上記の検討ではN=5の外接モデルを使いました.
しかし多数のボンディングワイヤを同時に電磁界解析するためには 実用的な解析時間で解析できるようにモデルを単純化する必要があります.
以下ではこの階段近似モデルの妥当性を検討し, より単純なモデルを求めます.
並列共振周波数
N=1,2,3,4,5で,しかも外接,内接の場合で解析してみると 並列共振周波数は最大±10%変動します.
特にNが小さいとviaの長さが波長に対して無視できない長さになるので Sonnetの解析モデルでは原理的に誤差が大きくなってしまいます.
インダクタンス
同様にインダクタンスも最大±10%変動しました. N=1で外接するモデルは最もインダクタンスが大きく, N=2で内接するモデルは最もインダクタンスが小さい結果になりました. Nが大きい時は概ねその中間に収束します.
右図の モデルbw_simple.zonのように ループ高さhLに対して, 高さ$$$0.85 \cdot h_L$$$,長さ$$$2 \cdot 0.85 \cdot h_L$$$ の形状で近似すると 解析負荷を大幅に小さくしながら, 階段近似による誤差を最小にすることができます.
円形ワイヤとリボンワイヤ
Sonnetで円形断面の導体モデルを扱うことは 解析負荷が非常に大きく現実的ではありません. もし断面が波長より非常に小さければ, 円形の断面と等価なリボン断面の導体に置き換えることができます.
例えば文献[2]によると 円形断面の導体のインダクタンス$$$L_\mathrm{round}$$$と リボン状断面の導体のインダクタンス$$$L_\mathrm{ribbon}$$$はそれぞれ特に簡単のために$$$k=2$$$とする場合もあります. これによる誤差は $$$d$$$=0.001cmのとき -1.8% $$$d$$$=0.01cmのとき-2.5%に過ぎません.
文献
- [1] 石飛,"波長による高周波の分類とそれぞれの性質--高周波に取り組む前に知っておくべき概念
- [2] "Transmission Line Design Handbook (Artech House Antennas and Propagation Library"
- [3] 石飛,"ボンディングワイヤの集中定数等価回路モデル抽出," 超高速高周波エレクトロニクス実装研究会 Vol.11,No.2. の発表に使った資料.