来源论文: https://arxiv.org/abs/2607.06511v1 生成时间: Jul 08, 2026 18:23

临界 J₁-J₂ 横场 Ising 链中的静态空位：作为参数化共形缺陷的深度解析

0. 执行摘要

一维量子自旋链中的杂质与缺陷效应是凝聚态物理和量子信息领域的核心课题。当体系统处于量子临界点（Quantum Critical Point, QCP）时，系统关联长度发散，局域杂质将演化为共形缺陷（Conformal Defects）。这类缺陷不仅打破了空间平移对称性，更在保持体系统共形对称性的同时，展现出独特且普适的反射与透射行为。

本研究聚焦于含有第一邻接耦合 $J_1$ 和第二邻接耦合 $J_2$ 的横场 Ising 链（TFIM），在量子临界线 $b = b_c(J_2)$ 上深入探讨了两个静态非磁性空位之间的相互作用，以及单个空位的传输特性和边界熵。通过采用高精度的密度矩阵重整化群（DMRG）算法（系统尺寸最大达到 $N = 300$ 格点），研究取得了以下突破性进展：

1. 确认了量子卡西米尔相互作用（Quantum Casimir Interaction）的代数衰减律：在整个临界线上，空位-空位结合能 $\Delta_b(r)$ 表现为幂律衰减 $|\Delta_b(r)| \sim r^{-\alpha}$。外推至热力学极限下，指数 $\alpha_\infty$ 紧邻共形场论（CFT）预测的普适值 $1.0$，并随 $J_2/J_1$ 的增加展现出从 $1.07$ 到 $1.15$ 的微弱但系统性的单调增长。
2. 刻画了部分透射共形缺陷（Partially Transmissive Conformal Defect）的连续演化：单空位结构中的自旋关联透射比 $T(k)$ 在远距离极限下饱和至非平凡的平台值 $T_\infty(J_2) \in [0.11, 0.33]$，直接证实了该缺陷并非简单的“断链”（对应 $T_\infty = 0$）或“平凡杂质”（对应 $T_\infty = 1$），而是处于一个受 $J_2$ 调控的 marginal 共形固定点上。
3. 精确提取了 Affleck-Ludwig 边界熵（Boundary Entropy）：通过引入一种新颖的共形坐标匹配减法方案，克服了开放边界条件（OBC）下的有限尺寸效应和伪尖峰佯攻，测得边界熵 $\log g_\infty \approx -0.073$。该值在整条临界线上近乎保持常数，显著高于 Ising 固定边界条件下的值 $-\ln \sqrt{2} \approx -0.347$，更接近自由边界。这表明即使是弱透射缺陷，也具备与自由边界极为接近的低纠缠熵特征。
4. 确立了体自旋关联指数 $\eta = 1/4$ 作为临界点标定的鲁棒判据：相较于极易受边界效应干扰的中心电荷 $c$ 拟合，利用 bulk 区域自旋相关函数的标度行为能够以极高精度（误差 $< 0.005$）标定临界场 $b_c$。

上述三类可观测物理量（$\alpha$, $T_\infty$, $\log g_\infty$）随第二邻接耦合 $J_2$ 呈现出平滑、单调的协同演化。这一物理图景完美契合由 $J_2$ 作为唯一控制参数、横跨强反射到部分透射的单参数共形缺陷族模型。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题

一维量子临界系统中的非磁性空位（Vacancy）究竟是作为系统的物理“割裂”（即断开成两个独立的开链），还是作为允许能量和量子相干性穿透的“半透膜”？

传统上，在系统处于有能隙（Gapped）的基态相区时，局域扰动会被快速屏蔽，空位之间的有效相互作用随着距离 $r$ 呈指数衰减。然而，在无能隙的临界点上，屏蔽效应消失，涨落处于长程关联状态。此时，空位将作为边界扰动，通过体临界涨落介导的长程量子力学力——即量子卡西米尔力——发生相互作用。理解这种卡西米尔效应的标度指数、缺陷处的透射率以及其对系统基态简并度（由 Affleck-Ludwig 边界熵表征）的贡献，是共形缺陷理论的关键所在。更进一步地，如何通过格点模型中的非微扰参数（如第二邻接耦合 $J_2$）去连续调控共形缺陷的普适类，是本研究试图回答的核心问题。

