该模型在物理学界引起广泛关注的原因在于：

1. 可解性与全息性：在大 $N$ 极限下，其由一类称为梅隆图（Melon Diagrams）的费曼图主导，可以通过 Schwinger-Dyson 方程解析求解。在低能区，该模型展现出涌现的共形对称性（$SL(2,\mathbb{R})$ 重新参数化对称性）及其显式破缺，其低能有效作用量为 Schwarzian 作用量，与二维近 AdS2 背景下的 Jackiw-Teitelboim (JT) 相对论引力完全对偶。
2. 最大混沌性：其低温条件下的 Lyapunov 指数达到马尔达西纳-沈克-斯坦福（Maldacena-Shenker-Stanford, MSS）绑定的上限 $\lambda_L = 2\pi k_B T / \hbar$，表现出强烈的量子混沌动力学。

1.4 理论基础之三：Kourkoulou-Maldacena (KM) 纯态与其引力对偶

为了构建具有明确几何解释的纯态，Kourkoulou 和 Maldacena 引入了一类特殊的边界态。首先，将 $N$ 个 Majorana 费米子配对，定义互不相干的自旋算符（对角算符）：

由于 $\hat{S}_k$ 彼此对易，可定义它们的共同本征态 $|B_{\vec{s}}\rangle$：

这些状态 $|B_{\vec{s}}\rangle$ 是完全去相干的费米子高斯态。接着，利用 SYK 哈密顿量对其进行虚时间（欧氏时间）演化，得到规范化的纯 KM 态：

全息引力解释：

在体区引力对偶中，$|B_{\vec{s}}(\beta)\rangle$ 对应于一个包含近 AdS2 视界的单边界永恒黑洞，但在其分叉视界（Bifurcation Horizon）后方，存在一个由边界状态 $\vec{s}$ 决定的**“世界末日粒子”（End-of-World Brane/Particle）**。参数 $\beta$ 对应于黑洞的倒数温度。当 $\beta \to 0$ 时，黑洞演化为极高温度，引力背景消失，退化为无相互作用的高斯态；而当 $\beta \to \infty$ 时，系统深入红外（IR）共形区，对应于接近极值黑洞的复杂量子态。

1.5 核心公式推导与解析方法细节

1.5.1 大 $N$ 极限下 KM 态 FAF 的显式简化

在大 $N$ 极限下，由于系统涌现出 $O(N)$ 连续对称性，该对称性包含一个离散子群——翻转群（Flip Group），其生成元可以独立翻转自旋序列 $\vec{s}$ 的分量。这一选择定则导致除了对角（配对内）关联外，所有非同名费米子之间的两点关联矩阵元在热力学极限下完全消失：

因此，协方差矩阵 $M_{mn}$ 简化为仅由两点“非对角”（Off-diagonal）关联构成。定义欧氏格林函数 $G_E(\tau) \equiv \frac{1}{N} \sum_{i=1}^N \langle T \chi_i(\tau) \chi_i(0)\rangle_{\beta}$，在大 $N$ 下利用对称性，FAF (Eq. 1) 的计算转化为：

该公式是本文最核心的解析成果，它将极其复杂的量子魔性直接与 SYK 模型的两点格林函数 $G_E(\beta/2)$ 关联起来！

1.5.2 实时演化的解析延续

若初始状态为高斯态 $|B_{\vec{s}}\rangle$，并在实时间 $t$ 下进行演化：

对应的 FAF 时间动力学可通过解析延续 $\tau = \beta/2 + i t$ 得到：

其中 $G^W(t) \equiv \frac{1}{N} \sum_{i=1}^N \text{tr} [ e^{-\beta H/2} \chi_i(t) e^{-\beta H/2} \chi_i(0) ] / Z_{\beta}$ 为正则化的 Wightman 关联函数。对于无损高斯初态（$\beta=0$ 极限），$G^W(t)$ 即为双时格林函数。

1.6 技术难点：指数级希尔伯特空间与无序平均的阵痛

