[论文解读] On the application of homotopy--perturbation and Adomian decomposition methods to the linear and nonlinear Schrödinger equations

Francisco M. Fernández|ArXiv.org|Aug 11, 2008

Fractional Differential Equations Solutions参考文献 3被引用 18

一句话总结

本文批判性地审视了同伦摄动法（HPM）和Adomian分解法（ADM）在线性和非线性Schrödinger方程中的应用，表明当应用于非物理的、可解析求解的情形时，这些方法实际上退化为指数函数的泰勒级数展开。主要贡献在于为抽象形式 $ u_t + i\hat{A}u = 0 $ 提出了HPM与ADM的广义框架，揭示了当这些方法应用于算符 $ \hat{A} $ 的本征态时，结果为平凡解，暴露了这些技术在量子力学背景下的方法论局限性。

ABSTRACT

I discuss a recent application of homotopy perturbation and Adomian decomposition methods to the linear and nonlinear Schrödinger equations. I propose a generalization of the procedure for the treatment of a wider class of problems.

研究动机与目标

• 分析并批判Sadighi与Ganji所提出的HPM与ADM在线性和非线性Schrödinger方程中的最新应用。
• 证明通过HPM与ADM获得的解在本质上等价于指数函数 $ e^{ibt} $ 的泰勒展开，尤其是在使用非物理初始条件时。
• 揭示所选测试案例缺乏物理相关性，例如违反平方可积性的初始函数以及模长为1的非线性解。
• 为HPM与ADM提出适用于抽象线性演化方程 $ u_t + i\hat{A}u = 0 $ 的广义形式化方法，强调当应用于 $ \hat{A} $ 的本征态时其结果的平凡性。
• 提醒研究人员警惕发表数学上正确但物理上无意义的结果，尤其是在偏爱同伦方法的期刊中。

提出的方法

• 作者将HPM与ADM应用于线性Schrödinger方程 $ u_t + i u_{xx} = 0 $，采用非物理的初始条件，如 $ u(x,0) = 1 + 2\cosh(2x) $，其不满足平方可积性。
• 对于非线性Schrödinger方程 $ i u_t + u_{xx} + \gamma |u|^2 u = 0 $，作者表明将解限制在 $ |u| = 1 $ 时，问题退化为线性方程 $ i u_t + u_{xx} + \gamma u = 0 $，该方程可解析求解。
• HPM与ADM被应用于这些简化情形，得到的幂级数解在数学上等价于 $ e^{ibt} $ 的泰勒级数，如 $ \sum_{n=0}^\infty \frac{(-it\alpha)^n}{n!} $。
• 作者将该方法推广至抽象方程 $ u_t + i\hat{A}u = 0 $，其形式解为 $ u = e^{-it\hat{A}} u_0 $，并证明HPM/ADM恢复了级数展开 $ \sum_{n=0}^\infty \frac{(-it\hat{A})^n}{n!} u_0 $。
• 该方法进一步扩展至初始态 $ u_0 = \sum_k c_k f_k $，其中 $ \hat{A} f_k = a_k f_k $，导出 $ u = \sum_k c_k e^{-ita_k} f_k $，其再次退化为关于 $ t $ 的幂级数。
• 作者认为，此类结果本质上是平凡的，仅因使用了复杂且非物理的初始条件，以及回避了标准解法，才显得非平凡。

实验结果

研究问题

• RQ1 HPM与ADM应用于线性Schrödinger方程时，其结果是否在本质上不同于指数函数的泰勒级数展开？
• RQ2 为何在HPM与ADM研究Schrödinger方程时使用非物理的初始条件（如 $ u(x,0) = 1 + 2\cosh(2x) $）？这对物理相关性有何影响？
• RQ3 将非线性Schrödinger方程限制在模长恒定 $ |u| = 1 $ 的解上，其程度在多大程度上使‘求解非线性问题’的宣称失效？
• RQ4 能否将HPM与ADM推广至抽象线性演化方程 $ u_t + i\hat{A}u = 0 $？此类推广有何含义？
• RQ5 为何基于平凡数学恒等式（如 $ e^{ibt} $ 的泰勒级数）的结果，仍被发表于同行评审期刊，并以新颖方法之名呈现？

主要发现

• HPM与ADM应用于线性Schrödinger方程并采用非物理初始条件（如 $ u(x,0) = 1 + 2\cosh(2x) $）时，所得解在数学上等价于 $ e^{-4it} $ 的泰勒级数，如 $ u(x,t) = 1 + 2\cosh(2x) \left[1 - 4it + \frac{(4it)^2}{2!} + \cdots \right] $，其收敛于已知精确解。
• 对于非线性Schrödinger方程，将解限制在 $ |u| = 1 $ 时，问题退化为线性方程 $ i u_t + u_{xx} + \gamma u = 0 $，HPM/ADM解恢复了 $ e^{i(\gamma - \alpha^2)t} $ 的泰勒级数，如 $ u(x,t) = e^{ix} \left[1 + it + \frac{(it)^2}{2!} + \cdots \right] $，其收敛于 $ e^{i(x+t)} $。
• HPM与ADM被证明等价于在 $ t = 0 $ 处计算 $ e^{ibt} $ 的泰勒级数，而非提供新颖或近似解法。
• 作者证明，该方法无法求解真正非线性或复杂的初值问题，因其仅适用于特定、可解析求解的非物理初始数据情形。
• 广义形式 $ u_t + i\hat{A}u = 0 $ 导出解 $ u = \sum_{n=0}^\infty \frac{(-it\hat{A})^n}{n!} u_0 $，其本质上即为指数算符的泰勒展开，确认了该方法缺乏新颖性。
• 本文结论认为，此类方法并非新颖或强大，而仅仅是用更繁琐的方式重述同一数学恒等式，其被接受仅因同伦方法导向期刊的编辑偏见。