1.2 理论基础：边界共形场论（BCFT）与 Affleck-Ludwig 熵

在一维量子临界系统（或等价的二维经典统计模型）中，体系统具有共形对称性，其标度行为由特定的共形场论（CFT）描述。对于处于 Ising 临界点的自旋链，体 CFT 的中心电荷为 $c = 1/2$，对应的标度算符（如能流算符 $\epsilon$、自旋算符 $\sigma$ 等）具有确定的标度维数。特别地，自旋关联函数的衰减指数为 $\eta = 1/4$。

当在临界系统中引入局域缺陷时，系统的行为需要用**边界共形场论（BCFT）**来刻画。一个保角的缺陷（Conformal Defect）在重整化群（RG）流入极限下构成一个不流动的缺陷固定点。该缺陷的物理性质可以通过以下三个基本物理量完全表征：

A. 卡西米尔相互作用与结合能

当系统存在两个缺陷，且距离为 $r$ 时，它们之间的基态能量变化（结合能）定义为：

其中 $E_0$ 为无缺陷洁净链（Clean Chain）的基态能量，$E_{1v}$ 为含单个中心空位的链基态能量，$E_{2v}(r)$ 为包含两个相距 $r$、且对称分布于链中心两侧的空位的链基态能量。在 BCFT 的预言下，该结合能表现为代数衰减：

其中 $\alpha$ 被称为卡西米尔指数。对于理想的共形墙，卡西米尔指数通常具有普适值 $\alpha = 1$。

B. 缺陷透射率 $T_\infty$

为了量化缺陷对量子相干性的阻碍程度，引入穿越缺陷的自旋关联函数比值 $T(k)$。设空位位于格点 $v = L/2$，则在空位两侧距离为 $k$ 的格点间自旋-自旋关联函数，与洁净链中对应格点间的关联函数之比定义为：

在共形缺陷图景中，当 $k$ 远大于格点间距但远小于系统尺寸时，$T(k)$ 将趋向于一个常数平台 $T_\infty$。根据 BCFT，若缺陷是“流向断裂”（Relevant Perturbation），则在长程极限下 $T_\infty \to 0$；若缺陷是“流向愈合”，则 $T_\infty \to 1$；若 $0 < T_\infty < 1$，则说明该缺陷由一个 marginal 算符描述，对应于部分透射固定点。

C. Affleck-Ludwig 边界熵 $\log g$

对于开放边界的一维量子链，其子系统（长度为 $\ell$）的 von Neumann 纠缠熵 $S(\ell)$ 遵循 Calabrese-Cardy 公式：

这里的 $\log g$ 即为著名的 Affleck-Ludwig 边界熵（其中 $g$ 为边界简并度）。在含有单个空位（位于 $v$）的链中，整个系统被视作一个整体。因此，空位在此扮演的是“体缺陷”而非“边界割裂”，但由于空位两端自旋的解耦，它对纠缠熵的修正等价于引入了一个特定的缺陷边界熵。若空位使链退化为两条完全独立的断链，则边界熵应呈现出两个独立自由边界的叠加；若空位保持部分连通，则其边界熵将偏离标准边界值，并作为反映缺陷反射强度的精细指标。

1.3 物理模型：含第二邻接相互作用的横场 Ising 链

本研究探讨的模型哈密顿量为：

其中 $\sigma^{x,z}_i$ 为 Pauli 矩阵。当 $J_2 = 0$ 时，模型退化为解析可解的标准横场 Ising 模型；当 $J_2 \neq 0$ 且为铁磁性（$J_2 > 0$）时，系统在临界场 $b = b_c(J_2)$ 处依旧属于 Ising 普适类（中心电荷 $c = 1/2$），而第二邻接项 $J_2$ 作为一个不改变 Bulk 普适类的“不相关算符”（Irrelevant Operator）存在。

