其中 $a$ 为玻尔半径。由于在原子单位制下 $a=1$，因此 $r_s = L/(2N)$，即平均粒子间距的一半。

关键部分在于相互作用势 $V(x_{ij})$ 的定义。为了处理周期性边界条件下的长程库仑势，必须采用一维埃瓦尔德（Ewald）求和技术。一维 Ewald 势描述了一个电子与另一个电子的所有周期性映像、以及与之相伴的背景电荷份额的相互作用，其数学表达式为：

此势能函数精细地消除了库仑力在一维 PBC 下的发散问题。常数项 $V_M$ 是所谓的马德隆（Madelung）能量，代表每个粒子与其自身映像及均匀背景电荷之间的净相互作用能。在实际数值模拟中，为了提高计算效率，通常使用三次样条函数（Cubic Spline Interpolation）对预先在极细密网格上计算好的 $V(x_{ij})$ 值进行在线插值。

1.3 技术难点：量子蒙特卡洛中的费米子负号问题在 1D 中的特殊消解

量子蒙特卡洛（QMC）模拟高维（2D 及 3D）费米子体系时，面临的最大技术障碍是费米子负号问题（Fermion Sign Problem）。当我们在虚时间路径积分表象下对费米配分函数进行展开时，粒子轨道的反对称交换会导致配分函数的加权项出现正负交替，从而导致统计方差随着粒子数 $N$ 和逆温度 $\beta = 1/k_B T$ 的增加呈指数级爆炸。为了规避此问题，传统高维 DMC 模拟通常采用“固定节点近似（Fixed-Node Approximation）”，但这会引入难以定量评估的变分偏差。

然而，在一维极限下，情况发生了根本性的变化。由于一维空间的拓扑局限性，粒子被严格限制在一条直线上。如果粒子之间存在强烈的短程库仑排斥（$V(r) \to \infty$ 当 $r \to 0$），那么粒子在直线上是绝对无法相互穿越的。这意味着，在虚时间演化过程中，任何粒子轨道的交叉都会因为极高的势垒而被完全抑制。在路径积分展开中，只有恒等排列（Identity Permutation）对配分函数有实质性贡献。

这一物理本质带来了两个重大方法学突破：

1. 负号问题彻底消失：一维 Fermi OCP 的 QMC 模拟不存在任何负号问题。我们可以实现真正数学意义上的“精确（Exact）”和“无偏（Unbiased）”计算。
2. 量子统计的独特性：由于恒等排列占绝对主导，体系在热力学量（如基态能量、对关联函数 $g(x)$、静态结构因子 $S(q)$）上表现得与玻色子或完全可分辨的“玻尔兹曼子（Boltzmannons）”完全相同。这些“对角（Diagonal）”性质完全独立于体系的自旋极化率 $\zeta$。然而，对于“非对角（Off-Diagonal）”性质，如单粒子密度矩阵 $n(x)$ 和动量分布 $n(q)$，量子统计效应和波函数的反对称性依然扮演着决定性角色，会展现出强烈的自旋极化依赖性。

1.4 方法细节：正则蠕虫算法与基态外推

为了获取完全无偏的物理基态性质，本工作采用了有限温度连续空间蠕虫算法（Continuous-space Worm Algorithm）的正则系综变体。相较于扩散蒙特卡洛（DMC），该方法有以下核心优势：

• 消除波函数偏差：DMC 通常需要引入人工设计的试探波函数（Trial Wave Function）进行引导，这会导致对角和非对角物理量的估计（特别是使用混合估计器 mixed estimators 时）带有人为偏置。蠕虫算法纯粹基于配分函数的虚时间路径积分展开，不依赖任何试探波函数。
• 精确的有限温到基态外推：通过在足够低的温度下进行模拟来逼近基态。在标准的 TLLT 理论中，当系统长度 $L \to \infty$ 且温度 $T \to 0$ 时，所有物理平均值都应该唯一取决于乘积 $LT$。虽然 1D OCP 偏离了 TLLT，但作者谨慎地利用 $LT = \text{const}$ 作为标度律来指导数值收敛，通过在固定 $LT$ 的前提下增加粒子数 $N$（由于密度恒定，增加 $N$ 等价于增大 $L$）和降低温度 $T$，成功将热力学量外推至零温极限。

