[论文解读] The modelling of a Josephson junction and Heun polynomials

S. I. Tertychniy|ArXiv.org|Jan 30, 2006

Quantum Mechanics and Non-Hermitian Physics参考文献 2被引用 17

一句话总结

本文通过谐波驱动的一阶非线性常微分方程（ODE）对约瑟夫森结的非线性动力学进行建模，将其转化为具有多项式系数的二阶线性ODE，该方程退化为双重合流惠特克方程（DCHE）。关键贡献在于识别出使多项式解存在的参数代数约束，这些多项式解在正实轴上构成正交且可归一化的系统，从而在超导量子器件的特定参数区域内实现精确解析解。

ABSTRACT

The first order nonlinear ODE \dot ϕ(t) + \sinϕ(t)=q(t),q(t)=B+A\cosωt, where A,B,ωare real constants, is considered, the transformation converting it to a second order linear homogeneous ODE with polynoimial coefficients is found. The latter is identified as a particular case of the double confluent Heun equation. The series of algebraic constraints on the constant parameters is found whose fulfillment leads to the existance of solutions representable through polynomials in explicit form. These polynomials are found to constitute the orthogonal normalizable system

研究动机与目标

• 分析在谐波偏置下的电阻性并联约瑟夫森结的非线性动力学。
• 将一阶非线性ODE转化为具有多项式系数的二阶线性ODE。
• 识别导致变换后方程存在多项式解的参数约束。
• 建立这些多项式解在正实轴上构成正交且可归一化系统的性质。
• 将问题与双重合流惠特克方程（DCHE）联系起来，并推导出显式的解结构。

提出的方法

• 使用复相位变量将约瑟夫森结的一阶非线性ODE转化为两个线性ODE的系统。
• 引入复变量 $ z = e^{i\omega t} $，将时变系数转换为 $ z $ 的有理函数。
• 通过变量变换 $ v, \check{v} $ 将系统重写为具有多项式系数的二阶线性ODE。
• 通过莫比乌斯变换将所得ODE识别为双重合流惠特克方程（DCHE）的特例。
• 通过假设有限次数的幂级数解形式，推导出多项式解的递推关系。
• 利用Wronskian型恒等式和权函数分析，建立多项式解的正交性与可归一性。

实验结果

研究问题

• RQ1 约瑟夫森结的非线性ODE在何种条件下存在精确多项式解？
• RQ2 如何将相位动力学的一阶非线性ODE转化为具有多项式系数的二阶线性ODE？
• RQ3 变换后ODE与双重合流惠特克方程（DCHE）之间的精确联系是什么？
• RQ4 参数 $ A, B, \omega $ 的何种代数约束会导致存在次数为 $ n $ 的多项式解？
• RQ5 所得多项式解在正实轴上是否正交且可归一化？

主要发现

• 非线性ODE $ \dot{\varphi} + \sin\varphi = B + A\cos\omega t $ 等价于一个具有多项式系数的二阶线性ODE，该方程被识别为双重合流惠特克方程（DCHE）的特例。
• 对参数 $ A, B, \omega $ 的代数约束（以非负整数 $ n $ 参数化）导致存在次数为 $ n $ 的多项式解 $ P_n(z) $，其由主方程（15）导出。
• 多项式解 $ P_n(z) $ 满足递推关系（66），并通过有限个 $ 2\times2 $ 矩阵的乘积确定，从而可通过（57）和（60）显式计算。
• 这些解在 $ \mathbb{R}^+ $ 上关于权函数 $ \Xi_{n_1,n_2} $ 构成正交系统，如恒等式（71）和积分条件（72）所示。
• 在包含因子 $ \exp(-\mu(z + z^{-1})) $ 的加权 $ L^2 $-范数下，解是可归一化的，其中 $ \mu = A/(4\omega) $，确保其物理意义。
• 在约束 $ 4\omega^2(\lambda + \mu^2) = 1 $ 下，谱方程（21）可分解为两个行列式：$ \det \mathbf{G}^{(+1)} = 0 $ 或 $ \det \mathbf{G}^{(-1)} = 0 $，从而简化数值计算。