基于分数作用量类变分方法的玻色-爱因斯坦凝聚体中的分数耗散：修正后的格罗斯-皮塔耶夫斯基动力学、临界温度及集体模式

《Physica C: Superconductivity and its Applications》：Fractional Dissipation in Bose–Einstein Condensates via Fractional Action-Like Variational Approach: Modified Gross–Pitaevskii Dynamics, Critical Temperature, and Collective Modes

　　摘要：我们在分数作用量变分方法框架下发展了玻色-爱因斯坦凝聚体动力学的分数扩展形式，其中经典作用量通过黎曼-刘维尔时间核得到修改。这一表述方式产生了一个广义的格罗斯-皮塔耶夫斯基方程，该方程包含一个内在的时间依赖耗散项，它打破了时间反演对称性，同时保留了理论的变分结构。由此产生的动力学在关键观测量上表现出普适的幂律缩放规律，这些观测量包括原子数、凝聚体密度以及托马斯-费米半径，且无需与外部热库耦合或引入现象学上的损失机制。通过结合卡普托分数时间导数与里兹分数空间算子，该模型进一步被推广，从而能够同时描述记忆效应和异常色散现象。激发谱仍保持类似博戈留波夫的结构，但出现了非零虚部，从而导致集体模式的动量无关阻尼。从热力学角度而言，该框架预测临界温度和凝聚体分数都会降低，这与平衡性质的分数重整化一致。总体而言，这种方法为量子流体中的非马尔可夫耗散、无尺度松弛以及异常输运提供了统一的描述，同时也为超冷原子系统提供了可实验验证的启示。

引言：1995年，人们在稀薄碱金属气体中首次实现了玻色-爱因斯坦凝聚，这一现象至今仍是研究宏观量子现象最重要的平台之一。玻色-爱因斯坦凝聚体为研究超流性、量子化涡旋、物质波干涉以及远离平衡态的量子动力学提供了高度可控的环境——而这些状态在强关联的凝聚态系统中往往难以实现[1,2]。在零温度下，弱相互作用凝聚体的动力学可以通过极值化与复序参量相关的平均场作用量来精确描述，由此得到的方程即为格罗斯-皮塔耶夫斯基方程[3,4]。

格罗斯-皮塔耶夫斯基方程在描述各种凝聚现象方面取得了显著成功。它能定量再现谐振陷阱中的平衡密度分布，准确预测呼吸模式和四极模式等集体激发频率，还为理解涡旋和孤子等拓扑激发提供了理论基础。尽管如此，传统的格罗斯-皮塔耶夫斯基方程存在一个固有局限：它描述的系统是严格守恒的。粒子数、能量以及凝聚相的辛结构都是完全守恒的，而耗散过程只能通过阻尼项、复杂势能或与热库及非凝聚部分的耦合以现象学方式纳入模型[5]。

这一局限引发了一个重要的概念性问题。实验研究中经常观察到原子逐渐丢失以及集体激发的阻尼现象，而这些现象并不总归因于技术缺陷。即便在超高真空条件下，背景压力低于10^-11毫巴，温度处于纳开尔文量级，长期存在的凝聚体也常常表现出集体模式的非指数衰减和光谱展宽现象。虽然三体复合以及与残留热原子的相互作用等因素会导致这些现象，但在接近零温度的情况下耗散现象依然存在，这就引发了这样一个问题：耗散是否可能是量子场的内在特性，而非仅仅源于环境耦合[[6], [7], [8], [9], [10]]。

正是由于理论预测与实验观察之间的这种矛盾，我们开始在分数作用量变分方法框架下研究玻色-爱因斯坦凝聚体动力学。这种形式主义通过用加权记忆核替换作用量中的标准时间积分，对传统的变分原理进行了推广[[11], [12], [13], [14], [15], [16]]。由此产生的理论引入了一个无量纲的分数参数α，该参数可以在α=1时的标准守恒动力学与α<1时的逐渐耗散行为之间连续过渡。与基于分数时间导数的表述方式不同，FALVA在运动方程中保留了常规的整数阶导数。该理论的分数特性则通过一个与(1?α)/(t?τ)成正比的普遍时间依赖摩擦项体现出来，这一项源自加权作用量的极值化过程。

这种形式主义具有很大吸引力，因为它是通过修改作用量本身来引入耗散的，而非在运动方程层面添加现象学上的耗散项。这样一来，它既保留了理论背后的拉格朗日量和哈密顿量结构，又自然地融入了不可逆动力学。分数作用量变分方法已经在包括具有线性阻尼的经典力学、具有有效宇宙摩擦的宇宙学模型以及量子振荡器在内的多种物理系统中展现出强大的预测能力，在量子振荡器中，它预测了更高的零点能量以及呈幂律而非指数形式的退相干现象。因此，将这一框架扩展到多体量子气体是一种自然而合理的进展。由于格罗斯-皮塔耶夫斯基方程本身源自U(1)对称作用量，因此在作用量层面直接嵌入分数记忆机制，为引入内在耗散、原子丢失以及相位扩散提供了理论上一致的途径，同时还能保持模型的变分基础。这样一来，该方法为描述玻色-爱因斯坦凝聚体实验中观察到的粒子缓慢丢失和集体模式阻尼现象提供了统一框架，同时也为理解孤立量子系统中不可逆性和非平衡行为的出现提供了新的见解。

本文的其余部分安排如下：第2节介绍该模型的基本设定和理论框架；第3节推导出广义的分数格罗斯-皮塔耶夫斯基方程，并探讨其主要的物理意义；最后，第4节总结主要结果并给出结论性意见。

基本设定：首先，我们令凝聚体序参量为ψ(r,t)，并将其标准化，使得N(t)=∫|ψ(r,t)|2dDr。传统的格罗斯-皮塔耶夫斯基理论源于作用量S=∫LdtdDr，其拉格朗日密度为：L(ψ,ψ*,?τψ,?τψ*,?ψ,?ψ*)=i?2(ψ*?tψ?ψ?tψ*)??22m|?ψ|2?V(r)|ψ|2?g2|ψ|4，其中m是构成凝聚体的单个玻色原子的质量，V(r)是限制凝聚体空间分布的外部陷阱势，g则是平均场近似下玻色子之间的有效相互作用强度。

广义分数格罗斯-皮塔耶夫斯基方程：虽然分数作用量变分方法通过核((t?τ)α?1/Γ(α))来捕捉时间记忆效应，但它仍然保留了整数阶导数，因此假设量子涨落在算子层面是在空间上局部传播、在时间上瞬间传播的。然而，许多超冷系统还存在着更强的非局域性：偶极玻色-爱因斯坦凝聚体具有长程的1/r3相互作用，处于能带边缘附近的光学晶格会产生有效的异常色散ks，其中s<2，此外还存在非马尔可夫性……

结论：总之，分数作用量变分方法为格罗斯-皮塔耶夫斯基框架提供了统一的扩展，通过非局域时间核，玻色-爱因斯坦凝聚体动力学获得了内在的记忆效应和不可逆性。由此得到的分数格罗斯-皮塔耶夫斯基方程保留了标准的平均场非线性结构，同时引入了一个普遍的时间依赖耗散项，这一项无需外部热库即可打破时间反演对称性。

作者贡献说明：Rami Ahmad El-Nabulsi：方法论、研究、形式分析、概念构建。Waranont Anukool：验证、资源协调、项目管理。