基于分数阶微分算子的分数阶薛定谔方程初边值问题的适定性分析与数值解研究

分数阶薛定谔方程是量子力学和分数阶微积分理论交叉融合产生的重要成果。在近些年，它在理论物理领域以及应用数学领域引起了较多的关注。该方程通过引入分数阶微分算子来对传统的薛定谔方程进行扩展，这样做能够更为精准地描述具有反常扩散特性的量子系统所呈现出的诸如非局域性、长程相关性以及记忆效应等复杂现象。分数阶微分算子的关键之处在于把整数阶导数推广到任意阶次，以此构建出更加灵活的数学模型。其数学表达式通常采用Riemann - Liouville定义或者Caputo定义，其中Riemann - Liouville定义更适合用于纯理论推导，而Caputo定义在处理包含初始条件的实际工程问题的时候更加便捷。

研究分数阶薛定谔方程的数值解，要先对其初边值问题的适定性进行分析。适定性主要涉及解的存在性、唯一性以及稳定性这三个关键的方面。在证明存在性的时候，一般会借助Galerkin方法或者不动点定理，把无限维问题转化成为有限维近似问题来求解；证明唯一性依靠的是能量估计方法，通过建立先验估计式来排除存在多个解的可能性；稳定性分析要求解对于初始条件和边界条件有连续依赖性，一般是通过构造合适的Lyapunov函数来实现。这些理论推导过程既需要有扎实的泛函分析基础，还需要结合分数阶算子所具有的特殊性质，例如非局部性和弱奇异性。

在数值求解方面，分数阶薛定谔方程的计算复杂度要明显高于传统的整数阶模型。目前经常使用的数值方法有有限差分法、有限元法和谱方法。有限差分法通过构造分数阶Grünwald - Letnikov近似格式，能够直接对分数阶导数进行离散，不过需要留意数值振荡问题；有限元法更加适合处理复杂几何区域，通过选取合适的分数阶Sobolev空间来建立变分形式；谱方法在处理周期边界条件问题的时候具有明显的优势，但是要谨慎处理分数阶导数的数值积分。在实际应用的时候，这些方法常常会和时间离散方案（例如Crank - Nicolson格式）相结合，从而形成完整的时空离散系统。

分数阶薛定谔方程的研究具有重要的应用价值。在材料科学领域，它能够更加准确地模拟半导体中载流子的反常输运行为；在生物物理领域，它为蛋白质折叠动力学研究提供了新的建模工具；在量子信息科学领域，它对开发新型量子计算器件具有指导作用。随着分数阶理论不断地完善以及计算技术持续地发展，这个方程会在更多的交叉学科领域发挥重要的作用，其数值解法研究也将为实际工程问题提供更为可靠的计算工具。

Riemann - Liouville分数阶微分算子$D_{a+}^{\alpha} f(t)$的定义式是这样的：

Caputo算子的定义和Riemann - Liouville分数阶微分算子不同。Caputo算子先对函数求导，之后再进行分数阶积分，其具体形式为：

和整数阶微分算子相比较，Riemann - Liouville分数阶微分算子与Caputo算子都有明显的非局部特性。函数在某一点的微分值不是由局部信息决定，而是和历史时刻的所有状态有关，这种和历史状态有关的特性也被叫做记忆性。可以举一个例子来说明，一阶导数$f'(t)$仅仅和$t$时刻的局部信息有关系，但是分数阶微分却需要用到积分区间$[a,t]$里的全部数据。

分数阶微分算子具有线性性质。对于任意常数$k$1、$k$2和函数$f(t)$、$g(t)$，有如下等式成立：

在逆运算关系方面，分数阶积分算子\(I^{\alpha}\)和微分算子\(D^{\alpha}\)是广义逆的关系，满足下面这个表达式：

其中$n - 1<\alpha<n$。复合运算的性质主要体现为半群性，其具体的形式是：

在分数阶薛定谔方程的实际应用当中，Caputo算子有一个突出的优势，那就是能够自然地处理初始条件。就拿时间分数阶薛定谔方程来说：

当使用Caputo算子时，初始条件$\psi(x,0)$能够直接代入进行计算，不会像使用RL算子那样遇到分数阶导数的初值问题。空间分数阶导数通常会用到Riesz分数阶拉普拉斯算子，其表达式为：

这种Riesz分数阶拉普拉斯算子的非局部特性，能够更加准确地描述量子系统里的长程相互作用。而这些Riemann - Liouville分数阶微分算子、Caputo算子以及Riesz分数阶拉普拉斯算子等算子所具有的特性，为后续开展分析解的存在性、唯一性等适定性问题的相关研究提供了非常扎实的数学基础。

### 2.2初边值问题解的存在性与唯一性

研究分数阶薛定谔方程初边值问题的适定性，是开展该方程数值解研究的理论根基。现在考虑这样一个情况，有一个时间分数阶薛定谔方程的初边值问题，这个问题具体呈现出下面这样的形式：

这里面的$\partial_t^\alpha$所代表的是Caputo型分数阶导数，其中参数$\alpha$的大小是处在满足$0<\alpha\leq1$这个范围里的；而$(-\Delta)^{\beta/2}$是Riesz分数阶拉普拉斯算子，关于$\beta$，它的取值范围是$1<\beta\leq2$。接下来看看这个问题的初边值条件设置，具体如下所示：

在上面这两个式子当中，$\Omega$是$\mathbb{R}^n$里的有界区域，$u_0(x)$是一开始就给定好的函数，$f(x,t)$则是外源项。

现在要去证明解的存在性，这里采用Banach不动点定理来完成这个证明。具体的做法是定义一个会作用在函数$v(x,t)$上的算子$A$，这个算子$A$的具体形式如下：

通过进行先验分析之后，能够得到这样的结果：

从这个结果能够知道，当$T$的值变得足够小的时候，算子$A$就会变成压缩映射。一旦算子$A$成为压缩映射，那么就会存在唯一的不动点，而这个唯一不动点的存在也就证明了解的存在性。

下面探讨解的唯一性。先做出一个假设，假设存在两个解分别为$u$1和$u$2，然后令$w = u$1 - u2。接着利用能量估计的方法去进行分析，经过这样的分析之后可以得到以下表达式：