1つのアブレーションパターンを手動で調整して結果を見るのではなく、千通り以上の組み合わせを一度に試し、その中から最良のものをツールに指し示させる——それがこのツールです。本ガイドでは、その原理と背後の計算、そして熱工学の専門家でなくても推奨結果を読み解ける方法を説明します。
熱異方性シミュレーターは「このパターンはどれくらい良いか？」という問いに答えます。追い込みには向いていますが、探索には不向きです——1回の試行ごとに設定・計算・読み取りが必要になるからです。パラメトリックスイープは問いを逆転させます。製造可能な溝形状のすべての中で、この材料に最適なのはどれか？
本ツールは一度の実行で1,008通りの組み合わせ（スリット幅12段階 × ピッチ21段階 × パターン4種）を評価し、ユーザーが設定する優先度スライダーに従って最適解を提示します。favored方向（x方向）の伝導を強めるか、逆方向（y方向）の遮断を強めるか、あるいはその間の妥協点を選べます。
具体的な用途は次の2つで、この順に使うのが基本です。
シミュレーション前の絞り込み。熱異方性シミュレーターを開く前に、どの範囲を探せばよいかを知りたいはずです——スリット幅は90 µmか200 µmか？ピッチは500 µmか1,500 µmか？ スイープは数秒でこの答えを与え、トレードオフ曲線によって最適点がなぜそこにあり他ではないのかを示します。
4つのパターンを対等な条件で比較。各デザイン（ストレート、テーパー、ヘリンボーン、六角形）にはそれぞれ固有の最適設定があり、デザイン間で共通ではありません。任意の共通設定で比較するのは不公平であるため、スイープは各デザインをそれぞれの最適点で比較します。
スイープの範囲は恣意的なものではありません。スリット幅の下限（90 µm）はレーザー加工プロセスの公称最小トレース幅であり、ピッチの範囲は49 × 49 mmのチップ上で実際に描画可能な範囲をカバーしています。つまり、このツールが提示するものはすべて製造可能です。
完全な数値ソルバーを1,008回解くには数分かかり、スライダーが実用に耐えなくなります。そこで本ツールは2段階の方式を採用しており、これが速度と信頼性両方の鍵になっています。
解析ベースは、材料とスリットが交互に並ぶ無限積層の厳密解です。帯に沿った方向では並列、横断方向では直列。すべてのスリット比率——最も細いものも含めて——で信頼性が保たれます。これはまさに、数値ソルバーがグリッドの細かさに敏感になる領域です。
デザイン係数は「同じ空隙率で、このパターンはストレートスリットよりどれだけ良い（あるいは悪い）か？」に答えます。これは比であり、同一グリッドで行われた2つの計算の比です。共通の誤差は打ち消し合います。そのため、絶対値ではソルバーの弱点を引き継ぐことなく、欲しい情報——デザイン間の順位——だけを保持できます。
すべてを支配するパラメータは空隙率φです——スリット幅をピッチで割った値。スリットを広げる（φが大きい）ほどy方向はより遮断され異方性は上がりますが、材料が減るためx方向の伝導は落ちます。スリットを狭めればその逆になります。絶対的な最適点は存在しません。トレードオフ曲線があり、その上でどこを選ぶかという選択があるだけです。
ここまでの数値は膜単体のものです。チップ上では、この膜ははるかに厚いEFFIBLUE基板の上に乗っており、基板はあらゆる方向に伝導するため、せっかく得た異方性を「希釈」してしまいます。右側パネルには重要な3つの数値——膜の異方性、システムの異方性、その間の希釈係数——が表示されます。デザインの実際の実現可能性は、パターンの選択そのものよりも、ここで決まることが非常に多いのです。
この希釈に効く2つのレバーがあり、いずれも調整可能です。基板を薄くすることと、膜のスリットに合わせて基板に溝を入れることです。どちらか一方がなければ、パターンをどう選んでもシステムの異方性は1に近いままです。
ここで扱うすべての量は、熱平衡状態における伝導率です。本ツールは熱が伝わる時間については何も語りません——それは熱拡散コンパレーターの役割です。直接的な帰結として、比熱容量はこの計算に一切関与しません（§09の注記を参照）。
熱異方性シミュレーターと同様、計算はパターンが同一のまま無限に繰り返されることを前提としています。チップ端の縁効果はモデル化されていません。
スイープは接合ブリッジをモデル化していません。しかし機械的強度に不可欠なこのブリッジは、遮断方向を短絡させます。完全版シミュレーターでは、膜の実際の異方性は典型的には ×100 前後まで下がります。したがってスイープの数値は、構成間の順位として読むべきであり、最終的な性能値としてではありません。
フォノン補正が有効な場合、用いられる特性長は膜厚のみです。これがここで正当化されるのは、スイープ全域にわたって2つのスリット間に残る壁が最小でも300 µmあり、典型的な膜厚の1,000倍に達するためです。一方、熱異方性シミュレーターは両者の最小値を採用します。その設定範囲では、はるかに薄い壁が許容されるからです。
ここで調整できるのは基板の厚みと溝加工率のみで、その伝導率はEFFIBLUEの基準値に固定されています。これを変化させたい場合は、熱異方性シミュレーターに移ってください。
結果（CSV） — スイープされた1,008通りの組み合わせすべてを1行ずつ：デザイン、スリット幅、ピッチ、φ、k_xx、k_yy、異方性。最適点だけでなく生のスイープデータそのものであり、独自の並べ替えや表計算ソフトでの別のトレードオフのプロット、計算そのままの実行結果のアーカイブに使えます。
ここには意図的に設定エクスポートも加工用ファイルエクスポートもありません。本ツールはジオメトリを選ぶためのものであり、それを生成するためのものではないからです。選定したパターンは熱異方性シミュレーター側で再現し、Gerber形式でエクスポートしてください。
このセクションは本ガイドの技術リファレンスです。実際に実行される計算式を順を追って示し、それぞれのモデル選択の理由を説明します。ツールの利用に必須ではありません。
サイズ効果補正が有効な場合、カタログの伝導率はまず次のように縮小されます。
これは熱輸送に適用したマティーセン則です。フォノンを妨げる機構は独立に2つあります——材料内部での衝突（平均自由行程Λ）と、膜の表面での反射（経路がt_fに制限される）です。衝突頻度は加算されるため、長さの逆数が加算されます：1/Λ_eff = 1/Λ + 1/t_f、すなわちΛ_eff = Λ·t_f/(Λ + t_f)。伝導率は平均自由行程に比例するため、k_eff/k = Λ_eff/Λ = t_f/(t_f + Λ) となり、これが式(1)です。
ここで用いる特性長は膜厚のみであり、熱異方性シミュレーターのようにスリット間に残る壁と比較することはしません。これはスイープ対象領域によって正当化されます——その壁は最小でも500 − 200 = 300 µmあり、典型的な膜厚より3桁大きい値です。したがって壁が律速要因になることは決してありません。
膜厚5 nm未満では、式(1)が前提とする拡散的な輸送という仮定はもはや擁護できません——輸送は二次元的あるいは弾道的になります。それでも本ツールは補正を適用しますが、そのことを明示します。別のモデルへ黙って置き換えるよりも、そのほうがましだからです。
与えられた空隙率φに対して、パターンはまず材料とスリットが交互に並ぶ無限積層として扱われます。
これらは近似ではありません。無限層状複合材に対しては厳密解です。帯に沿った方向では材料とスリットが並んで伝導し、コンダクタンスは幅の比率に応じて加算されます——算術平均。横断方向では熱が一方を通り、続いてもう一方を通らねばならず、加算されるのは抵抗です——調和平均。これらはワイナー限界でもあります——同じφを持つどのような微細構造も、k_xx⁰を超えることも、k_yy⁰を下回ることもできません。
このベースが数値ソルバーに代わってスイープの基盤として選ばれたのは、まさに小さいφでも厳密であり続けるためです。数値ソルバーはその領域で解像度に敏感になります——細いスリットはわずか数セルしか占めず、そのバリア効果は過小評価されます。スイープはφ = 0.06まで下がるため、これは致命的な問題になり得ました。
ベース式(2)は、ストレート、テーパー、ヘリンボーン、六角形アレイのどれであるかを区別しません——φしか知らないのです。この違いは、スイープ起動時に一度だけ完全な数値ソルバー（有限体積法＋過緩和法、熱異方性シミュレーターと同じエンジン）に対して行われる較正から生まれます。
すなわちφ5点 × デザイン4種 × 方向2軸 = 40回のソルバー実行で、サーバー側でキャッシュされ、（膜の伝導率、スリット媒質の伝導率）の組で索引されます。
ソルバーには既知の系統的な解像度バイアスがあります——わずか数セルしか占めないスリットのバリア効果を過小評価するのです。このバイアスは主にφとグリッドに依存します。つまり、テスト対象のデザインと、同じφ・同じグリッドで計算された「ストレート」の基準値との間で、大部分が共通しています。両者の比を取ることで、これは相殺されます。そのため、絶対値を擁護しづらくするような不正確さを引き継ぐことなく、求めている情報——このパターンが基準からどれだけ異なるか——だけを保持できます。「ストレート」が構成上1になるのもこのためです。それは他のデザインと並ぶ一つのデザインではなく、物差しなのです。
ソルバーが失敗した場合、本ツールはスイープを中断する代わりに、オフラインで較正された係数表にフォールバックします。5つの較正点の間では、係数はφに関して線形補間され、範囲外では一定値に保たれます。
この演算——3回の乗算——が1,008回繰り返されます。これがスイープの瞬時性の理由です。
この1,008点を、相反する2つの目的に基づいて順位づけする必要があります。それが本ツールが最大化する評価指標Fの役割です。
正規化された2つの項。それぞれがスイープ内の最大値でスケーリングされ、0から1の範囲に収まります。W/m·Kと純粋な比をそのまま足し合わせるのは意味を成さないため、2つの無次元量を比較するわけです。
和ではなく積。これは重み付き幾何平均です。決定的な性質として、どちらか一方の項がゼロに近づけば、積もゼロに近づきます。そのため、ほとんど何も伝導しない一方で異方性だけが際立って高い構成はゼロ点になります——これが望ましい挙動です。算術平均であれば、2番目の項の強さだけで首位に立ってしまうでしょう。
異方性への対数適用。スイープ対象領域全体で、k_xxはおよそ2倍程度しか変化しませんが、異方性は2～4桁にわたって変化します。対数がなければ、2番目の項が1番目を完全に圧倒し、優先度設定にかかわらず評価指標は「できるだけ広いスリット」の一角しか指し示さなくなってしまいます。
ガードレール。底の値はk_xxとAともに1で下限が設けられ、対数の分母も10で下限が設けられています。これにより、負の数に分数冪が適用されることも、対数がゼロになることもありません。
主推奨はストレートスリットの中からのみ探索されます。これは見落としではありません。伝導が等しい条件下では、テーパーやヘリンボーンパターンが異方性で得をするのは幾何学的効果によるものにすぎません——その不均一な形状が局所的にφを狭めているだけで、実質的な伝導の向上はないのです。これらを推奨として採用することは、実質的でない利点のために加工の複雑さを求めることになってしまいます。4つのデザインすべてはスイープされ続け、「デザイン別ベスト構成」の表にはそれぞれの最適点ですべて表示されます。
使われる式は熱異方性シミュレーターと全く同じもので、推奨された構成にのみ適用されます。
膜と基板は積み重なった2枚の伝導シートです。横方向にはそれぞれがk·tというコンダクタンスを持ち寄って並列に伝導します。この和を全厚みで割ることで、膜単体と直接比較できる等価伝導率に戻します。これは希釈現象の説明でもあります——t_fが数百ナノメートル程度であるのに対しt_sがミリメートルのオーダーであるため、基板側の項が数千倍も重く効き、膜が持つ美しい異方性は平均の中に埋もれてしまうのです。基板の溝加工は、まさに基板自体を異方的にするために存在します。それがK_xx/K_yyを引き上げる唯一の方法だからです。
システム異方性が×3を下回ると警告が表示されます。それを下回ると、膜の幾何形状はチップスケールではもはや見えなくなり、パターンの選択は基板の厚みに比べて二次的な要素になります。
成膜材料に関して繰り返し寄せられる2つの質問があります。比熱容量はどこに反映されているのか、そして界面熱抵抗はどう扱われているのか。3つの熱ツールで実際に採用されている前提とともに、率直にお答えします。
決め手となるのは体制（レジーム）です。定常状態では、比熱容量は一切効果を持ちません——温度場が落ち着いた後は、熱の分布を決めるのは伝導率のみです。逆に過渡状態では、伝導率とともに熱拡散率 α = k/(ρ·c_p) を決めるのは比熱容量であり、したがって温度が変化する速さを左右します。
3つのツールは次のように分かれます。
| ツール | 体制 | 膜のρ·c_p |
|---|---|---|
| パラメトリックスイープ | 定常状態 | 関与しない |
| 熱異方性シミュレーター | 定常状態 | 関与しない |
| 熱拡散コンパレーター | 過渡状態 | 明示的に関与する |
本ツールおよび熱異方性シミュレーターについては、これは見落としでも暗黙の組み込みでもありません。両者が解いている方程式——∇·(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）を持ち、これは面外方向の抵抗に直列に加わります。
これは膜／基板間のただ1つのグローバルな界面です。幾何モデルは2層しか持たないため、界面ごとに抵抗を持つ多層積層モデルも、結晶粒間やパターン要素間の横方向抵抗も存在しません。
参考となるオーダー：グラフェン（k_z = 6 W/m·K、厚さ1.675 nm）の場合、膜の項は約3×10⁻¹⁰ m²·K/Wであり、典型的な界面抵抗より5桁小さい値です。R_zにおいては、常に界面と基板が支配的であり、膜が支配することはありません。
無視されているもの、そしてそれは意図的な前提です。R_th,ifは、k_xx/k_yyの均質化にも、K_xxにも、したがって横方向拡散前線の拡散率α_effにも一切関与しません。コンパレーターは該当するスライダーの下にこのことを明示的に表示しています。本ガイドが扱うパラメトリックスイープに至っては、界面抵抗を一切持ちません——これは純粋に膜の面内のみを扱うツールです。
横方向輸送において、界面熱抵抗は熱の経路上に直列に入る抵抗ではありません——それは膜と基板という2つの並列経路の間の結合です。これを無視することは、z方向に完全な結合があると仮定すること——すなわち、あらゆる位置xにおいて膜と基板が同じ局所温度を持つという「等価平板」の前提――と厳密に等価です。これは対象としている厚みと観測時間に対しては整合的です。
この前提が破綻する体制は特定できます——観測時間が短い場合、あるいは伝導性の高い膜が伝導性の低い基板の上にある場合です。その場合、膜は基板を引き連れる代わりに先行してしまいます——このモデルが構造的に捉えられない効果です。もし該当するケースに近いと感じる場合は、R_th interfaceパラメータへの感度は熱異方性シミュレーターで検証できますが、それはあくまで面外方向のR_zに限られます。