空位（Vacancy）的微观微扰构造是本文的核心设计。如图 1 所示，空位通过以下方式实现：

• 彻底移除格点 $v$ 上的自旋，使得所有与格点 $v$ 相连的第一邻接耦合 $J_1$（即 $J_1 \sigma^x_{v-1}\sigma^x_v$ 和 $J_1 \sigma^x_v\sigma^x_{v+1}$）以及第二邻接耦合 $J_2$（即 $J_2 \sigma^x_{v-2}\sigma^x_v$ 和 $J_2 \sigma^x_v\sigma^x_{v+2}$）全部归零。
• 关键之处：保留跨越空位的第二邻接耦合 $J_2 \sigma^x_{v-1}\sigma^x_{v+1}$（在原始晶格中跨越空位，物理距离为 2）。

这种构造保证了即使移除了格点 $v$，空位两侧的子系统仍然通过保留的 $J_2$ 耦合保持物理上的连通。因此，整个体系仍然是一个单一的、连通的系统，包含了 $N-1$ 个自旋。通过调节 $J_2/J_1$，我们可以人为调控这唯一的传输通道的强度。这一精妙设计将 $J_2$ 从一个体不相关算符，转化为调控局域缺陷透射率的边缘相关算符（Marginal Defect Operator），从而生成了一个连续变化的共形缺陷族。

1.4 技术难点一：临界线 $b_c(J_2)$ 的高精度标定

在开放边界条件（OBC）下，利用常规方法（如磁化率峰值或有限尺寸纠缠熵标度拟合中心电荷 $c$）标定临界场 $b_c$ 会引入显著的边界伪差。例如，作者在 $J_2 = 0$ 的自由费米子极限下测试发现，直接利用中部纠缠熵拟合会给出 $c \approx 0.75$ 的错误高估值（真实值为 $c = 1/2$）。

解决方案：本研究独创性地采用体区域自旋-自旋关联函数作为标定判据。在 Ising 临界点上，系统体自旋关联函数必定遵循：

作者在开链的中心区域，选取对称分布的格点对，其间距限制在 $4 \le r \le 2N/5$ 范围内，以最大程度避开边界效应。通过二分法（Bisection）调节外加横场 $b$，直到在当前尺寸 $N$ 下拟合出的指数满足 $|\eta(b_c) - 0.25| < 0.005$，拟合判定系数 $R^2 > 0.97$。测试表明，此方案确定的 $b_c(N)$ 随尺寸 $N$ 的增大平滑收敛，且在 $N=250$ 与 $N=300$ 之间的变动小于 $2.5 \times 10^{-3}$。

1.5 技术难点二：空位系统边界熵的提取与匹配坐标减法方案

提取 Affleck-Ludwig 边界熵 $\log g$ 的标准做法是利用含有缺陷系统的纠缠熵减去无缺陷系统的纠缠熵：$\Delta S(\ell) = S_{\text{vac}}(\ell) - S_{\text{clean}}(\ell)$。然而，由于空位的存在，含空位系统的有效自旋数减少为 $N-1$，这导致在同一格点索引 $\ell$ 处，两者的共形尺寸发生了相对漂移：

• 对于空位左侧（$\ell < v$），真实的物理位置未变，匹配坐标为 $\ell_{\text{red}} = \ell$；
• 对于空位右侧（$\ell > v$），由于移除了一个点，右侧子系统的有效索引整体前移了一位，即匹配坐标应修正为 $\ell_{\text{red}} = \ell - 1$。

若不进行匹配，直接在相同的原始格点索引处做简单相减，将由于体对数项 $\frac{c}{6}\ln[\sin(...)]$ 中的自变量微小偏差，导致在空位中心处产生极其严重的“伪尖峰”（如图 6(a) 所示），使得边界熵的外推完全失效。

解决方案：引入共形匹配坐标减法方案，定义差值纠缠熵为：