在具体模拟实现中：

• 使用了高阶复合分解方案（第四阶传播子），确保虚时间步长误差得到极大抑制。对于所有关注的物理量，当虚时间步长 $\tau$ 逼近以下极限时：

热力学平均值达到完美的数值收敛。

2. 关键 Benchmark 体系、计算所得数据与性能数据分析

2.1 基态能量的验证（Benchmark vs. DMC）

为了验证数值方法的正确性并确立计算精度，作者首先在全极化（$\zeta=1$）体系下，将本方法得到的基态每粒子能量与经典的扩散蒙特卡洛（DMC）高精度外推数据（Lee & Drummond, Ref [28]）进行了直接对比。具体对比数据整理于下表：

表 1：一维 Fermi OCP 每粒子基态能量（单位：Hartree）对比

数据科学解读：

1. 惊人的一致性：在极宽的密度范围（$1 \le r_s \le 20$）内，两套独立方法得到的基态能量在统计误差范围内高度契合。这强有力地确立了正则蠕虫算法在处理一维长程库仑相互作用时的物理正确性。
2. 误差来源分析：DMC 计算由于采用了极好的试探波函数并进行了昂贵的系统尺寸外推（最高达 $N=100$ 粒子），其能量统计误差显著更小。而本工作并未刻意追求对能量的极限精度，而是将计算重心放在了无试探波函数偏置的结构与动量关联函数上，因为这些物理量在 DMC 中极易受到试探波函数的污染。

2.2 基态对关联函数 $g(x)$ 的尺寸坍缩与 Wigner 晶体行为

图 2 现象解析：标度律的检验

研究展示了在 $r_s = 2$ 处、两个不同系统尺寸（$N=256$ 与 $N=512$）和相应温度（满足 $LT = \text{const}$）下的对关联函数 $g(x)$。数据点在整条曲线上达成了完美的重合（Collapse）。这一现象极具物理启发性：即使在标准 TLLT 标度律理论失效的一维 OCP 体系中，在低至当前模拟的温度下，$LT = \text{const}$ 依然是一个极其优秀且实用的基态外推物理向导。

图 3 现象解析：拟晶态序的涌现

分析其峰值结构可以得出以下结论：

• 振荡幅度的增强：随着 $r_s$ 的增加（即体系密度降低，库仑排斥能相对于动能占绝对主导），$g(x)$ 的第一峰峰值急剧升高，且振荡向空间更远处延伸。在 $r_s = 20$ 时，振荡特征即使在 $x/r_s > 14$ 处依然极其显著，这明确预示着一维强关联极限下，体系高度趋向于形成拟晶态 Wigner 晶体。
• 高密度极限下的晶格倾向：即便在 $r_s=1$ 和 $2$ 的高密度（弱相互作用）极限下，虽然振荡幅度有所衰减，但 $g(x)$ 的特征周期性起伏依然清晰可见，这从数值上支持了 Schulz 提出的“一维 OCP 在任意密度下均存在拟晶态序”的预言。

2.3 静态结构因子 $S(q)$ 的发散特征

为了定量表征一维空间中的准长程序，作者直接在 QMC 模拟中计算了静态结构因子 $S(q)$。其关键物理现象展示在论文的图 4 中：

• 布拉格峰在 $q = 2q_F$ 处的涌现：在全极化系统中，当波矢 $q$ 对应于费米波矢的倍数 $2q_F$（其中 $q_F = \pi / (2r_s)$）时，$S(q)$ 展现出极其尖锐的峰值。
• 尺寸发散性：随着模拟体系尺寸从 $N=128$、 $N=256$ 递增到 $N=512$，$q = 2q_F$ 处的峰值高度单调递增，表现出明显的发散趋势。这说明该准晶态序在热力学极限下是切实存在的。由于一维不可避免的量子涨落，真实的“长程晶态序”被禁止，但这种“幂律衰减”的准长程序（Quasi-long-range order）得到了 QMC 的最直接确证。
• 高密度下的生存：即使在 $r_s = 2$ 的高密度下，虽然峰值生长的速度变慢，但在 $q=2q_F$ 处的峰值依然随系统尺寸的增加而清晰地发散（图 4 底部）。这表明一维长程库仑作用成功阻碍了量子熔化，将拟晶态稳定在了极高密度区域。

