[論文レビュー] 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）およびアドミアン分解法（ADM）が線形および非線形シュレーディンガー方程式に応用される場合を厳密に検討し、物理的に非現実的で解析的に解けるケースに適用した場合、これらの手法が指数関数のテイラー展開に帰着することを示している。主な貢献は、抽象的形 $ 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.

研究の動機と目的

• サジギとガンジが提示した、HPMおよびADMが線形および非線形シュレーディンガー方程式に応用された最近の研究の分析と批判を行う。
• 物理的に非現実的な初期条件が用いられた場合、HPMおよびADMによって得られる解が、特に $ e^{ibt} $ の指数関数のテイラー展開と自明に同等であることを示す。
• 初期関数が可積分性を満たさないなど、選ばれたテストケースに物理的意味の欠如が顕在されることを暴露する。
• 抽象的線形時間発展方程式 $ u_t + i\hat{A}u = 0 $ に適用可能なHPMおよびADMの一般化された形式を提案し、$ \hat{A} $ の固有状態に適用した場合の自明性を強調する。
• 数学的に正しいが物理的に意味のない結果が、ホモトピーに基づく手法を重視するジャーナルに掲載されることに対する研究者への警告を発する。

提案手法

• 著者は、物理的に非現実的な初期条件（例：$ u(x,0) = 1 + 2\cosh(2x) $）を用いて、線形シュレーディンガー方程式 $ u_t + i u_{xx} = 0 $ にHPMおよびADMを適用する。
• 非線形シュレーディンガー方程式 $ 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を線形シュレーディンガー方程式に適用した場合、標準的な指数関数のテイラー展開とは本質的に異なる結果が得られるか？
• RQ2 なぜ物理的に非現実的な初期条件（例：$ u(x,0) = 1 + 2\cosh(2x) $）が、シュレーディンガー方程式のHPMおよびADM研究で用いられるのか？その物理的意味への影響は何か？
• RQ3 非線形シュレーディンガー方程式を定数ノルム $ |u| = 1 $ に制限することで、非線形問題を解くという主張がどの程度無効化されるのか？
• RQ4 HPMおよびADMを抽象的線形時間発展方程式 $ u_t + i\hat{A}u = 0 $ に一般化できるか？その一般化にどのような意味があるか？
• RQ5 なぜ、$ e^{ibt} $ のテイラー展開のような自明な数学的恒等式に基づく結果が、新規手法であるが如く、査読付きジャーナルに掲載されるのか？

主な発見

• 物理的に非現実的な初期条件（例：$ u(x,0) = 1 + 2\cosh(2x) $）を用いた線形シュレーディンガー方程式へのHPMおよびADMの適用は、数学的に $ e^{-4it} $ のテイラー展開に等価な解を生成する。例として $ u(x,t) = 1 + 2\cosh(2x) \left[1 - 4it + \frac{(4it)^2}{2!} + \cdots \right] $ が得られ、これは既知の正確解に収束する。
• 非線形シュレーディンガー方程式において、$ |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は、$ e^{ibt} $ の $ t = 0 $ におけるテイラー展開を計算しているに過ぎず、新規または近似解法としての独自性を持たないことが示された。
• 著者は、これらの手法が真に非線形または複雑な初期値問題を解くことができないことを示し、非物理的初期データを有する特定の解析的に解けるケースにしか機能しないことを明らかにした。
• 一般化形式 $ u_t + i\hat{A}u = 0 $ は、解 $ u = \sum_{n=0}^\infty \frac{(-it\hat{A})^n}{n!} u_0 $ を与え、これは指数作用素のテイラー展開として自明に成立しており、手法の新規性の欠如を確認した。
• 本稿は、このようなアプローチが新規でも強力でもなく、単に同じ数学的恒等式をより煩雑な形で再表現しているに過ぎず、ホモトピー指向のジャーナルにおける編集的バイアスのため、受け入れられていると結論づけた。