其中 $S_{\text{vac}}$ 的对数项计算使用有效长度 $N-1$，而 $S_{\text{clean}}$ 使用 $N$。此消减法完美地消除了所有的体对数贡献和公共的自由边界常数，暴露出纯粹由缺陷引起的纠缠增量。最终，边界熵 $\log g$ 可通过在远离缺陷和物理边界的体区间（如 $0.2(N-1) \le \ell_{\text{red}} \le 0.8(N-1)$，且排除空位邻近的 $0.06N$ 个格点）取中位数来高精度确定。

2. 关键 Benchmark 体系与计算数据分析

2.1 标定数据列表（Table I 深度解析）

研究系统地计算了 5 个尺寸 $N \in \{100, 150, 200, 250, 300\}$ 在 10 个不同第二邻接耦合强度 $J_2/J_1 \in [0.1, 1.0]$ 下的物理参数。以下摘录部分极具代表性的计算数据：

数据揭示的关键物理规律包括：

• 临界场 $b_c$ 随着 $J_2$ 的增大而显著增加（从 1.1616 增至 2.4675），反映了铁磁性第二邻接项对无序相的更强抵抗力；
• 拟合得到的 $c_{\text{fit}}$ 稳定保持在 $0.75 \sim 0.80$ 左右，展现出高度一致的边界高估，从而证明了放弃使用中心电荷标定临界场的绝对正确性；
• 边界熵 $\log g$ 的尺寸波动极小（如在 $J_2=0.5$ 时，从 $N=100$ 的 $-0.0827$ 到 $N=300$ 的 $-0.0771$），说明匹配坐标消减法具有极高的数值稳定性。

2.2 结合能 $\Delta_b(r)$ 的代数衰减与卡西米尔指数的外推

图 2 展示了在 $N = 300$ 尺寸下结合能随距离 $r$ 的对数-对数图。在 $r \in [4, 32]$ 范围内，数据点展现出完美的线性标度行为（拟合优度 $R^2 \ge 0.996$），这直接宣示了量子卡西米尔力的无能隙、幂律衰减本征属性。对于不同 $J_2/J_1$，拟合得到的有效指数 $\alpha(N, J_2)$ 随系统尺寸的倒数 $1/N$ 呈现微弱的非单调漂移。通过对 $1/N \to 0$ 进行线性外推，研究得到了热力学极限下的卡西米尔指数：

这一结果极为关键：指数 $\alpha_\infty$ 极其接近 CFT 预言的理想共形墙相互作用普适指数 $\alpha = 1.0$。虽然随 $J_2$ 的增加，其值从 $1.07$ 缓慢单调漂移至 $1.15$，但由于这一增加斜率仅为 $0.09$，且 $N=300$ 处的有效指数已是所有尺寸中最低的（更接近 1.0），这表明该微小上漂在很大程度上属于未完全出清的有限尺寸修正。可以确信，该相互作用在本质上正是由 Ising 临界涨落介导的量子卡西米尔效应。

2.3 自旋关联透射比 $T_\infty$ 的演化

图 3 详尽展示了透射比 $T(k)$ 的空位穿透特性：

• 在经过极短的近域暂态（$k \lesssim 8$）后，$T(k)$ 迅速切入一个极其平坦的平原相区，该平台可一直稳定延伸至开链边界附近；
• 这一平原高度 $T_\infty(J_2)$ 在热力学极限下展现出强烈的 $J_2$ 依赖性。当 $J_2/J_1 = 0.1$ 时，$T_\infty \simeq 0.112$，表现出强烈的反射性；而当 $J_2/J_1 = 1.0$ 时，$T_\infty \simeq 0.328$，透射能力增强了近三倍。
• 物理判言：既然 $0 < T_\infty < 1$，则该缺陷成功避开了流向断裂（$T_\infty=0$）和流向完全愈合（$T_\infty=1$）这两个平凡的固定点，从而证实 $J_2$ 在格点上的引入成功诱导了一个由不相关体参数调控的、非平凡的部分透射共形缺陷族。

