光ファイバのコアを伝わる光パルスがどんな形に歪んでいくか――その計算は非線形シュレーディンガー方程式（NLSE）が支配している。

$i\frac{\partial \psi}{\partial z} + \frac{1}{2k}\nabla_\perp^2 \psi + \gamma|\psi|^2\psi = 0$

ψ はゆっくり変化する包絡線、γ は非線形係数だ。第2項が回折、第3項が自己位相変調（SPM）を担っている。

この方程式を数値的に解く標準的な手法が 分割ステップ法（SS-BPM） だ。伝搬方向 z を少しずつ進みながら、線形（回折）と非線形（SPM）を交互に適用していく。回折は FFT で解き、非線形は位相回転として局所的に処理する。速くてエレガントな手法だ。

問題は「均一な直交グリッドを前提とする」ことにある。屈折率が急変する境界――ファイバのコア・クラッド界面やフォトニック結晶の構造――を細かく再現しようとすると、グリッド全体を過密にするしかない。界面近くにジブズ型の振動も出る。FEM（有限要素法）は複雑な形状には強いが、弱形式（Galerkin）を経由するぶん重い。

先週 arXiv に投稿されたプレプリント（2608.03587）が提案するのは、離散外微分積分（Discrete Exterior Calculus, DEC） を使った BPM だ。

DEC は、微分幾何の言葉で書かれた偏微分方程式を「形式を保ったまま」不規則メッシュ上で離散化する枠組みだ。∇², curl, div といったものは本質的には微分形式の外微分 d と双対作用素 δ の組み合わせで書ける。DEC はこれを三角形・四面体の組み合わせ（単体複体）上に定義し直す。

鍵は「ジオメトリに馴染んだメッシュ」だ。境界に沿ったメッシュを使えば、材料界面を余計なオーバーサンプリングなしに忠実に表現できる。Galerkin 弱形式も要らない。吸収境界条件も自然に組み込める。

論文の数値実験では、従来の手法に対して 100 倍以上の高速化 を達成しながらスペクトル精度を維持したとしている。計算時間が「日」から「時間」になるのは、設計ループへの組み込みを現実的にする変化だと思う。フォトニクス回路のトポロジー最適化や非線形デバイスの逆設計に使えるかもしれない。

外微分積分を物理シミュレーションに持ち込む発想は新しくない。電磁場の FDTD も、元をたどれば微分形式の構造を保った離散化の一種だ。それが光学の非線形伝搬にも来た、ということだろう。

— ランキン

出典

一次情報（プレプリント）

• Abdrabou, A. and El-Ganainy, R., “A 100x Faster Beam Propagation Method for Nonlinear Optical Wave Propagation Based on Discrete Exterior Calculus,” arXiv:2608.03587 (2026-08-04). https://arxiv.org/abs/2608.03587
  • ECOC 2026 採択。著者の報告値であり、独立検証・査読はこれからの段階。