[论文解读] Dimer lambda_d Expansion Computer Computations

Paul Federbush|ArXiv.org|Apr 26, 2008

Theoretical and Computational Physics参考文献 3被引用 3

一句话总结

本文通过符号计算在 Maple 中实现了对 dimer $\lambda_d$ 展开式按 $1/d$ 的高精度计算机计算，将渐近级数扩展至 $1/d^3$ 项。作者计算了 $d=1,2,3$ 时的归一化系数 $\bar{J}_i$，验证了先前的手工计算结果，并推导出新的 $1/d^3$ 项；同时数值研究了 $\lambda_d$ 的两种近似序列，发现其在 $d=2$ 和 $d=3$ 时均无明显收敛性。

ABSTRACT

In a previous paper an asymptotic expansion for lambda_d in powers of 1/d was developed. The results of computer computations for some terms in the expansion, as well as various quantities associated to the expansion, are herein presented. The computations (in integer arithmetic) are actually done only for d=1, d=2, and d=3, but some results for arbitrary dimension follow from the general structure of the theory. In particular we obtain the next term in powers of 1/d, the 1/d^3 term in the asymptotic expansion for lambda_d, as well as checking the correctness of the terms previously calculated by hand. For d=2 and d=3 a second expansion for lambda_d is defined and numerically studied. Contrary to our hopes this second type of expansion appears not to be convergent. There is not enough evidence to argue if the expansions for lambda_d in powers of 1/d are convergent for either d=2 or d=3, but likely they are not. (Certainly convergence, but at a very slow rate, remains an interesting possibility in all our cases.)

研究动机与目标

• 通过计算机代数方法计算 $\lambda_d$ 在 $1/d$ 幂级数中的高阶项。
• 验证并扩展先前手工计算的 $\lambda_d$ 展开式中 $d=1,2,3$ 时的 $\bar{J}_i$ 系数。
• 研究 $d=2$ 和 $d=3$ 时 $\lambda_d$ 的两种不同近似序列的收敛行为。
• 开发并实现一个四步计算流程，以高效计算 $\bar{J}_i$ 至 $i=6$。

提出的方法

• 使用 Maple 进行符号计算，计算 $i=1$ 到 $6$ 以及 $d=1,2,3$ 时的 $\bar{J}_i = \frac{1}{N} \cdot \frac{1}{i!} \cdot J_i$。
• 计算过程分为四个步骤：生成具有重叠模式的 $s$-dimer 树（步骤1），计算连通瓷砖和（步骤2），根据瓷砖重叠计算 $J_s$（步骤3），并将结果代入主公式以获得近似值（步骤4）。
• 定义 $s$-dimer 树为一组有序的 $s$ 个定位瓷砖，使得每个前缀子序列均保持连通，从而实现高效的算法枚举。
• 该方法使用整数算术，避免浮点运算，以确保 $\bar{J}_i$ 符号求值的精度。
• 构建了两种近似序列：$A_r$ 基于直接的 $1/d^i$ 项，$B_r$ 则使用形式变量 $x$ 来建模生成函数。
• 计算 $d=2$ 和 $d=3$ 时的 $\bar{J}_6$ 耗时数周 CPU 时间，而 $\bar{J}_7$ 在当前方法下被认为不可行。

实验结果

研究问题

• RQ1 在 $\lambda_d$ 的渐近展开中，$1/d$ 的下一个项是什么，特别是 $1/d^3$ 的系数？
• RQ2 先前手工计算的 $\lambda_d$ 展开式中的 $\bar{J}_i$ 系数是否准确？能否通过计算机代数方法验证？
• RQ3 针对 $d=2$ 和 $d=3$ 的两种近似序列 $A_r$ 和 $B_r$ 是否收敛？
• RQ4 能否从 $d=1,2,3$ 时计算出的 $\bar{J}_i$ 值中推导出任意维度 $d$ 下 $\bar{J}_i$ 的通解形式？
• RQ5 将 $\bar{J}_i$ 的计算扩展至 $i=6$ 以上时，$d=2$ 和 $d=3$ 的计算复杂度如何？

主要发现

• $\lambda_d$ 的渐近展开被扩展至包含 $1/d^3$ 项：$\lambda_d \sim \frac{1}{2}\ln(2d) - \frac{1}{2} + \frac{1}{8d} + \frac{5}{96d^2} + \frac{5}{64d^3} + \cdots$
• $d=1,2,3$ 时的 $\bar{J}_i$ 系数被精确计算：$\bar{J}_1 = 0$，$\bar{J}_2 = \frac{1}{8d}$，$\bar{J}_3 = \frac{1}{12d^2}$，$\bar{J}_4 = -\frac{3}{32d^2} + \frac{3}{64d^3}$，$\bar{J}_5 = -\frac{1}{8d^3} - \frac{3}{80d^4}$，$\bar{J}_6 = \frac{7}{48d^3} - \frac{5}{64d^4} - \frac{1}{6d^5}$。
• 对于 $d=2$，$A_r$ 序列给出 $A_3 = 0.2784$，趋近于已知精确值 $\lambda_2 \approx 0.29156$，而 $B_r$ 序列在 $B_2 = 0.2921$ 处达到峰值后下降。
• 对于 $d=3$，$A_r$ 序列给出 $A_3 = 0.4462$，处于已知范围 $0.440075 \leq \lambda_3 \leq 0.457547$ 内，而 $B_r$ 序列在 $B_2 = 0.4538$ 处达到峰值后下降。
• $d=2$ 和 $d=3$ 的 $B_r$ 序列均未收敛，表现出振荡并最终下降，表明该形式展开可能不收敛。
• $d=2$ 和 $d=3$ 时 $\bar{J}_6$ 的计算耗时数周 CPU 时间，表明 $\bar{J}_7$ 的计算可能需要约一年时间，且需重大新理论洞察。