尽管大 $N$ 解析公式非常优雅，但在有限尺寸 $N$ 下验证这些解析结论、提取次领先修正（Subleading Corrections）并研究其在低能区（基态附近）的行为，面临着极高的数值计算挑战：

1. 维数灾难：$N$ 个 Majorana 费米子的希尔伯特空间维度为 $\mathcal{D} = 2^{N/2}$。当 $N=54$ 时，$\mathcal{D} = 2^{27} \approx 1.34 \times 10^8$。如此庞大的复数向量，仅存储一个态矢量就需要：

   $$1.34 \times 10^8 \times 16 \text{ Bytes} \approx 2.14 \text{ GB}$$

   进行完全精确对角化（ED）需要 $O(\mathcal{D}^3)$ 的算力，这在 $N > 30$ 时便绝对不可行。

2. 无序平均（Disorder Average）：SYK 模型是随机无序系统。为了消除样本间的涨落、提取干净的物理标度，特别是在大时间尺度（如 Heisenberg 时间）或低能基态附近，必须对数万乃至数百万个独立无序样本进行集成平均。例如，在计算 $N=18$ 的长时间演化时，论文使用了高达 $10^7$ 个无序实现！这一计算量对任何传统的 CPU 集群都是毁灭性的。

为了克服这一技术难点，作者团队开发了结合 GPU 硬件加速 与 Krylov 子空间（Lanczos 算法）投影技术 的高性能计算管线，实现了对最大 $N=54$ 规模系统的时间演化和状态制备，这是目前该领域的计算极限之一。

2. 关键 Benchmark 体系、计算数据与性能数据

在本节中，我们将逐一剖析论文中给出的关键 Benchmark 测试及其核心物理数据。

2.1 KM 纯态的欧氏演化与非平坦性调控 (Fig 1)

体系描述：

研究了 $q=4$ 的 SYK 模型中 KM 纯态 $|B_{\vec{s}}(\beta)\rangle$ 的 FAF。通过调控虚时间参数 $\beta$（即等效逆温度），观察量子魔性的变化。

数据分析与物理现象 (参考图 1)：

• 图 1(a)（解析计算） 展示了利用 Schwinger-Dyson (SD) 方程数值求解得到的 $2\mathcal{F}/N$ 随 $\beta$ 的变化曲线，并将其与共形极限下的解析逼近进行对比。当 $\beta \to 0$ 时，FAF $2\mathcal{F}/N \to 0$。随着 $\beta$ 增加，系统受到强关联相互作用的塑造，魔性单调上升。在低温极限下（$\beta J \gg 1$），FAF 逼近其理论最大值：

  $$\lim_{\beta \to \infty} \frac{2\mathcal{F}}{N} = 1$$

  此时，低温渐近线可通过共形微扰展开式（Eq. 7）极其精准地拟合：

  $$\mathcal{F}(\beta J) \approx N \left[ \frac{1}{2} - \frac{2\pi}{(\beta J)^2} + \frac{13.311}{(\beta J)^3} + \frac{9.291}{(\beta J)^4} + \dots \right]$$

• 图 1(b)（有限尺寸数值对比） 展示了从 $N=18$ 到 $N=54$ 的有限尺寸精确数值结果（彩色实线带符号）与大 $N$ 解析 SD 方程解（黑色虚线）的对比。可以看出：

  1. 数值结果与大 $N$ 理论符合得极好，随着 $N$ 增大，数值曲线迅速收敛至大 $N$ 解析曲线上。
  2. 即使在低温大 $\beta$ 区（$\beta \approx 50$），利用双对数插图（Inset log-log plot）可以清晰看到，$1 - 2\mathcal{F}/N$ 呈现出标准的幂律衰减，数值点与 SD 解析解完美重合，直接证实了大 $N$ 极限下选择定则（Eq. 5）的正确性。

2.2 实时动力学演化与饱和行为 (Fig 2)

体系描述：

初始状态为无关联高斯态 $|B_{\vec{s}}\rangle$（即 $\beta=0$），在实时间 $t$ 驱动下，研究 $\mathcal{F}(t)$ 的产生与饱和过程。

数据分析与物理现象 (参考图 2)：

