通过收缩快速混合格劳伯动力学直至达到独特性

Zongchen Chen Georgia Institute of Technology. Email: chenzongchen@gatech.edu. Research supported in part by NSF grant CCF-2007022.    Kuikui Liu University of Washington. Email: liukui17@cs.washington.edu. Research supported in part by NSF grant CCF-1907845 and ONR-YIP grant N00014-17-1-2429.    Eric Vigoda University of California, Santa Barbara. Email: vigoda@ucsb.edu. Research supported in part by NSF grant CCF-2007022.
摘要

对于一般的反铁磁2自旋系统,包括加权独立集上的硬核模型和反铁磁伊辛模型,当无限正则树位于李等人(2013)的唯一性区域时,最大阶数Δ图上的分割函数有一个𝖥𝖯𝖳𝖠𝖲 此外,在树非唯一性区域中,Sly(2010)表明,除非𝖭𝖯=𝖱𝖯,否则不存在𝖥𝖯𝖱𝖠𝖲来估计配分函数。 算法结果遵循 Weitz (2006) 提出的相关衰减方法或 Barvinok (2016) 开发的多项式插值方法。 然而,运行时间只是常数Δ的多项式。 对于硬核模型,Anari 等人(2020)最近的工作在树唯一性区域建立了简单单点马尔可夫链的快速混合,称为格劳伯动力学。 我们的工作通过考虑固定顶点 v 对其他顶点的总成对影响(而不是其他顶点对 v 的总影响)来简化他们对格劳伯动力学的分析,从而扩展他们对所有 2 旋转模型的工作并改进了混合时间。

更重要的是,我们的证明将三种不同的算法方法联系在一起:我们证明了所谓的树递归与合适的势函数的收缩,这是建立 Weitz 相关衰减方法和 Barvinok 多项式插值方法效率的主要技术,也建立了格劳伯动力学的快速混合。 我们强调,这种联系适用于所有 2 自旋模型(反铁磁和铁磁),并且相关衰减或多项式插值方法的现有证明立即意味着格劳伯动力学的快速混合。 我们的证明利用图配分函数是韦茨自回避游走树配分函数的除数。 这一事实催生了用于分析顶点影响的新工具,并且可能对复零的研究具有独立意义。

1简介

在最大度数 Δ 的一般图中的近似计数问题的计算复杂性与度数 Δ 的无限规则树(或向上)上的统计物理相变之间建立了显着的联系。在更一般的情况下为 Δ)。 这种联系适用于二态反铁磁自旋系统——独立集上的硬核模型和伊辛模型是此类系统最有趣的例子。

给定一个 n 顶点图 G=(V,E),2 个自旋模型的配置是对顶点的自旋 0,12n 分配。 2-自旋系统由三个参数定义:边权重β,γ>0和顶点权重λ>0 边缘参数β控制相邻1-自旋之间相互作用的(相对)强度,γ对应于相邻0-自旋,并且λ 是应用于具有 1-自旋的顶点的外部场。

每个自旋配置 σ{0,1}V 都分配有一个权重

wG(σ)=βm1(σ)γm0(σ)λn1(σ),

其中,对于自旋s{0,1},ms(σ)=#{uvE:σu=σv=s}是具有自旋s的单色边的数量,n1(σ)=#{vV:σv=1}是具有自旋的顶点的数量1(按照标准,参数已标准化,因此我们可以避免两个额外的参数)。 自旋构型的吉布斯分布由 μG(σ)=wG(σ)ZG(β,γ,λ), 给出,其中 ZG(β,γ,λ)=σ{0,1}Vβm1(σ)γm0(σ)λn1(σ) 是配分函数。

有两个特别有趣的例子:硬核模型和伊辛模型。 β=0γ=1时,唯一具有非零权重的配置是G的独立集和独立集σ的权重是w(σ)=λ|σ|;这个例子被称为硬核模型,其中参数λ对应于逸度。

β=γ的情况下,重要的量是单色边的总数m(σ)=m0(σ)+m1(σ),配置的权重σw(σ)=βm(σ)λn1(σ);这是经典的伊辛模型,其中参数β对应于逆温度,λ是外部场(λ=1表示没有外场)。 请注意,当β>1时,模型是铁磁,因为相邻顶点更喜欢具有相同的自旋,而β<1反铁磁伊辛模型。 在一般的2自旋系统中,当βγ>1时模型是铁磁的,当βγ<1时模型是反铁磁的。 (当 βγ=1 时,我们得到一个简单的产品分布。)

基本算法任务是从吉布斯分布中采样并估计配分函数。 对于近似采样问题,我们给出一个图 G 和一个 ϵ>0,我们的目标是从总变异距离内的分布 π 生成样本ϵ 的吉布斯分布 μG 在时间 poly(n,log(1/ϵ)) 上的分布。 高效的近似采样算法意味着近似计数问题[JVV86,ŠVV09]𝖥𝖯𝖱𝖠𝖲(完全多项式随机近似方案)。 回想一下,给定一个 n 顶点图 Gϵ,δ>0,𝖥𝖯𝖱𝖠𝖲 输出一个 (1±ϵ) 近似值ZG 的概率 1δ 在时间 poly(n,1/ϵ,log(1/δ)) 中,而 𝖥𝖯𝖳𝖠𝖲 是确定性模拟(即 δ=0 )。

近似采样问题的标准方法是马尔可夫链蒙特卡罗 (MCMC) 方法;事实上,有一个简单的马尔可夫链,称为格劳伯动力学 格劳伯动力学的工作原理如下:从时间 t 的配置 Xt 中,选择一个随机顶点 v,然后设置 Xt+1(w)=Xt(w) 为所有wv,最后我们从μ(σv|σw=Xt+1(w) for all wv)的条件分布中选择Xt+1(v) 对于硬核模型的情况,如果当前没有邻居被占用,则Xt+1(v)被设置为被占用(即自旋1),概率为λ/(1+λ),否则它被设置为空闲。

很容易验证格劳伯动力学是遍历的,吉布斯分布是唯一的平稳分布。 混合时间是保证从最差初始状态X0开始,Xt的分布在总变化距离1/4 的吉布斯分布。 目标是证明混合时间是n中的多项式,在这种情况下,链被称为快速混合

对于铁磁 Ising 模型(有或没有外场)的情况,Jerrum 和 Sinclair [JS93] 的经典结果通过 MCMC 方法为所有图给出了 𝖥𝖯𝖱𝖠𝖲 这是通用图的有效算法的唯一情况。 对于反铁磁二自旋模型,该图与常规树上的统计物理相变密切相关。

硬核模型的情况很好地说明了唯一性/非唯一性相变。 考虑以 r 为根的无限 Δ-正则树 T,并让 Th 表示在第一个 h 此相变捕获 Th 叶子处的配置是否在限制 h 内“影响”根。 对于硬核模型,我们可以考虑偶数高度树(对应于全偶数边界条件)与奇数高度树。 ph表示吉布斯分布μTh中根被占据的边际概率。 p𝖾𝗏𝖾𝗇=limhp2hp𝗈𝖽𝖽=limhp2h+1 如果成立,我们就说树唯一性成立;如果p𝖾𝗏𝖾𝗇=p𝗈𝖽𝖽不相等,我们就说树非唯一性成立。 对于所有Δ3,都存在临界逸度λc(Δ)=(Δ1)Δ1/(Δ2)Δ)[Kel85],其中树唯一性保持当且仅当λλc(Δ)

显着的联系是最大度 Δ 的一般图的算法相变发生在同一树临界点。 对于所有常量Δ、所有δ>0、所有λ<(1δ)λc(Δ)、所有最大度图Δ[Wei06] 提出了一个用于近似配分函数的𝖥𝖯𝖳𝖠𝖲 另一方面,对于所有δ>0,所有λ>(1+δ)λc(Δ),[Sly10,SS14,GŠV16]证明,除非𝖭𝖯=𝖱𝖯,没有用于估计配分函数的𝖥𝖯𝖱𝖠𝖲

一个重要的警告是,Weitz 算法的运行时间为 (n/ϵ)ClogΔ,其中近似因子为 (1±ϵ),常数 C 多项式取决于间隙 δ(回想一下,λ<(1δ)λc)。 Sinclair 等人[SST14]将Weitz的相关衰减算法扩展到树唯一性区域中的反铁磁Ising模型,以及相应树唯一性区域中的所有反铁磁2-自旋系统(如下详述) ) 作者:Li、Lu 和 Yin [LLY13]

Barvinok [Bar16] 提出了一种有趣的新算法方法,并由 Patel 和 Regts [PR17] 改进,利用复平面中配分函数零点的缺失来使用泰勒近似有效地近似配分函数对数的合适变换。 这种多项式插值方法被证明在与 Peters 和 Regts [PR19] 的 Weitz 结果相同的树唯一性区域中是有效的,尽管运行时间的指数取决于 Δ

长期以来,人们推测简单的格劳伯动力学正在树的独特性区域中迅速混合。 Anari、Liu 和 Oveis Gharan 最近证明了这一点[ALO20];他们证明,对于所有δ>0,只要λ<(1δ)λc(Δ),混合时间就是nO(exp(1/δ)) 我们改进了这个结果。 首先,我们将混合时间从 nO(exp(1/δ)) 改进到 nO(1/δ),如以下定理所述。

Theorem 1 (硬核模型)

Δ3 为整数,δ(0,1) 对于最大度Δ的每个n-顶点图G和每个0<λ(1δ)λc(Δ),硬核模型的格劳伯动力学的混合时间在 G 上,逸度 λO(n2+32/δ)

除非进一步改进 [AL20] 中的局部到全局参数,此界限是最佳的。 我们改进的结果源于更简单、更清晰的证明方法,该方法使我们能够将结果扩展到各种 2-spin 模型,与相关衰减算法的关键结果相匹配,并大大缩短运行时间。

我们的证明方法统一了近似计数的三种主要算法工具:相关衰减、多项式插值和 MCMC。 相关衰减和多项式插值方法的大多数已知结果都是通过在所谓的树递归上显示适当定义的势函数的收缩来证明的;树递归是韦茨的自回避游走树的结果,我们将在本文后面更详细地描述。 Shao 和 Sun 最近的一项工作[SS20] 统一了这两种方法,表明通常用于证明相关衰减算法效率的收缩也意味着(在一些额外的分析条件下)多项式插值方法是有效的。

在这里,我们证明了势函数的相同收缩也意味着格劳伯动力学的快速混合,以及我们改进的独立于 Δ 的运行时间;详细说明请参见45 我们的证明利用了有关韦茨自回避游走树的几个新工具,这些工具在部分 3中有详细介绍。 特别地,我们证明了图G的配分函数除以Weitz自回避游走树的配分函数;请参阅8 这一结果对于确定具有复杂参数的配分函数不存在零点具有潜在的独立意义,因为它使人们能够考虑自回避游走树。 这一结果还产生了一个新的、有用的等价物,用于根据自回避树来限制图中的影响,这加强了 Weitz [Wei06] 先前已知的联系;详情请参见8

作为一个简单的结果,我们在树唯一性区域中获得了反铁磁 Ising 模型的 Glauber 动力学的快速混合。 就边活动而言,Δ-正则树上 Ising 模型的两个临界点位于 βc(Δ)=Δ2Δβ¯c(Δ)=1βc(Δ)=ΔΔ2;第一个位于反铁磁区域,而第二个位于铁磁区域。 如果βc(Δ)<β<β¯c(Δ),则唯一性适用于Δ-常规树上的所有外部字段λ

如前所述,对于铁磁 Ising 模型,𝖥𝖯𝖱𝖠𝖲 被称为通用图 [JS93] 此外,Mossel和Sly[MS13]证明了铁磁Ising模型的Glauber动力学在1β<β¯c(Δ)时的O(nlogn)混合时间。 然而,反铁磁伊辛模型在树独特性区域的快速混合尚不清楚。

我们为情况 β>βc(Δ) 提供以下混合结果。 请注意,当 ββc 时,外部字段 λ 的某些值有一个额外的唯一性区域;该区域由 3 覆盖。

Theorem 2 (反铁磁伊辛模型)

Δ3 为整数,δ(0,1) 假设1>ββc(Δ)+δ(1βc(Δ))λ>0 然后对于最大度Δ的每个n-顶点图G,伊辛模型的格劳伯动力学在G上的混合时间边权重为β,外场λO(n2+1.5/δ)

我们的硬核模型和伊辛模型的结果适合一般反铁磁二自旋系统的更大框架。 回想一下,反铁磁情况是在 βγ<1 时。

对于一般的 2-自旋系统,适当的树相变更加复杂,因为存在树唯一性阈值在 Δ 中不是单调的模型。 因此,正如 [LLY13] 所考虑的那样,适当的概念是“最高Δ唯一性”。 粗略地说,如果对于每个整数 1,距离 处的所有顶点,我们说间隙 δ(0,1) 的唯一性在 d 正则树上成立> 从根开始对根的边缘有总的“影响”(1δ) 如果带有间隙 δ 的唯一性在所有 d 常规树上成立,则我们说最多 Δ 带有间隙 δ 的唯一性成立1dΔ;准确定义请参见部分 2

12都是以下一般快速混合结果的推论,该结果适用于整个树独特性区域中的一般反铁磁2自旋系统。

Theorem 3 (一般反铁磁2-自旋系统)

