非互易相互作用(non-reciprocal interactions)的哈密顿量(Hamiltonian)描述

《Nature Physics》：Hamiltonian description of non-reciprocal interactions

　　

摘　要：在许多体系中——包括沉降颗粒与鸟群——相互作用不源自势函数且通常为非互易的(non-reciprocal)，即不满足作用力与反作用力原理（牛顿第三定律）。因此无法定义常规的能量函数，也无法使用依赖于它的分析与数值工具。本文通过构造含辅助自由度(auxiliary degrees of freedom)的哈密顿量(Hamiltonian)，并在约束条件下使其生成原始的非互易动力学，以解决此局限。研究人员表明，基于该约束哈密顿量的蒙特卡洛(Monte Carlo)模拟可重现原始朗之万(Langevin)动力学的稳态与非稳态，并以具视锥(vision-cone)相互作用的耗散XY自旋为例进行了说明。该构造固有的辛结构(symplectic structure)使研究人员能应用哈密顿量工程(Hamiltonian engineering)中的既有思想，研究人员通过调节周期性（Floquet）驱动的振幅，将自旋相互作用在方格与点链式晶格几何之间调谐。总体而言，本构造为非互易体系推广统计力学与哈密顿动力学铺平了道路。

论文解读：非互易相互作用的哈密顿量描述

一、研究背景与意义

在经典与统计力学中，保守势场衍生出的相互作用满足牛顿第三定律，即互易(reciprocal)性，可写出哈密顿量(Hamiltonian)或能量函数，进而使用正则方程、正则变换(canonical transformation)、蒙特卡洛(Monte Carlo, MC)方法及平衡态统计力学进行分析。然而大量介观与宏观体系——如沉降胶体颗粒、活性胶体(droplets)、化学趋向性粒子及动物群体（鸟群、兽群）——其相互作用由非平衡介质（流体、化学场、温度场或视觉感知）介导，表现为非互易(non-reciprocal)相互作用，即粒子i对j的作用力F_j→i≠?F_i→j，无法从单一相互势函数导出，也不具备辛结构(symplectic structure)。这类体系本质上是非平衡态的，缺乏定义良好的能量泛函，导致：(1)无法直接应用基于Boltzmann因子的Metropolis或Glauber蒙特卡洛算法；(2)无法使用依赖于泊松括号(Poisson bracket)的哈密顿分析工具（如正则变换、Floquet理论等）。尽管近年对非互易相变、非互易Ising/XY模型有所研究，仍缺乏普适的哈密顿框架将平衡态分析工具迁移至非互易体系。本文（发表于Nature Physics）提出一种通用哈密顿量嵌入(Hamiltonian embedding)方案：为每个原始自由度引入辅助自由度，构造互易耦合的扩展哈密顿量，再施加初值约束(mirror constraint)使约化动力学严格还原原始非互易方程。研究人员证明该嵌入可用于定义约束Glauber MC动力学（等价于原始Langevin方程对应的Fokker–Planck方程），并借助辛结构做Floquet工程(Floquet engineering)调控非互易相互作用，为统计力学与哈密顿方法在非互易多体系统中的拓展奠定基础。

二、主要关键技术方法

研究人员针对一般两两非互易相互作用系统给出哈密顿嵌入通解证明，并以具视锥(vision-cone)相互作用的二维耗散XY自旋模型（及追逐-逃逸chase-and-run模型）为范例：(1)引入与原自旋θ_i共轭的辅助自旋φ_i（{θ_i, φ_j}=δ_ij），构造含原-原对称化相互作用H_SS与原-辅互易耦合H_Sa的总哈密顿量H=H_SS+H_Sa；(2)施加镜像约束θ_i?φ_i=π（等价于辅助自旋a_i=?S_i），证明该约束在哈密顿演化下守恒且约化出原始非互易Langevin方程（无噪声情形）；(3)基于约束哈密顿量定义Glauber转移概率ΔE=(1/2)[H(θ',φ)?H(θ,φ')]|_constr.，推导MC连续极限下的Fokker–Planck方程，证明其与原Langevin方程等价；(4)对周期驱动情形，将含时矢势项加入嵌入哈密顿量，用Floquet–Magnus展开得零阶有效哈密顿量H_F⁽⁰⁾，分析Bessel函数J₀(A_ij)导致的维度交叉(dimensionsal crossover)；(5)数值对比：约束Glauber MC与直接Langevin模拟比较序参量（磁化强度〈m〉、比热C_v、条纹序参量O(t)）及Floquet驱动下平均磁化随时间演化。

三、研究结果

XY spins with non-reciprocal vision-cone interactions（具非互易视锥相互作用的XY自旋）