2.4 单粒子密度矩阵 $n(x)$ 与动量分布 $n(q)$：自旋极化的独特烙印

作为非对角可观测物理量，$n(x)$ 和 $n(q)$ 直接暴露出量子统计效应。这也是本工作超越前人研究（Ref [28] 仅考虑了全极化 case）的最重要亮点。

图 5 与 图 6：全极化极限（$\zeta=1$）

• 单粒子密度矩阵 $n(x)$：在高密度 $r_s = 1$ 下，$n(x)$ 展现出非常清晰的、由费米面切断引起的 Friedel 振荡；而在极低密度 $r_s = 20$ 下，Friedel 振荡几乎完全消失，且 $n(x)$ 在 $x \gtrsim 2r_s$ 处开始转为负值，反映出强相互作用导致的局域化效应。
• 动量分布 $n(q)$：在 $r_s = 2$ 时，$n(q)$ 在 $q = q_F$ 处表现出一个极其陡峭的“台阶（Step）”，这非常类似于无相互作用的费米气体，带有鲜明的弱关联费米子统计特征。然而，当 $r_s$ 增大至 $20$ 时，由于库仑排斥导致的强关联作用，大量粒子被激发至高动量态，导致 $q = q_F$ 处的台阶被彻底抹平，动量分布呈现出一条平滑、宽广且具有长尾（Slowly decaying tail）的非费米液体曲线。令人惊叹的是，尽管空间排列上晶格序极度明显（Wigner 晶体），但 $n(q)$ 曲线本身平滑无奇，完全无法体现空间局域化的晶格结构。

图 7：不同自旋极化状态下的动量分布调控

极具启发性的物理现象发生在一维无极化（$\zeta=0$）和部分极化（$\zeta=0.5$）体系中：

1. 无极化与反铁磁对齐（$\zeta=0$）的“经典化”行为：
   • 在 $\zeta=0$ 且体系呈现完全反铁磁交替排列（即每个电子的近邻均为相反自旋）时，$n(q)$ 的量子统计特征（如费米台阶）几乎完全消失（图 7 顶部）。
   • 物理阐释：在一维强排斥体系中，相反自旋的费米子由于近邻无法交叉，它们在物理上可以被视为完全可分辨的粒子（Distinguishable particles / Boltzmannons）。这种“去相干”机制使得体系的动量分布退化为近乎经典的分布。这一结果完美证实了 Girardeau 的一维可分辨子映射理论。
2. 部分极化（$\zeta=0.5$）的多数-少数组分非对称性：
   • 图 7 底部给出了主组分（占 75% 粒子数）与次组分（占 25% 粒子数）的动量分布。
   • 物理阐释：由于主组分粒子数较多，在直线上不可避免地会经常遇到同自旋的近邻粒子，因此为了满足 Pauli 不相容原理，同自旋粒子间的量子交换（Exchange）极为频繁，导致主组分的 $n(q)$ 依然保留了强烈的量子费米子特征。相反，次组分的粒子在空间中极稀疏，其周围几乎全被相反自旋的主组分粒子包围。由于无法与近邻（相反自旋）发生量子交换，次组分的量子特征几乎被完全剥夺，其 $n(q)$ 曲线表现出近乎经典的宽缓平滑分布。这一发现深刻揭示了一维多体物理中量子统计与自旋构型的微妙共生关系。

3. 代码实现细节、复现指南与开源工具推荐

为了方便科研人员在自己的超级计算机或工作站上重现 Massimo Boninsegni 教授的物理发现，本节将详细梳理该工作的模拟算法架构，并提供一条可行的复现指南。由于本项工作采用的正则系综连续空间蠕虫算法（Canonical Worm Algorithm）多为课题组自研的高效 Fortran/C++ 代码，我们在此给出算法的数学逻辑实现细节及可替代的开源软件资源。

