Spectrally smooth broadband response via autocorrelation-constrained inverse design

技術要約：自己相関制約を用いた広帯域スペクトル平滑化のための逆設計

問題提起
フォトニクスにおける時間領域の逆設計は、広範な周波数領域にわたる光学デバイスの性能を最大化するための効率的な手法として広く認識されている。標準的なアプローチでは、ポインティングの定理から導かれる時間積分された場量（エネルギー、フラックス、または散逸電力など）を利用する。パーセバルの定理を通じて、これらの時間領域の目的関数は周波数積分された応答に対応する。しかし、著者らは決定的な限界を指摘している。すなわち、これらの積分された目的関数は計算効率には優れているものの、帯域内の全応答のみを測定するという点である。これらは、望ましくない帯域内のスペクトルリップルや、狭いスペクトル特徴（ディップなど）、あるいはサイドローブを捉えたり、ペナルティを与えたりすることができない。その結果、最適化アルゴリズムは全エネルギーの局所的極大値へと収束する可能性があるが、その過程で、目的関数からは実質的に不可視である鋭い欠陥を持つ、スペクトル的に「粗い」応答を生み出してしまう。既存の解決策（適応的なスペクトル重み付けなど）は、最適化中に不足している領域を明示的に特定し、励起パルスを動的に変更する必要があり、複雑さを増大させる。

手法
著者らは、パルスの適応を必要とせずにスペクトルの平滑性を定量化するために設計された、時間領域のメトリックを提案している。この手法の核心は、信号とそのエネルギースペクトル密度がフーリエ変換ペアを形成するというウィーナー・ヒンチン定理に基づいている。

• メトリック ($T_{AC}$): 著者らは、重み付けされた長ラグ（long-lag）自己相関エネルギーに基づくメトリックを定義している。広範で平滑なスペクトルは、急速に減衰する自己相関関数に対応する。逆に、狭いスペクトル特徴（ディップなど）は、自己相関関数において減衰が遅く、振動的なテールとして現れる。
• 定式化: メトリック $T_{AC}[s; s_0]$ は、応答信号 $s(t)$ の重み付けされた長ラグ自己相関エネルギーと、広帯域参照信号 $s_0(t)$ のそれとの比として定義される。重み付け関数はラグに対して二次関数的（$\tau^2$）に増加し、後半の時間領域の「リンギング」テールを強調するように設計されている。
  • $T_{AC} = 1$ は、参照信号と同様に滑らかな応答であることを意味する。
  • $T_{AC} > 1$ は、鋭いスペクトル特徴の存在を示す。
• 最適化フレームワーク: このメトリックは、有限差分時間領域法（FDTD）を用いた密度ベースのトポロジー最適化（TopOpt）フレームワークに組み込まれている。メトリックは以下の2つの方法で適用される：
  1. 制約条件として: 反射効率（$\eta_r$）を $T_{AC} \leq 1$ という条件下で最大化する。
  2. ペナルティ項として: 重み付けされた目的関数 $J = (1-\alpha)\eta_r - \alpha T_{AC}$ を最大化する。ここで $\alpha$ は効率とスペクトル平滑性のバランスをとる。

主な貢献

1. 新規の時間領域メトリック: 狭いスペクトル特徴の検出を時間領域の量（長ラグ自己相関エネルギー）へと変換し、標準的な時間領域の勾配ベース最適化器で扱えるようにするメトリックの導入。
2. ブラッグ回折格子への実証: 1次元誘電体ブラッグミラーの逆設計へのメトリックの適用。本研究は、全エネルギーの最大化だけではスペクトルの欠陥を防ぐには不十分であり、自己相関制約がこれらの欠陥を効果的に抑制することを実証している。
3. 物理的洞察: 本手法は、既知の物理的な設計戦略（アポダイズドおよびチャープ構造）を純粋に最適化を通じて再現できることを示しており、スペクトルの平滑な広帯域解へとソルバーを導くためのメトリックの能力を検証している。

結果
著者らは、300 THzを中心とする105 THzの帯域幅において、屈折率コントラストが弱い（$\Delta n \approx 0.22$）1次元誘電体回折格子を用いて手法を評価した。

• 無制約 vs 制約あり: 無制約の最適化では、結果として得られる反射スペクトルに顕著な狭いスペクトル特徴（$\pi$ 位相シフトされた回折共鳴に似たもの）が現れ、高い $T_{AC}$ 値（例：8.6）を示した。$T_{AC} \leq 1$ の制約を課すと、これらの特徴が一貫して除去され、滑らかでフラットトップな反射スペクトルが得られた。
• 統計的堅牢性: 100個のランダムな初期密度を用いた統計的評価により、制約が中央値の反射効率（$\eta_r$）を大幅に低下させないことが示された。中央値の $\eta_r$ はほぼ同一（無制約時 0.8096 vs 制約あり 0.8093）であり、これは制約が性能を制限するものではなく、正則化として機能していることを示している。
• パレート・トレードオフ: ペナルティ付き目的関数の重み付けパラメータ $\alpha$ を変化させることで、効率と平滑性の間のトレードオフをマッピングした。弱いペナルティ（$\alpha = 10^{-2}$）は、実際には効率の局所最適値を改善しつつ、$T_{AC}$ をほぼ1に達成した。より強いペナルティは、ピーク反射率を犠牲にして極端な平滑化をもたらすアポダイズド・プロファイルへと導いた。
• チャープ構造の発見: より長い設計領域（$L = 17.2 \mu m$）において、最適化は自然にチャープ回折格子構造へと収束した。得られたバイナリ設計のパラメトリック解析により、局所的な回折周期が線形に変化していることが判明し、これはターゲットとなるスペクトル境界のブラッグベクトルと一致した。このパターンに基づき再構成されたチャープ回折格子は、生のTopOpt設計と同等の性能を達成し、チャープ構造の物理的妥当性を裏付けた。

意義および主張
本論文は、提案されたメトリックが、局所的なスペクトル欠陥に対する直接的な感度を提供することで、統合された時間領域の目的関数が抱える中心的な限界に対処することを主張している。著者らは以下のように述べている：

• 当該メトリックは、励起パルスの適応を必要とせずに、広帯域のスペクトル変動を評価し抑制するための効率的な時間領域の尺度として機能する。
• これは、重み付き目的関数、制約、または事後評価ツールとして利用可能である。
• 特定のテンポラル・エンベロープを規定するのではなく、$T_{AC}$ をペナルティ化することは、より緩やかな制約を課すことになり、スペクトルの一様性を確保しつつ、詳細な時間的形状を決定するための自由度を最適化器に与える。
• このアプローチは、単純な回折格子を超えて、広帯域近接場結合、アクロマティック集光、広帯域吸収などの他の困難な逆設計問題にも適用可能である。

本研究はVILLUM FONDENの支援を受けており、基礎となるコードは一般に公開されている。