Δ3 为整数,δ(0,1) β,γ,λ为实数,使得0βγγ>0βγ<1λ>0 假设参数 (β,γ,λ) 最多 Δ 是唯一的,间隙为 δ 然后对于最大度Δ的每个n-顶点图G,反铁磁2-自旋的格劳伯动力学的混合时间G 上带有参数 (β,γ,λ) 的系统是 O(n2+72/δ)

我们还匹配了铁磁 2-自旋模型的现有相关衰减结果[GL18, SS20];结果请参见部分8,以及附录F 为证明。

1.1电位法混合

树递归在近似计数的研究中非常有用。 考虑一棵以 r 为根的树。假设 rd 子级,用 v1,,vd 表示。 对于 1iΔi,我们将 Tvi 定义为以 vi 为根的 T 的子树,其中包含 vi 的所有后代。 Rr=μT(σr = 1)/μT(σr = 0) 表示根的边际比率,并且 Rvi=μTvi(σvi = 1)/μTvi(σvi = 0) 对于每个子树。 由于Tvi的独立性,树递归是一个在给定Rv1,,Rvd的情况下计算Rr的公式。 更具体地说,我们可以编写 Rr=Fd(Rv1,,Rvd),其中 Fd:[0,+]d[0,+] 是一个多元函数,对于 (x1,,xd)[0,+]d

Fd(x1,,xd)=λi=1dβxi+1xi+γ.

然而,在本文中,我们特别关注边际比率的对数。 原因是我们要仔细研究[ALO20]中介绍的吉布斯分布μG成对影响矩阵 G并定义为每个 r,vV

G(r v)=μG(σv = 1σr = 1)μG(σv = 1σr = 0).

