[論文レビュー] Tangent equations of motion for nonlinear response functions

Atsushi Ono|arXiv (Cornell University)|Mar 23, 2026

Spectroscopy and Quantum Chemical Studies被引用数 0

ひとこと要約

本論文は Gateaux 微分に基づく接線方程式系 TEOM を用いて、実時間ダイナミクスから非線形応答関数を直接計算し、明示的な多点相関関数や不安定な有限差分法を回避する。

ABSTRACT

Nonlinear response functions, formulated as multipoint correlation functions or Volterra kernels, encode the dynamical and spectroscopic properties of physical systems and underpin a wide range of nonlinear transport and optical phenomena. However, their evaluation rapidly becomes prohibitive at high orders because of combinatorial (often factorial) scaling or severe numerical errors. Here, we establish a systematic and efficient framework to compute nonlinear response functions directly from real-time dynamics, without explicitly constructing multipoint correlators or relying on numerically unstable finite-difference methods for order-resolved extraction. Our approach is based on the Gateaux derivative with respect to the external field in function space, which yields a closed hierarchy of tangent equations of motion (TEOM). Propagating the TEOM alongside the original dynamics isolates each perturbative order with high accuracy, providing a term-by-term decomposition of physical contributions. The computational cost scales exponentially with response order in the fully general setting and reduces to polynomial complexity when all perturbation directions are identical; both regimes avoid the factorial scaling of explicit multipoint-correlator evaluations. We demonstrate the power of TEOM by computing frequency-resolved fifth-order response functions for a solid-state electron model and by obtaining nonlinear response functions up to the 49th order with controlled accuracy in a classical Duffing oscillator. We further show that our time-evolution formulation allows optical conductivities to be evaluated directly while remaining numerically stable even near zero frequency. TEOM can be incorporated seamlessly into existing real-time evolution methods, yielding a general framework for computing nonlinear response functions in quantum and classical dynamical systems.

研究の動機と目的

• 量子系と古典系の非線形応答関数の計算の重要性と課題を動機づける。
• 実時間ダイナミクスから次数分解された非線形応答を得るための Gateaux 微分ベースの EOM フレームワークを導入する。
• 場関数のプロパゲーションとともに状態の導関数を閉じた TEOM の階層を開発する。
• 多変量非線形応答カーネルを再構成する周波数分解抽出プロトコルを提供する。
• 量子モデルと古典モデルで手法を実演し、計算的スケーリングと安定性を解析する。

提案手法

• 非線形応答関数を n 次の遅延カーネルとその周波数領域表現として定義する。
• 外場に対する Gateaux 微分を導入して方向性機能微分を得る。
• 元のダイナミクスと共に状態導関数の閉じた TEOM 階層を導出・伝搬させる。
• 局所的な摂動場（コサイン成分とサイン成分）を用いて infinitesimal な変化から多変数カーネルを再構成する。
• infinitesimal Gateaux 微分と特定の摂動を用いた再構成手順を通じて周波数分解カーネルを抽出する。
• スケーリング、偶奇性の考慮、自動微分との接続といった実装上の側面を論じる。

実験結果

リサーチクエスチョン

• RQ1 実時間ダイナミクスから explicit な多点相関関数を構成せずに非線形応答関数を直接計算するにはどうすればよいか。
• RQ2 外場あるいは進化する状態に依存するダイナミカル発生器が使われても閉じた TEOM 階層は妥当か。
• RQ3 Gateaux 微分から周波数分解されたオフダイアゴナルな非線形カーネルをどう再構成するか。
• RQ4 量子系と古典系の高次応答に対する TEOM の計算スケーリング特性はどうなるか。
• RQ5 固体状態電子モデルや Duffing 発振子のようなベンチマークモデルでこのフレームワークはどう機能するか。

主な発見

• TEOM は有限差分の差分法を用いず、実時間ダイナミクスから非線形応答関数を安定に次数分解して再構成できる。
• TEOM 階層は元のダイナミクスとともに混合方向微分を伝搬させ、階数に対して指数的（階乗ではなく）スケーリングするため高次の扱いが改善される。
• 第五次の周波数分解カーネルが four-band 電子モデルで計算され、Duffing 発振子では非常に高次（最大 49 次）まで制御された精度で達成される。
• この方法はユニタリ、リリューボ空间、自己無条件平均場のようなダイナミクス発生器を含む、線形および非線形の発生器をサポートする。
• 時間発展形式は近 zero 周波数での数値安定性を保ちつつ、光学伝導係数を直接評価できる。