3.1 正则蠕虫算法的核心虚时间演化更新步骤

连续空间路径积分蒙特卡洛（PIMC）将配分函数 $Z = \text{Tr}(e^{-\beta \hat{H}})$ 写为沿虚时间轴 $\tau$ 离散化的路径积分。由于我们要模拟固定粒子数 $N$ 的正则系综，在处理非对角物理量（如 $n(x)$）时，必须允许路径在虚时间中临时“打开”（形成一条带有两个端点 $M$ 和 $I$ 的“蠕虫”路径）。

          [ 路径在虚时间 τ 中的演化示意图 ]

τ = β  o=====o=====o=====o=====o=====o=====o   (对角构型: 所有环闭合)
     |     |     |     |     |     |     |
τ = 0  o=====o=====o=====o=====o=====o=====o

                     ⬇ (Worm Open 更新)

τ = β  o=====o     o=====o=====o=====o=====o   (非对角构型: 产生端点 M 和 I)
     |     |     /   |     |     |     |
     |     |   M     |     |     |     |
     |     |         |     |     |     |
     |     |       I |     |     |     |
τ = 0  o=====o     o=====o=====o=====o=====o

复现该算法需要实现以下六种基础蒙特卡洛更新（Monte Carlo Updates）：

1. Open（打开）：随机选择一个闭合粒子路径环，在虚时间区间内将其切断，引入两个活跃端点：头部（Head, $M$）和尾部（Tail, $I$）。该步允许体系进入非对角（Worm）空间。
2. Close（闭合）：当头部 $M$ 和尾部 $I$ 在空间和虚时间上足够接近时，尝试通过桥接两个端点，将蠕虫路径重新闭合为正常的粒子环。
3. Insert（插入）：在空的物理空间中尝试插入一条全新的、带有两个端点的短蠕虫路径段，以增加物理采样几率。
4. Remove（移除）：Close 的逆过程，直接尝试擦除一条短蠕虫路径段。
5. Advance / Retreat（向前推进 / 向后收缩）：通过在虚时间上向前或向后对端点进行路径重构（通常结合热浴或 Levy 飞跃技术），从而在空间中快速移动端点。
6. Swap（交换）：当蠕虫端点移动到另一个粒子路径附近时，尝试重组路径连接，从而隐式地进行粒子轨道置换。在一维 OCP 中，此交换步骤由于巨大的两粒子排斥势垒，在热力学极限下成功率为 0。

3.2 1D Ewald 势的预计算与三次样条插值（Spline Interpolation）复现要点

一维 Ewald 势（Eq. 2）的无穷项求和是计算量的瓶颈，绝不能在蒙特卡洛每一步中实时计算。标准的工程复现路径如下：

• 第一步：建立一维超细密网格。网格范围为 $[0, L/2]$。网格步长建议设为 $\Delta x \le 10^{-4} a.u.$。
• 第二步：计算网格点上的 Ewald 实空间项和倒空间项之和。在网格点上精确求和（Eq. 2），计算时截断数 $n$ 至少取到使相邻项贡献低于双精度浮点数极限（$10^{-16}$）为止。
• 第三步：拟合三次样条函数。计算每个网格区间的样条系数 $A_i, B_i, C_i, D_i$，使得在模拟中对于任意粒子间距 $x_{ij}$，势能和受力均可通过下式极速获取：

• 第四步：扣除发散项。在 $x_{ij} \to 0$ 处，库仑发散必须正确截断。本工作在路径积分中自然通过极小的热德布罗意波长保障了粒子不重叠，但在样条插值最左端点仍需设置极小截断半径（例如 $r_{\text{cut}} = 10^{-6}$）。

3.3 零温极限（基态）外推方案复现

1. 标度约束：将模拟的长度 $L$（与粒子数 $N$ 成正比）与逆温度 $\beta = 1/T$ 绑定。设定一个常数标度参数，如：

2. 模拟序列：
   • 跑一组 $N=20, T=0.0125 \text{ Ha}$ 的基态前哨计算。
   • 跑一组 $N=40, T=0.00625 \text{ Ha}$ 的计算。
   • 跑一组 $N=80, T=0.003125 \text{ Ha}$ 的计算（如果算力充足，可挑战 $N=256, 512$）。