• 短时间行为（UV 极限，Fig 2(b)）：在极短时间内（$tJ < 3.56$），高阶 Padé 1000 逼近与短期 Taylor 展开（Eq. 9）能够完美复现格林函数的极点结构。系统 FAF 随时间呈二次方起步：$\mathcal{F}(t)/N \approx (Jt)^2 / 4 + \dots$
• 长时间饱和与 Ruelle-Pollicot 共振（Fig 2(a) & 2(d)）：
  • FAF 趋向饱和的时间尺度约为 $t \sim 1/J$，令人惊讶的是，这个饱和时间完全独立于系统尺寸 $N$。这与普通格点自旋链中魔性产生速率随 $N$ 增长的行为截然不同，体现了 SYK 模型的强非局域全对全关联特性。
  • 图 2(d) 展示了 $N=50$ 时，FAF 逼近其饱和值的差值 $|\mathcal{F}(t) - \mathcal{F}(\infty)|$ 随时间的指数衰减行为。数值拟合得到的指数衰减率约为 $\approx 2.11$。根据量子混沌理论，这一衰减率与正则化 Wightman 函数的衰减率直接相关，对应于领先的 Ruelle-Pollicot 共振 极点。解析预测的衰减率约为 $\approx 2.52$，数值与理论的接近证明了量子魔性的弛豫由混沌动力学的主导共振极点控制。
• 海森堡时间尺度的涨落 (Fig 2(c))：在极长时间（$t \sim e^{N/2}$，即海森堡时间，能级平均间距的倒数）下，对于 $N=18$ 体系进行 $10^7$ 次无序平均后，可以清晰地观察到类似谱甲板（Ramp-Plateau）的结构。这一现象与随机矩阵理论（RMT）的预测高度一致，表明量子魔性在极长时间下能够作为探测多体量子混沌谱统计特征的精细探针。

2.3 能级本征态的有限尺寸修正与稀疏化效应 (Fig 3)

为了评估全息引力描述在有限 $N$ 下的鲁棒性，必须深入研究 FAF 的有限尺寸修正（Finite-$N$ Corrections）。本节给出了本征态（Eigenstates）的系统标度行为。对于能量本征态 $|E\rangle$，其领先的 FAF 标度为 $\mathcal{F} \sim N/2$。

数据分析 (参考图 3)：

• 能谱中心（Center of Spectrum，高能区混沌态）：

  • 对于密集的经典 SYK 模型（Dense SYK），在能谱中心附近的本征态，其有限尺寸修正 $|2\mathcal{F}_{cen}/N - 1|$ 随 $N$ 呈指数衰减：

    $$\left|\frac{2\mathcal{F}_{cen}}{N} - 1\right| \propto N e^{-b N}, \qquad b \approx 0.33$$

    这一行为与 Haar 随机态极其相似，证实了光谱中心的本征态在量子魔性意义下已经高度接近极限混沌的 Haar 态。

  • 当对 SYK 进行稀疏化（Sparsification）处理（即随机将大部分耦合系数 $J_{ijkl}$ 设为零，保留比例 $p$ 满足 $k = p \binom{N}{4} \sim O(1)$），光谱中心的本征态修正依然保持指数级衰减（如 $N e^{-0.34 N}$ 等）。这表明稀疏化并没有从根本上破坏高能区的强混沌性质。

• 基态附近（Near Ground State，低能区全息对应区）：

  • 关键结论：在接近基态（例如第二激发态）的区域，有限尺寸修正行为发生了戏剧性的改变！

  • 在密集 SYK 中，虽然修正依然是指数级衰减，但衰减指数显著变小（$\approx 0.18$），导致在相同尺寸下，其偏差比能谱中心大出数个数量级。

  • 最重要的发现在于稀疏 SYK：在基态附近，次领先修正完全退化为幂律衰减（Power-law Decay）：

    $$\left|\frac{2\mathcal{F}_{cen}}{N} - 1\right| \approx \frac{1}{N^3}$$

    这一从指数到幂律的跃变具有深刻的全息启示：在低能区，稀疏化虽然保留了系统的某些大尺寸热力学性质，但它极大地增强了有限尺寸下的非高斯涨落，使得低能全息状态在有限 $N$ 下显著偏离了密集 SYK 对应的光滑引力双几何。这为在实验室芯片中重构全息黑洞提出了重要的尺寸与稀疏度限制限制条件。

