[论文解读] Approach and separation of quantum vortices with balanced cores

C. Rorai, Jack W. D. Skipper+2 more|arXiv (Cornell University)|Oct 6, 2014

Quantum, superfluid, helium dynamics参考文献 11被引用 3

一句话总结

本研究通过求解三维 Gross-Pitaevskii 方程的高精度数值模拟，研究了量子涡旋重连现象，引入改进的涡旋核心追踪与初始条件控制方法，揭示了反平行与正交涡旋构型下不同的标度律。关键发现为：正交涡旋在重连前呈现 δ ∼ |tₜ − t|¹ᐟ³ 的标度行为，重连后则为 δ ∼ |t − tₜ|²ᐟ³，与经典维度预测 δ ∼ |tₜ − t|¹ᐟ² 不同，这些差异可归因于早期持续存在的 Frenet-Serret 框架错位。

ABSTRACT

Using two innovations, smooth, but distinctly different, scaling laws for the numerical reconnection of pairs of initially orthogonal and anti-parallel quantum vortices are obtained using the three-dimensional Gross-Pitaevskii equations, the simplest mean-field non-linear Schrödinger equation for a quantum fluid. The first innovation suppresses temporal fluctuations by using an initial density profile that is slightly below the usual two-dimensional steady-state Padé approximate profiles. The second innovation is to find the trajectories of the quantum vortices from a pseudo-vorticity constructed on the three-dimensional grid from the gradients of the wave function. These trajectories then allow one to calculate the Frenet-Serret frames and the curvature of the vortex lines. For the anti-parallel case, the scaling laws just before and after reconnection obey the dimensional $δ\sim|t_r-t|^{1/2}$ prediction with temporal symmetry about the reconnection time $t_r$ and physical space symmetry about the $x_r$, the mid-point between the vortices, with extensions of the vortex lines formng the edges of an equilateral pyramid. For all of the orthogonal cases, before reconnection $δ_{in}\sim(t-t_r)^{1/3}$ and after reconnection $δ_{out}\sim(t-t_r)^{2/3}$, which are respectively slower and faster than the dimensional prediction. In these cases, the reconnection takes place in a plane defined by the directions of the curvature and vorticity. To define the structure further, lines are drawn that connect the four arms that extend from the reconnection plane, four angles $θ_i$ between these arms are found, then summed, giving $\sumθ_i>360^\circ$. This implies that the overall structure is convex or hyperbolic, as opposed to the acute angles of the anti-parallel pyramid.

研究动机与目标

• 解决量子流体中量子涡旋重连动力学缺乏直接微观洞察的问题。
• 应对在涡旋重连模拟中精确追踪涡旋核心位置并最小化数值波动的挑战。
• 探究不同初始构型（反平行与正交涡旋）是否导致重连过程中涡旋间距的标度律存在差异。
• 探索正交涡旋重连中观测到的异常标度律的几何与运动学起源。
• 将观测到的标度行为与 Frenet-Serret 框架对齐演化及曲率动力学联系起来。

提出的方法

• 使用分裂步谱方法对三维 Gross-Pitaevskii 方程进行数值积分，以模拟量子流体动力学。
• 通过波函数梯度构建伪涡度场，以在三维网格上精确确定涡旋核心位置。
• 采用定制的初始条件，使涡旋核心周围的密度分布保持平衡，以抑制虚假的时间密度波动并稳定标度行为。
• 计算涡旋线轨迹、曲率及 Frenet-Serret 框架，以分析重连过程中的局部几何与运动学特性。
• 通过在最近距离处曲率与涡度方向的平均值定义重连平面，并分析延伸涡旋臂之间的夹角。
• 比较反平行与正交初始构型下的标度律，重点关注最小间距 δ(t) 的时间演化。

实验结果

研究问题

• RQ1 不同的初始量子涡旋构型（特别是反平行与正交）是否会导致重连过程中最小间距的标度律存在差异？
• RQ2 为何正交涡旋重连表现出偏离经典预测 δ ∼ |tₜ − t|¹ᐟ² 的标度律？
• RQ3 Frenet-Serret 框架的几何对齐如何影响涡旋重连事件中观测到的标度行为？
• RQ4 异常标度律在不同初始间距下持续存在的程度如何？这些标度律是否可推广至宏观尺度的涡旋系统？
• RQ5 观测到的标度差异是否可与涡旋线的全局对齐结构相关联，特别是与从重连平面延伸出的四条臂所形成的总角度有关？

主要发现

• 对于初始为反平行的涡旋，重连前后分离距离 δ(t) 均遵循 δ ∼ |tₜ − t|¹ᐟ²，表明在 tₜ 处具有时间对称性，且在空间中点处具有物理对称性。
• 对于初始为正交的涡旋，重连前的标度为 δ ∼ |tₜ − t|¹ᐟ³（慢于经典预测），重连后的标度为 δ ∼ |t − tₜ|²ᐟ³（快于经典预测），表明存在异常行为。
• 正交涡旋的异常标度律并非瞬态现象，而是从早期即持续存在，且与 Frenet-Serret 框架的持续双曲对齐密切相关，延伸涡旋臂之间的夹角总和超过 360°。
• 正交涡旋的重连平面由最近距离处曲率与涡度方向的平均值定义，从该平面延伸出的四条臂构成一个凸的双曲结构。
• 观测到的标度律在不同初始间距下均具有鲁棒性，且并非数值噪声所致，这一点通过使用核心平衡初始条件抑制了非期望波动得到验证。
• 结果表明，正交涡旋的异常标度行为可能适用于宏观尺度，因为几何对齐在整个重连过程中持续存在，且与初始间距无关。