3. 结果判定：如果对于给定的物理量（如每粒子能量或 $g(x)$），这三组模拟的数据在统计误差范围内完全重合（图 1 与 图 2），则宣告该标度条件成功捕获了基态热力学极限性质。

3.4 推荐开源软件包与 GitHub 资源

虽然 Massimo Boninsegni 教授的个人代码没有完全开源，但社区中有几个极佳的开源 PIMC 和蠕虫算法框架可用于直接搭建本复现任务：

1. PIMC++ (Path Integral Monte Carlo in C++)

   • 项目简介：由 Urbana-Champaign 社区（Ceperley 组及合作者）维护的经典三维/一维连续空间路径积分蒙特卡洛程序包。支持 Ewald 求和以及大粒子数体系。
   • 适合任务：极适合用于重现 1D OCP 的对角性质（基态能量、$g(x)$、$S(q)$）。
   • Repo 链接检索：可在 GitHub 上检索 pimc 或访问 UIUC 凝聚态计算物理相关开源页面。

2. pi-qmc (Python Path Integral Quantum Monte Carlo)

   • 项目简介：一个轻量级、高度模块化的 Python / NumPy 连续空间 PIMC 实现。非常适合科研人员进行算法原理快速验证。
   • 适合任务：通过其 Worm 更新模块，可以极为直观地复现图 7 中自旋极化对 $n(q)$ 的调控行为。

4. 关键引用文献与评论

4.1 核心引用文献深度解读

本工作立足于一维多体物理数十年发展的坚实基础之上，以下四篇关键文献构成了其理论与方法学的核心支柱：

1. Ref [13] T. Giamarchi, Quantum Physics in One Dimension (Oxford University Press, 2003)

   • 地位：一维量子多体系统的圣经级著作。
   • 关联：详尽阐述了标准 Tomonaga-Luttinger 液体理论（TLLT）的数学构架。本工作研究一维 OCP，其根本目的就在于探索体系在何种程度上背离了 Giamarchi 书中经典的 TLLT 预言，特别是其独特的激发谱对数修正形式。

2. Ref [14] H. J. Schulz, Wigner crystal in one dimension, Phys. Rev. Lett. 71, 1864 (1993)

   • 地位：首次从解析上论证了一维无屏蔽库仑系统在任意密度下都具有维格纳晶格倾向的重要文献。
   • 关联：本工作模拟得到的对关联函数 $g(x)$ 以及静态结构因子 $S(q)$ 在极低温度和超大尺寸下的发散行为，是迄今为止对 Schulz 一维维格纳晶体假说最直接、最无偏的数值实证。

3. Ref [28] R. M. Lee and N. D. Drummond, Ground-state properties of the one-dimensional electron liquid, Phys. Rev. B 83, 245114 (2011)

   • 地位：前人最高精度的 1D 费米子 DMC 模拟工作。
   • 关联：该文利用传统的变分和扩散蒙特卡洛（VMC/DMC）获取了全极化一维电子气的能量及结构性质。本工作将其作为最重要的 Benchmark 对手。本工作的一大核心贡献在于：证明了即使不依赖 Lee & Drummond 文章中复杂的试探波函数，蠕虫算法也能在完全无偏的条件下，输出完全吻合且物理信息更为丰富的多体物理图像，并将研究范围推进到了非全极化区域。

4. Ref [38, 39] M. Boninsegni, N. Prokof’ev, and B. Svistunov, Worm algorithm for continuous-space path integral Monte Carlo simulations, Phys. Rev. Lett. 96, 070601 (2006)

   • 地位：连续空间蠕虫算法的奠基性文献。
   • 关联：该方法彻底攻克了传统 PIMC 无法高效采样大尺度宏观量子置换和非对角关联的缺陷。本工作正是该算法架构在一维费米长程系统中的一次成功延伸。

4.2 技术评论：本工作的局限性与潜在挑战

作为本领域的资深技术评论者，我们必须在赞赏本工作无偏、高精度特点的同时，冷静地指出其潜在的方法学局限性，以便同行科研人员后续攻关：