2.4 边界熵 $\log g_\infty$ 的微观物理图景

图 4 给出了最终外推得到的边界熵 $\log g_\infty$ 与 $J_2/J_1$ 的依赖关系：

• 在整个研究区间内，$\log g_\infty$ 几乎维持常数，数值落在极窄的范围内：$\log g_\infty \approx -0.073 \pm 0.003 \text{(stat)} \pm 0.005 \text{(window)}$；
• 这一常数高度远高于 Ising 的固定边界条件值（Fixed BC: $-\ln \sqrt{2} \approx -0.347$），而极其逼近自由边界条件值（Free BC: $0.0$）。

这一物理结论具有高度的普适性：在 BCFT 中，边界熵的大小本质上衡量了缺陷对系统局部相干性的破坏程度。一个“完全断裂”的缺陷会引入两个独立的自由边界，其边界熵的改变量在适当归一化后将对应显著的负值；而本模型中，即使在 $J_2/J_1 = 0.1$ 处透射率仅有 $T_\infty \approx 0.11$，其边界熵依旧极其接近自由开链边界。这深刻说明：只要缺陷两端存在哪怕极微弱的非零量子透射（即不发生绝对割裂），缺陷处的局部纠缠和状态密度便会更倾向于保持一个连通系统的特征，而非演变为两个分立边界。

2.5 协同演化的强证据：$\alpha_\infty$ 与 $T_\infty$ 的高度线性相关

图 5 给出了将外推得到的卡西米尔指数 $\alpha_\infty$ 直接对透射平台值 $T_\infty$ 作图的结果。令人瞩目的是，10 组独立计算的数据点在双物理量空间中画出了一条极其优美的直线：

其 Pearson 相关系数达到了惊人的 $r = 0.986$。这一极高线性相关性并非偶然，它强力印证了：这两个源于不同物理测度（一个是两点卡西米尔力，一个是单点关联穿透率）的物理量，在本质上是由同一个底层的共形缺陷参数（即部分透射固定点的边缘扰动耦合常数）单调控制的。这也为后续建立解析的共形缺陷映射提供了坚实的数值基石。

3. 代码实现细节与高精度复现指南

本节提供一份完整的基于 Julia 语言和开源张量网络库 ITensors.jl 的高精度 DMRG 计算复现方案，包含临界场 $b_c$ 的自动双分标定算法。

3.1 环境部署与软件包导入

首先确保安装了 Julia（推荐 1.9+）及 ITensors 库：

using Pkg
Pkg.add("ITensors")

3.2 核心哈密顿量构造函数

以下 Julia 代码定义了含有单个中央空位的 $J_1-J_2$ 横场 Ising 哈密顿量。请注意代码中对空位格点 $v$ 处耦合项的特殊清空逻辑，以及跨空位 $J_2$ 项的保留：

using ITensors

"""
构建含有单个中央空位的 J1-J2 横场 Ising 链哈密顿量
N: 洁净链的格点数 (空位链的有效自旋数为 N-1)
v: 空位对应的原始格点索引 (通常设为 N ÷ 2)
"""
function build_vacancy_hamiltonian(sites, J1::Float64, J2::Float64, b::Float64, v::Int)
    N = length(sites)
    os = OpSum()
    
    # 1. 单场项 (横场 b)
    for i in 1:N
        if i == v
            continue # 彻底移除空位格点的单场项
        end
        os += -b, "Sz", i
    end
    
    # 2. 第一邻接项 J1
    for i in 1:(N-1)
        if i == v || (i+1) == v
            continue # 任何触碰空位 v 的 J1 耦合全部归零
        end
        os += -J1, "Sx", i, "Sx", i+1
    end
    
    # 3. 第二邻接项 J2
    for i in 1:(N-2)
        # 如果第二邻接项直接触碰空位 v，则归零
        if i == v || (i+2) == v
            continue 
        end
        # 特殊保留：跨越空位 v 的 J2 项 (即自旋 v-1 与 v+1 之间的直接耦合)
        # 当 i == v-1 时，i+2 == v+1，该项得以自然执行并保留
        os += -J2, "Sx", i, "Sx", i+2
    end
    
    return MPO(os, sites)
