sft-wick：随机场理论中费曼图展开与数值评估的形式化框架及开源软件包

《Computer Physics Communications》：sft-wick: A formalism and package for Feynman-diagram expansion and evaluation in stochastic field theories

　　

当随机场动力学被表述为路径积分形式时，微扰理论变得系统化，但所得到的展开式会迅速产生组合爆炸式的增长。此处所针对的设定包括具有多个分量和多个维度的场、矩阵传播子、张量值耦合，以及由任意n阶累积量指定的非高斯驱动噪声。维克斯配对的数目呈阶乘增长，且分量指标必须经由张量值顶点进行路由。有用的输出并非原始的收缩列表，而是一张图表：每个拓扑对应一个条目，并附有多重性、耦合求和、符号和因果约束。研究人员介绍了sft-wick，一个开源Python包，用于构建这些图表并以数值方式计算其积分。给定作用量和可观测量，该包枚举拓扑不等价的费曼图，推导其代数系数，并根据用户提供的响应函数和累积量函数评估所得图表积分。核心算法先枚举空间拓扑，再路由分量指标，从而避免逐收缩的维克斯展开。响应场约束，包括响应-响应收缩为零、伊藤规定以及因果响应环的消除，均在枚举过程中强制执行。预测结果在两个工作示例中与直接朗之万模拟进行了验证，分别针对高斯和非高斯驱动噪声，并在规定的误差预算内与模拟结果一致。程序摘要程序名称：sft-wick；CPC库程序文件链接：（由编辑添加）；开发者仓库链接：https://github.com/StatFieldTheory/sft-wick；存档软件DOI：10.5281/zenodo.22284796（版本0.4.2，即本文所述版本；最新版本的概念DOI为10.5281/zenodo.20776358）；许可条款：BSD-3-Clause；编程语言：Python 3.10+；补充材料：仓库examples/目录中的工作示例（demo1、demo2、demo3和一个非局域顶点教程笔记本）、生成的测试套件验证目录（docs/verification/catalog.rst）以及完整文档（https://sft-wick.readthedocs.io）。问题性质：以马丁-西格利亚-罗斯（MSR）路径积分框架表述的随机场论中的微扰计算，需要枚举维克斯收缩，其数量随展开阶数呈阶乘增长，使得超越领头阶的手动评估不切实际。给定作用量（场、相互作用顶点、耦合）和可观测量，需要得到每一阶所有拓扑不等价费曼图的结构化列表及其代数系数，并且，在给定用户提供的驱动场统计量时，还需得到所需阶数下可观测量的数值预测。求解方法：sft-wick实现了一个空间拓扑引擎，用于枚举拓扑不等价的图结构而非单个维克斯配对，从而将缩放关系从阶乘降低为空间点数的多项式。每个拓扑的多重性通过附录中推导的组合公式解析计算。MSR约束（ψ–ψ关联子为零、伊藤规定R(x,x)=0、因果R环消除）在枚举过程中原生执行，以实现最大程度的剪枝。等价图通过规范的图形式进行分组，并可选配pynauty后端用于高阶计算。数值评估使用样条缓存的传播子、因果单纯形上的张量积高斯-勒让德求积（用于工作示例），或对于高维积分使用基于Sobol序列的准蒙特卡罗积分，并通过joblib实现令人尴尬的并行图评估。对于由可分离噪声（具有指数时间核）驱动的对角常数漂移，相关传播子由内置闭式表达式求值；其他核则使用高斯-勒让德表，其节点数由收敛检查确定，两条路径在代码切换处相互验证。YAML配置层配合命令行驱动器提供进度报告和试运行成本估算。附加说明（包括限制和不寻常特性）：该包支持多分量场论以及局域和非局域相互作用顶点。支持任意阶的驱动场累积量：高斯二点累积量K⁽²⁾具有内置的便捷接口，而更高阶累积量K^(m)（m?≥?3）被声明为携带任意秩m耦合张量（数值或可调用）的非局域仅ψ顶点，正如第二个工作示例（K⁽³⁾、K⁽⁴⁾和K⁽⁵⁾；examples/demo3）所实践的那样；与外部积分范围相比较窄的核可以以R-收缩形式提供。默认后端下，规范图形式的规模随积分变量数k按k!增长；可选的pynauty后端使用图规范化，建议用于微扰阶数?≥?6。无需外部计算机代数系统。