1. 标度外推假定的局限性： 作者在进行基态外推时，极度依赖于 $LT = \text{const}$ 这一标度关系。然而，正如作者本人在引言中所强调的，一维 OCP 并不遵循标准的 TLLT 理论。因此，在物理本质上，一维 OCP 在热力学极限（$L \to \infty, T \to 0$）下的严格零温物理量，可能存在偏离 $LT$ 标度的细微修正项。虽然在当前数值模拟精度（误差棒约在 $10^{-4}$ 数量级）内没有观察到这种偏离，但未来若要进行更高精度的量子化学交换关联能高阶拟合，这一标度假定的精确度仍需受到更严苛的检验。

2. 计算复杂度的温度惩罚： 由于需要模拟 $T \to 0$ 的物理基态，虚时间路径积分的“时间片”（Beads）数目 $M = \beta / \tau$ 会随着温度的降低而呈线性暴增。特别是在强关联极限 $r_s = 20$ 下，为了抑制时间片离散误差，$\tau$ 必须取得极小，这导致时间片数目轻松突破数万。此时，蠕虫算法路径更新的自关联时间会急剧变长，在大尺寸（如 $N=512$）时面临着严峻的计算资源瓶颈，难以向更低的强关联温度（如 $T < 10^{-4} \text{ Ha}$）继续推进。

3. 激发谱和谱函数的缺失： 由于虚时间路径积分 QMC 的物理天性，该方法无法直接获取实时间（Real-time）动力学性质。对于一维 OCP 最核心的理论元激发 $\omega_q \sim q\sqrt{\ln(1/q)}$，本工作无法直接给出动力学结构因子 $S(q, \omega)$ 的精确谱线。虽然可以通过诸如最大熵方法（Maximum Entropy Method）对虚时间虚关联函数进行解析延拓，但众所周知，解析延拓是一个极度不稳定的反问题，难以给出具有化学精度和尖锐线宽的激发谱。这一缺陷仍需结合精确的含时密度矩阵重整化群（t-DMRG）方法来弥补。

5. 补充与拓展探讨

为了帮助量子化学和凝聚态物理背景的读者建立更广阔的物理视野，我们在此对本工作涉及的内容进行三维维度的纵深拓展探讨。

5.1 维格纳结晶行为：1D 与 2D、3D 的本质差异

维格纳结晶（Wigner Crystallization）是电子气在低密度、强排斥极限下最令人着迷的相变现象。然而，在不同维度下，该现象背后的物理机制和表现形式有着天壤之别：

• 三维（3D）与二维（2D）电子气： 在 3D 和 2D 体系中，维格纳结晶是一个典型的一阶量子相变。当体系密度降低到某一临界点（例如 3D 中 $r_s \approx 106$，2D 中 $r_s \approx 37$）时，电子动能（自能）已不足以对抗强烈的库仑排斥，体系会自发对称性破缺，形成具有长程晶格序的体心立方（BCC）或三角点阵晶体。此时，静态结构因子 $S(q)$ 在倒格子矢量处会表现出不依赖于系统尺寸的、真正的 $\delta$ 函数布拉格尖峰。
• 一维（1D）电子气： 根据 Mermin-Wagner 定理，在一维多体系统中，任何连续的平移对称性在有限温度下都无法被自发破缺；即使在零温下，由于极强的量子涨落，真实的无限长程晶态序也是被禁止的。然而，一维 OCP 的特殊性在于，库仑相互作用一维无屏蔽特性过于强大，导致涨落受到极大抑制。因此，一维 OCP 发展出了一种拟晶态序（Quasi-crystalline order），其密度关联函数 $g(x)$ 表现出超平缓（幂律）衰减。本工作直接展示了，无需极低密度（如 3D 的 $r_s > 100$），哪怕在 $r_s=2$ 这样的高密度下，一维 OCP 就已经顽强地展现出了这种拟晶态特征。这说明在低维多体系统中，几何局限大大强化了局域化效应。

5.2 一维多体系统中的自旋-电荷分离（Spin-Charge Separation）

在一维强关联物理中，自旋-电荷分离是最具革命性的概念之一。对于满足标准 TLLT 的短程排斥系统：