end

3.3 洁净链哈密顿量与 DMRG 求解器

为了实现“共形匹配坐标减法”，我们同样需要高精度的洁净链基态求解器：

function build_clean_hamiltonian(sites, J1::Float64, J2::Float64, b::Float64)
    N = length(sites)
    os = OpSum()
    for i in 1:N
        os += -b, "Sz", i
    end
    for i in 1:(N-1)
        os += -J1, "Sx", i, "Sx", i+1
    end
    for i in 1:(N-2)
        os += -J2, "Sx", i, "Sx", i+2
    end
    return MPO(os, sites)
end

function run_dmrg_solve(H, sites; maxdim=600, noise_amp=1e-2)
    N = length(sites)
    # 采用 Neel 态作为初始态以规避临界区低纠缠局域极小值
    state = [isodd(n) ? "Up" : "Dn" for n in 1:N]
    psi0 = MPS(sites, state)
    
    # 极其关键的 Sweeps 调节器：在前半段维持强 Mixer 噪声以确保跃出亚稳态
    sweeps = Sweeps(30)
    setmaxdim!(sweeps, 10, 20, 50, 100, 200, 400, maxdim)
    setmindim!(sweeps, 1, 1, 10, 10, 10, 10, 10)
    setcutoff!(sweeps, 1e-11)
    # 在前 20 次扫掠中注入递减的强噪声，后 10 次不加噪声以达到最高精度收敛
    setnoise!(sweeps, noise_amp, noise_amp, noise_amp*0.1, noise_amp*0.01, 0.0)
    
    energy, psi = dmrg(H, psi0, sweeps; silent=true)
    return energy, psi
end

3.4 标定临界场的自动双分算法

下面的脚本展示了如何通过自旋关联指数 $\eta$ 的拟合，自动标定 $b_c$：

using LinearAlgebra

"""
对给定的 (N, J2) 自动寻找临界场 b_c
"""
function calibrate_critical_field(N::Int, J2::Float64; tol=0.005, max_iter=15)
    J1 = 1.0
    b_low = 1.0 + J2  # 物理下界
    b_high = 3.0 + 2*J2 # 物理上界
    b_mid = (b_low + b_high) / 2
    
    sites = siteinds("S=1/2", N)
    
    for iter in 1:max_iter
        H_clean = build_clean_hamiltonian(sites, J1, J2, b_mid)
        _, psi = run_dmrg_solve(H_clean, sites)
        
        # 计算体对称区自旋关联函数
        mid_point = N ÷ 2
        r_vals = Int[]
        corr_vals = Float64[]
        
        # 范围 4 <= r <= 2N/5
        for r in 4:floor(Int, 2*N/5)
            i = mid_point - r ÷ 2
            j = mid_point + r ÷ 2 + (isodd(r) ? 1 : 0)
            
            # 计算 <Sx_i Sx_j>
            # 注：ITensors 快速两点关联计算
            # 在高精度求解中需确保 Gauge 正确
            orthogonalize!(psi, i)
            op_i = op("Sx", sites[i])
            op_j = op("Sx", sites[j])
            
            # 构建夹持算符
            # 此处为简化示意，生产代码建议用更高效的 MPO 夹持或三张量收缩
            corr = inner(psi, OpSum() + ("Sx", i, "Sx", j), psi) # 示意逻辑
            
            push!(r_vals, r)
            push!(corr_vals, abs(corr))
        end
        
        # 线性拟合 log(corr) = -eta * log(r) + const
        X = log.(r_vals)
        Y = log.(corr_vals)
        mean_x = sum(X)/length(X)
        mean_y = sum(Y)/length(Y)
        slope = sum((X .- mean_x) .* (Y .- mean_y)) / sum((X .- mean_x).^2)
        eta = -slope
        
        println("Iter $iter: b_mid = $b_mid, fitted eta = $eta")
        
        if abs(eta - 0.25) < tol
            return b_mid
        elseif eta > 0.25
            # eta 偏大说明自旋关联衰减过快，处于无序相，需要减小横场 b
            b_high = b_mid
        else
            b_low = b_mid
        end
        b_mid = (b_low + b_high) / 2
    end
    return b_mid