本文解读围绕sft-wick软件包展开，该软件包面向基于Martin–Siggia–Rose（MSR）路径积分框架的随机场论微扰计算。研究背景方面，MSR形式主义将随机微分方程转化为场论语言，通过引入响应场将随机驱动场的统计性质系统性地传递给物理场，已在临界动力学、随机暴胀、宇宙结构形成、反应扩散系统及神经科学等领域得到广泛应用。然而，现有处理方式存在明显局限：教学性介绍通常局限于单方程或少量标量场，而多分量教科书处理则多置于高斯驱动噪声的Hohenberg–Halperin临界动力学分类之中。实际应用，如引力弱透镜统计场论，需要同时处理多分量随机偏微分方程、张量值耦合、矩阵传播子以及通过高阶累积量进入的非高斯驱动统计。这些要素虽各自见于文献，但从未被整合为统一的可操作方案。此外，即便在形式主义层面可行，微扰展开的实用计算也面临组合爆炸：2n个场算符的完全Wick配对数目为(2n?1)!!，在n=7时已超过10⁵，多分量情形下还需对分量指标经耦合张量路由，进一步加剧复杂度。通用量子场论工具如FeynArts、Cadabra等并不原生支持MSR特有的约束——响应-响应关联子为零、It?处方、因果响应环消除以及C与R双传播子结构，迫使每个阶次都需要大量手动干预。为此，研究人员开发了开源Python包sft-wick，以自动化完整微扰计算管线，并针对多分量、非高斯驱动场论提供系统化解决方案。

研究人员开展了系统的理论整合与算法实现工作。首先给出多分量非高斯MSR形式主义的自洽表述，明确自由传播子（因果响应传播子R与相关传播子C）的显式形式，并阐述Wick定理在MSR框架下的应用及其特有的因果约束。其次，实现了sft-wick软件包，其核心算法是先枚举空间拓扑，再路由分量指标，替代逐收缩的Wick展开，从而将复杂度从阶乘降为多项式。验证工作包括两个工作示例：第一个是双分量、一维空间、含二次非线性并受高斯噪声驱动的Langevin系统，sft-wick计算两点关联函数至四阶，并与Heun积分器直接数值模拟对比；第二个将驱动噪声替换为脉冲泊松过程（散粒噪声），其各阶累积量均有闭式，通过非局域ψ-only顶点处理非高斯累积量，计算交叉关联函数ξ₁₂，并与指数时间差分模拟对照。结论是，在误差预算内，微扰预测与模拟一致：高斯示例中四阶结果落在模拟的统计噪声（约0.6%）范围内；非高斯示例中，领先的FK⁽³⁾道正确重现奇数累积量产生的交叉关联，四阶修正F³K⁽³⁾和F³K⁽⁵⁾均被计算，未计算的O(F⁵)余项经几何外推估计为0.63%；在自由场极限下，sft-wick的数值结果与闭式符合至机器精度。重要意义在于，该软件首次将多分量、非高斯MSR形式主义的全套要素整合进一个自动化工具，使四阶及更高阶微扰计算对多分量场变得可行，为弱透镜统计场论等应用提供了可复用的计算平台。

关于主要关键技术方法：sft-wick采用空间拓扑引擎，在空间拓扑级别枚举拓扑不等价的图结构，并利用组合公式解析计算各拓扑的多重性，避免(2n?1)!!级别的Wick配对枚举。枚举过程中原生实施MSR剪枝规则：跳过ψ–ψ配对、在It?处方下禁止同点R配对、通过记忆化深度优先搜索拒绝形成有向R环的赋值。图等价归类通过规范边表（默认k!排列）或可选的pynauty图规范化实现。数值积分方面，时间积分在因果有序单纯形上采用张量积Gauss–Legendre求积（d≤5）或Sobol序列准蒙特卡罗方法（高维），传播子通过样条缓存，并使用闭式或双重数值积分评估C传播子。非局域累积量顶点可提供原始核或R-收缩形式，后者将ψ腿积分纳入用户可调用函数以降低积分维度。

研究结果部分按原文小标题归纳：