• 体系的集体激发可以彻底解耦为两个完全独立的通道：电荷密度波（等离激元）通道与自旋密度波（激子）通道。
• 这两套通道具有各自独立的激发速度 $v_c$ 和 $v_s$，通常 $v_c > v_s$。当一个物理电子注入系统时，它会瞬间“分裂”为一个携带电荷的 holon 和一个携带自旋的 spinon，并以不同的速度在空间传播。

在一维 OCP 中，由于长程库仑作用的介入：

• 电荷通道：由于排斥力极长，其激发能谱在一维极限下被强力“抬升”，形成了具有对角数修正的准能隙能谱 $\omega_q \sim q\sqrt{\ln(1/q)}$。电荷运动表现出等离子体振荡特征。
• 自旋通道：由于自旋交换属于电中性扰动，长程库仑相互作用在自旋通道中不扮演角色。自旋激发依然保持着完好的标准 TLLT 线性激发特征（即自旋波速度 $v_s$ 保持恒定，且满足 $K_s = 1$ 如果具有自旋旋转对称性）。

本工作在图 7 底部展现的部分极化体系 $n(q)$，实际上就是这种自旋-电荷分离在基态动量分布中的隐性物理投射。主组分和次组分展现出截然不同的量子统计完整性，正是因为它们在自旋通道的局域相互作用与电荷通道的长程主导排斥之间，达成了一种奇特的动力学不平衡。

5.3 实验现实化：碳纳米管、量子线与冷原子光晶格的桥梁

尽管一维 OCP 听起来像是一个理想化的纯理论模型，但在当今纳米科技和超冷原子物理的发展下，它已经拥有了非常坚实的实验对应平台：

1. 半导体量子线（GaAs/AlGaAs Quantum Wires）： 通过高精度的分子束外延（MBE）和纳米刻蚀技术，可以将电子运动严格限制在一条一维沟道中。由于衬底介电常数的屏蔽作用，电子间的实际相互作用形式介于 $1/r$ 势与短程势之间，是检验本工作拟晶态理论的理想固体物理平台。
2. 单壁碳纳米管（SWCNTs）： 由于其超高的长径比，电子在碳纳米管表面的轴向运动是极其完美的一维费米气体行为。实验上通过输运测量和拉曼光谱观测到的电荷不均匀重构，与本工作在强关联下计算得到的 $g(x)$ 密度振荡特征高度吻合。
3. 冷原子一维光晶格（1D Optical Lattices with Cold Polar Molecules）： 物理学家现在可以通过调整激光干涉图案，将带有永久偶极矩的极性分子（如 KRb）锁定在一维管道中。利用 Feshbach 共振技术，可以自由调节粒子间的相互作用强度，从而在一维极限下人工模拟从弱排斥到极强排斥的完整 $r_s$ 区间。本工作给出的无偏动量分布 $n(q)$ 数据，可作为冷原子飞行时间质谱（Time-of-Flight Spectroscopy）实验进行基态量子统计检验的标准理论参考线。

5.4 总结：对量子化学交换关联能（XC Functional）开发的启示

对于量子化学和密度泛函理论（DFT）研究人员而言，一维等离子体多体能量计算并非纸上谈兵。当代 DFT 的核心成功在于本地密度近似（LDA）和广义梯度近似（GGA），而这些近似泛函中对于三维和二维电子气交换关联能（$\epsilon_{xc}(r_s)$）的拟合，完全依赖于 Ceperley & Alder 1980 年著名的多体 QMC 计算数据。

当前，一维微纳器件（如分子导线、单分子晶体管）的量子化学计算需求日益暴增，但市面上大部分主流商用量子化学软件直接套用三维 LDA 泛函处理一维体系会产生灾难性的系统误差（如高估离域、低估能隙）。Massimo Boninsegni 教授在这项工作中给出的精确到小数点后第 4 位的无偏基态能量、静态结构因子和非对角密度矩阵数据，正是未来开发**专用一维交换关联泛函（1D Density Functional Theory）**梦寐以求的精准数据库，将直接助力低维纳米材料电子结构计算精度的质的飞跃。