研究人员在正方格上定义经典XY自旋S_i=(cos θ_i, sin θ_i)，近邻相互作用强度J_ij(θ_i)依赖施力自旋自身取向与视锥半角Ψ——仅当邻点落入自旋i前方视锥内时J_ij=J否则0。由于J_ij(θ_i)≠J_ji(θ_j)，相互作用非互易，Langevin方程为\dot{θ}_i=?Σ_jJ_ij(θ_i)sin(θ_i?θ_j)+√(2T)η_i(t)，无法从势函数导出。

Constrained Hamiltonian embedding（约束哈密顿量嵌入）

研究人员引入辅助自旋层φ_i，构造哈密顿量H=H_SS+H_Sa，其中H_SS=?Σ_?ij?[J_ij(θ_i)+J_ji(θ_j)]cos(θ_i?θ_j)为对称化原-原耦合，H_Sa=Σ_?ij?[?J_ij(θ_i)cos(θ_j?φ_i)?J_ji(θ_j)cos(θ_i?φ_j)]为原-辅互易耦合。取约束θ_i(0)?φ_i(0)=π（mirror constraint），可证该约束在Hamilton方程下保持不变，代入后\dot{θ}_i=Σ_jJ_ij(θ_i)sin(θ_i?θ_j)，严格还原原始非互易动力学。约束子流形测度为零但不与外围轨道相交，非互易分布可独立嵌入更大哈密顿相空间。H|_constr.≡0，仅为运动生成元而非系统能量。该构造对惯性与过阻尼(overdamped)非互易两两相互作用均适用。

Equivalence between Hamiltonian Monte-Carlo and Langevin dynamics subject to white noise（哈密顿蒙特卡洛与受白噪声朗之万动力学的等价性）

研究人员基于嵌入哈密顿量定义约束Glauber转移率w(θ_i→θ_i') = 1/2[1?tanh(ΔE(θ_i,θ_i')/(2T))]，ΔE=(1/2)[H(θ',φ)?H(θ,φ')]|_constr.，辅助变量在计算中消去。通过Master方程连续极限推导得Fokker–Planck方程?_tρ=αΣ_i?_θi[T?_θiρ?F_i(θ)ρ]，与原始Langevin方程所得Fokker–Planck方程一致（α=Δ2/(12T)为MC-物理时间标度因子），证明约束Glauber MC采样等价于原始非互易Langevin过程。

Non-equilibrium steady states（非平衡稳态）

对视锥XY模型（Ψ=0.85π，L=100），约束Glauber MC与Langevin模拟给出相同临界温度T_c/J=0.79±0.04（磁化强度?m?与比热C_v峰位一致）；而基于selfish energy的MC给出T_c/J≈0.86且有偏差。表明约束哈密顿嵌入MC正确复现非平衡稳态。

Non-stationary states（非稳态——追逐-逃逸chase-and-run相互作用）

对x方向具非互易J_→≠J_←、y方向互易J_rec的XY模型，初始竖直条带图案在J_←≈?J_→下产生持续左移振荡（时间平移对称性破缺）。约束Glauber MC与Langevin模拟的复序参量O(t)=L^?1Σ_ye^iθ_1,y(t)实部/虚部演化全程吻合，证实嵌入MC可捕捉非稳态动力学。

Hamiltonian engineering for non-reciprocal systems（非互易体系的哈密顿量工程）

利用嵌入哈密顿量的辛结构，研究人员对具周期驱动（矢势A_ij=aA₀cos ψ_ijcos ωt）的视锥XY模型作Floquet–Magnus展开，零阶有效哈密顿量H_F⁽⁰⁾含Bessel因子J₀(A_ij)调制各键耦合。调驱动幅A₀至J₀(z_α)=0处（z_α为J₀零点），x方向相互作用被抑制，系统发生二维正方格→一维独立链维度交叉，长程有序消失（与Mermin–Wagner定理一致）。Langevin模拟显示随A₀增大T_c降低，高频(ω?J)下有效方程与精确驱动动力学符合良好。证明哈密顿嵌入使Floquet工程可用于非互易多体系统相互作用的精细调控。

四、讨论与结论总结（翻译浓缩）

既往个别非互易系统（如Lotka–Volterra模型）曾有哈密顿描述，本文给出适用于任意两两非互易相互作用的普适哈密顿嵌入：加倍自由度并加互易原-辅耦合，初值约束还原原始动力学且在演化下保持。与Galley的双变分原理(final-time equality)不同，本文采用初值约束(initial-time constraint)自动沿轨线保持。该嵌入本身不引入新物理，但使非互易系统可用平衡/非平衡物理的哈密顿工具——快速MC采样绕过Langevin长瞬态、正则变换求慢变量、Floquet工程调控相图、未来可拓展至量子化嵌入及非互易热力学描述。结论：研究人员构建了含辅助自由度的约束哈密顿嵌入，证明约束Glauber MC等价于原始Langevin过程（含非稳态），并利用辛结构实现Floquet调控非互易相互作用，为将统计力学与哈密顿动力学推广至非互易体系提供通用框架。