**Worked example I: Gaussian driving field**——研究人员考虑双分量Langevin方程，线性算符A_ab=?δ_ab，高斯噪声具有因子化时空协方差（指数核，参数λ=0.05, σ_t=0.3, σ_x=1.0）。直接数值模拟采用Heun积分器（Δt=0.02）和自回归精确生成时间相关噪声，共10⁵个独立实现。微扰展开至四阶，sft-wick在二阶识别出6个拓扑不等价图，四阶为64个；代入具体耦合后许多图系数为零。Gauss–Legendre求积在因果时间单纯形上给出确定性结果。结果显示，对于ξ₁₁和ξ₂₂，四阶预测在全部时间范围内落在模拟的1σ误差带内；空间依赖在t=1,15,30三个时刻均与模拟定量一致；逐阶贡献呈几何抑制，每增加二阶约抑制因子λ≈0.05。残差主要来源为模拟噪声（约0.6%），截断误差估计约0.3%，求积误差小于0.1%（最大时间处约2%的相应道），Heun时间步偏差约0.04%可忽略。

**Worked example II: non-Gaussian driving field and a non-local vertex**——将高斯噪声替换为过滤泊松过程，各阶累积量由Campbell定理给出闭式，均因子化通过单一源点。参数取ν=2, h=1, σ_t=0.5, σ_x=1，并定义n≡νσ_tσ_x，n=1时偏度为0.63。非高斯部分在MSR作用量中产生仅含ψ的非局域顶点。通过腿数计数规则K^(m)在阶数m?1首次进入。由于对称性，ξ₁₂在Gaussian示例中恒为零，因此直接度量非高斯贡献。结果显示：在F=0自由场极限下，sft-wick复现闭式至相对误差1.4×10^?16（m=3）和6.6×10^?16（m=4）；模拟验证了偏度随1/n变化的规律（因子4与预期4.04一致）。在s=0.2时，ξ₁₂的三个领先道被计算：二阶FK⁽³⁾、四阶F³K⁽³⁾（30个图）和F³K⁽⁵⁾（6个图），t=3时两个修正分别为领先道的7.9%和0.094%，未计算的O(F⁵)余项经几何外推估计为0.63%。模拟采用指数时间差分法（6个种子各4×10⁵实现），ξ₁₂(t)各时刻均在0.90σ内一致，ξ₁₂(r)在1.35σ内，对角ξ₁₁在0.69σ内，7.3×10⁶条轨迹无一发散。

**Algorithm**——算法分为五个阶段：作用量展开（多顶点多重集组合）、Wick收缩（空间拓扑级枚举并应用三条剪枝规则，多重性公式μ=∏_vm_v!/(2^N_self·∏_ek_e!)）、图收集（规范形式去重，并联合查找简化传播子指标）、耦合评估（代入耦合张量、求和耦合仅指标，应用组合因子和MSR相位）、数值积分（方向等价类分解、时间排序DAG、Gauss–Legendre或拟蒙特卡洛、并行）。复杂度对比显示，直接Wick枚举在n=1时需15个配对，sft-wick仅需2.96×10^?4秒并返回6个图；n=3时配对数为1.0×10⁴，sft-wick仅返回80个图并耗时1.22×10^?1秒。

总结讨论部分指出，sft-wick作为开源包自动化了MSR框架下从作用量和可观测量到符号图表达式及数值预测的完整流程。核心创新是空间拓扑引擎结合MSR特定剪枝和图同构检测，避免了全量Wick配对枚举，使多分量场的四阶及更高阶计算可行。验证显示高斯驱动下逐阶几何抑制（λ=0.05时每偶数阶约λ倍），四阶与模拟在统计噪声内一致；非高斯情形下，仅由奇数累积量产生的交叉关联被正确重现，且四阶修正被完整计算。限制包括规范图形式随积分变量数k按k!缩放（推荐高阶使用pynauty后端），以及高阶累积量核若在相对时间上较窄则应以R-收缩形式提供。未来方向包括自动微分集成、GPU高维积分、平移不变系统的动量空间表示及重求和技术。sft-wick的开发由弱引力透镜统计场论应用所驱动，该应用详见姊妹论文[16]，研究结论是：该软件包为多分量、非高斯MSR微扰计算提供了通用且验证过的工具，可将预测推进至手动或通用QFT工具无法企及的阶数，并在规定误差内与直接数值模拟相符。