[ALO20]中,作者表明,如果G的最大特征值有适当的界限,那么格劳伯动力学就是快速混合。 我们在本文中做出的一个重要观察是,影响 G(r v) v 上的 r 可以被视为 logRr 相对于 v 处的日志外部字段的导数(参见 12 因此,我们使用对数比率更方便。 为此,我们将树递归重写为 logRv=Hd(logRv1,,logRvd),其中 Hd:[,+]d[,+] 是一个函数,对于 (y1,,yd)[,+]d

Hd(y1,,yd)=logλ+i=1dlog(βeyi+1eyi+γ).

观察H=logFexp 此外,我们定义

h(y)=(1βγ)ey(βey+1)(ey+γ)

对于 y[,+],以便 yiHd(y1,,yd)=h(yi) 对于每个 i

为了证明我们的主要结果,我们使用了已广泛用于建立相关性衰减的势方法。 通过为对数比率选择合适的势函数,我们表明给定顶点的总影响随着距离呈指数衰减,从而建立了格劳伯动力学的快速混合。 让我们首先明确我们对潜力的要求。 对于每个整数 d0,我们定义一个有界区间 Jd,其中包含 d 度顶点处的所有对数比率。更具体地说,当 βγ<1 时,我们让 Jd=[log(λβd),log(λ/γd)] ;当 βγ>1 时,我们让 Jd=[log(λ/γd),log(λβd)] 此外,将J=d=0Δ1Jd定义为包含次数小于Δ的所有对数比率的区间。

Definition 4 ((α,c)-势函数).

Δ3 为整数。 β,γ,λ为实数,使得0βγγ>0λ>0 Ψ:[,+](,+) 为具有图像 S=Ψ[,+] 和导数 ψ=Ψ 的可微且递增函数。 对于任意 α(0,1)c>0,如果 Ψ 满足以下条件,我们就说它是关于 Δ(β,γ,λ)(α,c) 势函数

  1. 1.

    (收缩)对于每个整数 d 使得 1d<Δ 和每个 (y~1,,y~d)Sd,我​​们有

    HdΨ(y~1,,y~d)1=i=1dψ(y)ψ(yi)|h(yi)|1α

    其中 HdΨ=ΨHdΨ1yi=Ψ1(y~i)1id,以及y=Hd(y1,,yd)

  2. 2.

    (有界)对于每个y1,y2J,我们有

    ψ(y2)ψ(y1)|h(y1)|cΔ.

(α,c)-势的定义中,应该将y视为顶点处的对数边际比,势函数为logR 以下定理在给定 (α,c)-势函数的情况下建立了格劳伯动力学的快速混合。

Theorem 5

Δ3 为整数。 β,γ,λ为实数,使得0βγγ>0λ>0 假设对于某些 α(0,1)c>0,存在相对于 Δ(β,γ,λ)(α,c) 势。 然后对于最大度Δ的每个n-顶点图G,2-自旋系统的格劳伯动力学的混合时间带有参数 (β,γ,λ)GO(n2+c/α)

我们在部分3中概述了我们的证明。 请注意,在 45 中,常数 c 允许依赖于最大度数 Δ 和参数 (β,γ,λ) 总的来说。 例如,对 [LLY13] 中的势进行简单的黑盒应用将为 有界性 条件提供 c=Θ(Δ),从而产生 nΘ(Δ) 混合。 然而,这对于具有潜在无限度数的图来说是不可取的。 我们的贡献之一在于,我们证明了 Boundedness 条件在一个通用常数 c Δ(β,γ,λ) 无关的情况下成立。 因此,我们的混合时间是 O(n2+c/δ),除了 1/δ 之外,指数中没有任何参数。

部分7中,我们给出了(α,c)-势的稍微更一般的定义,它放宽了有界性 条件,并且对于我们分析具有 0β<1<γ 的反铁磁 2 自旋系统是必要的。 5 对于这一更大类的势仍然成立。

我们注意到,在潜在方法的所有先前工作中,结果和证明总是以 FdR 的树递归和 Φ 的形式呈现, R 的潜在函数。事实上,我们的结果也可以翻译成(Fd,Φ)的语言。 要看到这一点,从Hd=logFdexp开始,如果我们选择Φ=Ψlog,则可以直接检查HdΨ=ΨHdΨ1=ΦFdΦ1=FdΦ,从而检查HdΨ=FdΦ 这意味着当且仅当相应的收缩条件对于(Fd,Φ)成立时,4中的收缩条件对于(Hd,Ψ)成立。 有界性条件也可以等效地表述为(Fd,Φ) 尽管如此,在本文中我们选择使用 (Hd,Ψ) 出于以下两个原因。 首先,正如前面提到的,事实是 G(r v) logRr 的导数,因此可以很自然地考虑对数比率的树递归。 事实上,直接使用 (Hd,Ψ) 呈现我们的结果和证明比切换到 (Fd,Φ) 更容易、更清晰。 其次,我们将使用的势函数 Ψ 是通过转换 Ψ=Φexp[LLY13] 中的精确势 Φ 获得的。111更准确的说,我们还乘了一个常数因子,它只是简化了我们的计算,并没有多大关系;另请注意,[LLY13] 表示 φ 的势函数及其 Φ=φ 的导数。 有趣的是,这个势能的导数就是ψ=|h| 那么收缩条件有一个很好的形式:i=1dh(y)h(yi)1αBoundedness 条件仅涉及 h(y) 的上限。 这似乎揭示了 [LLY13] 中神秘的势函数 Φ,并且还表明 Hd 是树递归的一个有意义的变体考虑。 再添加一个证据,对于很多情况(例如,Δ2Δ<βγ<ΔΔ2),选择了潜在的Φ=log,这仅仅意味着我们可以选择Ψ恒等函数和 Hd 本身正在收缩,没有任何不平凡的潜力。

2021 年 7 月修订。

这篇论文在 FOCS 2020 上发表后,在[LLY13]中发现了关于反铁磁二自旋系统唯一性区域描述的一个小错误。 该错误已在最新版本的[LLY13]中修复。 在本次修订中,我们更新了部分7附录E中相应的结果和证明 依赖于 [LLY13] 中的更改;特别地,36根据当前唯一性区域的描述进行调整。 我们指出,这些变化纯粹是技术性的,不会影响我们主要结果(例如5)的有效性。

致谢。

我们要感谢 Shayan Oveis Gharan 和 Nima Anari 激发了讨论。 我们也感谢匿名审稿人提出的有益的意见和建议。 我们感谢 Yitong Yin 与我们沟通[LLY13]的最新更新,并提供有关修改附录中的陈述和结果证明的有用说明E,特别是36

2 预赛

混合时间和光谱间隙

P 为具有平稳分布 μ 的有限状态空间 Ω 上遍历(即不可约且非周期)马尔可夫链的转移矩阵。 Pt(x0,)表示从x0Ω开始的t步之后链的分布。 P混合时间定义为

Tmix(P)=maxx0Ωmin{t0:Pt(x0,)μ()TV14}.

如果μ(x)P(x,y)=μ(y)P(y,x)对于所有x,yΩ,我们说P可逆的 如果P是可逆的,则P只有实特征值,可用1=λ1λ|Ω|1表示。 P光谱间隙定义为1λ2,P绝对光谱间隙为定义为λ(P)=1max{|λ2|,|λ|Ω||} 如果 P 相对于内积 ,μ 也是半正定的,则​​ P 的所有特征值都是非负的,因此 λ(P)=1λ2 也是非负的。 最后,混合时间和绝对光谱间隙的关系为

Tmix(P)1λ(P)log(4minxΩμ(x)). (1)

独特性

Δ3 为整数或 Δ= β,γ,λ为实数,使得0βγγ>0βγ<1λ>0 对于 1d<Δ,定义

fd(R)=λ(βR+1R+γ)d

并用Rd表示fd的唯一不动点。 对于δ(0,1),我们说参数(β,γ,λ)up-to-Δ唯一的,间隙δ如果|fd(Rd)|<1δ 对于所有 1d<Δ

比例及影响力

考虑图 G=(V,E) 上的 2-自旋系统。 ΛVσΛ{0,1}Λ 对于所有vV\Λ,我们将v处的边际比率定义为

RGσΛ(v)=μG(σv = 1σΛ)μG(σv = 0σΛ).

对于所有u,vV\Λ,我们定义uv(成对)影响

GσΛ(u v)=μG(σv = 1σu = 1,σΛ)μG(σv = 1σu = 0,σΛ).

(成对)影响矩阵编写GσΛ,其条目由下式给出 GσΛ(u v).

韦茨的自回避行走树

G=(V,E) 为连通图,rVG 的顶点。 自回避行走(SAW)树定义如下。 假设顶点集V存在全序。从 r 出发的自回避步行是一条路径 r=v0v1v,使得所有 0i<j 都具有 vivj SAW树Tsaw(G,r)是一棵以r为根的树,由所有带有deg(v)=1的自回避行走r=v0v1v以及附加的关闭循环的又一个顶点(即 r=v0v1vvi 对于某些 0i2 ,例如 {v,vi}E)。 请注意,G 的顶点可能在 SAW 树中具有多个副本,并且除了叶子之外,顶点的度数将被保留。 请参阅1了解示例。

我们可以使用相同的参数(β,γ,λ)Tsaw(G,r)上定义一个2-自旋系统,其中一些叶子被固定到特定的自旋。 更具体地说,对于附加了 vi 的自回避行走 r=v0v1v,我们将 vi 修复为旋转 1 如果 vi+1<v 相对于 V 上的总排序,如果 vi+1>v 则旋转 0 对于每个 vV,我们用 𝒞v 表示 Tsaw(G,r)v 的所有自由(未固定)副本的集合。 对于 ΛV 和部分配置 σΛ{0,1}Λ,我们通过将自旋 σv 分配给 v 的每个副本 v^𝒞v 开始,并删除 v^ 的所有后代,为每个 vΛ 定义具有条件 σΛ 的 SAW 树。 请注意,一般来说,来自 𝒞vv 的不同副本可以接收不同的自旋分配。 最后,如果每个顶点 v 都有不同的字段 λv,则来自 𝒞vv 的所有副本都将具有相同的字段SAW 树中的 λv

Refer to caption
图1 G 和以 r 为根的自回避游走树 Tsaw(G,r)Tsaw(G,r) 中具有相同标签的顶点是 G 中相同顶点的副本。 (/:固定为旋转1/0。)

3 主要结果的证明大纲

步骤 1 ([ALO20]):光谱独立意味着快速混合。

我们的证明建立在[ALO20]的基础上,他证明了从最大程度图上的硬核分布中采样的格劳伯动力学最多Δ混合在O(nexp(O(1/δ)))步中每当 λ(1δ)λc(Δ) 时。 他们证明的关键要素之一是他们称之为谱独立性的概念。 [ALO20]表明光谱独立性意味着快速混合。 请注意,GσΛ 的对角线条目是 1,而不是 [ALO20] 中原始定义中的 0

Definition 6 (谱独立性[ALO20])

我们说 n 顶点图 G 上的吉布斯分布 μG(η0,,ηn2) 谱无关,如果对于大小为 kσΛ{0,1}Λ 的每个 0kn2ΛV,都有 λmax(GσΛ)1ηk

Theorem 7 ([ALO20]).

如果 μ 是一个 (η0,,ηn2) 谱独立分布,则从 μ 采样的格劳伯动力学至少具有谱间隙

1ni=0n2(1ηini1).

我们现在的主要目标是限制 GσΛ 的最大特征值。

步骤 2:自我回避行走树保留影响。

从标准线性代数,我们知道 GσΛ 的最大特征值的上限是 1-范数 GσΛ10>=1>maxrεVΣvεV|GσΛ(v r)|,对应于顶点 r 的总影响,以及无穷范数 GσΛ∞0>=1>maxrεVΣvεV|GσΛ(r v)|,对应于r的总影响。在[ALO20]中,作者使用GσΛ1作为λmax(GσΛ)的上限。 粗略地说,他们证明了在所有边界条件下,固定顶点 r 上的绝对影响总和的上限是以 r 为根的自避让行走树中 r 上的最大绝对影响。 在本文中,我们将使用 GσΛ 来代替上限 λmax(GσΛ) 事实上,如果我们看看自回避树中 r 的影响,则事实更多。 我们证明对于每个顶点vV,影响 GσΛ(r v) G中的内容以影响力总和的形式保存在以r为根的自回避游走树T=Tsaw(G,r) TσΛ(r v^) v 的所有副本 v^

我们建立这一事实的方法是将配分函数视为 λ 中的多项式。 事实上,考虑每个 vV 具有任意外部字段 λv 的更一般情况将很有用。 𝝀={λv:vV} 表示字段。 对于ΛVσΛ{0,1}Λ,以σΛ为条件的σ{0,1}V\Λ的权重定义为wG(σσΛ)=βm1(σσΛ)γm0(σσΛ)vV\Λλvσv,其中mi(σΛ)i=0,1V\Λ 中至少有一个端点的 i-i 边的数量。 此外,ZGσΛ=σ{0,1}V\ΛwG(σσΛ)是以σΛ为条件的配分函数。 我们将βγ视为一些固定常量,并将𝝀视为n=|V|变量。 从这个意义上说,我们将权重wG(σσΛ)视为𝝀中的单项式,将配分函数ZGσΛ视为𝝀中的多项式。 此外,边际比率RGσΛ(v)及其影响 GσΛ(r v) for r,vV 都是 𝝀 中的函数。 我们的主要结果是 G 的配分函数除以每个 rVTsaw(G,r) 的配分函数。 由此,我们表明 SAW 树保留了根的影响,并重新建立了 Weitz 的著名结果[Wei06],请参阅13

Lemma 8.

G=(V,E) 为连通图,rV 为顶点,ΛV\{r} 使得 G\Λ 连通。 T=Tsaw(G,r)为以r为根的G的自回避游走树。然后对于每个σΛ{0,1}Λ,ZGσΛ除以ZTσΛ 更准确地说,存在一个独立于 λr 的多项式 PG,rσΛ=PG,rσΛ(𝛌),使得

ZTσΛ=ZGσΛPG,rσΛ. (2)

作为推论,对于每个顶点 vV

GσΛ(r v)=v^𝒞vTσΛ(r v^), (3)

其中 𝒞vTv 的所有自由(未固定)副本的集合。

Remark 1

我们强调,为了限定 G 中顶点的总影响,只需要 Eq. 38,这可以用纯粹的组合方式证明。 然而,我们相信G的多元配分函数的整除性Eq. 2及其自回避游走树可能具有独立的兴趣。

我们注意到,整除性语句 Eq. 2 的单变量版本已经出现在 [Ben18] 中配分函数复数根研究中的硬核模型和零场伊辛模型的[LSS19] 8,我们可以得到 ΣvεV|GσΛ(r v)|ΣvεVT |TσΛ(r v)| 对于任何固定的r。这意味着,我们只需要对树的所有影响力之和设置上限,即可获得 λmax(GσΛ) 的上限。

第三步:给予良好潜力的影响力的衰减。

树递归为我们提供了一个很好的工具来递归计算树的顶点(对数)比率。 正如我们在 12 中所示,影响 GσΛ(r v) 实际上是 r 处对数边际比的导数版本。因此,可以自然地使用树递归来关联这些影响。 然后,我们应用势方法(该方法已在文献中广泛使用)来建立相关性的衰减(强空间混合)。 下面的引理表明,对距离 k 的绝对影响之和随 k 呈指数衰减,可以将其视为成对影响的衰减。

Lemma 9.

如果存在关于Δ(β,γ,λ)(α,c)-潜在函数Ψ,其中α(0,1)c>0,然后对于每个 ΛVT\{r}σΛ{0,1}Λ 和所有整数 k1

vLr(k)|TσΛ(r v)|c(1α)k1

其中Lr(k)表示距r距离为k的所有自由顶点的集合。

然后通过组合789来证明5 我们将其证明留给附录A

第四步:找到好的潜力。

作为最后一步,我们需要找到 4 中定义的 (α,c) 势函数。 我们选择的潜在Ψ正是来自[LLY13]的那个,适应对数边际比率和树递归H(参见部分6了解更多详情)。 我们证明,如果参数 (β,γ,λ) 最多 Δ 是唯一的,且有间隙 δ(0,1) 以及 βγ>Δ2Δγ1,则 Ψ(α,c) 势。

Lemma 10.

Δ3 为整数。 β,γ,λ为实数,使得0βγγ>0βγ<1λ>0 假设 (β,γ,λ) 最多 Δ 是唯一的,间隙为 δ(0,1) 隐式定义函数 Ψ

Ψ(y)=ψ(y)=(1βγ)ey(βey+1)(ey+γ)=|h(y)|,Ψ(0)=0. (4)

如果βγ>Δ2Δ,则Ψ是具有αδ/2c1.5(α,c)势函数。 如果 βγΔ2Δγ1,那么 Ψαδ/2c18 是一个 (α,c) 电位;如果 β=0,我们可以进一步取 c4

我们从 510 推断出 βγ>Δ2Δγ1 情况的 3 其证明可参见附录A βγΔ2Δγ>1 的情况比较棘手。 正如[LLY13]第5节所讨论的,当βγΔ2Δγ>1时,对于某些λ>0,自旋系统在于唯一性任意图的区域,即使具有无界度(即,最多 唯一)。 因此,在这种情况下,顶点的总影响可能与 Θ(Δ/δ) 一样大,从而导致 nΘ(Δ/δ) 混合时间。 为了解决这个问题,我们考虑固定顶点的绝对影响的适当加权总和,这也是影响矩阵的最大特征值的上限。 然后将 45 修改为稍强的版本。 本案的陈述及证据载于部分7附录D

本文的其余部分安排如下。 部分4中,我们证明了8关于SAW树的属性。 部分 5中,我们建立了关于势法影响衰减的9 我们验证部分6中的收缩条件以选择势能。 部分 7专门讨论βγΔ2Δγ>1的情况,其中更通用的版本需要 45;缺失的证明可以在附录0>D1>中找到。在附录3>E4>2>中,我们验证了有界性5>条件及其在所有情况下对我们潜力的概括。 我们在部分 8中考虑铁磁自旋系统,证明留给附录 F 。我们在附录A中证明了我们的所有主要结果。

4 保留自回避行走树的影响

在本节中,我们将展示[Wei06]中介绍的自回避行走(SAW)树(另见[SS05]),保留了根的所有影响,从而建立8 为此,我们将 G 的配分函数(视为外部字段 𝝀 的多项式)除以 SAW 树的配分函数。 由此我们证明根顶点 rG 中另一个顶点 v 的影响与对 v 使用我们的证明方法,我们表明根的边缘保留在 SAW 树中,重新建立了 Weitz 著名的结果 [Wei06],以及与 v 相关的所有成对协方差> 被保留。

Theorem 11.

G=(V,E) 为连通图,rV 为顶点,ΛV\{r} 使得 G\Λ 连通。 T=Tsaw(G,r)为以r为根的G的自回避游走树。然后对于每个σΛ{0,1}Λ,ZGσΛ除以ZTσΛ 更准确地说,存在一个多项式 PG,rσΛ=PG,rσΛ(𝛌) 使得

ZTσΛ=ZGσΛPG,rσΛ.

此外,多项式PG,rσΛ独立于λr

Remark 2

11 的证明可以进行调整,以给出 Eq. 38 的纯组合证明。 就像在[Wei06,定理3.1]的证明中一样,我们可以通过顶点分裂和伸缩来进行,其中不是伸缩边际比率的乘积,而是伸缩单顶点影响的总和。

我们注意到,[Ben18] 证明了硬核模型的 11 的单变量版本,而 [LSS19] 对于零模型显示了类似的结果。具有均匀边缘权重的伊辛场模型。 我们的结果适用于所有 2 自旋系统和每个顶点的任意字段。 我们还可以以简单的方式将其推广到每条边的任意边权重。 至关重要的是,商多项式PG,rσΛ与根处的场λr无关,由此我们可以立即推导出根的边际保留和影响。

在证明11之前,我们首先给出它的一些结果。 对于所有 u,vV\Λ,我们将 v 处的 marginal 定义为 MGσΛ(v0>)1>=2>μG(v = 1σΛ) (为了方便起见,我们将事件 σv=i 写为 v=i),以及 uv协方差 > 作为

KGσΛ(u,v)=μG(u = v = 1σΛ)μG(u = 1σΛ)μG(v = 1σΛ).

以下引理将我们感兴趣的量与(对数)配分函数的适当导数联系起来。 引理的第一部分和第二部分是民间传说。

Lemma 12.

对于每个图 G=(V,E)ΛVσΛ{0,1}Λ,以下内容成立:

  1. 1.

    对于所有vV,

    (λvλv)logZGσΛ=MGσΛ(v);
  2. 2.

    对于所有u,vV,

    (λvλv)(λuλu)logZGσΛ=(λvλv)MGσΛ(u)=KGσΛ(u,v);
  3. 3.

    对于所有u,vV,

    (λvλv)logRGσΛ(u)=GσΛ(u v).
证明。

前两部分是标准的。 为了完整性,我们将这两个事实的证明包含在附录B中。 对于第 3 部分,我们从第 2 部分推断出

(λvλv)logRGσΛ(u)=(λvλv)log(MGσΛ(u)1MGσΛ(u))=(λvλv)MGσΛ(u)MGσΛ(u)(1MGσΛ(u))=KGσΛ(u,v)KGσΛ(u,u).

仍有待证明的是

GσΛ(u v)=KGσΛ(u,v)KGσΛ(u,u),

这实际上适用于任何两个二元随机变量。 为了看到这一点,我们首先计算 KGσΛ(0>u1>,2>u3>)4>⋅5>ℐ7>G8>σ0>Λ1>9>6>​2>(u v) 根据定义:

KGσΛ(u,u)GσΛ(u v)
= μG(u = 1σΛ)μG(u = 0σΛ)[μG(v = 1u = 1,σΛ)μG(v = 1u = 0,σΛ)]
= μG(u = 1,v = 1σΛ)μG(u = 0σΛ)μG(u = 1σΛ)μG(u = 0,v = 1σΛ)
= μG(u = 1,v = 1σΛ)μG(u = 0,v = 0σΛ)μG(u = 1,v = 0σΛ)μG(u = 0,v = 1σΛ).

同时,协方差可以写为

KGσΛ(u,v) =μG(u = 1,v = 1σΛ)μG(u = 1σΛ)μG(v = 1σΛ)
=μG(u = 1,v = 1σΛ)μG(u = 0,v = 0σΛ)μG(u = 1,v = 0σΛ)μG(u = 0,v = 1σΛ).

这表明 GσΛ(u v)=KGσ Λ0>(2>u3>,4>v5>)6>1>/7>K9>G0>σ2>Λ3>1>8>​4>(6>u7>,8>u9>)0>5> 从而确立了第 3 部分。

我们从 11 和以下引理的第二项推导出 8 11的证明在4.1中给出。

Lemma 13.

G=(V,E) 为连通图,rV 为顶点,ΛV\{r} 使得 G\Λ 连通。 T=Tsaw(G,r)为以r为根的G的自回避游走树。那么对于每个 σΛ{0,1}Λ 我们有:

  1. 1.

    ([Wei06,定理3.1])根边缘的保留 r:

    MGσΛ(r)=MTσΛ(r)andRGσΛ(r)=RTσΛ(r);
  2. 2.

    保留协方差和影响r:对于每个vV,

    KGσΛ(r,v)=v^𝒞vKTσΛ(r,v^)andGσΛ(r v)=v^𝒞vTσΛ(r v^).

    其中,𝒞vvT.中的所有自由(非固定)副本的集合。

证明。

11,存在多项式PG,rσΛ=PG,rσΛ(𝝀),使得ZTσΛ=ZGσΛPG,rσΛPG,rσΛ独立于λr 那么从 12 可以得出

MTσΛ(r)=(λrλr)logZTσΛ=(λrλr)(logZGσΛ+logPG,rσΛ)=(λrλr)logZGσΛ=MGσΛ(r),

因此RTσΛ(r)=RGσΛ(r) 对于第二项,再次从 12 我们得到

KGσΛ(r,v)=(λvλv)MGσΛ(r)=(λvλv)MTσΛ(r).

回想一下,对于 SAW 树 T 上的自旋系统,来自 𝒞vv 的每个自由副本 v^ 都具有相同的外部字段 λv^=λv 然后,根据导数的链式法则和12,我们推导出

KGσΛ(r,v)=v^𝒞v(λv^λv^)MTσΛ(r)λv^λvλvλv^=v^𝒞vKTσΛ(r,v^).

最后,我们有

GσΛ(r v)=(λvλv)logRGσΛ(r)=(λvλv)logRTσΛ(r)=v^𝒞vTσΛ(r v^),

其中最后一个等式如下。

4.1 11的证明

在提出证明之前,让我们首先回顾一下之前介绍的符号和定义。 𝝀={λv:vV} 表示所有顶点处的字段集。 对于ΛVσΛ{0,1}Λ,以σΛ为条件的σ{0,1}V\Λ的权重由下式给出

wG(σσΛ)=βm1(σσΛ)γm0(σσΛ)vV\Λλvσv,

其中对于 i=0,1,mi(σΛ) 表示两个端点都接收自旋 i 且至少其中一个位于 V\Λ 中的边数>。 σΛ为条件的配分函数定义为ZGσΛ=σ{0,1}V\ΛwG(σσΛ) 对于 SAW 树,我们以相同的方式定义条件权重和配分函数。 特别是,回想一下,当我们在 SAW 树上修复条件 σΛ 时,我们还会删除每个 vΛv^𝒞v 的所有后代。

对于每个 vV\Λi{0,1},我们应编写 v=i 来表示配置集,使得 σv=i (即 {σ{0,1}V\Λ:σv=i}) 并让 ZGσΛ(v = ) 是所有配置的权重之和 v = . 我们进一步扩展这个符号并写成 ZGσΛ(U = σU) 对于每个 UV\ΛσU{0,1}U 对于 SAW 树,我们也采用相同的符号。

11 的证明。

我们将证明存在一个独立于 λr 的多项式 PG,rσΛ=PG,rσΛ(𝝀),使得

ZTσΛ(r = 1)=ZGσΛ(r = 1)PG,rσΛandZTσΛ(r = 0)=ZGσΛ(r = 0)PG,rσΛ. (5)

方程 5的高级证明思路与[Wei06,定理3.1]中的相应结果类似>。 mV\Λ 中至少有一个端点的边数。 我们对 m 使用归纳法。当 m=0 时,该语句自 T=G 起就很简单。 假设 Eq. 5 对于所有图形和所有小于 m 边的条件成立。 假设根 rd 个邻居 v1,,vd 定义G为将顶点r替换为d顶点r1,,rd,然后连接{ri,di}得到的图对于1id

首先考虑(G\{r})\Λ仍然连接的情况。 对于每个i,让Gi=Gri 使用相同的参数 (β,γ,𝝀)Gi 上定义 2-自旋系统,加上顶点 r1,,ri1 固定为自旋的附加条件0ri+1,,rd固定为自旋1;我们用 σUiUi={v1,,vd}\{vi} 来表示这种条件。 然后,可以通过以下递归过程生成T=Tsaw(G,r) 另请参阅2了解说明。

Refer to caption
图2 自回避行走(SAW)树的递归构造。 这里TiGi的SAW树,其根为i=1,2,3vi (/:固定为旋转1/0。)
算法:Tsaw(G,r)
  1. 1.

    对于每个i,让Ti=Tsaw(Gi,vi)加上条件σUi

  2. 2.

    通过将 {r,vi} 连接为 1id,使 T=Tsaw(G,r) 成为 rT1,,Td 的并集;输出T

为了证明的目的,我们还考虑 G 上具有相同参数 (β,γ,𝝀)2-自旋系统,但我们让顶点 r1,,rd 没有字段(即为 1id 设置 λri=1 而不是 λr)。 然后我们观察到

ZGσΛ(r = 1)=λrZGσΛ(r1 = 1,,rd = 1),

将自旋 1 替换为 0 也是如此。 对于1id,让σΛi表示条件σΛσUi的并集,其中Λi=ΛUi 然后对于每个 1id 我们有

ZGσΛ(r1 = 0,,ri1 = 0,ri = 1,,rd = 1)=βZGiσΛi(vi = 1)+ZGiσΛi(vi = 0).

请注意,两侧都独立于字段λr:对于左侧,所有riG上都没有用于旋转系统的字段;对于右侧,请记住,我们不计算每个 Gi 的条件分区函数的固定顶点的权重。 现在定义 QG,rσΛ=QG,rσΛ(𝝀) 通过

QG,rσΛ=i=2dZGσΛ(r1 = 0,,ri1 = 0,ri = 1,,rd = 1),

它独立于 λr 然后我们得到

ZGσΛ(r = 1)QG,rσΛ =λri=1dZGσΛ(r1 = 0,,ri1 = 0,ri = 1,,rd = 1)
=λri=1d(βZGiσΛi(vi = 1)+ZGiσΛi(vi = 0)).

使用类似的论证,我们也有

ZGσΛ(r = 0)QG,rσΛ =i=1dZGσΛ(r1 = 0,,ri = 0,ri+1 = 1,,rd = 1)
=i=1d(ZGiσΛi(vi = 1)+γZGiσΛi(vi = 0)).

由于我们假设 (G\{r})\Λ 是连通的,因此对于每个 i 来说,图 Gi\Λ 也是连通的。 然后,根据归纳假设,对于每个 i 都存在一个多项式 PGi,viσΛi=PGi,viσΛi(𝝀) 使得

ZTiσΛi(r = 1)=ZGiσΛi(r = 1)PGi,viσΛiandZTiσΛi(r = 0)=ZGiσΛi(r = 0)PGi,viσΛi;

这些多项式独立于 λr,因为 Gi 的条件分配函数不涉及 λr 现在如果我们让

PG,rσΛ=QG,rσΛi=1dPGi,viσΛi,

那么从树递归可以得出

ZTσΛ(r = 1) =λri=1d(βZTiσΛi(vi = 1)+ZTiσΛi(vi = 0))
=λri=1d(βZGiσΛi(vi = 1)PGi,viσΛi+ZGiσΛi(vi = 0)PGi,viσΛi)
=ZGσΛ(r = 1)QG,rσΛi=1dPGi,viσΛi
=ZGσΛ(r = 1)PG,rσΛ.

另一个平等 ZTσΛ(r = 0)=ZGσΛ(r = 0)·PG,r σΛ 以同样的方式建立。 这样就完成了(G\{r})\Λ已连接的情况的证明。

如果(G\{r})\Λ有两个或多个连通分量,那么我们可以通过每个分量的SAW树构造Tsaw(G,r) 回想一下,G 是通过将顶点 r 拆分为图中 G 中的 d 个副本来定义的。假设 G\Λ 具有整数 k2k 连接分量。 G(1),,G(k) 为每个组件导出的子图,以及来自与其相邻的 Λ 的顶点。 对于每个 j,让 G(j) 是通过将 r 的所有副本收缩为一个顶点 r(j)G(j) 得到的图形,并让 T(j)=Tsaw(G(j),r(j)). 观察一下,一旦我们收缩 T(1),,T(k) 的根 r(1),,r(k),生成的树就是 Tsaw(G,r)

我们在每个 G(j) 上使用相同的参数 (β,γ,𝝀) 定义 2-spin 系统,只是顶点 r(j) 没有字段(即 λr(j)=1 而不是 λr)。 对于1jk,令Λ(j)=ΛV(G(j))σΛ(j)为限制在Λ(j)上的配置σΛ 然后 G(j)\Λ(j) 与每个 j 连接,并且由于 k2,每个具有条件 σΛ(j)G(j) 的数量少于m 边缘。 因此,我们可以应用归纳假设;即,对于1jk,存在一个独立于λr的多项式PG(i),r(i)σΛ(j)=PG(i),r(i)σΛ(j)(𝝀),使得

ZT(j)σΛ(j)(r(j) = 1)=ZG(j)σΛ(j)(r(j) = 1)PG(j),r(j)σΛ(j)andZT(j)σΛ(j)(r(j) = 0)=ZG(j)σΛ(j)(r(j) = 0)PG(j),r(j)σΛ(j).

我们将多项式 PG,rσΛ=PG,rσΛ(𝝀) 定义为

PG,rσΛ=j=1kPG(j),r(j)σΛ(j).

然后很容易检查

ZTσΛ(r = 1) =λrj=1kZT(j)σΛ(j)(r(j) = 1)=λrj=1k(ZG(j)σΛ(j)(r(j) = 1)PG(j),r(j)σΛ(j))
=ZGσΛ(r = 1)j=1kPG(j),r(j)σΛ(j)=ZGσΛ(r = 1)PG,rσΛ,

和类似地 ZTσΛ(r = 0)=ZGσΛ(r = 0)·PG,r σΛ. 定理如下。

5 对树木的影响范围

在本节中,我们研究根对树中其他顶点的影响。 我们给出根对固定距离处所有顶点的总影响的上限。 为此,我们应用了势方法,该方法已用于建立相关衰减属性(例如,参见[LLY12,LLY13,GL18])。 给定任意势函数 Ψ,我们的上限是根据 Ψ 的属性,涉及 HdΨ1|ψ| 上的界限,其中 ψ=Ψ 然后,在 Ψ(α,c) 势的情况下,我们推导出 9

假设T=(VT,ET)是一棵以r为根、最大度数至多Δ的树。 ΛVT\{r}σΛ{0,1}Λ 是任意且固定的。 考虑 T 上的 2-自旋系统,参数为 (β,γ,λ),以 σΛ 为条件。 我们需要限制影响力 TσΛ(r v) 从根r到另一个顶点vVT 请注意,如果在删除 Λvr 断开连接,则 TσΛ(r v)=0 由自旋系统的马尔可夫性质。 因此,我们可以假设,通过删除所有此类顶点,Λ 只包含 T 的叶子。

对于顶点 vVT,令 Tv=(VTv,ETv) 为以 v 为根的 T 的子树,其中包含 v 的所有后代>;请注意Tr=T 我们将 Lv(k)VT\Λ 写为 Tv 中距 v 距离为 k 的所有自由顶点的集合。 我们特别关注子树Tvv处的边际比率,为了简单起见,写成Rv=RTvσΛ(v) logRv 通过树递归 H 关联。如果顶点 vd 子节点,用 u1,,ud 表示,则树递归由下式给出

logRv=Hd(logRu1,,logRud),

其中对于 1dΔ(y1,,yd)[,+]d

Hd(y1,,yd)=logλ+i=1dlog(βeyi+1eyi+γ).

另请记住,对于 y[,+],我们定义

h(y)=(1βγ)ey(βey+1)(ey+γ)

和所有 1idΔyiHd(y1,,yd)=h(yi)

下面的引理允许我们使用任意势函数将所有影响的总和从根限制到距离 k

Lemma 14.

Ψ:[,+](,+) 为具有图像 S=Ψ[,+] 和导数 ψ=Ψ 的可微且递增(势)函数。 Δr表示根r的次数。 然后对于每个整数 k1

vLr(k)|TσΛ(r v)|ΔrAΨBΨ(max1d<Δsup𝒚~SdHdΨ(𝒚~)1)k1

在哪里

AΨ=maxuLr(1){|h(logRu)|ψ(logRu)}andBΨ=maxvLr(k){ψ(logRv)}.

在证明14之前,我们首先提出对树的影响的两个有用的属性。 首先,在[ALO20]中表明,影响满足以下形式的树链式法则。

Lemma 15 ([ALO20,引理B.2])

假设 u,v,wVT 是三个不同的顶点,使得 u 位于从 vw 的唯一路径上。然后

TσΛ(v w)=TσΛ(v u)TσΛ(u w).

其次,对于树上的两个相邻顶点,一个对另一个的影响由函数h给出。

Lemma 16.

vVTu 为子树 Tvv 的子级。 然后

TσΛ(v u)=h(logRu).
证明。

引理可以通过影响力的显式计算来证明。 在这里,我们利用 12 提出了一个更微妙的证明,它给出了影响力与函数 h 之间关系的一些见解。我们假设v在子树Tv中有d个子节点,分别用u1=uu2,,ud表示。 我们还假设,作为比统一字段更通用的设置,每个顶点 w 都附加到其自己的字段 λw 上。 那么 12 和树递归意味着

TσΛ(v u) =TvσΛ(v u)=(λuλu)logRv
=(λuλu)Hd(logRu1,,logRud)
=i=1dlogRuiHd(logRu1,,logRud)(λuλu)logRui
=i=1dh(logRui)TuiσΛ(ui u)=h(logRu),

最后一个相等是因为 TuiσΛ​0>(ui u)=0 对于 uiu TuσΛ(u u)=1.

我们现在准备证明14

14 的证明。

对于顶点vVT,用dv表示其子节点的数量;请注意dr=Δr u1,,uΔr 为根 r 的子级。我们可以假设 r 的所有这些子元素都是自由的,因为如果 ui 是固定的,那么 TσΛ(r ui)=0 根据定义。 然后通过1516,我们得到

vLr(k)|TσΛ(r v)| =i=1Δr|TσΛ(r ui)|vLui(k1)|TσΛ(ui v)|
=i=1Δr|h(logRui)|vLui(k1)|TσΛ(ui v)|
=i=1Δr|h(logRui)|ψ(logRui)vLui(k1)ψ(logRui)|TσΛ(ui v)|.

因此,我们得到

vLr(k)|TσΛ(r v)|Δrmax1iΔr{|h(logRui)|ψ(logRui)}max1iΔr{vLui(k1)ψ(logRui)|TσΛ(ui v)|}. (6)

接下来,我们通过归纳法证明,对于每个顶点 uVT\{r} 和每个整数 k0 我们有

vLu(k)ψ(logRu)|TσΛ(u v)|maxvLu(k){ψ(logRv)}(maxwVTusup𝒚~SdwHdwΨ(𝒚~)1)k. (7)

观察一下,一旦我们建立 Eq. 7,引理就会立即通过插入 Eq. 7 代入方程 6 我们将使用k上的归纳法来证明方程7 k=0时,如果uΛ固定,则Lu(0)=没有任何内容可显示;否则,Eq. 7 变为

ψ(logRu)|TσΛ(u u)|ψ(logRu),

自从 TσΛ(u u)=1. 现在假设 Eq. 7 对于某个整数 k10 成立(对于每个顶点 uVT\{r} )。 uVT\{r} 为任意值,并用 w1,,wd 表示 u 的子级,其中 1d<Δ (如果 d=0Lu(k)=Eq. 7 基本成立)。 再次到 1516 我们有

vLu(k)ψ(logRu)|TσΛ(u v)| =i=1dψ(logRu)|TσΛ(u wi)|vLwi(k1)|TσΛ(wi v)|
=i=1dψ(logRu)ψ(logRwi)|h(logRwi)|vLwi(k1)ψ(logRwi)|TσΛ(wi v)|.

使用归纳假设,我们得到

vLu(k)ψ(logRu)|TσΛ(u v)|
i=1dψ(logRu)ψ(logRwi)|h(logRwi)|maxvLwi(k1){ψ(logRv)}(maxwVTwisup𝒚~SdwHdwΨ(𝒚~)1)k1
maxvLu(k){ψ(logRv)}(maxwVTu\{u}sup𝒚~SdwHdwΨ(𝒚~)1)k1i=1dψ(logRu)ψ(logRwi)|h(logRwi)|
maxvLu(k){ψ(logRv)}(maxwVTusup𝒚~SdwHdwΨ(𝒚~)1)k,

最后的不等式由此得出

i=1dψ(logRu)ψ(logRwi)|h(logRwi)| =i=1d|Ψ(logRwi)HdΨ(Ψ(logRw1),,Ψ(logRwd))|
=HdΨ(Ψ(logRw1),,Ψ(logRwd))1.

至此方程 7成立,从而完成了引理的证明。

然后我们推导出 9 作为推论。

9 的证明。

由于 Ψ 是一个 (α,c) 势,因此 Contraction 条件意味着

max1d<Δsup𝒚~SdHdΨ(𝒚~)11α.

同时,由于子树Tv中的顶点vVT\{r}的度小于Δ,所以我们有logRvJ 那么 Boundedness 条件意味着对于所有 uLr(1)vLr(k)

ψ(logRv)ψ(logRu)|h(logRu)|cΔ.

因此,我们得到

ΔrAΨBΨ=ΔrmaxuLr(1){|h(logRu)|ψ(logRu)}maxvLr(k){ψ(logRv)}c.

引理紧随 14 而来。

6 验证良好的潜力:收缩

在本节中,我们第一步来证明10 Δ3 为整数。 β,γ,λ为实数,使得0βγγ>0βγ<1λ>0 回想一下,我们通过导数定义了我们的势函数 Ψ:[,+](,+)

Ψ(y)=ψ(y)=(1βγ)ey(βey+1)(ey+γ),Ψ(0)=0. (1)

我们在附录C中包含一个简短的证明,以表明Ψ是明确定义的。 If (β,γ,λ) is up-to-Δ unique with gap δ(0,1), then we show that Ψ satisfies the Contraction condition for α=δ/2. 这适用于唯一性区域中的所有参数(β,γ,λ),无需要求γ1 在后面的 附录 E中,我们建立了 Ψγ1 时的 有界性 条件,完成了 10 的证明。 γ>1的情况比较复杂,留给Section7

在给出证明之前,我们首先指出势函数 Ψ 本质上与 [LLY13] 中使用的势函数 Φ 相同(请注意,[LLY13] 使用 φ 作为势函数的符号,使用 Φ=φ 作为其导数)。 回想一下,边际比率的树递归由函数 Fd:[0,+]d[0,+] 给出,其中 1dΔ 使得对于所有 (x1,,xd)[0,+]d

Fd(x1,,xd)=λi=1dβxi+1xi+γ.

来自 [LLY13] 的势函数 Φ:[0,+](,+) 通过其导数隐式定义为

Φ(x)=φ(x)=1x(βx+1)(x+γ),Φ(1)=0.

以下引理解释了我们如何从 Φ 获得潜在的 Ψ

Lemma 17.

我们有Ψ=1βγ(Φexp);即,Ψ(y)=1βγΦ(ey) 代表所有 y[,+]

证明。

检查起来很简单

ψ(y)=(1βγ)ey(βey+1)(ey+γ)=1βγey1ey(βey+1)(ey+γ)=1βγeyφ(ey).

所以,

Ψ(y)=0yψ(t)dt=1βγ0yetφ(et)dt=1βγ1eyφ(s)ds=1βγΦ(ey).

结合[LLY13]引理12、13和14的结果,我们得到当(β,γ,λ)在时,势函数Φ满足以下梯度界限独特性区域。 请注意,这可以被视为收缩条件,但对于ΦFd而言。

Theorem 18 ([LLY13]).

SΦ=Φ[0,+]Φ 的图像。 如果参数 (β,γ,λ) 最多为 Δ 唯一,且间隔 δ(0,1),则对于每个整数 d 使得 1d<Δ 和每个 (x~1,,x~d)SΦd

FdΦ(x~1,,x~d)11δ

其中FdΦ=ΦFdΦ1

回想一下我们在部分1.1中的定义。 就对数边际比而言,树递归由函数 Hd:[,+]d[,+] 描述,其中 1dΔ 使得对于每个 (y1,,yd)[,+]d

Hd(y1,,yd)=logλ+i=1dlog(βeyi+1eyi+γ).

观察 Hd=logFdexp,因为我们从比率转向对数比率。 我们现在准备好为 Ψ 建立收缩条件。

Lemma 19.

SΨ=Ψ[,+]Ψ 的图像。 如果参数 (β,γ,λ) 最多为 Δ 唯一,且间隔 δ(0,1),则对于每个整数 d 使得 1d<Δ 和每个 (y~1,,y~d)SΨd

HdΨ(y~1,,y~d)11δ

其中HdΨ=ΨHdΨ1

证明。

x 的线性函数 a: 定义为 a(x)=1βγx 然后17给出Ψ=aΦexp,从而给出Ψlog=aΦ 由此可见,对于每个 1d<Δ

HdΨ=ΨHdΨ1=ΨlogFdexpΨ1=aΦFdΦ1a1=aFdΦa1.

这意味着,对于每个 (y~1,,y~d)SΨd 我们有

HdΨ(y~1,,y~d)=1βγFdΦ(x~1,,x~d)

其中 x~i=y~i/1βγ 代表 1id 然后,对于每个i

y~iHdΨ(y~1,,y~d)=1βγx~iFdΦ(x~1,,x~d)dx~idy~i=x~iFdΦ(x~1,,x~d).

这意味着 HdΨ(y~1,,y~d)=FdΦ(x~1,,x~d) 对于所有 (y~1,,y~d)SΨd 而言,引理则来自 18

7 剩余的反铁磁情况:βγΔ2Δγ>1

在本节中,我们讨论 βγΔ2Δγ>1 的情况。 正如[LLY13]中所研究的,在这种情况下,唯一性区域更加复杂。 例如,存在一个临界λc>0,使得具有λ<λc2自旋系统处于任意图的唯一性区域中;也就是说,(β,γ,λ) 最多是 是唯一的。 为了处理大度数,我们需要放宽4中的有界性条件,并定义更通用版本的(α,c)势。 我们将看到 5 对于这个一般的 (α,c) 潜力仍然成立。 其背后的原因是,为了限制影响力矩阵的最大特征值,只需考虑具有大度数的顶点的绝对影响力的顶点加权和。

Remark 20

我们在部分E.1中提供了有关唯一性区域的更多背景信息。 请注意,在最近的 [LLY13] 修订版中,作者更新了案例 βγΔ2Δγ>1 的唯一性区域的描述,修复了一个小错误在以前的版本中。 本文本节及附录E中的陈述和证明也根据新版本[LLY13]<进行了相应调整/t3>.

回想一下,我们的目标是限制矩阵 GσΛ 的最大特征值。 我们可以通过限制绝对行总和的上限来做到这一点 ΣvεV\Λ |GσΛ(r v)| 固定的 r,从而为我们提供了 λmax(GσΛ) 的有效上限。 然而,当βγΔ2Δγ>1时,这种方法不起作用。 在这种情况下,势 Ψ 无法成为独立于 Δ 的通用常数 c(α,c) 势。 事实上,不存在这样的 (α,c) 势作为绝对行和 ΣvεV\Λ |GσΛ(r v)| 可以与Θ(Δ)一样大。 特别是,如果参数 (β,γ,λ) 最多为 唯一,这意味着自旋系统对于任意图都具有唯一性,则绝对行和 ΣvεV\Λ |GσΛ(r v)| 可以是 Θ(n),其中 n=|V| 我们举一个具体的例子来说明这种情况。

Example 21.

考虑由参数 β=0γ>1λ>0 在以 r 为中心、有 Δ 叶的星形图上指定的反铁磁性 2 自旋系统。 一个简单的计算表明 |G(r v)|=λλ+γ 对于任何叶顶点vr 因此, Σvr|G(r v)|=Δλλ+γ. 现在,从 γ>1 开始,我们有

λc=λc(γ,Δ)=min1<d<Δγd+1dd(d1)d+1=Θγ(1),

强迫 Σvr|G(r v)|=θγ(Δ) 即使 λ<λc 位于唯一性区域。 然而,我们仍然有 λmax(G)=O(1) 因为 Σvr|G(v r)|=O(1) .

为了解决这个问题,人们可能需要考虑绝对列总和,涉及对固定顶点的绝对影响的总和。 然而,这不允许我们使用图和 SAW 树之间的美丽连接,如 8 所示。 相反,我们在这里考虑 GσΛ 绝对行和的顶点加权版本,它也是最大特征值的上限。

Lemma 22.

ρ:V+ 为顶点的正权重函数。 如果存在一个常量 ξ>0,对于每个 rV 我们有

vV\Λρv|GσΛ(r v)|ξρr, (8)

然后λmax(GσΛ)ξ

证明。

𝒫=diag{ρv:vV\Λ} 那么假设就等价于𝒫1GσΛ𝒫ξ 由此可见λmax(GσΛ)=λmax(𝒫1GσΛ𝒫)ξ

然后,我们修改 4(α,c)-势的定义,以允许较弱的 有界性 条件。 我们注意到 234 之间唯一的两个区别是:我们允许 Δ=;并且有界性条件被放宽到我们所说的一般有界性 回想一下,对于每个 0d<Δ,我们在 βγ<1 时设置 Jd=[log(λβd),log(λ/γd)],在 βγ>1 时设置 Jd=[log(λ/γd),log(λβd)]

Definition 23 (一般(α,c)-势函数)

Δ3 为整数或 Δ= β,γ,λ为实数,使得0βγγ>0λ>0 Ψ:[,+](,+) 为具有图像 S=Ψ[,+] 和导数 ψ=Ψ 的可微且递增函数。 对于任何α(0,1)c>0,我们说Ψ是一个一般(α,c)势函数 Δ(β,γ,λ)如果满足以下条件:

  1. 1.

    (收缩)对于每个整数 d 使得 1d<Δ 和每个 (y~1,,y~d)Sd,我​​们有

    HdΨ(y~1,,y~d)1=i=1dψ(y)ψ(yi)|h(yi)|1α

    其中 HdΨ=ΨHdΨ1yi=Ψ1(y~i)1id,以及y=Hd(y1,,yd)

  2. 2.

    (一般有界性)对于所有整数 d1,d2 使得 0d1,d2<Δ 和所有实数 y1Jd1,y2Jd2,我​​们有

    ψ(y2)ψ(y1)|h(y1)|2cd1+d2+2.

请注意,一般有界性 是比有界性 更弱的条件。 要看到这一点,如果潜在函数 Ψ 满足带有参数 c有界性,则对于每个 0di<Δ 和每个 yiJdi 其中 i=1,2 我们有

ψ(y2)ψ(y1)|h(y1)|cΔ2cd1+d2+2.

以下定理概括了 5 并表明一般的 (α,c) 势函数足以建立格劳伯动力学的快速混合。

Theorem 24.

Δ3 为整数或 Δ=+ β,γ,λ为实数,使得0βγγ>0λ>0 假设对于某些 α(0,1)c>0 存在相对于 Δ(β,γ,λ) 的一般 (α,c) 势。 然后对于最大度Δ的每个n-顶点图G,2-自旋系统的格劳伯动力学的混合时间带有参数 (β,γ,λ)GO(n2+2c/α)

然后,我们给出了 10 的对应关系,表明当 βγΔ2Δγ>1 时,Ψ 是一个一般的 (α,c) 电位。 然后从2425获得本例的3

Lemma 25.

Δ3 为整数。 β,γ,λ为实数,使得0β<1<γβγΔ2Δ 假设 (β,γ,λ) 最多 Δ 是唯一的,间隙为 δ(0,1) 那么由Eq.4隐式定义的函数Ψ就是一个通用的(α,c)-势函数与αδ/2c18;如果β=0,我们可以进一步取c4

24的证明可以在附录D中找到。对于25,Ψ收缩条件来自19,并且一般有界性1>是与所有其他案例一起在附录3>E4>2>中得到证明。

8 铁磁案例

在铁磁情况下,[GL18, SS20] 中给出了最著名的相关衰减结果。 使用[GL18][SS20]中的势函数,我们显示以下两个结果,与已知的相关衰减结果相匹配。 事实上,[SS20] 的势函数是常量 α=Θ(δ)cO(1)(α,c) 势函数。

Theorem 26.

固定一个整数Δ3、正实数β,γ,λ0<δ<1,并假设(β,γ,λ)满足以下三个条件之一:

  1. 1.

    Δ2+δΔδβγΔδΔ2+δ,λ是任意的;

  2. 2.

    βγΔΔ2λ(1δ)γmax{1,βΔ1}((Δ2)βγΔ);

  3. 3.

    βγΔΔ2λ11δ(Δ2)βγΔβmin{1,1/γΔ1}

那么恒等函数Ψ(y)=y(基于[SS20]中给出的势)是(α,c)-α=Θ(δ)的势函数,并且cO(1) 此外,对于每个最大度数至多Δn-顶点图G,G 与参数 (β,γ,λ)O(n2+c/δ),对于通用常量 c>0

Remark 3

条件1包括铁磁情况1<βγΔδΔ2+δ和反铁磁情况Δ2+δΔδβγ<1 请注意,在这两种情况下,(β,γ,λ) 最多是 Δ 唯一,且有间隙 δ 对于反铁磁性情况,与 公式 410 给出的电势 Ψ 的约束 αδ/2 相比,同一函数 Ψ 是一个 (α,c) 电位,具有 c1.5 和更好的收缩率 αδ 对于具有β=γ>1(伊辛模型)的铁磁情况,已知[MS13]的更强结果,这给出了O(nlogn)混合。

[GL18] 的势函数确实是一个 (α,c) 势函数,但不幸的是,c 必须依赖于 Δ 我们得到以下结果,该结果弱于无界度图[GL18]中的相关衰减算法。

Theorem 27.

固定一个整数Δ3和满足β1γβγΔΔ2λ<(γβ)βγβγ1的非负实数β,γ,λ 然后对于每个n-顶点图G,其最大度数至多为Δ,铁磁2-自旋系统的格劳伯动力学在G 与参数 (β,γ,λ)O(nC),对于常量 C 仅取决于 β,γ,λ,Δ,但不取决于 n.

这些定理的证明在附录F中提供。

参考

  • [AL20] Vedat Levi Alev and Lap Chi Lau “Improved analysis of higher order random walks and applications” In Proceedings of the 52nd Annual ACM SIGACT Symposium on Theory of Computing (STOC), 2020, pp. 1198–1211
  • [ALO20] Nima Anari, Kuikui Liu and Shayan Oveis Gharan “Spectral Independence in High-Dimensional Expanders and Applications to the Hardcore Model” In arXiv preprint arXiv:2001.00303, 2020
  • [Bar16] Alexander Barvinok “Combinatorics and Complexity of Partition Functions” Springer AlgorithmsCombinatorics, 2016
  • [Ben18] Ferenc Bencs “On trees with real-rooted independence polynomial” In Discrete Mathematics 341.12 Elsevier, 2018, pp. 3321–3330
  • [GL18] Heng Guo and Pinyan Lu “Uniqueness, Spatial Mixing, and Approximation for Ferromagnetic 2-Spin Systems” In ACM Transactions on Computation Theory 10.4, 2018
  • [GŠV16] Andreas Galanis, Daniel Štefankovič and Eric Vigoda “Inapproximability of the Partition Function for the Antiferromagnetic Ising and Hard-Core Models” In Combinatorics, Probability and Computing 25.4 Cambridge University Press, 2016, pp. 500–559
  • [JS93] Mark Jerrum and Alistair Sinclair “Polynomial-time approximation algorithms for the Ising model” In SIAM Journal on Computing 22.5 SIAM, 1993, pp. 1087–1116
  • [JVV86] Mark R Jerrum, Leslie G Valiant and Vijay V Vazirani “Random generation of combinatorial structures from a uniform distribution” In Theoretical Computer Science 43 Elsevier, 1986, pp. 169–188
  • [Kel85] F.. Kelly “Stochastic Models of Computer Communication Systems” In Journal of the Royal Statistical Society. Series B (Methodological) 47.3, 1985, pp. 379–395
  • [LLY12] Liang Li, Pinyan Lu and Yitong Yin “Approximate Counting via Correlation Decay in Spin Systems” In Proceedings of the 23rd Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), 2012, pp. 922–940
  • [LLY13] Liang Li, Pinyan Lu and Yitong Yin “Correlation Decay Up to Uniqueness in Spin Systems” In Proceedings of the 24th Annual ACM-SIAM Symposium on Discrete Algorithms (SODA), 2013, pp. 67–84
  • [LSS19] Jingcheng Liu, Alistair Sinclair and Piyush Srivastava “Fisher zeros and correlation decay in the Ising model” In Journal of Mathematical Physics 60.10 AIP Publishing LLC, 2019, pp. 103304
  • [MS13] Elchanan Mossel and Allan Sly “Exact thresholds for Ising-Gibbs samplers on general graphs” In Annals of Probability 41.1, 2013, pp. 294–328
  • [PR17] Viresh Patel and Guus Regts “Deterministic Polynomial-Time Approximation Algorithms for Partition Functions and Graph Polynomials” In SIAM Journal on Computing 46, 2017, pp. 1893–1919
  • [PR19] Han Peters and Guus Regts “On a conjecture of Sokal concerning roots of the independence polynomial” In The Michigan Mathematical Journal 68.1, 2019, pp. 33–55
  • [Sly10] Allan Sly “Computational Transition at the Uniqueness Threshold” In Proceedings of the 51st Annual IEEE Symposium on Foundations of Computer Science (FOCS), 2010, pp. 287–296
  • [SS05] Alexander D. Scott and Alan D. Sokal “The Repulsive Lattice Gas, the Independent-Set Polynomial, and the Lovász Local Lemma” In Journal of Statistical Physics 118.5, 2005, pp. 1151–1261
  • [SS14] Allan Sly and Nike Sun “The Computational Hardness of Counting in Two-Spin Models on d-Regular Graphs” In The Annals of Probability 42.6, 2014, pp. 2383–2416
  • [SS20] Shuai Shao and Yuxin Sun “Contraction: A Unified Perspective of Correlation Decay and Zero-Freeness of 2-Spin Systems” In Proceedings of the 47th International Colloquium on Automata, Languages, and Programming (ICALP), 2020, pp. 96:1–15
  • [SST14] Alistair Sinclair, Piyush Srivastava and Marc Thurley “Approximation Algorithms for Two-State Anti-Ferromagnetic Spin Systems on Bounded Degree Graphs” In Journal of Statistical Physics 155.4, 2014, pp. 666–686
  • [ŠVV09] Daniel Štefankovič, Santosh Vempala and Eric Vigoda “Adaptive simulated annealing: A near-optimal connection between sampling and counting” In Journal of the ACM 56.3, 2009, pp. 1–36
  • [Wei06] Dror Weitz “Counting Independent Sets Up to the Tree Threshold” In Proceedings of the 38th Annual ACM Symposium on Theory of Computing (STOC), 2006, pp. 140–149

