[论文解读] Numerical investigation of the second harmonic effects in the LJJ

P. Kh. Atanasova, T. L. Boyadjiev+2 more|arXiv (Cornell University)|May 31, 2010

Physics of Superconductivity and Magnetism参考文献 9被引用 4

一句话总结

本文通过连续牛顿法与四阶数值积分方法，对具有二次谐波电流-相位关系的长约瑟夫森结（LJJ）中的双正弦-高斯方程进行了数值研究。研究揭示了新的通量子解，并表明二次谐波的符号与大小显著影响通量子的稳定性，包括在负二次谐波贡献下出现共存的稳定单通量子态。

ABSTRACT

We study the long Josephson junction (LJJ) model which takes into account the second harmonic of the Fourier expansion of Josephson current. The sign of second harmonic is important for many physical applications. The influence of the sign and value of the second harmonic on the magnetic flux distributions is investigated. At each step of numerical continuation in parameters of the model, the corresponding nonlinear boundary problem is solved on the basis of the continuous analog of Newton's method with the 4th order Numerov discretization scheme. New solutions which do not exist in the traditional model have been found. The influence of the second harmonic on stability of magnetic flux distributions for main solutions is investigated.

研究动机与目标

• 研究约瑟夫森结电流-相位关系中的二次谐波对长约瑟夫森结中磁通分布的影响。
• 利用数值分岔分析，研究在不同二次谐波参数下通量子解的稳定性。
• 识别在传统正弦型电流-相位模型中不存在的新静态通量子解。
• 考察二次谐波符号与振幅在决定通量子态稳定性及共存性方面的作用。

提出的方法

• 采用包含电流-相位关系中一阶与二阶谐波项的双正弦-高斯方程对LJJ进行建模。
• 应用连续牛顿法（CANM）求解非线性两点边值问题，实现高阶精度。
• 采用四阶数值积分方法（Numerov有限差分格式）对微分方程进行离散化，以确保数值稳定性和精度。
• 求解相关的斯特鲁姆-刘维尔特征值问题，通过最小特征值λ₀评估解的稳定性。
• 在参数（a₁, a₂, hₑ, γ, l）上执行数值延拓，追踪解分支并检测分岔点。
• 通过网格细化与收敛性检验（采用多种空间步长h = 0.15625, 0.078125, 0.0390625）验证结果可靠性。

实验结果

研究问题

• RQ1 在约瑟夫森结电流-相位关系中引入二次谐波，如何影响长约瑟夫森结中静态磁通分布的存在性与结构？
• RQ2 二次谐波的符号与大小对通量子解稳定性的影响是什么，特别是对分岔点与稳定区间的影响？
• RQ3 当二次谐波非零时，是否会出现新的静态解？它们与传统单正弦-高斯模型中的解有何不同？
• RQ4 在何种条件下，多个稳定的单通量子态会共存？这如何依赖于二次谐波参数a₂？
• RQ5 对于不同a₂值，斯特鲁姆-刘维尔问题的最小特征值λ₀随外磁场hₑ的变化如何演化，以指示稳定性转变？

主要发现

• 当二次谐波为负（a₂ < 0）时，在外磁场一定范围内，存在两个不同的稳定单通量子态（Φ¹与Φ¹*），在a₂ ≈ -0.8处被观测到。
• 当a₂ ∈ [0, 1]时，单通量子态Φ¹在更高的hₑ下变得稳定，且随着a₂超过0.7，分岔点先向左移动后向右移动。
• 当a₂ > 0.7时，不稳定临界场h_cr增大，表明正二次谐波贡献越大，系统稳定性越强。
• 当a₂ ∈ [-1, 0]时，随着a₂从负值增加，λ₀(hₑ)曲线向右移动；在a₂ ≈ -0.8处，曲线分裂为两条相交分支，表明存在两个稳定单通量子态的共存。
• 当a₂ → 0时，单通量子态Φ¹的能量趋近于1，与标准无限结模型中单个通量子的已知能量一致。
• 对于所有a₂，通量子数N[Φ¹]保持为1，因为∂N/∂a₂ = 0，证实了单通量子态的拓扑稳定性。