2.4 有限尺寸 KM 态的 FAF 标度性质 (Fig 4)

体系描述：

研究了不同温度参数 $\beta$ 下，纯 KM 态的有限尺寸偏差 $|1 - \mathcal{F}_{pure}/(b_{\beta} N)|$ 随 $N$ 的变化关系，其中 $b_{\beta}$ 是大 $N$ 极限下的渐近线性斜率。

数据分析 (参考图 4)：

• 在大 $N$ 下，$\mathcal{F}_{pure} \approx b_{\beta} N$。其中斜率 $b_{\beta}$ 随着 $\beta$ 增加，从 $0$（当 $\beta \to 0$）连续单调增长到 $1/2$（当 $\beta \to \infty$），提供了一个连续调控量子魔性的物理旋钮。
• 对于领先修正的有限 $N$ 行为，图 4 显示其修正大小随着 $\beta$ 的增加而急剧上升（从 $\beta=0.1$ 到 $\beta=50$，修正增大超过两个数量级）。
• 尽管对于中等大小的 $N$，该修正呈现指数衰减趋势，但在最大尺寸 $N$ 附近，曲线出现了明显的偏离。作者指出，这可能是由于在进行大 $N$ 线性拟合以提取 $b_{\beta}$ 时，引入了微小的系统残差误差，这暗示了在超低温极限下，多体状态的魔性具有更加复杂的量子相干结构。

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

为了使量子化学与凝聚态计算领域的科研人员能够快速复现本文的数值计算结果，本节提供一份基于 Julia（推荐，在处理算子代数和稀疏矩阵演化上性能极佳）与 Python 的核心实现方案及算法流程指南。

3.1 Majorana 费米子代数的数值构建

为了在计算机中表示满足 $\{\chi_i, \chi_j\} = \delta_{ij}$ 的 Majorana 费米子，我们必须通过 Jordan-Wigner 变换将它们映射到标准的 Pauli 矩阵上。对于 $N$ 个 Majorana 费米子，需要 $M = N/2$ 个量子比特。定义如下：

对于 $k = 1, \dots, M$：

Julia 核心复现代码：构建 Majorana 算符

using LinearAlgebra
using SparseArrays

# 定义基本 Pauli 矩阵
const Ix = sparse([1 0; 0 1] .+ 0im)
const Xx = sparse([0 1; 1 0] .+ 0im)
const Yy = sparse([0 -1im; 1im 0])
const Zz = sparse([1 0; 0 -1] .+ 0im)

"""
构建第 i 个 Majorana 费米子算符，共 N 个，需要 M = N/2 个 qubits
"""
function construct_majorana(i::Int, N::Int)
    @assert N % 2 == 0 "N 必须为偶数"
    M = div(N, 2)
    k = div(i + 1, 2)  # 对应的 qubit 指标
    is_odd = (i % 2 != 0)
    
    op_list = Vector{SparseMatrixCSC{ComplexF64, Int}}(undef, M)
    for p in 1:M
        if p < k
            op_list[p] = Zz
        elseif p == k
            op_list[p] = is_odd ? Xx : Yy
        else
            op_list[p] = Ix
        end
    end
    
    # Kronecker 积合成大算符
    op = op_list[1]
    for p in 2:M
        op = kron(op, op_list[p])
    end
    return op / √2
end

3.2 SYK 哈密顿量与 KM 态制备

利用上述 Majorana 算符，我们可以极易地构造出 $q=4$ 的 SYK 哈密顿量：