附录A主要结果证明

本节我们给出1235的证明。

5 的证明。

请注意,由于格劳伯动力学的转移矩阵 P 具有所有非负特征值,因此我们有 λ(P)=1λ2(P) ,因此为了推导出混合,它足以满足下界 1λ2(P) 我们通过使用7来做到这一点。 对于足够小的ηi,足以显示(η0,,ηn2)的光谱独立性。

要绑定ηi,只需绑定 ΣvεV\{r}0>|GσΛ(r v)| 对于具有 n=|V| 顶点的所有图 G=(V,E) 以及 i 顶点的子集 Λ 上的所有边界条件 σΛ 我们声明如下:

vV\{r}|GσΛ(r v)|min{cα,C(ni1)} (9)

其中 C(0,1) 是仅取决于 β,γ,λ,Δ 的常数。 第一个上限 cδ 推导如下

vV\{r}|GσΛ(r v)| vVT\{r}|TσΛ(r v)| (8; T=Tsaw(G,r))
=k=1vLr(k)|TσΛΛ(r v)| (split the sum by levels)
ck=1(1α)k1 (9)
=cα.

第二个上限 C(ni1) 更简单。 直观上,这意味着每个绝对成对影响 |GσΛ(r v)| 至多是某个常数C,因此绝对影响力的总和上限为C(ni1) 下面的两个主张给出了更精确的陈述,其证明在部分A.2中提供。

