[論文レビュー] Transformation of the Eilenberger Equations of Superconductivity to a Scalar Riccati Equation

N. Schopohl|arXiv (Cornell University)|Apr 6, 1998

Control and Stability of Dynamical Systems参考文献 3被引用数 26

ひとこと要約

本稿では、超伝導のEilenberger方程式を、スカラーのリッカティ方程式に変換することで、新しいパrametrizationを導入した。これにより、線形化されたBogoliubov-de Gennes方程式の初期値問題を用いた効率的な数値的解法が可能となり、固有関数や固有値の明示的計算を経ずに、準古典的伝搬関数および局所状態密度を再構成できる。この手法により、安定的かつ高速な計算的手法が得られる。

ABSTRACT

A new parametrisation of the Eilenberger equations of superconductivity in terms of the solutions to a scalar differential equation of the Riccati type is introduced. It is shown that the quasiclassical propagator, and in particular the local density of states, may be reconstructed, without explicit knowledge of any eigenfunctions and eigenvalues, by solving a simple initial value problem for the linearised Bogoliubov-de Gennes equations. The Riccati parametrisation of the quasiclassical propagator leads to a stable and fast numerical method to solve the Eilenberger equations. For some spatially varying model pair potentials exact solutions to the Eilenberger Equations are found.

研究の動機と目的

• 超伝導におけるEilenberger方程式をより効率的かつ数値的に安定した方法で解くための手法を開発すること。
• 準古典的伝搬関数の再構成にあたり、固有関数や固有値の明示的計算を回避すること。
• 簡略化された初期値問題を用いて、局所状態密度などの物理的観測量を正確に計算できること。
• 将来的には、磁気的効果や非平衡状態を含む一般化を可能とする。

提案手法

• 準古典的伝搬関数の構造に基づいて導出されたスカラーのリッカティ方程式を用いて、Eilenberger方程式を再定式化する。
• 線形化されたBogoliubov-de Gennes方程式の初期値問題の解を用いて、準古典的伝搬関数を再構成する。
• リッカティ解の解空間を制約するために、正規化条件 $\widehat{g}^2 = -\pi^2 \hat{1}$ を利用する。
• 粒子-hole対称性およびスピン回転対称性を活用し、スピンシンギュレート状態では $4\times4$ 行列問題を $2\times2$ 形式に簡略化する。
• リッカティパラメータ化により、フェルミ面に沿った軌道に沿った安定な数値的積分が可能になる。
• 放物型円筒関数や超幾何関数などの特殊関数を用いて、特定のモデル（例：渦核、ドメイン壁）に対して解析的解を導出する。

実験結果

リサーチクエスチョン

• RQ1 Eilenberger方程式を、固有値や固有関数の明示的計算を避ける形に再定式化できるか？
• RQ2 スカラーのリッカティ方程式から、準古典的伝搬関数を効率的に再構成できるか？
• RQ3 リッカティパラメータ化が、超伝導輸送問題の数値的安定性および高速性を向上させるか？
• RQ4 渦核やドメイン壁のような非自明な対称性を持つペア関数に対して、正確な解析的解を得られるか？
• RQ5 この手法は、磁気的効果や非平衡状態を含む系へどのように一般化できるか？

主な発見

• Eilenberger方程式がスカラーのリッカティ方程式に成功裏に変換され、準古典的伝搬関数の新しいパラメータ化が可能になった。
• 固有関数や固有値の計算を経ずに、線形化されたBogoliubov-de Gennes方程式の初期値問題からの単純な解から、準古典的伝搬関数を再構成できる。
• Eilenberger方程式を解くにあたり、数値的に安定かつ計算的に効率的なアルゴリズムが得られた。
• 放物型円筒関数および超幾何関数を用いて、渦核（式106）およびドメイン壁（式108）モデルに対して正確な解析的解が導出された。
• パラメータ $\xi = \hbar v_F / |\Delta_\infty|$ を持つドメイン壁の解は、$\widehat{g}(x)$ に対して明示的な表現を提供し、ゼロエネルギーに中間ギャップ状態が存在することを示した。
• ギャップ端縁（$E = |\Delta_\infty|$）では、伝搬関数が代数的減衰を示し、束縁状態の存在と整合的である。