end

4. 关键引用文献与方法论局限性批判

4.1 关键参考文献分析

1. Affleck & Ludwig, PRL 67, 161 (1991)：
   • 理论地位：提出了边界共形场论中的关键可观测物理量——边界简并度（或称 Affleck-Ludwig 边界熵 $\log g$），并证明了著名的边界重整化群流 $g$-定理。本工作提取的 $\log g$ 物理概念即源于此。
2. Oshikawa & Affleck, Nucl. Phys. B 495, 533 (1997)：
   • 理论地位：建立了 Ising CFT 中共形缺陷线的边界 CFT 描述，明确指出局部自旋耦合微扰在 Ising 临界点属于 marginal 算符，为本文将 $J_2$ 视为缺陷透射率连续调节源的物理逻辑奠定了坚实的场论基础。
3. Calabrese & Cardy, J. Stat. Mech. P06002 (2004)：
   • 理论地位：给出了一维临界链子系统纠缠熵的标准对数标度律。本文用于分离边界熵的基准 Calabrese-Cardy 对数背景项均基于该文献中的标准形式。
4. Silva et al., Solid State Commun. 134, 313 (2005)：
   • 理论地位：本研究的前期探索性工作。当时受限于计算手段，仅在小尺寸有能隙相区内讨论了空位排斥，并发现相互作用呈指数衰减。本工作则成功将其推进到无能隙临界线这一极具挑战性的卡西米尔标度相区。

4.2 局限性评论与非微扰缺陷物理的未解之谜

尽管本工作在数值模拟和算法创新方面堪称典范，但从前沿量子多体理论和量子化学的严苛视角审视，仍存在以下三点显著的局限性：

1. 缺乏解析的格点到场论映射（Analytic Mapping Gap）

作者提出了一套极其优美的唯象线性规律（如 $\alpha_\infty \propto T_\infty$ ），并断言 $J_2$ 参数化地驱动了整个共形缺陷族。然而，文章并未能给出格点哈密顿量微扰算符与共形缺陷场论算符（例如著名的 Roy-Saleur 缺陷模型中的连续参数）之间的严格数学映射。这使得文中的结论在很大程度上停留于“唯象观测”层面，未能给出透射率 $T_\infty$ 随格点微扰参数 $J_2/J_1$ 演化的首要原理（Ab-initio）解析解析式。

2. 有限系统尺寸（$N=300$）对卡西米尔指数的外推限制

尽管 DMRG 计算在 $N=300$ 处已极度精确，但量子卡西米尔相互作用极其微弱（在 $r=30$ 时，能量量级已跌落至 $10^{-2}$ 甚至更低）。如作者在图 2(b) 中所展示的，有效标度指数 $\alpha(N)$ 随尺寸的演化不仅微弱，而且呈现出令人困惑的非单调行为。这表明在 $N=300$ 尺寸下，由系统真实边界（OBC 的两个端点）产生的长程有限尺寸效应对空位结合能依然存留显著的残留干扰。若要实现彻底的指数去混叠（De-aliasing），在未来的研究中必须借助更强大的张量网络算法（如含杂质的无限尺寸 MPS，即 iMPS 算法）直接在热力学极限下开展计算。

3. 拟合窗口 $r \in [4, 32]$ 选择的主观性

卡西米尔力的拟合区间上限被严格限制在 $r=32$（仅占系统尺寸的 10% 左右）。这一选择虽然成功规避了空位靠近物理边界时发生的有限尺寸崩塌，但极窄的拟合跨度（仅仅半个数量级）在严苛的标度理论看来，极易受到非渐近标度修正（Sub-leading corrections to scaling）的污染。这在一定程度上解释了为何外推得出的 $\alpha_\infty \simeq 1.07$ 始终略高于理论值 $1.00$。未来的计算需要将系统尺寸推进到 $N \sim 1000$ 以上，使拟合窗口能够平稳跨越至少两个十倍频程（Decades）。