Claim 28 (反铁磁案例)

固定整数Δ3和实数β,γ,λ,并假设0βγγ>0βγ<1λ>0 那么,对于每个最大阶数至多为 Δn 顶点图 G,参数为 (β,γ,λ)G 上的反铁磁 2 自旋系统是 Cn 光谱独立的,其常数 0<C<1 仅取决于 β,γ,λ,Δ 此外,如果 (β,γ,Δ) 最多是 Δ 唯一,那么我们可以放弃对 Δ 的依赖。

Claim 29 (铁磁外壳)

固定一个整数Δ3和正实数β,γ,λ,并假设βγβγ>1 那么,对于每个最大阶数至多为 Δn 顶点图 G,参数为 (β,γ,λ)G 上的铁磁 2 自旋系统是 Cn 光谱独立的,常数 0<C<1 仅取决于 β,γ,λ,Δ

有了Eq.9,我们立即看到到7

1λ2(P)1ni=0n2(1ηini1)1n(1C)2c/α1i=0n2c/α1(1cα1ni1).

利用 1xexp(xx2) 对于所有 0x12 的事实(可以通过微积分直接证明),我们得到

i=0n2c/α1(1cα1ni1)=j=2c/αn1(1cα1j)exp(cαj=2c/αn11jc2α2j=2c/αn11j2).

