[论文解读] Superharmonic double-well systems with zero-energy ground states: Relevance for diffusive relaxation scenarios

Piotr Garbaczewski, V. A. Stephanovich|arXiv (Cornell University)|Apr 24, 2021

Quantum chaos and dynamical systems参考文献 53被引用 1

一句话总结

该论文提出一种计算机辅助谱方法，用于近似具有零能基态的超谐双阱薛定谔算符的低能谱，该算符与受限势阱中的扩散弛豫过程相关。该方法首次系统地恢复了这些拟精确可解系统的能隙和能级，揭示了当势阱阶数增加时，谱收敛于有界区域上的诺伊曼拉普拉斯算符。

ABSTRACT

Relaxation properties (specifically time-rates) of the Smoluchowski diffusion process on a line, in a confining potential $ U(x) \sim x^m$, $m=2n \geq 2$, can be spectrally quantified by means of the affiliated Schr\"{o}dinger semigroup $\exp (-t\hat{H})$, $t\geq 0$. The inferred (dimensionally rescaled) motion generator $\hat{H}= - \Delta + {\cal{V}}(x)$ involves a potential function ${\cal{V}}(x)= ax^{2m-2} - bx^{m-2}$, $a=a(m), b=b(m) >0$, which for $m>2$ has a conspicuous higher degree (superharmonic) double-well form. For each value of $m>2$, $ \hat{H}$ has the zero-energy ground state eigenfunction $ ho _*^{1/2}(x)$, where $ ho _*(x) \sim \exp -[U(x)]$ stands for the Boltzmann equilibrium pdf of the diffusion process. A peculiarity of $\hat{H}$ is that it refers to a family of quasi-exactly solvable Schr\"{o}dinger-type systems, whose spectral data are either residual or analytically unavailable. As well, no numerically assisted procedures have been developed to this end. Except for the ground state zero eigenvalue and incidental trial-error outcomes, lowest positive energy levels (and energy gaps) of $\hat{H}$ are unknown. To overcome this obstacle, we develop a computer-assisted procedure to recover an approximate spectral solution of $\hat{H}$ for $m>2$. This task is accomplished for the relaxation-relevant low part of the spectrum. By admitting larger values of $m$ (up to $m=104$), we examine the spectral "closeness" of $\hat{H}$, $m\gg 2$ on $R$ and the Neumann Laplacian $\Delta _{\cal{N}}$ in the interval $[-1,1]$, known to generate the Brownian motion with two-sided reflection.

研究动机与目标

• 解决在具有零能基态的超谐双阱系统中计算能谱的分析与数值方法缺失的问题。
• 开发一种计算机辅助程序，用于近似薛定谔算符 $ \hat{H} = -\Delta + V(x) $ 的低能谱，其中 $ V(x) = ax^{2m-2} - bx^{m-2} $，$ m > 2 $。
• 研究当 $ m \to \infty $ 时 $ \hat{H} $ 的谱行为，特别是其在 $[-1,1]$ 上收敛于诺伊曼拉普拉斯算符的性质。
• 为研究由此类势阱控制的扩散过程中弛豫速率提供一个框架，其中标准的WKB方法和瞬子方法失效。

提出的方法

• 通过分解 $ \rho(x,t) = \Psi(x,t) \rho^{1/2}_*(x) $，将斯莫卢奇夫斯基-福克-普朗克动力学转化为 $ L^2(\mathbb{R}) $ 中的薛定谔型问题，得到方程 $ \partial_t \Psi = \Delta \Psi - V(x)\Psi $。
• 定义有效势 $ V(x) = \frac{1}{2} \left( \frac{b^2}{2} + \nabla b \right) $，其中 $ b(x) = -\nabla U(x) $，$ U(x) = \kappa_m x^m / m $，当 $ m > 2 $ 时，该势呈现超谐双阱形式。
• 采用数值射击法结合谱配点法，近似 $ \hat{H} $ 的特征值与特征函数，重点关注低能谱。
• 通过无量纲化处理，设定扩散系数 $ D = 1 $，将谱问题简化为 $ \hat{H} = -\Delta + V(x) $。
• 将 $ \hat{H} $ 在大 $ m $ 时的谱与 $[-1,1]$ 上的诺伊曼拉普拉斯算符 $ \Delta_N $ 进行比较，建立在 $ m \to \infty $ 极限下的谱收敛性。
• 通过收敛性检验及与已知精确解（当 $ m = 4 $ 时可用）的对比，验证结果的可靠性。

实验结果

研究问题

• RQ1 对于 $ m > 2 $ 的超谐双阱势，是否存在一种计算机辅助方法可恢复 $ \hat{H} $ 的低能谱，其中解析解不可用？
• RQ2 当势阱阶数 $ m $ 增加时，$ \hat{H} $ 的能隙行为如何，特别是在 $ m \to \infty $ 的极限下？
• RQ3 大 $ m $ 时 $ \hat{H} $ 的谱与 $[-1,1]$ 上诺伊曼拉普拉斯算符的谱有何异同？
• RQ4 为何标准方法如WKB和瞬子微积分在计算此类系统中的非零特征值时失效？
• RQ5 在扩散弛豫动力学背景下，零能基态的物理意义是什么？

主要发现

• 所提出的计算机辅助方法成功计算出 $ m $ 高达 $ 10^4 $ 时 $ \hat{H} $ 的前几个正特征值与能隙，此前在这些参数范围内尚无解析或数值工具可用。
• 随着 $ m $ 增大，$ \hat{H} $ 的最低非零特征值减小，表明高阶双阱势中的弛豫速率更慢。
• 当 $ m $ 较大时，$ \hat{H} $ 的谱收敛于 $[-1,1]$ 上的诺伊曼拉普拉斯算符 $ \Delta_N $，表明在超谐极限下存在普遍的弛豫行为。
• 零能基态 $ \rho^{1/2}_* $ 为单峰，且在 $ V(x) $ 的不稳定极大值处达到峰值，与标准双阱量子系统的直观预期相反。
• 控制指数弛豫速率的能隙 $ \lambda_1 $ 在大 $ m $ 时近似满足 $ \lambda_1 \sim 1/m^2 $，表明弛豫速度呈幂律衰减。
• 该方法证实 $ \hat{H} $ 为拟精确可解系统，仅基态特征值 $ \lambda = 0 $ 可解析获得，而更高能级通过标准技术仍不可达。