[论文解读] 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

一句话总结

该论文通过将埃伦伯格方程转化为标量里卡蒂方程，提出了一种超导性的新参数化方法，从而能够通过线性化博戈留波夫-德戈夫近似的初值问题实现高效数值求解。该方法无需显式计算本征函数或本征值，即可重构准经典传播器和局域态密度，提供了一种稳定且快速的计算方法。

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.

研究动机与目标

• 开发一种更高效且数值稳定的求解超导性中埃伦伯格方程的方法。
• 消除在准经典传播器重构过程中显式计算本征函数和本征值的需求。
• 通过简化的初值问题实现物理可观测量（如局域态密度）的精确计算。
• 在未来扩展中推广该方法以包含参数磁效应和非平衡条件。

提出的方法

• 利用准经典传播器结构导出的标量里卡蒂方程，重新表述埃伦伯格方程。
• 通过线性化博戈留波夫-德戈夫方程的初值问题的解来重构准经典传播器。
• 利用归一化条件 $\widehat{g}^2 = -\pi^2 \hat{1}$ 来约束里卡蒂解空间。
• 利用粒子-空穴对称性和自旋旋转对称性，在自旋单重态条件下将 $4\times4$ 矩阵问题简化为 $2\times2$ 形式。
• 里卡蒂参数化使得沿费米面轨迹的数值积分具有稳定性。
• 针对特定模型（如涡旋核心、畴壁）推导出解析解，使用抛物柱面函数和合流超几何函数等特殊函数。

实验结果

研究问题

• RQ1 埃伦伯格方程能否被重新表述为避免显式计算本征值和本征函数的形式？
• RQ2 能否从标量里卡蒂方程高效重构准经典传播器？
• RQ3 里卡蒂参数化是否能带来一种数值稳定且快速求解超导输运问题的方法？
• RQ4 能否为非平凡配对势（如涡旋核心或畴壁）获得精确的解析解？
• RQ5 该方法在包含参数磁效应或非平衡条件的系统中如何推广？

主要发现

• 成功地将埃伦伯格方程转化为标量里卡蒂方程，从而实现了准经典传播器的新参数化。
• 无需本征函数或本征值，即可通过线性化博戈留波夫-德戈夫方程的简单初值问题重构准经典传播器。
• 该方法产生了一种数值稳定且计算高效的求解埃伦伯格方程的算法。
• 利用抛物柱面函数和合流超几何函数，推导出涡旋核心（公式106）和畴壁（公式108）模型的精确解析解。
• 当 $\xi = \hbar v_F / |\Delta_\infty|$ 时，畴壁的解给出了 $\widehat{g}(x)$ 的显式表达式，揭示了零能处的零能束缚态。
• 在能隙边缘（$E = |\Delta_\infty|$），传播器表现出代数衰减，与束缚态的存在一致。