现在自从

j=2c/αn11jj=2n1j1ndxx=logn

j=2c/αn11j2j=21j(j1)=1,

我们推断

1λ2(P)(1C)2c/α1e(c/α)2n(1+c/α).

然后从 Eq. 1 得出定理。

3 的证明。

我们利用 524,只要存在 (α,c) 电位,就会出现 O(n2+cα) 混合;如果存在一般的 (α,c) 电位,就会出现 O(n2+2cα) 混合。 我们使用 Eq. 4 给出的势,它是 [LLY13] 中势函数的改编到对数边际比率。 (β,γ,λ) 最多为 Δ 唯一且间隙 δ(0,1) 时,它是一个 (α,c) 势或一般 (α,c)-由1025计算的势,其中αδ/2c 定理如下。

1 的证明。

30之后的部分A.1中,λ(1δ)λc(Δ)意味着最多Δ 具有间隙δ/4 的唯一性。 γ1开始,我们可以再次诉诸10来获得(α,c)αδ/8c4的势。 1 然后是 5O(n2+32/δ) 混合。

2 的证明。

31之后的部分A.1中,ββc(Δ)+δ(1βc(Δ))意味着最多Δ 具有间隙δ 的唯一性。 再次呼吁 10,我们获得 (α,c)αδ/2c1.5 的势。 2 然后是 5O(n2+3/δ) 混合。

虽然我们在技术上通过使用 [LLY13] 势能得到 O(n2+3/δ),但我们可以通过使用平凡的恒等函数作为势能将其改进为 O(n2+1.5/δ) 混合。 参见26的第一种情况(在 F.1中证明)和备注 3

A.1 参数 paps 方面的唯一性差距

在本节中,我们陈述并证明3031,它们将参数差距与唯一性差距联系起来。

Claim 30 (硬核模型;来自 [ALO20] 的引理 C.1)

固定整数 Δ30<δ<1β=0,γ>0 如果λ(1δ)λc(γ,Δ),则(β,γ,λ)最多Δ是唯一的,间隔δ/4

Claim 31 (大βγ)

固定整数Δ30<δ<1 如果 βγΔ2Δ+δ(1Δ2Δ)=Δ2(1δ)Δ,则 (β,γ,λ) 最多 Δ 是唯一的,对于所有 λ 来说,间隔为 0<δ<1 注意如果β=γ,这正是条件ββc(Δ)+δ(1βc(Δ))

证明。

考虑 d<Δfd(R)=λ(βR+1R+γ)d 边际比率的单变量递归。求微分,我们有

fd(R) =dλ(βR+1R+γ)d1(βR+γβR+1(R+γ)2)=d(1βγ)λ(βR+1R+γ)d1(βR+1)(R+γ)
=d(1βγ)fd(R)(βR+1)(R+γ).

在唯一的固定点Rd,我们有fd(Rd)=Rd,所以

|fd(Rd)|=d(1βγ)Rd(βRd+1)(Rd+γ).

37 时,我们得到了上限

|fd(Rd)|d1βγ(1+βγ)2=d1βγ1+βγ.

由于我们假设 βγΔ2(1δ)Δ,我们得到

d1βγ1+βγdΔ(Δ2(1δ))Δ+(Δ2(1δ))=d1δΔ1+δ(1δ)dΔ1.

由于所有 d<Δ 最多为 1δ,因此我们最多有 Δ 个唯一性,间隙为 δ

A.2 恒定大小图的谱独立界限

在本节中,我们证明具有少于 O(c/α) 多个顶点的图的谱独立边界,因为对于具有如此少顶点的图,基于树递归收缩的边界变得微不足道。

28 的证明。

如果Rv表示顶点vG的边际比率,则RvλβΔ γ1的情况下,我们也有Rvλ/γΔ;如果γ>1,我们就有Rvλ 由此可见,我们立即就有了界限