"""
构造无序平均的 SYK 哈密顿量
"""
function construct_syk_hamiltonian(N::Int, J::Float64=1.0)
    H = spzeros(ComplexF64, 2^(div(N, 2)), 2^(div(N, 2)))
    majoranas = [construct_majorana(i, N) for i in 1:N]
    
    # 生成高斯随机耦合
    variance = 6.0 * J^2 / N^3
    std_dev = √variance
    
    for i in 1:N, j in (i+1):N, k in (j+1):N, l in (k+1):N
        J_ijkl = randn() * std_dev
        term = majoranas[i] * majoranas[j] * majoranas[k] * majoranas[l]
        H += J_ijkl * term
    end
    return H
end

对于 Krylov 子空间演化（例如计算虚时间演化 $|B_{\vec{s}}(\beta)\rangle \propto e^{-\beta H/2} |B_{\vec{s}}\rangle$），我们可以使用 KrylovKit.jl 中的指数传播子算法 exponentiate。这避免了对大矩阵进行显式对角化。

using KrylovKit

"""
制备 KM 纯态
"""
function prepare_km_state(H, s::Vector{Int}, beta::Float64, N::Int)
    M = div(N, 2)
    # 1. 构造初始配对高斯态 |B_s>
    # 满足 S_k |B_s> = s_k |B_s>
    # 我们可以通过对角化共同本征算符，或者直接在计算基底中选择对应的自旋状态。
    # 最简单的方法是直接从全 1 基态出发，通过投影算子 P_k = (I + s_k S_k)/2 作用在任意非零态上并归一化。
    dim = 2^M
    state = randn(ComplexF64, dim)
    state /= norm(state)
    
    # 施加投影算符
    for k in 1:M
        Sk = -2im * construct_majorana(2k-1, N) * construct_majorana(2k, N)
        proj = (I + s[k] * Sk) / 2
        state = proj * state
        state /= norm(state)
    end
    
    # 2. 虚时间演化 e^{-beta * H / 2} |B_s>
    # 使用 Krylov 指数工具进行精确虚时传播
    km_state, info = exponentiate(H, -beta / 2, state; ishermitian=true, tol=1e-10)
    km_state /= norm(km_state)
    
    return km_state
end

3.3 FAF 的高效数值计算管线

在获得状态矢量 $|\Psi\rangle$（例如 km_state）后，我们需要计算两点协方差矩阵 $M_{mn}$，并据此计算 FAF。注意到对于大规模系统，$M$ 的维度仅为 $N \times N$（远小于希尔伯特空间维度），因此在 $M$ 空间内的迹和矩阵乘法计算量极小，真正的瓶颈在于计算 $N(N-1)/2$ 个期望值 $\langle\Psi| \chi_m \chi_n |\Psi\rangle$。

"""
计算任意纯态的 FAF (k=1)
"""
function calculate_faf(state::Vector{ComplexF64}, N::Int)
    M_mat = zeros(ComplexF64, N, N)
    majoranas = [construct_majorana(i, N) for i in 1:N]
    
    # 仅需计算上三角，利用反对称性填充
    for m in 1:N
        for n in (m+1):N
            # 计算算符 [chi_m, chi_n] = 2 * chi_m * chi_n (m < n)
            # 注意到 chi_m 乘以 chi_n 在 m != n 时是反对易的
            op = 2.0 * majoranas[m] * majoranas[n]
            val = dot(state, op * state)  # <Psi| [chi_m, chi_n] |Psi>
            M_mat[m, n] = -1im * val
            M_mat[n, m] = 1im * val
        end
    end
    
    # M^T * M = - M^2 (由于 M 是反对称实矩阵，当期望值为实数时)
    # 这里直接按照定义计算 tr(M^T * M)
    # M_mat 的伴随矩阵为 M_mat^H. 因为 M 是纯虚反对称，M^T = -M = M^H
    tr_MtM = sum(abs2, M_mat)
    
    faf = N/2 - tr_MtM
    return faf
end

GPU 加速说明 (PyTorch / CUDA.jl)：

当 $N > 30$ 时，上述稀疏矩阵与密集向量的乘法 op * state 应完全转移至 GPU。在 Julia 中，只需导入 CUDA.jl 并将状态矢量及哈密顿量包装为 CuArray 即可：