5. 临界共形缺陷理论的深度拓展与前沿应用

5.1 边缘扰动算符的重整化群（RG）流：为什么空位不流向“断裂”？

从重整化群的视角来看，在一个一维 Ising 临界链中引入局域缺陷耦合，其对应的 RG 方程极为特殊。如果我们在链中引入一个一维自旋耦合杂质 $J_{\text{impurity}} \sigma^x_v \sigma^x_{v+1}$，由于自旋算符 $\sigma^x$ 的标度维数为 $d_\sigma = 1/8$，该耦合微扰在体临界点上对应的标度维数为：

这表明普通的单键耦合扰动是一个强相关算符（Relevant Perturbation）。在重整化群的流动下，任何微小的杂质减弱都会迅速被放大，最终系统会流向断裂固定点（即 $T_\infty = 0$ 的物理切割状态）。

然而，在本文构造的空位几何中，情况发生了彻底的改变。我们移除了一整点，但保留了跨越该点的第二邻接相互作用。在临界线上，跨越空位的能量流和自旋流耦合，其主导的微扰算符对应于体系统的能量流算符 $\epsilon$（标度维数为 $d_\epsilon = 1$）。由于缺陷处的局域标度维数恰好等于空间维度 1，该扰动算符变为了一个边缘算符（Marginal Perturbation）。正是这种特殊的对称性设计，使得算符的 RG beta 函数在缺陷耦合不为零时完全消失，从而在参数空间中产生了一条连续的临界固定点线（Line of Critical Fixed Points）。这正是本文数据展现出连续、平滑演化图景的深层物理机理。

5.2 在准一代磁性材料中的实验观测可能性

在实际的固体化学和凝聚态材料中，横场 Ising 链体系拥有极佳的实验物理对应物。例如：

• 铌酸钴（$\text{CoNb}_2\text{O}_6$）：被公认为实现一维横场 Ising 模型的经典材料。通过施加临界横向磁场（约 5.5 Tesla），可以将其完美驱动至 Ising 临界点。
• 非磁性杂质掺杂：如果在 $\text{CoNb}_2\text{O}_6$ 的 Cobalt 链中，利用少量的非磁性元素（如锌离子 $\text{Zn}^{2+}$ 或镁离子 $\text{Mg}^{2+}$）替代磁性 $\text{Co}^{2+}$ 离子，即可在晶格上完美制造出静态非磁性空位。

实验上，可以通过以下手段探测本文预言的共形缺陷效应：

1. 非弹性中子散射（INS）：通过测量低能动态结构因子 $S(q, \omega)$ 在空位掺杂后的演化，可以直接提取穿越空位的透射比 $T_\infty$；
2. 极低温比热测量：根据 Affleck-Ludwig 理论，边界熵 $\log g$ 会直接贡献一个极低温下的残余无序熵。通过对比洁净样片与极稀释空位掺杂样片的低温比热系数 $\gamma = C/T$，有望直接在实验上定量测出这一不平凡的边界熵变 $\log g_\infty \approx -0.073$。

5.3 总结与物理展望

本研究通过高精度的 DMRG 数值计算与精妙的边界坐标匹配消除方案，彻底理清了临界 $J_1-J_2$ 横场 Ising 链中静态空位的物理本质。研究确证了：

• 量子临界涨落介导的量子卡西米尔力遵循接近于 $1.0$ 的幂律衰减指数；
• 跨空位保留的第二邻接耦合成功扮演了调控缺陷透射普适类的边缘算符角色；
• 缺陷边界熵展现出极强的“不全断即连通”的拓扑非敏感性。

这一成果为后续在强关联量子化学体系中研究多重杂质的长程协同效应（如临界介导的杂质自组装）、乃至设计基于共形缺陷的一维量子调控器件，奠定了极其扎实的理论与数值基石。