|G(u v)|{|λλ+γΔλβΔ1+λβΔ|=λ(1βΔγΔ)(λ+γΔ)(1+λβΔ),if γ1|λ1+λλβΔ1+λβΔ|=λ(1βΔ)(λ+1)(1+λβΔ),o.w.

对于所有u,vG 请注意,这些常量小于 1,并且仅取决于 β,γ,λ,Δ,从而产生第一个声明。

现在,当最多 Δ 唯一性成立时,我们继续消除对 Δ 的依赖。 我们有以下案例:

  1. 1.

    如果γ>1,我们立即获得λ1+λ的边界,它独立于Δ

  2. 2.

    如果β=0γ1,则λ(1βΔγΔ)(λ+γΔ)(1+λβΔ)=λλ+γΔλγΔ 由于 (β,γ,λ) 的唯一性最多为 Δ,因此我们必须有 λλc(γ,Δ)=min1<d<Δγd+1dd(d1)d+1γΔ(Δ1)Δ1(Δ2)ΔγΔO(1/Δ) 由此可见λγΔO(1/Δ)

  3. 3.

    如果βγ>Δ2Δγ1,那么

    λ(1βΔγΔ)(λ+γΔ)(1+λβΔ)1βΔγΔ1e2.
  4. 4.

    如果 βγΔ2Δ,那么让 Δ0 成为 1<d<Δ 中的最大值,使得 βγ>d2d. 如果λλc(β,γ,Δ),那么到35,我们有

    λ(1βΔγΔ)(λ+γΔ)(1+λβΔ)λγΔO(Δ0/Δ).

    如果λλ¯c(β,γ,Δ),那么再次到35,我们有

    λ(1βΔγΔ)(λ+γΔ)(1+λβΔ)1λβΔO(Δ0/Δ).
29 的证明。

证明与反铁磁情况相同,我们在此省略。

附录 B 12 的证明(第 1 部分和第 2 部分)

12 的证明(第 1 部分和第 2 部分)。

为了看到第一个等式,我们直接计算并得到

(λvλv)logZGσΛ =1ZGσΛ(λvλv)ZGσΛ
=1ZGσΛσ{0,1}V\Λ(λvλv)(βm1(σ)γm0(σ)wVλwσw)
=1ZGσΛσ{0,1}V\Λσv(βm1(σ)γm0(σ)wVλwσw)
=σ{0,1}V\ΛσvμG(σσΛ)=MGσΛ(v).

对于第 2 部分,使用上面的结果,我们还可以得到

(λvλv)(λuλu)logZGσΛ
= (λvλv)(1ZGσΛ(λuλu)ZGσΛ)
= 1ZGσΛ(λvλv)(λuλu)ZGσΛ1(ZGσΛ)2(λvλv)ZGσΛ(λuλu)ZGσΛ
= 1ZGσΛ(λvλv)(σ{0,1}V\Λσu(βm1(σ)γm0(σ)wVλwσw))MGσΛ(u)MGσΛ(v)
= 1ZGσΛσ{0,1}V\Λσu(λvλv)(βm1(σ)γm0(σ)wVλwσw)MGσΛ(u)MGσΛ(v)
= 1ZGσΛσ{0,1}V\Λσuσv(βm1(σ)γm0(σ)wVλwσw)MGσΛ(u)MGσΛ(v)
= σ{0,1}V\ΛσuσvμG(σσΛ)MGσΛ(u)MGσΛ(v)
= KGσΛ(u,v).

附录 C Ψ 的技术引理

以下引理意味着由 Eq. 4 给出的势 Ψ 是明确定义的。

Lemma 32.

对于所有 β,γ>0 这样的 βγ<1,我们有

+(1βγ)ey(βey+1)(ey+γ)<+.
证明。

对于 + 一侧,我们有

0+(1βγ)ey(βey+1)(ey+γ)=0+1βγβey+γey+βγ+1<0+1βey<+.

类似地,对于 一侧,我们有

0(1βγ)ey(βey+1)(ey+γ)<01γey<+.

附录D势法混合:24的证明

在本节中,我们以与 5 相同的方式证明 24,如 Section 3 中所述 这里的主要区别是我们考虑绝对影响的加权和 ΣvεV\Λ ρv|GσΛ(r v)| 其中 ρ:V+ 是权重函数。 这足以让我们限制影响矩阵的特征值,如 22 所示。 我们将选择顶点v的权重为ρv=Δv,度为v。给定一般的(α,c)势,以下引理为我们提供了对距离k的绝对影响的加权和的上限。 特别是,它概括了9

Lemma 33.

如果存在关于Δ(β,γ,λ)的一般(α,c)-潜在函数Ψ,其中α(0,1)c>0,然后对于每个 ΛVT\{r}σΛ{0,1}Λ 和所有整数 k1

vLr(k)Δv|TσΛ(r v)|2c(1α)k1Δr

其中Lr(k)表示距r距离为k的所有自由顶点的集合。

为了证明33,我们首先声明对于任何权重函数ρ14的以下推广。 34 的证明与 14 相同,此处省略。

Lemma 34.

Ψ:[,+](,+) 为具有图像 S=Ψ[,+] 和导数 ψ=Ψ 的可微且递增(势)函数。 Δr表示根r的次数。 然后对于每个整数 k1

vLr(k)ρv|TσΛ(r v)|ΔrAΨBΨρ(max1d<Δsup𝒚~SdHdΨ(𝒚~)1)k1

在哪里

AΨ=maxuLr(1){|h(logRu)|ψ(logRu)}andBΨρ=maxvLr(k){ρvψ(logRv)}.

然后我们证明3324

33 的证明。

Δv表示顶点vVT\{r}的度数,用dv=Δv1表示子树Tvv的度数>。 为所有 vVT 选择顶点权重为 ρv=Δv 由于 Ψ 是一般的 (α,c) 势,因此 收缩 条件意味着

max1d<Δsup𝒚~SdHdΨ(𝒚~)11α.

由于 logRvJdv 根据 Jd 的定义,一般有界性条件意味着对于所有 uLr(1)vLr(k) ,

ψ(logRv)ψ(logRu)|h(logRu)|2cΔu+Δv.

因此,我们得到

ΔrAΨBΨρ=ΔrmaxuLr(1){|h(logRu)|ψ(logRu)}maxvLr(k){Δvψ(logRv)}2cΔr.

引理紧随 34 而来。

24 的证明。

24 的证明几乎与 5 相同。 我们指出,这里唯一的区别是我们考虑给定顶点的绝对影响的加权和。 由于 SAW 树保留了顶点度数,因此我们仍然可以应用 8 然后,结合722833,我​​们完成了定理的证明。

附录E验证良好的潜力:有界性

在本小节中,我们将展示由 Eq. 定义的势函数 Ψ有界性一般有界性 条件。 > 4 在不同的参数范围内。 结合19,我们完成了1025的证明。

E.1中,我们以[LLY13]的工作为基础,介绍了参数(β,γ,λ)唯一性区域的背景。 然后,我们在部分E.2中显示有界性一般有界性 技术引理的证明留给部分E.3

E.1 唯一性区域的预备知识

在本节中,我们简要描述参数(β,γ,λ)的唯一性区域。 这里的所有结果及其证明可以在最新版本的[LLY13]的引理21中找到。

Δ3 为整数,β,γ,λ 为实数。 我们假设 0βγγ>0βγ<1λ>0 对于 1dΔ 定义

fd(R)=λ(βR+1R+γ)d

并用Rd表示fd的唯一不动点。 回想一下,如果 |fd(Rd)|<1δ 对于所有 1d<Δ 都是唯一的,那么参数 (β,γ,λ) 直至 Δ 都是唯一的,间隙为 δ(0,1)

β=0时,自旋系统称为硬约束模型 在这种情况下,存在一个外部场的临界阈值,定义为

λc=λc(γ,Δ)=min1<d<Δγd+1dd(d1)d+1,

当且仅当 λ<λc 时,参数 (0,γ,λ) 最多为 Δ 唯一。 特别是,当 γ1 临界字段由下式给出时

λc=λc(γ,Δ)=γΔ(Δ1)Δ1(Δ2)Δ.

β>0时,自旋系统称为软约束模型 如果 βγ>Δ2Δ,则 (β,γ,λ) 对于所有 λ>0 来说最多是唯一的 Δ 如果βγΔ2Δ唯一性区域更加复杂,我们现在描述。

Δ¯=1+βγ1βγ,

这样对于每个 1d<Δ¯ 我们有 d1βγ1+βγ<1,对于每个 dΔ¯ 我们有 d1βγ1+βγ1 对于每个 Δ¯d<Δ,我们将 x1(d)x2(d) 定义为二次方程的两个正根

d(1βγ)x(βx+1)(x+γ)=1.

更具体地说,x1(d)x2(d) 由下式给出

x1(d)=θ(d)θ(d)24βγ2βandx2(d)=θ(d)+θ(d)24βγ2β

在哪里

θ(d)=d(1βγ)(1+βγ).

请注意,θ(d)2βγ 适用于所有 dΔ¯ 对于 i=1,2 我们让

λi(d)=xi(d)(xi(d)+γβxi(d)+1)d.

然后,当且仅当 λ 属于以下状态时,参数 (β,γ,λ) 最多是唯一的Δ

𝒜=Δ¯d<Δ[(0,λ1(d))(λ2(d),)]. (10)

特别是,当 γ1 存在两个关键阈值 0<λc<λ¯c 时,参数 (β,γ,λ) 最多可达 Δ 唯一,当且仅当如果 λ<λcλ>λ¯c(即 𝒜=(0,λc)(λ¯c,)),其中

λc=λc(β,γ,Δ)=minΔ¯d<Δλ1(d)andλ¯c=λ¯c(β,γ,Δ)=maxΔ¯d<Δλ2(d)=λ2(Δ1).

以下关键字段的界限对我们稍后的证明很有帮助。

Lemma 35.
  1. 1.

    如果 β=0,则对于每个整数 d 使得 1<d<Δ 我们有

    λc4γd+1d1.
  2. 2.

    如果 β>0βγΔ2Δ,则对于每个整数 d这样Δ¯d<Δ我们有

    λ1(d)18γd+1θ(d)andλ2(d)θ(d)18βd+1

    其中θ(d)=d(1βγ)(1+βγ)

35的证明被推迟到SectionE.3

E.2 有界性证明

本节我们通过在相应范围内建立有界性一般有界性来完成1025的证明参数。

Δ3 为整数。 β,γ,λ为实数,使得0βγγ>0βγ<1λ>0 回想一下,势函数 Ψ 定义为

Ψ(y)=ψ(y)=(1βγ)ey(βey+1)(ey+γ)=|h(y)|,Ψ(0)=0. (1)

令人惊讶的是发现ψ=|h|,因为潜在的Ψ正是[LLY13]中的那个,如17所示>。 这似乎不是巧合,它提供了一些直觉,说明为什么 [LLY13] 的潜力有效。 更重要的是,ψ=|h|这一事实有助于我们证明有界性一般有界性 回想一下,对于 0d<Δβγ<1,我们让 Jd=[log(λβd),log(λ/γd)] 为具有 d 子节点的顶点的对数边际比范围。 然后对于每个 0di<ΔyiJdi 其中 i=1,2,我们有

ψ(y2)ψ(y1)|h(y1)|=|h(y1)||h(y2)|. (11)

下面的引理给出了 |h(y1)||h(y2)| 的上限,从中和 Eq. 11 我们推导出 有界性一般有界性 立即。 引理中的括号指示边界应用于哪个引理。

Lemma 36.

Δ3 为整数。 β,γ,λ为实数,使得0βγγ>0βγ<1λ>0 假设参数 (β,γ,λ) 最多 Δ 是唯一的,间隙为 δ(0,1) 那么对于所有整数 d1,d2 使得 0d1,d2<Δ 以及所有实数 yiJdi 其中 i=1,2,以下成立:

  1. H。

    硬约束模型: β=0λ<λc

    1. H.1。

      (10) 如果 γ1,则

      |h(y1)|4Δ.
    2. H.2。

      (25) 如果 γ>1,则

      |h(y1)||h(y2)|8d1+d2+2.
  2. S。

    软约束模型: β>0λ𝒜

    1. S.1。

      (10) 如果 βγ>Δ2Δ,则

      |h(y1)|1.5Δ.
    2. S.2。

      (10) 如果 βγΔ2Δγ1 ,那么

      |h(y1)|18Δ.
    3. S.3。

      (25) 如果 βγΔ2Δγ>1 ,那么

      |h(y1)||h(y2)|36d1+d2+2.

下面的引理很有帮助,其证明可以在 Section E.3 中找到。

Lemma 37.

对于每个 y[,+] 我们有

|h(y)|=|1βγ|ey(βey+1)(ey+γ)|1βγ|1+βγ.

我们在这里给出36的证明。

36 的证明。

我们使用部分E.1中的符号和结果。

H. 硬约束模型:β=0λ<λc

H.1。 γ1

对于每个 y1Jd1,我们从 35 推断出

ey1λγd1λcγΔ14γΔ2.

因此,

|h(y1)|=ey1ey1+γ4γΔ24γΔ2+γ=4Δ+24Δ.

H.2。 γ>1

y¯=y1+y22d¯=d1+d22 然后我们得到

|h(y1)||h(y2)|=ey1ey1+γey2ey2+γ=1(1+γey1)(1+γey2)11+γey¯,

其中最后一个不等式由 AM-GM 不等式得出

(1+γey1)(1+γey2)=1+γ(ey1+ey2)+γ2e2y¯1+2γey¯+γ2e2y¯=(1+γey¯)2.

由于 yiJdi 对于 i=1,2,我们有

ey¯=ey1ey2λγd1λγd2=λγd¯.

如果d¯2,那么我们从35γ>1推断出

ey¯λcγd¯4γd¯1.

它遵循

|h(y1)||h(y2)|11+γey¯11+d¯14=4d¯+38d1+d2+2.

如果d¯<2,那么很容易看出

|h(y1)||h(y2)|18d1+d2+2.

S. 软约束模型:β>0λ𝒜

S.1。 βγ>Δ2Δ

对于每个 y1J,我们从 37 推断出

|h(y1)|1βγ1+βγ1Δ11.5Δ.

S.2。 βγΔ2Δγ1

在本例中,我们有 λ<λcλ>λ¯c,其中 λc,λ¯c 是两个关键字段。 首先考虑λ>λ¯c 对于每个 y1Jd1,我们从 35β<1 推断出

ey1λβd1λ¯cβΔ1θ(Δ1)18β

其中θ(d)=d(1βγ)(1+βγ) 因此,

|h(y1)|=(1βγ)ey1(βey1+1)(ey1+γ) =1βγβey1+γey1+(1+βγ)
1βγθ(Δ1)18+(1+βγ)=18(1βγ)(Δ1)(1βγ)+17(1+βγ)18Δ.

接下来我们考虑λ<λc 对于每个 y1Jd1,我们从 35γ1 推断出

ey1λγd1λcγΔ118γθ(Δ1).

因此,

|h(y1)|=1βγβey1+γey1+(1+βγ)1βγθ(Δ1)18+(1+βγ)18Δ.

S.3。 βγΔ2Δγ>1

y¯=y1+y22d¯=d1+d22dL=d¯dR=d¯ 我们首先考虑一些简单的情况。 如果d¯2那么很容易看出

|h(y1)||h(y2)|16d1+d2+2.

如果 d¯>2dLΔ¯,那么我们从 37 推断出

|h(y1)||h(y2)|1βγ1+βγ=1Δ¯2d1+d226d1+d2+2.

因此,下面我们可以假设d¯>2dL>Δ¯

由于参数 (β,γ,λ) 直至 Δ 都是唯一的,因此我们有 λ𝒜 其中的制度 𝒜Eq. 10 给出。 观察一下

𝒜(0,λ1(dL))(λ2(dR),)(λ2(dL),λ1(dR))

其中最后一个间隔仅在 λ2(dL)<λ1(dR) 时才为非空。 这意味着 λ 至少包含在三个间隔之一中。 我们通过分别考虑这三种情况来建立界限。

情况 1:λ<λ1(dL) 根据柯西-施瓦茨不等式,我们有

|h(y1)||h(y2)| =1βγβey1+γey1+(1+βγ)1βγβey2+γey2+(1+βγ)
1βγ(βey1+γey1)(βey2+γey2)+(1+βγ). (12)

因此,我们得到

|h(y1)||h(y2)|1βγγey¯+(1+βγ).

由于yiJdi对于i=1,2γ>1,我们从35推断出

ey¯λγd¯λ1(dL)γdL18γθ(dL),

其中θ(dL)=dL(1βγ)(1+βγ) 它遵循

|h(y1)||h(y2)|1βγγey¯+(1+βγ)1βγθ(dL)18+(1+βγ)36d1+d2+2.

情况 2:λ>λ2(dR) 类似地,我们从 Eq. 12 得到

|h(y1)||h(y2)|1βγβey¯+(1+βγ).

由于yiJdi对于i=1,2β<1,我们从35推断出

ey¯λβd¯λ2(dR)βdRθ(dR)18β,

其中θ(dR)=dR(1βγ)(1+βγ) 它遵循

|h(y1)||h(y2)|1βγβey¯+(1+βγ)1βγθ(dR)18+(1+βγ)36d1+d2+2.

情况 3:λ2(dL)<λ<λ1(dR) 我们可以假设d1d2 通过 Eq. 12,我们得到

|h(y1)||h(y2)|1βγβγey2y12+(1+βγ).

由于 yiJdi 对于 i=1,2β<1<γ,我们有

ey2y1βd2γd1βdLγdR.

同时,我们从 35 推断出

θ(dL)18βdL+1λ2(dL)<λ<λ1(dR)18γdR+1θ(dR),

这意味着

βγey2y12βdL+1γdR+1θ(dL)θ(dR)18θ(dL)18.

它遵循

|h(y1)||h(y2)|1βγβγey2y12+(1+βγ)1βγθ(dL)18+(1+βγ)36d1+d2+2.

E.3 技术引理的证明

35 的证明。

1. 对于每个 1<d<Δ 我们有

λcγd+1dd(d1)d+1=γd+1d1(dd1)d4γd+1d1,

对于所有整数 d>1,最后一个不等式是从 (dd1)d4 得出的。

2. 对于每个 Δ¯d<Δ 我们有

x1(d)=2γθ(d)+θ(d)24βγ2γθ(d).

观察函数 x+γβx+1βγ<1 时在 x 中单调递增,因此我们推断

x1(d)+γβx1(d)+12γθ(d)+γ2βγθ(d)+1=γ2+d(1βγ)(1+βγ)2βγ+d(1βγ)(1+βγ)=γd+1d1.

所以,

λ1(d)=x1(d)(x1(d)+γβx1(d)+1)d2γθ(d)γd(d+1d1)d18γd+1θ(d)

对于所有整数 d>1,最后一个不等式是从 (d+1d1)d9 得出的。

第二部分可类似证明。 对于每个 Δ¯d<Δ 我们有

x2(d)=θ(d)+θ(d)24βγ2βθ(d)2β,

因此,

x2(d)+γβx2(d)+1θ(d)2β+γθ(d)2+1=1βd(1βγ)(1+βγ)+2βγd(1βγ)(1+βγ)+2=1βd1d+1.

然后我们得出结论:

λ2(d)=x2(d)(x2(d)+γβx2(d)+1)dθ(d)2β1βd(d1d+1)dθ(d)18βd+1,

对于所有整数 d>1,最后一个不等式再次从 (d+1d1)d9 得出。

37 的证明。

我们从 AM-GM 不等式推断出

|h(y)|=|1βγ|βey+γey+1+β|1βγ|2βγ+1+β=|1βγ|1+βγ.

附录F铁磁情况的证明

F.1 26的证明

26 的证明。

自始至终,我们都使用“微不足道的势”函数Ψ(y)=y 请注意,ψ(y)=1 是一个常数函数。 现在,我们证明收缩有界性 我们分为三种情况。

  1. 1.

    我们首先证明Contraction部分。 37,对于所有y[,+],我们有

    |h(y)||1βγ|1+βγ1δΔ1.

    现在让我们证明有界性条件。 由上面的不等式我们有

    |h(y)|1Δ11.5Δ

    对于Δ3

  2. 2.

    对于 Contraction 部分,从 log(λmax{1,1/γΔ1})yilog(λmax{1,βΔ1}) 开始,我们有

    |Hd(𝒚)yi| =|h(yi)|=βγ11+βγ+γeyi+βeyiβγ11+βγ+γeyi
    βγ11+βγ+γλmax{1,βΔ1}.

    由于我们假设 λ(1δ)γmax{1,βΔ1}((Δ2)βγΔ),因此我们有上限

    βγ11+βγ+(Δ2)βγΔ1δ =(1δ)βγ1(Δ1δ)βγ(Δ1+δ)
    =(1δ)βγ1(Δ1δ)(βγ1)+2δ
    1δΔ1δ(1Θ(δ))1Δ1.

    现在,我们证明有界性条件。 请注意,从 λγmax{1,βΔ1}((Δ2)βγΔ) 开始,就得出 ylog(λmax{1,βΔ1})log(γ(Δ2)βγΔ) 简单的计算表明 γ(Δ2)βγΔγβ ,因此到 37 时,我们有

    |h(y)| |h(log(γ(Δ2)βγΔ))|(βγ1)elog(γ(Δ2)βγΔ)elog(γ(Δ2)βγΔ)+γ
    =(βγ1)11+(Δ2)βγΔ=βγ1(Δ2)(βγ1)1O(1/Δ).
  3. 3.

    对于 Contraction 部分,从 log(λmax{1,1/γΔ1})yilog(λmax{1,βΔ1}) 开始,我们有

    |Hd(𝒚)yi| =|h(yi)|=βγ11+βγ+γeyi+βeyiβγ11+βγ+βeyi
    βγ11+βγ+βλmax{1,1/γΔ1}.

    由于我们假设 λ11δ(Δ2)βγΔβmin{1,1/γΔ1},因此我们有上限

    βγ11+βγ+(Δ2)βγΔ1δ

    正如我们在上面的情况 2 中计算的那样,它的上限再次为 (1Θ(δ))1Δ1

    现在,我们证明有界性条件。 请注意,从 λ(Δ2)βγΔβmin{1,1/γΔ2 开始,就得出 ylog(λmin{1,1/γΔ1}log((Δ2)βγΔβ) 简单的计算表明 (Δ2)βγΔβγβ ,因此到 37 时,我们有

    |h(y)| |h(log((Δ2)βγΔβ))|(βγ1)1β(Δ2)βγΔβ+1
    =βγ1(Δ2)(βγ1)1O(1/Δ).

F.2 27的证明

在本小节中,我们使用 [GL18] 的结果来证明 27 它们的潜在函数由其边际比率的导数隐式定义为

Φ(R)=ϕ(R)=min{βγ1αγlogλ+γβλ+1,1RlogλR}

对于仅依赖于 β,γ,λ 的常量 0α1(有关精确定义,请参阅 [GL18])。 在我们的背景下,对数比率的相应潜力是

Ψ(y)=ψ(y)=eyϕ(ey)=min{βγ1αγlogλ+γβλ+1ey,1logλey}

并且受常数限制,具体取决于 log(λ/γΔ1)ylogλβ,γ,λ,Δ

[GL18] 中的主要技术成果之一是表明树递归与势函数 Φ 收缩,导数 ϕ 的边界为存在仅依赖于 β,γ,λ 的正常数 C1,C2 的感觉,使得 C1ϕ(R)C2 对于所有 0Rλ 而言。 [GL18]将此类函数称为通用势函数

在我们的上下文中,我们可以得到 Ψ 是一个 (α,c) 势函数,它满足 4,但有一个取决于 γ,Δ 的常数 c 事实上,最坏的情况,我们有

maxy1,y2ψ(y2)ψ(y1)ψ(logλ)ψ(log(λ/γΔ1))=λβγ1αγlogλ+γβλ+1βγ1αlogλ+γβλ+1λγΔ=γΔ1.

更准确地说,我们从 [GL18] 得到以下结果,以对数边际比率表示。

Theorem 38.

假设β,γ,λ是满足β1γβγ1λ<(γβ)βγβγ1的非负实数。 那么对于取决于 β,γ,λ 的常数 0<α<1 和取决于 β,γ,λ,Δ 的常数 c>0 来说,函数 Ψ 是一个 (α,c) 势函数。

5结合,这使得O(nC)与仅取决于β,γ,λ,Δ的常量C混合。 我们注意到这比[GL18]中的相关性衰减结果弱,因为C不依赖于Δ,因此对于任意图表。

附录 G 混合速度稍快

在本节中,我们通过更仔细地考虑我们基于收缩证明的(非平凡)谱独立性界限与我们获得的(平凡)谱独立性界限之间的权衡,稍微优化某些反铁磁2自旋系统的混合时间结果在 Section A.2 中用于处理恒定大小的图形。

Proposition 39

假设 [n] 子集上的分布 μ 对于 ηimin{a,(ni1)b} 来说是 (η0,,ηn2) 谱独立的,对于某些 a00b1 那么从μ采样的格劳伯动力学具有至少1nΩ(abn)a的谱间隙

证明。

假设我们已经以 c 为条件,将元素的一部分“输入/输出”。 所得分布既是b(1c)n-光谱独立的,又是a-光谱独立的。 边界b(1c)n优于a的确切阈值c由下式给出

c=1abn

我们注意到这样的 c 仅在 01abn1 或等效的 bna 时才有意义。 现在,我们基于固定至多 c 顶点分数,对所有条件分布应用 a 谱独立界限。 否则我们应用(ni1)b-谱独立性。 我们获得最终的光谱间隙下界

1n(1b)(1c)nk=0cn(1ank1)

观察一下

(1b)(1c)n=(1b)abexp(a)

我们还有

k=0cn(1ank1) exp(ak=0cn1nk1)
exp(a(k=0n21nk1lognk=cn+1n21nk1log(1c)n))
exp(alog11c)
exp(alogbna)
(abn)a

将它们放在一起,我们就得到了所需的下界。

有了这个结果,我们可以将其应用于具有 βγΔ2Δ,γ1β=0,γ1 的反铁磁模型,因为查看 28 的证明,我们有这样的系统 Cn-大致与CO(1/Δ)光谱无关。

Corollary 40 (软约束)

修复整数Δ31<Δ¯<Δ β,γ,λ0为满足Δ¯2Δ¯βγΔ¯1Δ¯+1γ1的非负实数。 进一步假设 (β,γ,λ) 最多可达 Δ 且具有间隙 0<δ<1 然后,对于每个 n 顶点图 G,其最大度数至多为 Δ,从具有参数 (β,γ,λ)O(Δ¯nΔ)O(1/δ) 步骤混合。

Corollary 41 (硬约束)

修复整数 Δ3,修复 β=0,并让 0γ1,λ0 最多为 Δ 唯一,间隙为 0<δ<1 然后,对于每个 n 顶点图 G,其最大度数至多为 Δ,从具有参数 (β,γ,λ)-以O(nΔ)O(1/δ)步骤混合。