using CUDA
# 将哈密顿量与初始状态转移至 GPU
H_gpu = cu(H)
state_gpu = cu(state)
# 使用 CUDA 优化的 KrylovKit 进行演化，速度可提升 10-100 倍！

对于开源库，可以直接调用经典的 SYK 数值求解包，如 GitHub 上的 syk-clean 或相关的全息引力研究 repo，它们提供了经过深度优化的 Majorana 费米子 GPU 算子库。

4. 关键引用文献与局限性评论

4.1 关键引用文献及原文中的作用

1. [40] I. Kourkoulou and J. Maldacena, “Pure states in the syk model and nearly-$AdS_2$ gravity” (2017), arXiv:1707.02325.
   • 作用：这是本项工作的奠基性文献，首次定义了 KM 纯态，并指明了其对偶于包含世界末日分岔时空粒子的 AdS2 黑洞。本文所有关于 $|B_{\vec{s}}(\beta)\rangle$ 的解析与数值探讨均建立在其涌现翻转对称性（Flip Symmetry）的基础之上。
2. [21] P. Sierant, P. Stornati, and Y. Turkeshi, “Fermionic magic resources of quantum many-body systems”, PRX Quantum (2026).
   • 作用：本文的核心工具——费米子抗平坦性（FAF）由此文献首次提出。它提供了一个极其简洁高效的代数框架，避开了直接计算指数级复杂的稳定子里尼熵（SRE），使大系统尺寸的量子魔性定量分析成为可能。
3. [31] J. Maldacena and D. Stanford, “Remarks on the sachdev-ye-kitaev model”, Phys. Rev. D (2016).
   • 作用：提供了大 $N$ 极限下 SYK 模型的 Schwinger-Dyson 解析求解细节、低温共形微扰展开式（Eq. 7）以及大 $q$ 极限下的格林函数行为。是本文进行解析逼近和 SD 方程数值验证的直接理论支撑。
4. [34] S. Bera and M. Schirò, “Non-stabilizerness of sachdev-ye-kitaev model”, SciPost Phys. (2025).
   • 作用：作为对比文献，利用传统的 SRE 研究了 SYK 的魔性。本文的研究进一步将其拓展到具引力对偶的 KM 纯态和实时演化共振领域，并揭示了稀疏化的关键效应。

4.2 对本工作局限性与开放性问题的学术评论

作为一篇发表在多体物理与高能引力前沿的杰出工作，该论文成功搭建了量子魔性与全息黑洞时空之间的桥梁。然而，从严格的量子化学与多体物理计算视角来看，该研究仍存在以下局限性与尚未解决的开放性问题：

1. FAF 作为量子魔性度量的操作性定义问题：
   • 局限性：尽管 FAF 在计算上极为方便，且在一系列典型状态下与稳定子 Rényi 熵（SRE）表现出强烈的单调相关性，但它在严格的量子信息论中并非一个完美的魔性单调量（Magic Monotone）。对于某些特定的非高斯操作，FAF 可能会发生非物理的非单调跃变。因此，基于 FAF 得到的全息物理结论，其精细结构是否在真正的 SRE 下完全成立，仍需更严谨的数学确证。
2. 非共形区与过渡区的全息引力对偶缺乏精确对应：
   • 局限性：论文指出，可以通过调节 $\beta$ 使得 FAF 在 $0$ 到 $N/2$ 之间连续调。然而，极低温共形区（$\beta J \gg 1$）才是 JT 引力对偶成立的唯一合法区间。当 $\beta \to 0$（高温区，魔性趋于0）或处于中等温度过渡区时，体区引力描述中会出现强烈的量子引力涨落，此时“世界末日粒子”和光滑度规时空的几何图像将不再有效。因此，在这部分区间“调谐量子魔性”在引力侧缺乏清晰、干净的物理对应物。
3. 稀疏化本征态基态附近幂律行为的引力解释空白：
   • 局限性：论文数值发现了稀疏 SYK 模型中，接近基态的能级本征态的有限尺寸修正从指数衰减退化为幂律衰减 $\propto 1/N^3$。这是一个非常迷人的物理发现，但作者并未能在体区引力侧给出一个明确的物理图像。究竟何种 Bulk 效应（是虫洞涨落、非局域关联的丧失，还是时空本征涨落）导致了魔性次领先修正的幂律退化？这仍是一个亟待解决的理论悬案。

