Ginestium™膜にレーザーで刻まれたスリットパターンが、どのようにして熱を一方向には自在に走らせ、もう一方向には遮断させるのか——そして熱工学の専門家でなくても、本ツールが計算する結果をどう読めばよいのか。本ガイドではその仕組みを説明します。
EFFIBLUEチップ上では、発熱点からの熱を一方向（放熱器や筐体端の方向）へ素早く逃がしつつ、もう一方向（隣接する敏感な部品の方向）には広がらせたくない、という要求がよくあります。中実の金属膜はあらゆる方向に同じように伝導するため、このような方向選別には向きません。
解決策は、Ginestium™膜に極めて細いレーザースリットを規則的なパターンで刻むことです。うまく選ばれたパターンは、一方向にはほぼ自由に熱を流し、もう一方向には強く食い止めます。本シミュレーターは、与えられたパターンと設定に対して、この方向選別がどの程度機能するかを、実際にアブレーション加工する前に計算します。
本ツールは具体的な問いに答えます——「このスリットパターンを、この膜厚で、この基板の上に刻んだ場合、favored方向（x方向）ではどの程度の伝導が得られ、遮断したい方向（y方向）ではどうなるか？」。さらに2つの実用的な指標も提供します——このデザインは機械的に持ちこたえるか、そして膜と基板の熱膨張差によって熱サイクルを重ねるうちに疲労するおそれはないか。
本ツールは加工前の探索用です。パターンを設定し、結果を見て、調整し、満足のいく結果が得られたら、レーザーアブレーション装置が使用する形式で図面をエクスポートします。
膜を上から見たところを想像してください——規則的な細いスリット（空気、真空、その他の絶縁体で満たされている）で区切られた、中実の帯状構造です。帯に沿った方向では、熱はほぼ障害物なしに流れます——これがfavored方向であるx方向です。逆方向へ進むには、必ずスリット（絶縁体）を横切らねばならず、これが熱を強く食い止めます——これが遮断したいy方向です。
この計算の結果として得られるのは2つの数値です——k_xx（favored方向の実効伝導率）とk_yy（遮断方向の実効伝導率）。両者の比 k_xx / k_yy が得られる異方性であり、大きいほど方向選別がうまく機能していることを意味します。
膜は真空中に浮かんでいるわけではありません。はるかに厚いEFFIBLUE基板の上に載っています。したがって熱は2つの異なる方向へ出て行くことができ、本ツールはこれを2つの異なるモデルで扱います。
具体的には、システムk_xxとシステムk_yyはこの基板による希釈を考慮した値であり、一方面外方向R_thは「熱が積層構造を垂直方向に通過しなければならないとき、どれだけの総抵抗に出会うか」という問いに答えます。
非常に薄い膜では、熱の運ばれ方はもはやバルク材料と全く同じではありません。熱を運ぶ担い手（フォノン）が膜の端に近いことで妨げられるからです。膜が彼らの「平均自由行程」（材料ごとに固有の特性距離）に対して薄いほど、実際の伝導率はバルク材料の値を下回っていきます。「閉じ込め効果でkを補正」チェックボックスをオンにするとこの補正が適用されます——これは既定では無効になっており、非常に薄い膜に対して意図的に有効化するものです。
モデルは何を前提としているかを知って初めて役立ちます。以下は各計算の背後にある前提です。
計算は、スリットパターンが見渡す限り同一のまま繰り返されることを前提としています（1つの「セル」を解いて、それを頭の中で複製する）。実際のチップ端ではパターンが途切れますが、本ツールはこの縁効果をモデル化していません。
k_xxとk_yyは熱平衡状態（システムが落ち着くまでの時間が十分に経過した状態）における伝導率です。本シミュレーターは、熱が時間とともに広がる速さについては何も語りません——それは併設ツールである熱拡散コンパレーターの役割です。直接の帰結として（詳細は§10の注記）、比熱容量はこの計算に一切関与しません。
接合ブリッジはしばしば数マイクロメートルしかなく、計算グリッドには細かすぎます。その局所的な熱短絡効果は、較正済みの解析式（高解像度ソルバーに対して5%未満の誤差で検証済み）によって追加されるのであり、パターンの他の部分のようにセルごとに解かれるわけではありません。
閉じ込め補正は「拡散的」な熱輸送を前提としています——これは十分に厚い内部を持つ膜に対して成立します。5 nmを下回ると、この前提は疑わしくなります（輸送はむしろ二次元的・弾道的になります）。本ツールは補正を黙って無効化するのではなく、適用範囲から外れていることを明示します。
「機械的脆弱性」と「CTEリスク」は単純な定性的スコア（除去された材料の割合、ブリッジ密度、熱膨張差×サイクル振幅）であり、リスクのある構成を素早く見つけるためのものです——有限要素法による機械シミュレーションや実際の疲労試験の結果ではありません。
界面熱抵抗（R_th interface）はユーザーが与える単一の数値であり、物理的な接触モデルではありません——感度を探るには便利ですが、最終結論を出す前には実測値に照らして較正すべきです。これは垂直方向経路にのみ作用します。詳細は§10の注記を参照してください。
本ツールは限界事例を隠しません——計算は続けつつ、警告を出します。
k_yyの信頼性が低下します。グリッド解像度を上げるか、ピッチを広げてください。設定（JSON） — 本ツールのすべての設定（材料、幾何形状、基板、ブリッジなど）をエクスポートします。計算結果は含まれません。セッションを後で再開したり、同僚に渡してそのまま再インポートしてもらったりする際に便利です。
Gerber（.gbr） — レーザーアブレーション装置が要求する形式で、アブレーションパターンの実際のベクター図面をエクスポートします。これはこのプロセスの実用的な到達点であり、パターンが仮想的に検証されたら、このファイルはそのまま製造工程へ送られます。
このセクションは本ガイドの技術リファレンスです。サーバー側エンジンが実際に実行する計算式を実行順に示し、それぞれのモデル選択の理由を説明します。ツールの利用に必須ではありません。
サイズ効果補正が有効な場合、入力された伝導率はまず次のように縮小されます。
これは熱輸送に適用したマティーセン則です。フォノンを妨げる機構は独立に2つあります——構造内部での衝突（平均自由行程Λ）と、構造の壁での反射（経路がL_cに制限される）です。衝突頻度は加算されるため、長さの逆数が加算されます：1/Λ_eff = 1/Λ + 1/L_c、すなわちΛ_eff = Λ·L_c/(Λ + L_c)。伝導率は平均自由行程に比例するため、k_eff/k = Λ_eff/Λ = L_c/(L_c + Λ) となります。
L_cにおける最小値は、2つの閉じ込め効果のうちどちらが効いているかを示します——膜厚か、あるいは2つのスリット間に残る壁の厚みかです。本シミュレーターでは、パラメトリックスイープとは異なり、どちらも起こり得ます。だからこそ最小値を明示的に取っているのです。
同じ係数が面外方向の伝導率k_zにも適用されます——閉じ込め効果は両方向の輸送を抑制するためです。
比 k_x/k_y は確かにAに等しくなり、求められている通りです。しかし重要なのは、積 k_x·k_y が k_eff² のまま保たれることです——幾何平均が保存されるのです。したがって異方性のつまみを回しても、材料が運べる総熱量がひそかに変化することはなく、2方向への配分だけが変わります。A = 1のとき、両方向でちょうどk_effが得られます。
部分深度の加工では、スリットの下に材料が残ります。膜厚のうち、上部の割合pはスリット媒質で満たされ、下部の割合(1 − p)はまだ中実な膜のままです。横方向輸送から見ると、この2つの層は並んで——並列に——伝導するため、調和平均ではなく厚み加重の算術平均になります。加工が膜を貫通する場合（p = 1）、k_gapに帰着します。
パターンは繰り返し単位セルに縮約されます——幅はブリッジ周期（六角形の場合はピッチ×アスペクト比）、高さはピッチの2倍。各セルには、スリットに該当すればk_slit、そうでなければk_xとk_yが割り当てられます。次にこのセル上で純伝導問題を2回解きます——一度はx方向、一度はy方向です。
隣接する2つのセルの間では、熱は一方の半分を通過し、続いてもう一方の半分を通過します。これは直列の2つの抵抗であり、等しい2つの半分の等価伝導率はその調和平均です。算術平均を取ると、熱が材料／スリットの境界を「漏れ出る」ことになり、k_yyを系統的に過大評価してしまいます——これはまさに本ツールが小さくすべき値です。もう一つ有用な性質として、調和平均は2つの伝導率のどちらか一方がゼロになった瞬間にゼロになるため、高真空のセルは真の障壁として機能します。
(5)は有限体積法における流束収支です——平衡状態では、あるセルに4つの面から流入するものはそのまま流出するため、そのセルの温度は隣接セルの温度を面コンダクタンスで重み付けした平均になります。すでに更新された値を再利用しながら全セルを走査するのがガウス=ザイデル法です。これは収束しますが、細かいグリッドでは遅くなります。係数ω > 1（逐次過緩和法）は計算された補正を超えて外挿し、反復回数をおよそ1桁削減します。ωは]1 ; 2[の範囲内に収める必要があり、方式が安定するために不可欠です。ここでは1.9が採用されています。
境界条件：テスト対象の方向に垂直な2つの面には温度1と0が課され、他の2つの面には周期境界条件が課されます——これは「無限パターン」という前提を数値的に翻訳したものです。総流束は50反復ごとに測定され、その相対変化が2×10⁻³を下回った時点で計算は停止します。上限は6,000反復です。
課された温度差はセルの全長にわたってちょうど1であるため、この流束はさらなる正規化なしにそのまま実効伝導率になります。メッシュは、スリットを横切る方向に複数のセルが入るよう設計されており（高さ方向に44～220セル、幅方向に24～220セル）、全体で最大46,000セルに制限されています。それを超える場合、グリッドは相似的に縮小され、「解像度不足のスリット」警告が発生します。
交互に並ぶ帯の無限積層に対しては、これらは近似ではなく厳密解です——帯に沿った方向ではコンダクタンスが加算され、横断方向では抵抗が加算されます。これらはワイナー限界でもあります——同じ空隙率を持つどのような微細構造も、k_xx^idealを超えることも、k_yy^idealを下回ることもできません。だからこそ二重の用途があります——ソルバーを実行せずに済む即座の表示と、常にこの2つの間に収まるべき数値解の妥当性チェックです。
ステップ4のソルバーはブリッジなしで実行されます。その効果は後から解析的に再導入されます。
ブリッジは数マイクロメートル程度ですが、グリッドの各セルは数十マイクロメートルあります。これを直接解くには、あらゆる箇所でグリッドを細かくする必要があり——局所的な1つの詳細のためにセル数がおよそ×100になるコストがかかります。補正(8)はより少ない計算でより多くを実現します——ソルバーが実際に計算したスリット帯の伝導率を取り出し、そこにブリッジをその帯の中だけの並列短絡として加え、再構成します。0.88という係数はレンガ状オフセットを反映しています——互い違いに配置されたブリッジは1本の連続した経路として整列しないため、短絡効果がやや弱まります。この一連の補正は、高解像度ソルバーに対して5%未満の誤差で較正されています。
k_xxは影響を受けません——ブリッジは帯を横切る方向に材料を追加するだけで、帯に沿った伝導には何も変化を与えません。
膜と基板は積み重なった2枚の伝導シートです。横方向にはそれぞれがk·tというコンダクタンスを持ち寄って並列に伝導します。この和を全厚みで割ることで、膜単体と直接比較できる等価伝導率に戻します。これは希釈現象の機構的な説明でもあります——t_fが数百ナノメートル程度であるのに対しt_sがミリメートルのオーダーであるため、基板側の項が数千倍も重く効き、膜の異方性は平均の中に埋もれてしまうのです。基板の溝加工は、基板自体を異方的にするために存在します——それがK_xx/K_yyを引き上げる唯一の方法だからです。
垂直方向では、熱は3つの層を一つずつ順番に通過します——面積あたり抵抗が重みなしで加算されます。これがモデル中で界面熱抵抗が関与する唯一の箇所です——§10の注記を参照してください。k_zは材料の面内伝導率とは独立した、真の面外方向伝導率であることに注意してください。グラフェンのような層状材料では、その差は100倍に達し、面内の値を使えばR_zはばかげたほど楽観的な値になってしまいます。
(12)には明確な物理的意味があります——それは熱サイクル1回あたりに膜と基板の間に蓄積される差動ひずみです。2,500 ppmは膨張差0.25%に相当し、薄膜の剥離が文献で報告されているオーダーです。このスコアは粗いものですが、根拠があります。
(11)にはそのような根拠はありません——これは、除去された材料（φ、脆弱化させる要因）とブリッジ密度（b、補強する要因）を対比させた経験的なランク付けであり、より頑丈な形状に対して2つの割増係数が加わります。係数2と3には単位がなく、いかなる機械モデルからも導かれていません。これは信号機のようなものとして扱うべきであり、寸法設計の基準として使うべきではありません。
成膜材料に関して繰り返し寄せられる2つの質問があります。比熱容量はどこに反映されているのか、そして界面熱抵抗はどう扱われているのか。3つの熱ツールで実際に採用されている前提とともに、率直にお答えします。
決め手となるのは体制（レジーム）です。定常状態では、比熱容量は一切効果を持ちません——温度場が落ち着いた後は、熱の分布を決めるのは伝導率のみです。逆に過渡状態では、伝導率とともに熱拡散率 α = k/(ρ·c_p) を決めるのは比熱容量であり、したがって温度が変化する速さを左右します。
3つのツールは次のように分かれます。
| ツール | 体制 | 膜のρ·c_p |
|---|---|---|
| 熱異方性シミュレーター | 定常状態 | 関与しない |
| パラメトリックスイープ | 定常状態 | 関与しない |
| 熱拡散コンパレーター | 過渡状態 | 明示的に関与する |
本ツールについては、これは見落としでも暗黙の組み込みでもありません。ステップ4で解かれる方程式——∇·(k∇T) = 0——にはρ·c_pが現れず、そもそも本シミュレーターはパラメータとしてこれらの項を持っていません。3つのツールが共有する材料カタログには各材料のρとc_pがきちんと記録されていますが、それを消費するのは熱拡散コンパレーターだけです。
熱拡散コンパレーターでは、比熱容量は膜＋基板システムの実効拡散率を通じて関与します。
成膜層についてここで明示的に置かれている前提は2つあります。まず、膜のアブレーションされた部分は何も蓄熱しません——スリット媒質（空気、真空、エアロゲル）はこの尺度では有意な蓄熱能力を持たないとみなされ、それが (1 − ψ) の係数です。次に、ρとc_pは300 Kの値として一定に扱われます——温度依存性はモデル化されていません。
数値は共有サーバーカタログに由来します。
| 材料 | ρ (kg/m³) | c_p (J/kg·K) | データの位置づけ |
|---|---|---|---|
| Ginestium™ | 2,100 | 710 | 代替値 — 300Kにおける配向性合成熱分解炭素。Ginestium™自体の実測値ではない |
| グラフェンCVD few-layer | 2,200 | 700 | 文献値 |
| ダイヤモンドCVDナノ結晶 | 3,400 | 536 | NCD膜での実測値（302K） |
| 銅（対照） | 8,960 | 385 | 文献値、バルク銅 |
| 基板 | ユーザー入力（既定値 3,900 / 880） | — | |
「カスタム」材料にはρもc_pもありません。この場合、コンパレーターは値を推測する代わりに「拡散率なし」と明示的に返します。
厚みによる重み付けは膜の寄与を押し潰します。コンパレーターの既定設定——Ginestium™ 400 nm、基板10 µm——では、膜が占めるのはシステムの熱容量のわずか約1.6%にすぎません。基板を2 mmにすると、その割合は0.01%未満まで下がります。実効拡散率は実質的に「システムの横方向伝導率を基板の熱容量で割ったもの」になります。
直接の帰結として、拡散前線は膜のc_pに対してほとんど鈍感です。無視されているからではなく——式には確かに含まれています——積層構造の幾何学がその重みをほとんど与えないからです。膜が効くのは、式の分母にある比熱容量よりも、はるかに分子にあるk·tを通してです。
無視されているわけではありませんが、その扱いは部分的であり、面外方向に限られます。
実際に存在するもの。本シミュレーターのR_th interfaceパラメータ（0～100 ×10⁻⁶ m²·K/W、既定値10。熱拡散コンパレーターでも同じスライダー）は、面外方向の抵抗に直列に加わります——これは§09の式(10)そのものであり、ここに再掲します。
これは膜／基板間のただ1つのグローバルな界面です。幾何モデルは2層しか持たないため、界面ごとに抵抗を持つ多層積層モデルも、結晶粒間やパターン要素間の横方向抵抗も存在しません。
参考となるオーダー：グラフェン（k_z = 6 W/m·K、厚さ1.675 nm）の場合、膜の項は約3×10⁻¹⁰ m²·K/Wであり、典型的な界面抵抗より5桁小さい値です。R_zにおいては、常に界面と基板が支配的であり、膜が支配することはありません——これはGinestium™についても同様で、その高いk_z（1,000 W/m·K）により膜の項は無視できるほど小さくなります。
無視されているもの、そしてそれは意図的な前提です。R_th,ifは、ステップ4のk_xx/k_yy均質化にも、ステップ7のK_xx/K_yyにも、コンパレーターの拡散前線の拡散率α_effにも一切関与しません。コンパレーターは該当するスライダーの下にこのことを明示的に表示しています。パラメトリックスイープに至っては、界面抵抗を一切持ちません——これは純粋に膜の面内のみを扱うツールです。
横方向輸送において、界面熱抵抗は熱の経路上に直列に入る抵抗ではありません——それは膜と基板という2つの並列経路の間の結合です。これを無視することは、z方向に完全な結合があると仮定すること——すなわち、あらゆる位置xにおいて膜と基板が同じ局所温度を持つという「等価平板」の前提――と厳密に等価です。これは式(9)にすでに暗黙のうちに含まれています。式(9)は膜＋基板の組に単一の伝導率しか割り当てていないからです。
この前提が破綻する体制は特定できます——観測時間が短い場合、あるいは伝導性の高い膜が伝導性の低い基板の上にある場合です。その場合、膜は基板を引き連れる代わりに先行してしまいます——このモデルが捉えられない効果です。R_th interfaceパラメータへの感度はここで検証できますが、それは面外方向のR_thにしか現れません。
これらの単位を扱えなくても本ツールは使えます——スライダーと数値入力欄がほとんどの作業を担ってくれます。この早見表は、表示される結果を読み解く際の補助として使ってください。