5. 补充探讨：大 $q$ 极限与黑洞准正规模式 (QNM)

为了使该研究的物理解析更为完备，论文在 End Matter 部分将 $q=4$ 的 SYK 推广到了更为普适的 大 $q$ 极限（Large-$q$ Limit）。这一部分极具物理美感，因为它允许我们在不需要借助任何 SD 方程数值求解的情况下，得到完全闭合形式（Closed-form）的解析解。

5.1 大 $q$ 极限下的 FAF 解析求解

对于一般的 $q$ 体相互作用 SYK 哈密顿量（Eq. 10），定义重新标度后的耦合常数 $\mathcal{J} = \sqrt{q} J / 2^{(q-1)/2}$。在大 $q$ 极限下，欧氏格林函数具有如下精确解析解：

其中参数 $v$ 由超越方程 $\beta \mathcal{J} = \frac{\pi v}{\cos(\pi v / 2)}$ 唯一确定。将该式代入 FAF 表达式中，我们立刻得到了大 $q$ 极限下 KM 纯态量子魔性的完全解析显式解：

高低温极限下的解析行为 (参考图 5(a))：

• 低温极限（$\beta \mathcal{J} \gg 1$）：通过对 Eq. 17 进行系统性级数展开，可以得到完全显式的微扰项（Eq. 18）：

  $$\frac{1}{N} \mathcal{F}(\beta) \approx \frac{1}{2} - \frac{1}{2} \left[ 1 + \frac{2}{q} \left( \log \frac{\pi}{\beta \mathcal{J}} - \frac{2}{\beta \mathcal{J}} + \frac{2}{(\beta \mathcal{J})^2} - \frac{8+\pi^2}{3(\beta \mathcal{J})^3} \right) \right]^4$$

• 高温极限（$\beta \mathcal{J} \ll 1$）：展开式表现为标准的渐近级数形式（Eq. 19）：

  $$\frac{1}{N} \mathcal{F}(\beta) \approx \frac{(\beta \mathcal{J})^2}{2q} - \frac{5q+9}{48q^2} (\beta \mathcal{J})^4 + \frac{94q^2 + 225q + 90}{2880q^3} (\beta \mathcal{J})^6 + \dots$$

  从图 5(a) 可以看到，这些大 $q$ 解析线（灰色虚线）与通过求解 Schwinger-Dyson 方程得到的不同 $q$ 值的精确数值曲线完美吻合。随着 $q$ 增大（从 $q=2$ 到 $q=20$），数值曲线平滑且迅速地逼近大 $q$ 解析预测，展现了这一解析方法的强大威力。

5.2 准正规模式与量子魔性的衰减率

对于实时间演化，利用大 $q$ 极限下的 Wightman 函数（Eq. 20）：

将其直接代入 FAF 的演化公式（Eq. 21），得到：

对其进行长时间渐近展开（$t \to \infty$），我们发现量子魔性逼近其饱和值的规律为：

物理启示：

这一指数级衰减对应的衰减率显然为 $\Gamma = 8\mathcal{J}/q$。在广义相对论和黑洞物理中，黑洞受到外界微扰后的无辐射衰减响应由**准正规模式（Quasi-normal Modes, QNM）**控制。论文指出，在大 $q$ 极限下，其 FAF 的极点结构（准正规频率）为：

其最低阶（主导极点，即 $n=0$）对应的虚部正是：

这是一个极其令人振奋的发现！它直接在数学上证明了：在边界上，强关联混沌系统其量子魔性的晚期弛豫衰减率，恰好完美等于其体区全息对偶黑洞的领先准正规模式衰减率。 量子魔性的流失与时空扰动的消散，在全息动力学中达成了完美的代数统一。这为未来基于量子模拟设备观测黑洞准正规模式提供了除传统的两点格林函数外，更为敏感、蕴含更多高阶信息的全新非线性量子相干探测手段。