自旋系统中相关性衰减至唯一性

Liang Li
Shandong University
li.liang@sdu.edu.cn
This work was done when this author visited Microsoft Research Asia.
   Pinyan Lu
SHUFE
lu.pinyan@mail.shufe.edu.cn
   Yitong Yin
Nanjing University
yinyt@nju.edu.cn
摘要

我们给出了在一般图上具有强空间混合的二态反铁磁自旋系统的完整表征。 我们证明,当且仅当系统在无限正则树上具有唯一的吉布斯测度时,二态反铁磁自旋系统在所有最大度数最多 Δ 的图上具有强空间混合到 Δ,其中 Δ 可以是有界的,也可以是无界的。 因此,当在度数高达 的无限正则树上满足唯一性条件时,在最大度数最多 Δ 的图上,存在一个二态反铁磁自旋系统的配分函数 FPTAS Δ 特别是,如果所有无限正则树都满足唯一性,则任意图都存在 FPTAS。 作为特殊情况,这涵盖了一般结构图上两态反铁磁系统的所有先前算法结果。

结合 Jerrum-Sinclair 和 Goldberg-Jerrum-Paterson 的二态铁磁自旋系统的 FPRAS,以及 Sly-Sun 和独立于 Galanis-S̆tefankovic̆-Vigoda 的最新硬度结果,这给出了完整的分类,除了相变边界,所有二态自旋系统的近似性,在有度限制的图族或所有图族上。

1简介

自旋系统在统计物理、应用概率和计算机科学领域得到了深入研究,作为捕捉局部相互作用和约束如何影响粒子系统宏观特性的本质的通用框架。 系统通常用图来描述,每个顶点处于固定数量的状态之一(称为自旋),边指定系统的邻域关系。

G(V,E) 为图表,q 为自旋状态数。 系统的配置是q|V|可能的分配σ:V[q]之一。 每个配置σ都有一个能量H(σ)作为所有边和顶点的总和,使得边(u,v)E的贡献由对称函数确定自旋状态σ(u)σ(v),顶点vV的贡献由其自旋状态σ(v)的函数确定。 配置σ的权重是w(σ)=exp(H(σ)T),其中T是温度。 我们专注于两种状态的自旋系统。 在归一化之前,二态自旋系统完全由三个参数 (β,γ,λ) 捕获,其中 βγ 确定边缘贡献的对称函数,λ,也称为外部场,决定顶点贡献的函数。 吉布斯测度是所有配置的自然概率分布,配置 σ 的概率为 w(σ)Z,其中归一化因子 Z=w(σ) 称为配分函数。 配分函数编码了有关自旋系统宏观行为的丰富信息。 然而,对于几乎所有重要的设置来说,计算配分函数的精确值是非常困难的。

自旋系统最重要的特性之一是相关性衰减,它表示两个顶点的边缘吉布斯分布之间的相关性随着两个顶点之间的距离而快速衰减。 此属性也称为弱空间混合 [26] 具有更大算法意义的是强烈的空间混合,这表明在其他顶点存在任意固定自旋的情况下相关性会衰减。 对于二态自旋系统,强空间混合可能意味着配分函数的有效近似算法。 那么,通过系统的参数来表征在图的任意实例上表现出强空间混合的系统是一个重要的问题。

如果相邻顶点有利于自旋一致,则二态自旋系统是铁磁的,否则是反铁磁的。 对于所有二态铁磁自旋系统,配分函数可以有效地近似[12, 10] 在反铁磁区域,相关衰减在配分函数的近似性中起着核心作用。 人们认为(参见[19]),对于此类模型,配分函数的逼近性由无限正则树上吉布斯测度的唯一性来表征,这相当于无限正则树上的弱空间混合普通树。 该条件称为唯一性条件

在开创性的著作[26]中,Weitz 表明,对于硬核模型来说,强空间混合的特点是唯一性条件,而在 Sinclair、Srivastava 和 Thurley 最近的著作中[21] 反铁磁伊辛模型也证明了相同的特征。 这两个模型都是重要的特殊双态反铁磁自旋系统。 在硬度方面,Sly [23] 和 Galanis 等人 [6] 为硬核模型证明了这一点,并且最近对于 Sly 和 Sun [24] 以及 Galanis、S̆tefankovic̆ 和 Vigoda [7] 独立提出的一般二态反铁磁自旋系统,违反唯一性条件意味着配分函数的不可近似性。 对于我们对二态自旋系统中的相关衰减和计算的完整理解,有两个问题仍然悬而未决:强空间混合的表征和一般二态自旋系统的近似性。

1.1我们的结果

我们通过无限正则树上吉布斯测度的唯一性来表征二态反铁磁自旋系统,该系统在一般有度有界图或任意图上表现出强烈的空间混合。

Theorem 1

对于任何有限Δ2Δ=,二态反铁磁自旋系统在最大度图上最多Δ具有强空间混合当且仅如果系统在所有dΔ的无限d-正则树上表现出唯一性。

由于Weitz的自回避游走树构造[26],有度有界图上的二态自旋系统的强空间混合立即意味着配分函数的FPTAS。 事实上,我们展示了之前的工作[15]中引入的更强大的相关性衰减概念,即计算有效的相关性衰减,它不仅为有度有界图提供 FPTAS,还为任意图提供 FPTAS当满足相应的唯一性条件时,具有无界度。

Theorem 2.

对于任意有限的Δ3Δ=,在最多Δ的图上,二态反铁磁自旋系统的配分函数存在一个FPTAS > 如果对于所有dΔ,系统参数位于无限d-正则树的唯一性区域的内部。

在上述两个定理中,Δ=情况代表无界度图。

由于 Sly 和 Sun [24] 针对一般二态自旋系统的最新硬度结果以及 Galanis、S̆tefankovic̆ 和 Vigoda [7] 针对一般二态自旋系统的独立结果不太通用的设置,违反唯一性条件意味着配分函数的不可近似性。

Theorem 3 (由于[24][7])

对于任何有限的Δ3Δ=,除非NP=RP,不存在用于二态反铁磁自旋系统在图上的配分函数的FPRAS如果对于某些dΔ,系统参数位于无限d-正则树的非唯一性区域的内部,则最大度至多Δ

[24] 中的原始定理适用于 d - 具有固定 d 的正则图,这立即意味着有度有界图或任意图的硬度。 [24][7] 中的硬度条件要求在具有 d3d 不规则树上的非唯一性。 但无限2-正则树(即无限路径)上的唯一性对于任何二态反铁磁自旋系统始终成立,因此定理3中的条件就足够了。

对于最大次数为 2 或更小的图,可以在多项式时间内精确计算配分函数。 因此,定理2和定理3与两态铁磁自旋系统[12, 10]的FPRAS一起给出了几乎完整的(除了在相变边界)所有二态自旋系统的配分函数的近似性分类,在有度有界图族或所有图族上。

1.1.1 规律性和单调性

在统计物理学中,相关性衰减通常在正则图甚至结构对称图(例如 Bethe 晶格)上研究。 尽管对于仅考虑正则图的硬度结果将加强下界,但从算法的角度来看,考虑通用图上的自旋系统更为普遍。 [26]中硬核模型的近似算法和反铁磁Ising模型的[21]中的近似算法均适用于具有有限最大度的一般图。 正如[21, 24]中所讨论的,对于具有固定dd规则图上的一般双态自旋系统,在参数平移之前,硬核模型和伊辛模型是完整的。然而,这并不包括最一般的情况,即一般图上有界或无界度的一般双态自旋系统。 造成这种差距的根本原因是一般自旋系统的非单调性。

充分研究的硬核和反铁磁 Ising 模型(以及所有具有 β,γ1 的双态自旋系统)都是单调自旋系统,从某种意义上说,无限 d 上的唯一性 -正则树意味着所有具有较小度数的无限正则树的唯一性。 这种单调性在一般的二态自旋系统中并不一定成立。

[26]中,Weitz 在硬核模型中建立了以下含义:

Claim ([26]中的定理2.3)

d-正则树上的强空间混合意味着最大程度至多d的所有图上的强空间混合。

[26]中,Weitz在没有证明的情况下也指出这一推论对于所有二态自旋系统都成立(事实上,这在[21]中的反铁磁伊辛模型中得到了严格证明>)。 有了这一点,为有度有界图上的二态自旋系统设计近似算法就可以简化为验证 d 正则树上的强空间混合。 这已被接受为关于二态自旋系统的事实,并已成为此类系统的近似算法的构建块(参见[21]中的定理2.8以及[22中的讨论) ,23])。 [22]中还提出了该主张是否适用于一般多态自旋系统的猜想。

作为我们分析的副产品(参见第 5 节),我们发现 d 正则树和最大度图上的强空间混合之间的这种可信的含义大多数 d 仅适用于单调自旋系统(包括硬核和反铁磁 Ising 模型)。 对于一般的二态自旋系统,唯一性以及所有有度有界图之间的强空间混合的最坏情况确实是常规树,但可能不再是最高度的树。 这种复杂性的一个好处是,较高的度数可能会产生更快的相关性衰减,从而使 FPTAS 能够用于具有无界度数的图。

这些新现象表明,一般的二态自旋系统比经过充分研究的单调自旋系统(例如硬核和反铁磁伊辛模型)具有更丰富的结构。 前一种方法通过d正则树上的强空间混合,在一般图上的单调自旋系统和正则图上的一般自旋系统上取得了成功,但在处理一般图上的一般自旋系统时遇到了障碍。 我们通过任意树而不是d-常规树上的强空间混合,给出了一般二态自旋系统中相关性衰减的统一方法。 这种方法是在我们之前处理无界度图的工作[15]中发起的。 在本文中,我们设计了一种统一的基于势的分析,它适应任意树的不规则性和一般二态自旋系统的非单调性,并给出所有二态反铁磁自旋系统的紧密相关衰减结果关于图的有度族和所有图的族。

1.1.2 主要结果的意义

通过解决唯一性条件,我们可以以各种阈值形式重申我们的主要结果。 定理2作为特殊情况涵盖了一般结构图上二态反铁磁自旋系统的所有先前算法结果,并清除了以前未发现的情况。

在互动方面。

我们可以固定外场λ并讨论(β,γ)的可处理区域。 由于βγ的作用是对称的,我们可以进一步固定其中一个并讨论另一个的可处理范围。 [10, 15]中使用了该公式。

我们的主要结果可以重述如下:对于任何Δ,对于无限正则树上的唯一性有一个临界阈值γc(β,λ,Δ),直到程度Δ,这样如果γc(β,λ,Δ)<γ<1β存在最大度数最多为Δ的图的FPTAS;特别地,γc=γc(β,λ,)>1是绝对常数,因此如果γ(γc,1β)则任意图都存在FPTAS。

作为特殊情况,这涵盖了[10]中关于反铁磁系统的所有算法结果以及[15]中的所有结果,扩展了这些先前作品中的可处理区域,并且也考虑有度有界图。

在外场方面。

在硬核伊辛模型和反铁磁伊辛模型研究的推动下,我们可以固定(β,γ)并讨论外场的可处理范围。

由于βγ的对称作用,我们可以假设βγ 我们的主要结果可以在特定设置下重述如下:

  • 硬约束(当 β=0 时):对于任何 Δ,λc(γ,Δ)=min1<d<Δγd+1dd(d1)d+1 是无限正则树上唯一性的关键阈值,最高可达 Δ 度> 这样,如果 λ<λc(γ,Δ) 存在最大度数最多 Δ 的 FPTAS。

    对于γ1,临界阈值等于λc(γ,Δ)=γΔ(Δ1)Δ1(Δ2)Δ 在无界度的无限正则树上不存在满足唯一性的外场λ>0 这与没有度数限制的硬核模型的硬度结果[4, 23]一致。 一个特别有趣的情况是当γ=1时,在这种情况下,模型是具有逸度λ的硬核模型,并且λc(1,Δ)=(Δ1)Δ1(Δ2)Δ是临界阈值。 这涵盖了[26]的结果。

    对于 γ>1,除了度受限图的结果外,还存在一个绝对正常数 λc(γ)=mind>1γd+1dd(d1)d+1,它可以降低所有 Δλc(γ,Δ) 的边界,这样,如果 λ<λc(γ),就存在一个适用于任意图的 FPTAS。

  • 软约束(当 β>0 时):如果 βγ>Δ2Δ 的唯一性条件对于任何外部域 λ>0 都成立,并且对于最大度最多为 Δ 的图,总是存在一个 FPTAS。

    现在假设βγΔ2Δ Δ¯为满足βγd1d+1的最小d 对于每个整数 Δ¯d<Δ,对于无限 (d+1) 正则树上的唯一性,存在两个阈值 λ1(β,γ,d)λ2(β,γ,d),这样如果λ<λ1(β,γ,d)λ>λ2(β,γ,d)对于所有整数Δ¯d<Δ,最大度数最多为Δ的图存在一个FPTAS。

    此外,如果γ1,则上述唯一性条件可以简化为λ<λcλ>λ¯c,其中λc=λc(β,γ,Δ)minΔ¯d<Δλ1(β,γ,d)λ¯c=λ¯c(β,γ,Δ)maxΔ¯d<Δλ2(d) 特别是,对于反铁磁伊辛模型β=γ,我们有λ¯c=1/λc,并且当|logλ|>logλ¯c时唯一性成立。 这涵盖了[21]中反铁伊辛模型的结果。

    对于无界最大度Δ,如果γ1,则在所有度的无限正则树上不存在满足唯一性的λ>0,这与硬度结果一致[10];如果γ>1,则当λ<λ1(β,γ,d)λ>λ2(β,γ,d)对于所有dΔ¯时,所有无限正则树的唯一性成立,在这种情况下,存在用于任意图的 FPTAS。

1.2相关作品

自旋系统配分函数的近似已得到广泛研究[12,13,9,3,11,25,4,5,17] 在一项开创性的工作[12]中,Jerrum 和 Sinclair 为铁磁 Ising 模型设计了完全多项式时间随机逼近方案 (FPRAS)。 后来在[10]中,通过将参数转换为铁磁伊辛模型,将 FPRAS 扩展到所有二态铁磁自旋系统。 同样在[10]中,Goldberg、Jerrum 和 Paterson 给出了任意图上两态反铁磁自旋系统的 FPRAS 和不可逼近性结果。 Dyer、Frieze 和 Jerrum 在 [4] 中提出了一种基于随机正则二部图的小工具,Mossel、Weitz 和 Wormald 在 [19] 中也对此进行了分析研究有界图上的不可逼近性。 人们普遍认为,反铁磁自旋系统的近似性转变是由无限树上的唯一性相变捕获的。 这是在[19]中作为猜想公开提出的。 Sly在[23]中针对硬核模型证明了该猜想。 Galanis 等人[6]中将该结果改进为多种参数。 最近,Sly 和 Sun [24] 证明了无限正则树上所有非唯一二态反铁磁自旋系统的硬度。 Galanis、S̆tefankovic̆ 和 Vigoda 在[7]中独立证明了无外场反铁磁 Ising 模型的相同结果。

由 Weitz [26] 以及 Bandyopadhyay 和 Gamarnik [1] 独立开发的相关衰减技术是设计分区的确定性完全多项式时间近似方案 (FPTAS) 的强大工具函数(其他重要示例包括 [8, 2])。 [26]中,Weitz引入了强空间混合的概念,并用它为硬核模型设计了达到唯一性阈值的FPTAS。 Sinclair、Srivastava 和 Thurley 最近在[21]中研究了另一个最重要的二态反铁磁自旋系统,即反铁磁伊辛模型,其中提出了一种更强大的消息衰减方法引入来分析强烈的空间混合并使 FPTAS 达到唯一性阈值。 Restrepo 等人[20]中开发了一种强大的技术,它利用图的特定结构进行强空间混合。 通过利用图的结构,在网格上实现了比唯一性区域更广泛的可处理区域。 在之前的工作[15]中,我们给出了在具有无界度的任意图上无外场的二态反铁磁自旋系统的FPTAS,直到唯一性阈值的连续松弛。 本文使用的方法是在[15]中提出的,但是[15]中的分析无法在一定程度上分离唯一性,因此无法处理程度有界图族。

2 定义和基本知识

两态自旋系统由图G=(V,E)描述。 系统的配置是顶点状态的2|V|可能分配σ:V{0,1}之一。 我们还使用蓝色绿色两种颜色来表示这两种状态。 配置的权重可以描述为各个边和顶点的贡献的乘积。 A=[A0,0A0,1A1,0A1,1]b=(b0,b1) 配置σ:V{0,1}的权重由下式给出

w(σ)=(u,v)EAσ(u),σ(v)vVbσ(v).

吉布斯测度ρ(σ)=w(σ)Z(G)定义的所有配置的概率分布。 归一化因子Z(G)=σw(σ)称为配分函数

我们可以将 {𝚋𝚕𝚞𝚎,𝚐𝚛𝚎𝚎𝚗} 边和绿色顶点的贡献标准化为 1。 因此 A=[β11γ] 代表某些 β,γ0,b=(b0,b1)=(λ,1) 代表某些 λ>0 由于蓝色和绿色的角色是对称的,我们可以假设βγ而不失一般性。 三个参数(β,γ,λ)0βγλ>0完全指定了一个二态自旋系统。 具有β=γ的二态自旋系统是伊辛模型,具有β=0,γ=1或对称β=1,γ=0的二态自旋系统是硬核模型

如果相邻顶点有利于不同的自旋,即如果βγ<1,则二态自旋系统称为反铁磁 不失一般性,我们重点关注βγ的情况。

Definition 4

(β,γ,λ) 如果0βγγ>0βγ<1λ>0,则为反铁磁

通过 βγ 的对称性以及情况 β=γ=0 的平凡性,对于所有非平凡的二态反铁磁系统,该定义是完整的

2.1相关性衰减

吉布斯测度定义了每个顶点的状态边缘分布。 pv 表示顶点 v 被着色为蓝色的概率。 σΛ{0,1}ΛΛV 的配置。 我们将顶点vΛ 固定顶点称为vΛ自由顶点。 我们使用pvσΛ表示在Λ的配置固定为σΛ的条件下v被着色为蓝色的边际概率。

Definition 5

如果对于图族中的任意图 G=(V,E)、任意 vV,ΛVσΛ,τΛ{0,1}Λ 来说,图族上的自旋系统具有 强空间混合性

|pvσΛpvτΛ|exp(Ω(dist(v,S))),

其中 SΛσΛτΛ 不同的子集,dist(v,S) 是距 v 的最短距离到 S 中的任意顶点。

弱空间混合可以通过测量相对于dist(v,Λ)而不是dist(v,S)的衰减来定义。 空间混合特性在统计物理学中也称为相关衰减。

T 为一棵以 v 为根的树。将 RTσΛ=pvσΛ/(1pvσΛ) 定义为在施加条件 σΛ 时根 v 为蓝色和绿色的概率之间的比率(约定为 RTσΛ=pvσΛ=1 时)。 假设 vd 个子级,并且 Ti 是以第 i 个子级为根的子树。 由于子树的独立性,我们有一个简单的递归来计算RTσΛ

RTσΛ =λi=1dβRTiσΛ+1RTiσΛ+γ. (1)

G(V,E) 为图表。 类似地定义RG,vσΛ=pvσΛ/(1pvσΛ) 与树的情况相反,由于循环引起的依赖性,对于一般图G来说,计算RG,vσΛ并不容易。 [26]中,引入了一种称为自回避游走(SAW)树的结构,它将一般图中边缘分布的计算减少到树中的计算。 具体来说,给定一个图 G(V,E) 和一个顶点 vV SAW 树 TSAW(G,v) 是一棵以 v 为根的树,具有新的顶点集 VSAW(它有效地枚举了源自 v 的所有路径) G 并且可能包括固定叶子)。 此外,任何顶点集SΛV都映射到相应的SSAWΛSAWVSAW,任何配置σΛ{0,1}Λ都映射到σΛSAW{0,1}ΛSAW 如果不会引起歧义,我们会滥用这种表示法并编写 S=SSAWσΛ=σΛSAW 给定图形 G(V,E)vVSV,令 distG(v,S)G 中距 vS 中的任意顶点。

Theorem 6 (Weitz定理3.1[26])

G(V,E) 为图形,vVσΛ{0,1}ΛΛV 上的配置,以及 SV T=TSAW(G,v) 其认为T的最大阶数等于GdistG(v,S)=distT(v,S)RG,vσΛ=RTσΛ的最大阶数。 此外,Tv 的任何邻域都可以在与邻域大小成正比的时间内构建。

2.2 唯一性条件

我们考虑无限 (d+1) 不规则树 𝕋^d 上吉布斯度量的唯一性,其中由于 𝕋^d 的对称结构,递归由 fd(x)λ(βx+1x+γ)d 给出。

x^dfd(x)的正不动点,即x^d=f(x^d) 已知[14, 18],𝕋^d上的二态反铁磁自旋系统在|fd(x^d)|=1处发生相变,并且当|fd(x^d)|1 这引发了以下定义。

Definition 7

(β,γ,λ) 是反铁磁性的。 x^d为函数fd(x)=λ(βx+1x+γ)d的正不动点。我们说 (β,γ,λ)up-to-Δ 唯一的,如果对于所有整数 1d<Δ

|fd(x^d)|=λd(1βγ)(βx^d+1)d1(x^d+γ)d+1=d(1βγ)x^d(βx^d+1)(x^d+γ)<1.

特别是,如果(β,γ,λ)最多唯一,那么它就是普遍唯一

达到Δ的唯一性相当于系统在无限正则树上的弱空间混合达到Δ度。 唯一性条件可以用各种阈值形式来描述,在附录A中给出。

唯一性是通过要求 |fd(x^d)|<1 来定义的。 以下引理指出,只要唯一性条件成立,|fd(x^d)| 就受到绝对常数的限制。

Lemma 8.

(β,γ,λ) 是反铁磁性的。 如果 (β,γ,λ) 最多为 Δ 唯一,则存在一个绝对常数 c<1,它仅取决于 βγλΔ,例如 |fd(x^d)|c 适用于所有 1d<Δ

证明。

对于有限的Δ,引理是微不足道的。 接下来需要证明的是,在普遍唯一性的情况下,随着 d 增长到无穷大,|fd(x^d)| 不能任意接近 1。 如果(β,γ,λ)是普遍唯一的,由于引理21.2,我们必须有γ>1 对于反铁磁(β,γ,λ),β1γ,因此固定点x^d=λ(βx^d+1x^d+γ)dλγd,因此|fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ)dλγd 引理如下。

3 一般图上的强空间混合

本节我们证明定理1 唯一性条件的必要性是微不足道的,因为一般图上的强空间混合意味着常规树上的弱空间混合。 接下来就剩下证明下面的定理了。

Theorem 9

(β,γ,λ) 是反铁磁性的。 对于任何有限的Δ2Δ=,如果(β,γ,λ)最多是Δ唯一,则参数的二态自旋系统(β,γ,λ)在最多Δ的最大度图上具有强空间混合。

我们的方法是证明任意树上最大度至多Δ的强空间混合,这通过自回避游走树构造隐含了定理。 请注意,与 [26][21] 不同,我们分析任意树而不是常规树的衰减。 这是因为对于一般的二态反铁磁自旋系统,最大程度d的所有图之间强空间混合的最坏情况可能不再是d-正则树。 我们将在5节中详细解释这一点。

给定最大度最多为 Δ 的任意图 G(V,E)ΛV 上的任意配置 σΛ{0,1}Λ 和任意 SΛ,固定任意顶点 vV、根据定理 6,可以构造一棵自避走树 T=TSAW(G,v),使得 T 的最大度受 ΔdistG(v,S)=distT(v,S)RG,vσΛ=RTσΛ 的约束。 回想一下 RG,vσΛ=pvσΛ/(1pvσΛ) 因此 pvσΛ=RG,vσΛ/(1+RG,vσΛ) 对于任何 σΛτΛ

|pvσΛpvτΛ|=|RG,vσΛ1+RG,vσΛRG,vτΛ1+RG,vτΛ||RG,vσΛRG,vτΛ|=|RTσΛRTτΛ|.

这引发了以下定义。

Definition 10

T 为以顶点 v 为根的树,τΛ{0,1}ΛΛV 上的配置,SΛ 为顶点集。 Rvδv 定义为:对于所有 σΛ{0,1}ΛRvRTτΛRv+δv 只在 S 上与 τΛ 不同。

通过为 T=TSAW(G,v) 构造这样的 Rvδv 并证明 δvexp(Ω(dist(v,S))) 就足以证明定理 9 >。

T 为一棵以 v 为根的树,它有 d 个子节点 v1,,vd,Ti 为子树以 vi 为根。 它认为

RTσΛ =f(RT1σΛ,,RTdσΛ)λi=1dβRTiσΛ+1RTiσΛ+γ.

下限和上限RvRv+δv可以递归地构造如下。 基本情况是: (1) vS,在这种情况下Rv=0δv=; (2) vΛS,即v固定为所有σΛ中相同的颜色,此时δv=0Rv=(或Rv=0)如果v固定为蓝色(或绿色)。 对于vΛ,由于f(R1,,Rd)对于反铁磁(β,γ,λ)是单调递减的,

Rv =f(Rv1+δv1,,Rvd+δvd),
Rv+δv =f(Rv1,,Rvd),

其中 RviRvi+δviRTiσΛ1id 相应的下限和上限。 特别是,当d=0时,按照惯例,空乘积等于1,因此Rv=Rv+δv=λ,这与v是没有子节点的自由顶点的情况一致。

通过f(R1,,Rd)的单调性,很容易检验上面构造的RvRv+δv是否满足定义10的要求。 我们的目标是证明,当唯一性成立时,δv 的递归深度呈指数衰减。 一种简单的方法是尝试证明 δ 在每个递归步骤中以恒定速率收缩。 但这不一定是真的才能保证指数衰减。 事实上,在某些情况下,误差不会单步衰减,而是长期衰减。 为了克服这个问题,我们使用潜力 Φ来摊销收缩并表明δΦ以恒定速率收缩。 我们选择势函数为

Φ(R)=1R(βR+1)(R+γ).

我们正在分析任意树上具有不规则程度的腐烂。 为了适应这种不规则性,势函数不能以d作为输入,而只能携带有关当前顶点的分布信息,但它必须能够为步骤提供正确的补偿 -在 R 的任何状态下以及对于满足足够唯一性的所有自旋系统而言,明智的衰减。 附录C中讨论了引导我们得到这个良好势函数的启发式过程。

φ(R)为满足φ(R)=Φ(R)的单调函数。 我们定义

ϵvφ(Rv+δv)φ(Rv),

相应地,ϵviφ(Rvi+δvi)φ(Rvi)1id

我们定义一个函数α(d;x1,,xd)如下:

α(d;x1,,xd) (1βγ)(λi=1dβxi+1xi+γ)12(βλi=1dβxi+1xi+γ+1)12(λi=1dβxi+1xi+γ+γ)12i=1dxi12(βxi+1)12(xi+γ)12.

以下引理是通过应用中值定理获得的。 之前在 [15, 20] 中使用了类似的例程。 引理的证明在附录B中。

Lemma 11.

以下内容适用于ϵv

  1. 1.

    (与δv相关) ϵv=δvΦ(R~)对于某些R~[Rv,Rv+δv].

  2. 2.

    (绝对界限) 假设 γ>1T 的最大度以常数为界,如果 vΛRv+δv=O(1) ,如果 viS 对于所有 1idϵv=O(1)

  3. 3.

    (逐步收缩)存在Ri~[Rvi,Rvi+δvi],1id,使得

    ϵvα(d;R~1,,R~d)max1id{ϵvi}.

To prove the strong spatial mixing, we first relate δv to ϵv by Item 1 of Lemma 11, and then apply induction on the depth in T, with Item 2 of Lemma 11 as basis and Item 3 of Lemma 11 as induction step. 然后我们需要限制收缩率α(d;x1,,xd) 请注意,z(βz+1)(z+γ)(1+βγ)21 对于 z[0,),因此它无条件地适用于所有 xi[0,)1id,即

α(d;x1,,xd) d, (2)
α(d;x1,,xd) dλγd+1. (3)

凭借唯一性,可以证明以下更紧的收缩界。

Lemma 12.

(β,γ,λ) 是反铁磁性的。 如果 (β,γ,λ) 最多为 Δ 唯一,则存在一个常量 α<1,使得对于任何整数 1d<Δ 和任何 x1,,xd[0,+)1id,则成立α(d;x1,,xd)α

这个引理是我们分析的技术核心。 它关键取决于势函数的选择。 在深入研究引理12的形式证明之前,我们注意到该引理可以隐含定理9

定理证明9

T=TSAW(G,v)G,其最大度数至多为 Δ 那么 T 的最大度最多为 Δ,因此根 vT 中最多有 Δ 个子顶点,而 T 中的每个其他顶点都有少于 Δ 个子顶点。 我们为 T 中的每个子树递归构造 Ruδuϵu

t=dist(v,S) 在不损失一般性的前提下,反复应用例题 113 项,我们在 T 中有一条路径 u1u2ut2u1=v,使得 ϵujα(dj;x1,,xdj)ϵuj+1 对于 j=1,2,,t3,其中 djujxi[0,) 的子节点数,1idj

请注意,d1Δdj<Δ 代表所有其他 j 假设 (β,γ,λ) 最多为 Δ 是唯一的。 如果 Δ 是有界的,那么根据定理 12,存在一个常数 α<1,使得 ϵujαϵuj+1 对于 2jt3,并且 ϵvd1ϵu2Δϵ2 由于 (2);如果 Δ=,那么根据定理 12,对于所有 1jt3 ϵujαϵuj+1 在这两种情况下我们都有ϵv=O(αtϵut2)

由于引理 11 的第 1 项,δv=ϵvΦ(R~)=O(1Φ(R~)αtϵut2) 对于某些 R~[Rv,Rv+δv] 而言。 然后我们绑定Φ(R~)ϵut2 由于引理 21 的项 2,(β,γ,λ) 最多 Δ 唯一这一事实意味着 Δ 是有界的或 γ>1 请注意,v 必须是自由的,否则该定理的证明很简单,并且 ut2 的子级都不在 S 中,因为 dist(v,S)=t 因此,根据引理11ϵut2=O(1)R~Rv+δv=O(1)的项2,这意味着Φ(R~)=1R~(βR~+1)(R~+γ)=Ω(1)

总之,如果(β,γ,λ)最多是Δ唯一,则存在一个常量α<1,使得δv=O(αt) 正如本节开头所讨论的,这证明了定理 9

引理12的证明。

本节的其余部分致力于引理 12 的证明。 鉴于 (β,γ,λ) 最多 Δ 是唯一的,因此存在绝对常数 α<1,使得 α(d;x1,,xd)α 对于任何 1d<Δx1,,xd0

我们定义 α(d;x1,,xd) 的对称版本:

αd(x) α(d;x,,xd)=d(1βγ)(xλ(βx+1x+γ)d)12(βx+1)12(x+γ)12(βλ(βx+1x+γ)d+1)12(λ(βx+1x+γ)d+γ)12.

以下引理表明,通过使用 Cauchy-Schwarz 不等式以及算术和几何平均值,对称情况主导 α(d;x1,,xd) 的最大值。

Lemma 13.

(β,γ,λ) 是反铁磁性的。 那么对于任何整数d和任何x1,,xd[0,+),都存在一个x¯[0,+)使得α(d;x1,,xd)α(d,x¯)

证明。

zi=βxi+1xi+γ 然后是zi(β,1γ]xi=1γziziβ zi表达α(d;x1,,xd)

α(d;x1,,xd) =(λi=1dzi)12(βλi=1dzi+1)12(λi=1dzi+γ)12i=1d(zi1γ)12(ziβ)12.

由于柯西-施瓦茨不等式,

i=1d(zi1γ)12(ziβ)12d(1di=1d(zi1γ)(ziβ))12=d(1+βγ1di=1d(ziγ+βzi1))12.

由于算术平均值和几何平均值不等式,

d(1+βγ1di=1d(ziγ+βzi1))12 d(1+βγγ(i=1dzi)1dβ(i=1dzi)1d)12.

z¯=(i=1dzi)1d 然后结合上面的计算,

α(d;x1,,xd) (λz¯d)12d(1+βγγz¯βz¯1)12(βλz¯d+1)12(λz¯d+γ)12=dλz¯d(z¯1γ)(z¯β)(βλz¯d+1)(λz¯d+γ).

假设 x¯βx¯+1x¯+γ=z¯ 然后 x¯[0,+) 并用 βx¯+1x¯+γ 替换 z¯,我们有

α(d;x1,,xd) d(1βγ)(x¯λ(βx¯+1x¯+γ)d)12(βx¯+1)12(x¯+γ)12(βλ(βx¯+1x¯+γ)d+1)12(λ(βx¯+1x¯+γ)d+γ)12=αd(x¯).

Lemma 14.

(β,γ,λ) 是反铁磁性的。 如果 (β,γ,λ) 最多为 Δ 唯一,则存在一个常量 α<1,这样对于任何整数 1d<Δ,它都满足αd(x)α 对于所有 x0

证明。

d 固定为任意正整数。 我们描述 αd(x) 达到最大值时 x 的值。 我们表示fd(x)=λ(βx+1x+γ)d。对 αd(x)x 求导,我们得到

αd(x)=d(1βγ)G(x)2G(x),

其中G(x)=xfd(x)(βfd(x)+1)(fd(x)+γ)(βx+1)(x+γ),其导数为

G(x) =fd(x)d(1βγ)x(βfd(x)+1)(fd(x)+γ)(βx+1)2(x+γ)2(γβx2d(1βγ)xγβfd(x)2(βfd(x)+1)(fd(x)+γ)).

由于 x 的范围超过 [0,),因此函数 γβx2d(1βγ)xx 中严格递减,范围从 +,函数γβfd(x)2(βfd(x)+1)(fd(x)+γ)x中严格递增,并且有一个有界范围。 因此,方程

γβx2d(1βγ)x =γβfd(x)2(βfd(x)+1)(fd(x)+γ). (4)

(0,)有唯一解,记为xd 此外,它认为

G(x){>0if 0x<xd,=0if x=xd,<0if x>xd. (5)

对于αd(x)也是如此。 因此,对于任何固定的d,αd(x)x=xd时达到最大值。

因此,对于所有x0

αd(x)αd(xd) =d(1βγ)(xdfd(xd)(βxd+1)(xd+γ)(βfd(xd)+1)(fd(xd)+γ))12
=(d(1βγ)(γβxd2)(βxd+1)(xd+γ)fd(xd)(γβfd(xd)2))12 (6)
α~d(xd).

(βfd(xd)+1)(fd(xd)+γ)代入(4),得到式(6)。

x^dfd(x)的正不动点,即x^d=fd(x^d) 然后我们声称,如果 (β,γ,λ) 最多是 Δ 唯一,则对于任何整数 1d<Δ,α~d(xd)α~d(x^d) 是唯一的,要看到这个声明足以暗示引理,请注意,替换 x^d=fd(x^d) 后,我们有 α~d(x^d)=d(1βγ)x^d(βx^d+1)(x^d+γ)=|fd(x^d)| 并且由于引理 8,如果 (β,γ,λ) 最多是 Δ 唯一,则存在一个常量 c<1 使得 |fd(x^d)|<c 对于所有整数 1d<Δ

然后我们证明这一主张。 假设 (β,γ,λ) 最多为 Δ 是唯一的且 1d<Δ 那么足以证明,如果 x^dxd,则 正在减少,如果 x^d>xd,则 α~d(x) 正在增加。

情况 1:x^dxd 由于 (5),G(x^d)0 注意

G(x^d)=d(1βγ)(γβx^d2)x^d2(βx^d+1)3(x^d+γ)3(1d(1βγ)x^d1(βx^d+1)(x^d+γ)).

由于唯一性,|fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ)<1,因此1d(1βγ)x^d1(βx^d+1)(x^d+γ)>0 G(x)0 结合,我们有 γβx^d20 由于函数fd(x)是单调递减的,并且x^d是其不动点,因此γβfd(xd)2γβfd(x^d)2=γβx^d20 由于xd满足(4),γβxd2γβfd(xd)2必须同时为正或负,因此也成立γβxd20 那么(γβx2)(βx+1)(x+γ)fd(x)(γβfd(x)2)都是正数,并且在x[x^d,xd]中单调递减。 因此,α~d(xd)α~d(x^d)

情况 2:x^d>xd 根据与上面相同的论证,它成立γβfd(x^d)2=γβx^d2<0γβfd(xd)2<0γβxd2<0 因此,(γβx2)(βx+1)(x+γ)fd(x)(γβfd(x)2)都是负数,并且在x[xd,x^d]中单调递减,因此它们的乘积是正数,并且在x[xd,x^d]中递增。 因此,α~d(xd)α~d(x^d)

引理12是通过组合引理1314来证明的。 这样就完成了定理9的证明。

普通树上的强烈空间混合。

作为我们分析的副产品,我们证明了规则树的强空间混合定理。 当图 G 本身是一棵常规树时。 所有顶点(根除外)都具有相同的数量。 并且证明中出现的所有d(摘录根的一个)都相同且等于该元数。 那么度数达到Δ的所有无限正则树上唯一性成立的条件就可以用无限Δ正则树上的唯一性来代替。

Theorem 15

对于二态反铁磁自旋系统,在任何无限Δ-正则树上,唯一性意味着强烈的空间混合。

结合硬核模型[26]和反铁磁Ising模型[21]的相同定理,并转换一般二态参数,可以得到相同的结果这些模型的反铁磁自旋系统(如[24, 21]中讨论)。 然而,与硬核模型和反铁磁 Ising 模型不同,对于一般二态反铁磁自旋系统,定理 15 本身不足以暗示最大程度图上的强空间混合 Δ 这将在 5 节中讨论。

4 算法含义

本节我们证明定理2 也就是说,如果反铁磁性 (β,γ,λ)Δ 是唯一的,那么对于最大度最多为 Δ 的任何图 G,存在分割函数 Z(G) 的 FPTAS,特别是普遍唯一性意味着任意图 G 的 FPTAS。

众所周知,Z(G) 可以通过以下标准过程从 pvσΛ 计算得出。 v1,,vn 枚举 G 中的顶点。对于 0in,令 σi 为固定第一个 i 顶点 v1,,vi 的配置,如下所示: σi(vj)=σi1(vj) for 1ji1σi(vi)被固定,使得piPr[σi(vi)σi1]1/3 特别是,σn{0,1}VV 的配置。对于 σn 的吉布斯度量,ρ(σn)=p1p2pn 以及 ρ(σn)=w(σn)Z(G)都成立,因此 Z(G)=w(σn)p1p2pn,其中权重 w(σn)=(u,v)EAσn(u),σn(v)vVbσn(v) 可以针对任何特定的 σnn 的多项式时间精确计算。请注意,pi 等于 pviσi11pviσi1 因此,如果pvσΛ可以在n1ϵ中时间多项式的加性误差ϵ内近似,则配置σi 可以有效地构造,使得所有 pi 都远离 0,因此乘积 p1p2pn 可以在时间上在 (1±nϵ) 因子内近似n1ϵ 中的多项式,这意味着 Z(G)

有界度图。

G 为最大度数至多为 Δ 的图,v 为任意顶点。 可以构造一棵自回避游走树T=TSAW(G,v),使得RG,vσΛ=RTσΛ 我们可以使用3节中描述的递归过程来计算RTσΛ的上下界,设置为对于所有大于t的顶点距离根 v 几步远,使用简单边界 0RTσΛ 那么定理9的证明表明,递归过程返回R0R1,使得R0RTσΛR1R1R0=O(αt) 对于一些常量 α<1 假设 (β,γ,λ) 最多是 Δ 唯一。 请注意RTσΛ=RG,vσΛ=pvσΛ1pvσΛ p0=R0R0+1p1=R1R1+1 然后 p0pvσΛp1

p1p0=R1R1+1R0R0+1R1R0=O(αt). (7)

递归过程在时间 O(Δt) 内运行,因为它只需要构造自回避游走树的前 t 层。 如果Δ有界,通过设置t=logαϵ,这给出了一个算法,该算法在时间多项式的加性误差O(ϵ)内逼近pvσΛn1ϵ,这意味着 Z(G) 的 FPTAS。

任意图形。

G 为任意图,v 为任意顶点。 T=TSAW(G,v) 我们使用[15]中引入的计算高效相关性衰减方法来处理无界度。 直观上,使用这种方法,我们观察到精细度量中的相关性衰减,而不是图距离,这样在这个新度量中,中等大小的邻域足以保证所需的相关性衰减。

类似地,我们使用3节中描述的递归过程来计算RTσΛ的上限和下限,但这次终止条件依赖于定义如下的新深度。

Definition 16

T 为有根树,M>1 为常量。 对于 T 中的任意顶点 v,定义 v 的基于 M 的深度,表示为 M(v):如果 v 是根节点,则定义为 M(v)=0;如果 vud 子节点之一,则定义为 M(v)=M(u)+logM(d+1)

M>1 被修复。 B() 表示具有基于 M 的深度 < 的所有顶点及其在 T 中的子代和孙子的集合。归纳起来可以验证|B()|n2M 当当前vB()不在内时,递归将用于估计RTσΛ,直到v不再在B()内,在这种情况下,将使用微分边界0RTσΛ

ϵv 的定义如 3 节中所述。 重复应用定理 113 项,在不违背一般性的前提下,我们在 T 中有一条路径 u1u2uk 从根 u1=vukM(uk)M(uk1)<、使得ϵujα(dj;x1,,xdj)ϵuj+1j=1,2,,k,其中djujxi[0,)1idj的子节点数。 克服无界度引起的爆炸的关键是观察收缩α(d;x1,,xd)随着度d的增加而急剧减小。

Lemma 17.

(β,γ,λ) 是反铁磁性的。 如果(β,γ,λ)满足通用唯一性,则存在常量α<1M>1,使得对于任何整数d1和任何xi[0,),1id,它成立α(d;x1,,xd)αlogM(d+1)

证明。

假设(β,γ,λ)具有普遍唯一性。 由于引理12,存在一个常数α<1使得α(d;x1,,xd)α 根据引理21的项2,通用唯一性意味着γ>1,因此存在一个常数M>1使得dλγd+1αlogM(d+1) 对于所有dM 由于 (3),它对于所有 dM 都成立 α(d;x1,,xd)dλγd+1αlogM(d+1) 请注意,α(d;x1,,xd)α 表示α(d;x1,,xd)αlogM(d+1) 代表d<M。因此,α(d;x1,,xd)αlogM(d+1)适用于所有d。 ∎

根据引理17,存在常量α<1M>1,使得

ϵv ϵukj=1kαlogM(dj+1)ϵukαj=1klogM(dj+1)=ϵukαM(uk)ϵukα.

对于 3 中使用的符号,SB() 的补集。 请注意,uk 的所有子代都在 B() 中,因此它们都不在 S 中,而根据 Lemma 212 项,普遍唯一性意味着 γ>1 因此,根据引理11的项2,它支持ϵuk=O(1) 因此,ϵvϵukα=O(α)

δv=R1R0,其中R0R1是递归过程返回的边界,使得R0RTσΛR1 通过与3δv=ϵvΦ(R~)=O(ϵv)=O(α)节中相同的分析。 然后通过(7),在O(α)的加性误差内近似边际概率pvσΛ 递归的运行时间为O(nB())=O(n3M) 通过设置 =logαϵ,我们得到了一种算法,它可以在 O(ϵ) 的加法误差范围内逼近 pvσΛ,所用时间为 n1ϵ 的多项式,这意味着对于任意图 GZ(G) 的 FPTAS。

异质自旋系统。

我们在上一节和本节中的分析实际上适用于异构自旋系统,该系统允许每个顶点 v 具有不同的恒定外场 λv>0

Theorem 18.

对于每个顶点v处具有参数βγ和外场λv的二态反铁磁异质自旋系统,对于任何有限的Δ2Δ=,如果对于所有v,(β,γ,λv)最多是Δ唯一则自旋系统具有强空间混合性,并且在最大度图上具有至多Δ的FPTAS。

5 一般二态自旋系统的非单调性

在本节中,我们证明以下定理。

Theorem 19.

存在二态反铁磁自旋系统,其在无限d-正则树上表现出强空间混合,但在至多d的所有最大度图上不表现出强空间混合。

在开创性的工作[26]中,Weitz证明了对于硬核模型,无限d正则树上的强空间混合意味着最大度图上的强空间混合最d([26]中的定理2.3)。 他进一步指出,同样的含义适用于所有两种状态的自旋系统。 那是,

对于任何二态自旋系统,无限d正则树上的强空间混合意味着最大度数最多d的图上的强空间混合。

这一主张对于当前对二态自旋系统相关衰减的理解以及为此类系统设计 FPTAS 发挥了重要作用。 [21] 中引用了该主张的算法形式作为所有二态自旋系统的定理([21] 中的定理 2.8),并被证明用于反-铁磁伊辛模型。 [22]中提出了这个主张是否适用于多状态自旋系统的猜想。

在这里我们澄清这一说法仅适用于有限设置下的二态自旋系统,但并不适用于所有一般的二态自旋系统。 这反驳了[22]中的猜想,并表明了人们普遍认为d-正则树代表了所有最大度图之间强空间混合的最坏情况,最多d 不能推广到一般的二态自旋系统。 我们首先描述一个区域,证明该说法是正确的。

Lemma 20.

对于0β,γ1,无限d-正则树上的强空间混合意味着最大度数最多为d的树上的强空间混合。

证明。

给定最大度数为 d 的有根树,对于具有 k<d1k 子代的每个顶点,我们可以附加 d1k 虚拟子代具有固定的自旋态(分布)。 这是[26][21]中使用的方法:对于硬核模型,虚拟子项被固定为空闲,对于反铁磁Ising模型,虚拟子项被固定为空闲子级在自旋态上呈均匀分布。 在这两种情况下,虚拟孩子对其父母没有影响。 一般来说,我们将每个虚拟子节点的分布固定为 (p0,p1) 满足

p0+p1 =1,
βp0+p1 =p0+γp1.

0β,γ1时,系统在0p0,p11中有解。 i 个子级的比率 Ri1ik 和虚拟子级 k<id1Ri=p0/p1 的比率,父级的比率由递归给出

λi=1d1βRi+1Ri+γ=λi=1kβRi+1Ri+γ,

与原始数量相同。

由于自回避游走树构造(定理6),它认为对于0β,γ1,无限d-正则树上的强空间混合意味着强空间混合和最大度至多d图上的FPTAS。

对于0β,γ1,自旋系统表现出以下单调性质:无限d正则树上的唯一性意味着所有更小度的无限正则树上的唯一性。 这可以通过以下推理来验证:根据定理15,在无限d-正则树上,唯一性意味着强烈的空间混合,对于0β,γ1 ,意味着所有较小度的无限正则树上的强空间混合(包括唯一性)。

存在非单调的二态反铁磁自旋系统。 我们可以选择满足γ>1λ>λc𝗎𝗉𝗉𝖾𝗋(β,γ)的反铁磁(β,γ,λ),其中λc𝗎𝗉𝗉𝖾𝗋(β,γ)是Lemma 218项中给出的普遍唯一性的上临界值。 由于 Lemma 213 项,对于 γ>1 而言,在所有足够大的 d 中,唯一性在 d 不规则树上成立。另一方面,由于 Lemma 218 项,当 λ>λc𝗎𝗉𝗉𝖾𝗋(β,γ) 时,(β,γ,λ) 不可能是普遍唯一的,这意味着存在一个有限的 d,使得系统在 d 不规则树上是非唯一的。

对于这样的非单调系统,由于定理15,唯一性意味着对于足够大的d,d正则树上存在强空间混合,但是强空间混合不适用于较小程度的规则树(因为较小树上的非唯一性)。 因此,d正则树上的强空间混合与最大度至多d图上的强空间混合之间的含义对于一般的二态自旋系统并不成立。

致谢。

我们要感谢 Alistair Sinclair 和 Piyush Srivastava 对本文早期版本的有益评论。

致谢 2021 年修订。

我们要感谢Xiaoyu Chen、Weiming Feng 和Xinyuan 张在论文原始(2013)版本中找到引理21 中的错误并提出修复建议。

参考

  • [1] A. Bandyopadhyay and D. Gamarnik. Counting without sampling: Asymptotics of the log-partition function for certain statistical physics models. Random Struct. Algorithms, 33:452–479, December 2008.
  • [2] M. Bayati, D. Gamarnik, D. Katz, C. Nair, and P. Tetali. Simple deterministic approximation algorithms for counting matchings. In Proceedings of the thirty-ninth annual ACM Symposium on Theory of Computing, STOC ’07, pages 122–127, New York, NY, USA, 2007. ACM.
  • [3] M. Dyer, M. Jerrum, and E. Vigoda. Rapidly mixing markov chains for dismantleable constraint graphs. In Proceedings of a DIMACS/DIMATIA Workshop on Graphs, Morphisms and Statistical Physics, 2001.
  • [4] M. E. Dyer, A. M. Frieze, and M. Jerrum. On counting independent sets in sparse graphs. SIAM Jounal on Computing, 31(5):1527–1541, 2002.
  • [5] M. E. Dyer and C. S. Greenhill. On markov chains for independent sets. Journal of Algorithms, 35(1):17–49, 2000.
  • [6] A. Galanis, Q. Ge, D. Štefankovič, E. Vigoda, and L. Yang. Improved inapproximability results for counting independent sets in the hard-core model. In The 15th International Workshop on Randomization and Computation, RANDOM ’11, 2011.
  • [7] A. Galanis, D. Stefankovic, and E. Vigoda. Inapproximability of the partition function for the antiferromagnetic ising and hard-core models. Arxiv preprint arXiv:1203.2226, 2012.
  • [8] D. Gamarnik and D. Katz. Correlation decay and deterministic fptas for counting list-colorings of a graph. In Proceedings of the eighteenth annual ACM-SIAM symposium on Discrete algorithms, SODA ’07, pages 1245–1254, Philadelphia, PA, USA, 2007. Society for Industrial and Applied Mathematics.
  • [9] L. A. Goldberg and M. Jerrum. A polynomial-time algorithm for estimating the partition function of the ferromagnetic ising model on a regular matroid. In L. Aceto, M. Henzinger, and J. Sgall, editors, Proceedings of the 38th International Colloquium on Automata, Languages and Programming, volume 6755 of ICALP ’11, pages 521–532. Springer, 2011.
  • [10] L. A. Goldberg, M. Jerrum, and M. Paterson. The computational complexity of two-state spin systems. Random Struct. Algorithms, 23(2):133–154, 2003.
  • [11] M. Jerrum. A very simple algorithm for estimating the number of k-colorings of a low-degree graph. Random Struct. Algorithms, 7(2):157–166, 1995.
  • [12] M. Jerrum and A. Sinclair. Polynomial-time approximation algorithms for the ising model. SIAM Journal on Computing, 22(5):1087–1116, 1993.
  • [13] M. Jerrum, A. Sinclair, and E. Vigoda. A polynomial-time approximation algorithm for the permanent of a matrix with nonnegative entries. Journal of the ACM, 51:671–697, July 2004.
  • [14] F. Kelly. Stochastic models of computer communication systems. Journal of the Royal Statistical Society. Series B (Methodological), 47(3):379–395, 1985.
  • [15] L. Li, P. Lu, and Y. Yin. Approximate counting via correlation decay in spin systems. In Proceedings of the Twenty-Third Annual ACM-SIAM Symposium on Discrete Algorithms, SODA ’12, pages 922–940, Kyoto, Japan, 2012. SIAM.
  • [16] L. Li, P. Lu, and Y. Yin. Correlation Decay up to Uniqueness in Spin Systems. In Proceedings of the Twenty-Fourth Annual ACM-SIAM Symposium on Discrete Algorithms, SODA ’13, pages 67-84, New Orleans, Louisiana, USA, 2013. SIAM.
  • [17] M. Luby and E. Vigoda. Approximately counting up to four (extended abstract). In Proceedings of the twenty-ninth annual ACM symposium on Theory of computing, STOC ’97, pages 682–687, New York, NY, USA, 1997. ACM.
  • [18] F. Martinelli, A. Sinclair, and D. Weitz. Fast mixing for independent sets, colorings, and other models on trees. Random Struct. Algorithms, 31(2):134–172, 2007.
  • [19] E. Mossel, D. Weitz, and N. Wormald. On the hardness of sampling independent sets beyond the tree threshold. Probability Theory and Related Fields, 143:401–439, 2009.
  • [20] R. Restrepo, J. Shin, P. Tetali, E. Vigoda, and L. Yang. Improved mixing condition on the grid for counting and sampling independent sets. In Proceedings of the 52nd Annual IEEE Symposium on Foundations of Computer Science, FOCS ’11, 2011.
  • [21] A. Sinclair, P. Srivastava, and M. Thurley. Approximation algorithms for two-state anti-ferromagnetic spin systems on bounded degree graphs. In Proceedings of the Twenty-Third Annual ACM-SIAM Symposium on Discrete Algorithms, pages 941–953. SIAM, 2012.
  • [22] A. Sly. Uniqueness thresholds on trees versus graphs. The Annals of Applied Probability, 18(5):1897–1909, 2008.
  • [23] A. Sly. Computational transition at the uniqueness threshold. In Proceedings of the 2010 IEEE 51st Annual Symposium on Foundations of Computer Science, FOCS ’10, pages 287–296, Washington, DC, USA, 2010. IEEE Computer Society.
  • [24] A. Sly and N. Sun. The computational hardness of counting in two-spin models on d-regular graphs. Arxiv preprint arXiv:1203.2602, 2012.
  • [25] E. Vigoda. Improved bounds for sampling coloring. In Proceedings of the 40th Annual Symposium on Foundations of Computer Science, FOCS ’99, pages 51–, Washington, DC, USA, 1999. IEEE Computer Society.
  • [26] D. Weitz. Counting independent sets up to the tree threshold. In Proceedings of the thirty-eighth annual ACM symposium on Theory of computing, STOC ’06, pages 140–149, New York, NY, USA, 2006. ACM.

附录 A 唯一性阈值

以下引理将唯一性条件转换为其各种阈值形式。 引理(2除外)不用于主要结果的证明,而是用于主要结果的解释以及与先前结果的比较,这些结果大多以阈值形式表述。

Lemma 21.

(β,γ,λ) 是反铁磁性的。

  1. 1.

    (β,γ,λ) 最多为 2 个唯一。

  2. 2.

    如果γ1,那么唯一性不成立于无限d-所有足够大的正则树d

  3. 3.

    如果γ>1,则唯一性保持无限d-所有足够大的正则树 d

  4. 4.

    对于任何Δ(包括Δ=),存在一个临界阈值γc=γc(β,λ,Δ),当且仅当γ(γc,1β)时,(β,γ,λ)直到Δ是唯一的。特别是,γc(β,λ,)>1γc(β,λ,)=γc(β,λ,Δ) 对于某个有限的Δ.

  5. 5.

    如果β=0,对于任何Δ(包括Δ=),存在一个临界阈值λc=λc(γ,Δ)=min1<d<Δγd+1dd(d1)d+1,使得(β,γ,λ)直到Δ是唯一的,当且仅当λ(0,λc).

  6. 6.

    如果βγ>Δ2Δ,则对于任何外部字段λ(β,γ,λ) 最多为-Δ 唯一。

  7. 7.

    如果β>0,对于任何Δ(包括Δ=) βγΔ2Δ,λ的制度满足最多-0>Δ 唯一性如下。

    Δ¯1+βγ1βγ。对于任意整数Δ¯d<Δ,令x1(d)x2(d)为方程d(1βγ)x(βx+1)(x+γ)=1的两个正根>,更具体地说,

    x1(d)=1βγ+d(1βγ)(1βγ+d(1βγ))24βγ2β,
    x2(d)=1βγ+d(1βγ)+(1βγ+d(1βγ))24βγ2β.

    λi(d)xi(d)(xi(d)+γβxi(d)+1)d(其中i=1,2)被定义为整数Δ¯d<Δ .

    那么 (β,γ,λ) 最多 - Δ 唯一当且仅当 λ属于以下政权

    Δ¯d<Δ((0,λ1(d))(λ2(d),)). (8)

    如果此外γ1,则(β,γ,λ)可达-Δ 唯一当且仅当 λ(0,λc)(λ¯c,),其中

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

    特别是,当β=γ时,它成立λcλ¯c=1,因此(β,β,λ) 最多 - Δ 唯一当且仅当 |logλ|>logλ¯c。0>

  8. 8.

    只有当 γ>1 时,(β,γ,λ) 才可能是普遍唯一的。如果γ>1,则存在有限正常数λc𝗅𝗈𝗐𝖾𝗋=λc𝗅𝗈𝗐𝖾𝗋(β,γ)λc𝗎𝗉𝗉𝖾𝗋=λc𝗎𝗉𝗉𝖾𝗋(β,γ),其中λc𝗅𝗈𝗐𝖾𝗋λc𝗎𝗉𝗉𝖾𝗋、使得(β,γ,λ)λ<λc𝗅𝗈𝗐𝖾𝗋条件下是普遍唯一的,而(β,γ,λ)只有在λ<λc𝗎𝗉𝗉𝖾𝗋条件下是普遍唯一的。

Remark (SODA’13版本错误)

在论文的早期版本中,引理21的第7项有错误(这是SODA'13会议版本中的引理3.1[16] )。 在该版本中,唯一性机制 (8) 被无条件且错误地简化为双区间简单形式 (0,λc)(λ¯c,),实际上,当 β,γ1 然而,当γ>1时,存在这样的βγ<1,使得唯一性机制(8)可能变得比最多由两个区间组成更复杂。 例如,当 β=1122γ=100Δ=23 时,λ 的机制使得 (β,γ,λ) 为up-to-Δ 唯一,由三个间隔 (0,102.664)(104.339,106.967)(109.444,) 组成。

当前版本纠正了错误并更改了引理21的项目7和项目8的陈述。 我们进一步指出,该错误不会影响论文的主要结果,即从 up-to-Δ 唯一性到强空间混合和 FPTAS 的含义,因为引理 21不用于主要结果的证明(除了第2项),而仅用于解释各种阈值形式的唯一性机制以及与其他已知结果的比较。

引理21的证明。

fd(x)=λ(βxd+1xd+γ)dx^d=fd(x^d)fd(x)的正不动点。 然后

|fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ).

(β,γ,λ) 是反铁磁性的。 0βγγ>0βγ<1,即 β<1

  1. 1.

    很容易验证对于βγ<1,(βx+1)(x+γ)(1βγ)x>0对于任何正x>0 因此,当d=1时,我们有|f1(x^1)|=(1βγ)x^1(βx^1+1)(x^1+γ)<1,这意味着(β,γ,λ)是最多2个唯一的。

  2. 2.

    对于所有足够大的 d,它认为

    λβdexp(d(1βγ)d3) <d(1βγ)3β, and
    λexp(γdd(1βγ)3) >γd(1βγ)3.

    根据矛盾,假设|fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ)1 然后,

    1d(1βγ)x^d(βx^d+1)(x^d+γ)=d(1βγ)βx^d+γx^d+(1+βγ)d(1βγ)βx^d+γx^d+2.
    • 情况1:

      x^dγ 然后γx^d1 因此,

      1d(1βγ)βx^d+γx^d+2d(1βγ)βx^d+3,

      这意味着x^dd(1βγ)3β 然而,它认为

      x^d =λ(βx^d+1x^d+γ)dλ(βx^d+1x^d)dλ(β+βd(1βγ)3)d
      λβdexp(d(1βγ)d3)<d(1βγ)3β,

      矛盾。

    • 案例2:

      x^d<γ 然后βx^dβγ<1 因此,

      1d(1βγ)βx^d+γx^d+2d(1βγ)γx^d+3,

      这意味着x^dγd(1βγ)3 然而,它认为

      x^d =λ(βx^d+1x^d+γ)dλ(x^d+1)dλ(1+γd(1βγ)3)d
      λexp(γdd(1βγ)3)>γd(1βγ)3,

      矛盾。

  3. 3.

    γ>1 固定点x^d=λ(βx^d+1x^d+γ)dλγd,以及

    |fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ)dλγd,

    对于所有足够大的 d 来说,它严格小于 1。因此,对于所有足够大的d,无限d-正则树具有唯一性。

  4. 4.

    我们首先表明,存在一个临界阈值 γc=γc(β,λ,Δ),使得 (β,γ,λ) 最多 Δ 是唯一的,当且仅当 γ(γc,1β) 时。 足以证明,如果反铁磁 (β,γ,λ) 的唯一性可达 Δ,则 (β,γ,λ) 的唯一性可达 Δ 对于任何 γ>γβγ<1 都是唯一的。

    回想一下,x^dfd(x)=λ(βx+1x+γ)d 的正不动点。还让 x^d 表示 x=λ(βx+1x+γ)d 的正解。

    我们首先展示x^d<x^d 根据矛盾,假设x^dx^d 由于对于反铁磁(β,γ,λ),fd(x)是单调递减的,我们有

    x^d=λ(βx^d+1x^d+γ)dλ(βx^d+1x^d+γ)d>λ(βx^d+1x^d+γ)d=x^d,

    矛盾。

    自从

    λ(β+(1βγ)x^d+γ)d=x^d<x^d=λ(β+(1βγ)x^d+γ)d,

    我们有

    (1βγ)x^d+γ<(1βγ)x^d+γ.

    对于x^d<x^d,也成立

    x^dβx^d+1=1β+1x^d<1β+1x^d=x^dβx^d+1.

    将以上两个不等式相乘,我们有

    d(1βγ)x^d(βx^d+1)(x^d+γ)<d(1βγ)x^d(βx^d+1)(x^d+γ).

    请注意,这些是当参数为(β,γ,λ)(β,γ,λ)时在各自固定点处的绝对导数。 因此,如果 (β,γ,λ) 最多为 Δ 唯一,则 (β,γ,λ) 最多为 Δ 唯一。

    由于引理的 2 部分,如果 γ1|fd(x^d)|>1 对于所有足够大的 d,因此 γc(β,λ,)>1 并且由于引理的 3 部分,对于任何 γγc(β,λ,)>1,随着 d 增长到无穷大,|fd(x^d)| 任意接近 0,因此γc(β,λ,)=γc(β,λ,Δ)对于有限的Δ

  5. 5.

    β=0,|fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ)=dx^dx^d+γ时,唯一性条件|fd(x^d)|<1相当于x^d<γd1(这里我们假设d>1,因为由于引理的 1 部分,对于 d=1 ,唯一性始终成立)。 回想一下x^d=λ(1x^d+γ)d。那么 |fd(x^d)|<1 当且仅当

    λ=x^d(x^d+γ)d<γd+1dd(d1)d+1.

    λc=λc(γ,Δ)=min1<d<Δγd+1dd(d1)d+1 它认为当且仅当 λ(0,λc)(β,γ,λ) 最多 Δ 是唯一的。

  6. 6.

    我们注意到 |fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ)x^d 中不是单调的。 它在x^d=βγ处达到最大值。 因此,如果对于任何1d<Δ,d(1βγ)x(βx+1)(x+γ)<1对于x=βγ,则(β,γ,λ)最多是Δ唯一对于任何λ βγ>d1d+1 对于所有 1d<Δ 时,即当 βγ>Δ2Δ 时,此条件成立。

  7. 7.

    θ(d)1βγ+d(1βγ) 它对所有 dΔ¯=1+βγ1βγ 都适用 θ(d)2βγ 因此,x1(d),x2(d) 是明确定义的,可以表示为:

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

    回想一下 λi(d)xi(d)(xi(d)+γβxi(d)+1)d 代表 i=1,2

    下面的λ2(d)的单调性和无界性很容易验证。

    Claim 22

    λ2(d) d 中单调递增,并随着 d 的增长而趋于无穷大。

    证明。

    观察到对于0βγβγ<1dΔ¯,它认为ββγd1d+1<1x+γβx+1x 而且,

    x2(d)+γβx2(d)+1 = 1β(d1)(1βγ)+θ(d)24βγ(d+1)(1βγ)+θ(d)24βγ
    d+1d1(d1)(1βγ)(d+1)(1βγ)=1.

    很容易验证 x2(d)d 中不断增加,并且随着 d 的增长而无限制。 因此,我们可以得出结论,λ2(d)=x2(d)(x2(d)+γβx2(d)+1)dd 中不断增加,并且随着 d 的增长而无界。

    我们继续分析唯一性机制。 对于 1d<Δ¯,它保持 βγ>d1d+1,在这种情况下,由于此引理的第 6 部分,它始终保持 |fd(x^d)|<1

    回想一下βγΔ2Δ 仍然需要分析 |fd(x^d)| 中的那些 Δ¯d<Δ 对于这样的 d,我们有 βγd1d+1,并且由于 x1(d)x2(d) 是方程 d(1βγ)x(βx+1)(x+γ)=1 的根>,当且仅当 x^d<x1(d)x^d>x2(d) 时,才成立 |fd(x^d)|=d(1βγ)x^d(βx^d+1)(x^d+γ)<1 请注意,x(x+γβx+1)dx 中单调递增。 因此:

    x^d<x1(d) λ=x^d(x^d+γβx^d+1)d<x1(d)(x1(d)+γβx1(d)+1)d=λ1(d),
    x^d>x2(d) λ=x^d(x^d+γβx^d+1)d>x2(d)(x2(d)+γβx2(d)+1)d=λ2(d),

    这意味着对于任何Δ¯d<Δ,|fd(x^d)|<1当且仅当λ(0,λ1(d))(λ2(d),)

    因此,当且仅当 λ 属于该政权时,(β,γ,λ) 最多 Δ 是唯一的

    Δ¯d<Δ((0,λ1(d))(λ2(d),)),

    这正是唯一性机制 (8)。 回想起那个

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

    其中最后一个方程是由于 λ2(d) 的单调性造成的。

    然后我们证明,对于 γ1,唯一性机制 (8) 变为 (0,λc)(λ¯c,) 我们需要以下断言,其证明推迟到本节末尾。

    Claim 23

    对于Δ¯d0<d1,它总是成立

    (0,mind0d<d1λ1(d))(maxd0d<d1λ2(d),)d0d<d1((0,λ1(d))(λ2(d),));

    此外,如果 λ1(d+1)λ2(d) 对于所有整数 d0d<d11,则

    (0,mind0d<d1λ1(d))(maxd0d<d1λ2(d),)=d0d<d1((0,λ1(d))(λ2(d),)).
    Claim 24

    以下内容成立:

    • 如果 γ1,则 λ1(d+1)λ2(d) 对于所有整数 dΔ¯;

    • 如果γ>1,则存在有限的d0=d0(β,γ)Δ¯使得λ1(d+1)λ2(d) 对于所有dd0

    现在假设γ1 根据声明24,对于所有整数Δ¯d<Δ1,λ1(d+1)λ2(d) 然后由于声明 23,唯一性机制 (8) 变为

    Δ¯d<Δ((0,λ1(d))(λ2(d),))=(0,λc)(λ¯c,).

    特别是考虑反铁磁伊辛模型β=γ 很容易验证x1(d)x2(d)=γβ=1,因此

    λ1(d)λ2(d)=x1(d)x2(d)(x1(d)+γβx1(d)+1x2(d)+γβx2(d)+1)d=1.

    因此,由于权利要求22所保证的λ2(d)的单调性,λ1(d)=1/λ2(d)是单调递减的。 结果是 λc=λ1(Δ1)λ¯c=λ2(Δ1)λcλ¯c=λ1(Δ1)λ2(Δ1)=1 那么当且仅当 |logλ|>logλ¯c(β,β,λ) 最多是 Δ 是唯一的。

  8. 8.

    β=0 的情况在该引理的 5 部分中处理。 在这种情况下,由于引理的 5 部分,(0,γ,λ) 最多是 Δ 唯一,当且仅当 λ(0,λc(γ,Δ))其中λc(γ,Δ)=min1<d<Δγd+1dd(d1)d+1 特别是,如果 γ1,λc(γ,)=mind>1γd+1dd(d1)d+1 变为 0;如果 γ>1, 变为正且有限的 λc(γ)>0 因此,仅当 γ>1 时,(0,γ,λ) 才是普遍唯一的,并且当 γ>1 时,(0,γ,λ) 是普遍唯一的,当且仅当 λ(0,λc(γ)). 我们可以令 λc𝗅𝗈𝗐𝖾𝗋(γ)=λc𝗎𝗉𝗉𝖾𝗋(γ)=λc(γ) 和引理如下。

    对于 β>0 的情况,由于引理的 7 部分,对于所有 ΔΔ¯1+βγ1βγ,元组 (β,γ,λ) 是 up-to -Δ 唯一当且仅当 λ 属于以下机制:

    Δ¯d<Δ((0,λ1(d))(λ2(d),)).

    γ1 时,此状态变为 (0,λc(β,γ,Δ))(λ¯(β,γ,Δ),),其中 λc(β,γ,Δ)=minΔ¯d<Δλ1(d)λ¯c(β,γ,Δ)=maxΔ¯d<Δλ2(d)=λ2(Δ1)

    Claim 25

    以下内容成立:

    • 如果 γ1,则 λ1(d) 接近 0,如 d 增长到无穷大;

    • 如果 γ>1,则 λ1(d) 无界,如 d 增长到无穷大。

    证明。

    对于整数 dΔ¯,让 θ(d)1βγ+d(1βγ),并让

    λ^1(d)2γd+1θ(d)(d+1d1)d. (9)

    γ1时,很容易验证随着d增长到无穷大,λ^1(d)接近0。

    此外,λ1(d)的上限以λ^1(d)为界。 具体来说,对于dΔ¯

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

    由于当βγ<1x+γβx+1x中增加,我们有

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

    因此,对于所有 dΔ¯

    λ1(d) =x1(d)(x1(d)+γβx1(d)+1)d2γd+1θ(d)(d+1d1)d=λ^1(d). (10)

    因此,当γ1时,随着d增长到无穷大,λ1(d)趋近于0。

    另一方面,当γ>1时,对于dΔ¯

    λ1(d)=x1(d)(x1(d)+γβx1(d)+1)dx1(d)γd=2γd+1θ(d)+θ(d)24βγ,

    随着 d 的增长,这显然是无限的。

    λ1时,由于上述说法,随着d增长到无穷大,λ1(d)趋近于0,这意味着λc(β,γ,)=0;根据声明22,随着d增长到无穷大,λ2(d)也趋于无穷大,这意味着λ¯c(β,γ,)= 因此,只有当γ>1时,(β,γ,λ)才可以是全局唯一的。

    现在假设γ>1 根据声明 25,随着 d 的增长,λ1(d) 是无界的。 因此,λ1(d) 在有限的d=d(β,γ)Δ¯ 处达到其最小值。 显然,对于任何有限的 dΔ¯,λ1(d) 总是有限的且为正值。 因此,λc(β,γ,)=λc(β,γ,d(β,γ))是有限且正的。 λc𝗅𝗈𝗐𝖾𝗋(β,γ)=λc(β,γ,)>0 表示这个有限的正阈值。

    另一方面,根据权利要求 24,对于 γ>1,存在一个有限的 d0=d0(β,γ)Δ¯,使得 λ1(d+1)λ2(d) 适用于所有 dd0 λc𝗎𝗉𝗉𝖾𝗋(β,γ)=mindd0λ1(d) 显然,λc𝗎𝗉𝗉𝖾𝗋(β,γ)λc𝗅𝗈𝗐𝖾𝗋(β,γ)>0 并且 λc𝗎𝗉𝗉𝖾𝗋(β,γ)λ1(d0(β,γ)) 也是有限的。

    由于声明23以及声明22λ2(d)的单调性和无界性,-唯一性的制度,即dΔ¯((0,λ1(d))(λ2(d),)),夹层如下:

    (0,λc𝗅𝗈𝗐𝖾𝗋(β,γ)) dΔ¯((0,λ1(d))(λ2(d),))
    dd0((0,λ1(d))(λ2(d),))
    =(0,λc𝗎𝗉𝗉𝖾𝗋(β,γ))(λ2(),)
    =(0,λc𝗎𝗉𝗉𝖾𝗋(β,γ)).

    因此,当 λ<λc𝗅𝗈𝗐𝖾𝗋(β,γ) 且仅当 λ<λc𝗎𝗉𝗉𝖾𝗋(β,γ) 时,(β,γ,λ) 是普遍唯一的。

现在还有待证明索赔23和索赔24

索赔证明23

首先,很容易验证对于Δ¯d0<d1,它总是成立

(0,mind0d<d1λ1(d))(maxd0d<d1λ2(d),)d0d<d1((0,λ1(d))(λ2(d),)),

因为对于所有d0d<d1,它总是成立

(0,mind0d<d1λ1(d))(0,λ1(d)) and (maxd0d<d1λ2(d),)(λ2(d),).

接下来,我们证明如果对于所有整数 d0d<d11λ1(d+1)λ2(d),则两种机制相等。 这是通过归纳法证明的。

d0<d1d0+1时,权利要求中的相等关系是微不足道的,因为在这种情况下,最多有一个整数 d 满足 d0d<d1 那么等式两边就变成了

((0,λ1(d))(λ2(d),)).

这就奠定了归纳基础。

对于归纳假设,假设 d1>d0 等式成立:

d0d<d1((0,λ1(d))(λ2(d),))=(0,λ1<d1)(λ2<d1,), (11)

我们在哪里表示

λ1<d1mind0d<d1λ1(d) and λ2<d1maxd0d<d1λ2(d).

对于归纳步​​骤,我们希望在 d1 加 1 时建立相等关系。 因此另外,我们假设λ1(d+1)λ2(d)为唯一整数d[d11,d1),即

λ1(d1)λ2(d11). (12)

为了完成归纳,只需验证

d0d<d1+1((0,λ1(d))(λ2(d),))=(0,λ1<d1+1)(λ2<d1+1,),

假设归纳假设 (11)。 因此,只要验证一下就足够了

((0,λ1<d1)(λ2<d1,))((0,λ1(d1))(λ2(d1),))
= (0,λ1<d1+1)(λ2<d1+1,). (13)

请注意,d1[d1,d1+1) 中的唯一整数。 由于λ2(d)的单调性,我们有λ2(d1)λ2<d1=λ2(d11) 因此,(λ2(d1),)(λ2<d1,) 以及 (13) 与 λ1(d1)λ2<d1=λ2(d11) 一样长,这是由假设 (12) 保证的。

索赔证明24

回想一下 (9) 中定义的函数 λ^1(d):对于整数 dΔ¯,设 λ^1(d)2γd+1θ(d)(d+1d1)d,其中 θ(d)1βγ+d(1βγ) 到 (10),对于所有 dΔ¯,我们有 λ1(d)λ^1(d)

它还成立λ^1(d)λ2(d) 对于所有 dΔ¯,我们有 x2(d)=θ(d)+θ(d)24βγ2βθ(d)2β 由于x+γβx+1βγ<1时在x中单调递增,因此成立x1(d)+γβx1(d)+1θ(d)2β+γθ(d)2+1=1βd1d+1 同时,对于所有dΔ¯=1+βγ1βγ,βγ(d1d+1)2θ(d)24βγ成立。 所以,

λ2(d) =x2(d)(x2(d)+γβx2(d)+1)dθ(d)2βd+1(d1d+1)d2γd+1θ(d)(d+1d1)d=λ^1(d). (14)

假设γ1 也不难观察到 λ^1(d)dΔ¯>1 中正在减少。 事实上,如果γ1,那么2γd+1θ(d)>0dΔ¯中单调递减;我们可以验证函数h(d)(d+1d1)d对于dΔ¯>1的单调性,如下所示:

(lnh(d)) =2dd21+ln(d+1d1)0 as d,
(lnh(d))′′ =4d21>0,

这保证了 h(d)d 中单调递减,因此对于 dΔ¯,λ^1(d)=2γd+1θ(d)h(d)d 中单调递减,因为 2γd+1θ(d)h(d) 都是正数,并且 d 中是递减的。

因此,当γ1时,由于(10)、λ^1(d)的单调性,以及(14),对于所有dΔ¯,

λ1(d+1)λ^1(d+1)λ^1(d)λ2(d).

这就证明了当γ1时,λ1(d+1)λ2(d)对于所有dΔ¯

最后,当γ>1时,通过矛盾假设λ1(d+1)>λ2(d)对于某些dΔ¯ 然后我们对这种糟糕的 d 推导出有限上限 d0(β,γ)。到(10),我们得到λ1(d+1)λ^1(d+1)=2γd+2θ(d+1)h(d+1)2γd+2θ(d)h(d),其中最后一个不等式是由于θ(d)增加而h(d)减少;到 (14),我们得到 λ2(d)θ(d)2βd+11h(d) 因此 λ1(d+1)>λ2(d) 意味着

(1βγ)dθ(d)2h(d)2<4βγ2. (15)

回想一下,对于 dΔ¯=1+βγ1βγ,它始终保持 θ(d)24βγ 在不失一般性的情况下,我们只考虑足够大的 d2,因此 h(d)h(2)=9 因为 h(d) 正在减小。 因此不等式 (15) 意味着

(1βγ)d<81γ,

这为我们提供了 λ1(d+1)>λ2(d) 的必要条件 d<log1/βγ(81γ),其中 ddΔ¯d2

d0(β,γ)max{2,Δ¯,log1βγ(81γ)}.

然后它必须为所有 dd0(β,γ) 保留该 λ1(d+1)λ2(d)

附录 B 引理证明 11(中值定理方法)

我们证明引理11 这些符号在 3 节中定义。

  1. 1.

    根据中值定理,存在 R~[Rv,Rv+δv] 使得

    ϵv =φ(Rv+δv)φ(Rv)=φ(R~)δv=Φ(R~)δv.
  2. 2.

    假设每个顶点最多有 k 个子节点。 然后由于假设 k 有界或 γ>1

    如果 vΛ,则 δvRv+δv=f(Rv1,,Rvd)f(0,,0)=λγd,其中 0dk 如果γ>1,则Rv+δvλ=O(1);如果k是有限的,则Rv+δvmax{λγk,λ}=O(1)

    根据中值定理,存在Ri~[Rvi,Rvi+δvi],1id,使得

    ϵv =φ(f(Rv1,,Rvd))φ(f(Rv1+δv1,,Rvd+δvd))
    =φ(f(R1~,,Rd~))(δv1,,δvd)
    =(1βγ)(f(R1~,,Rd~))12(βf(R1~,,Rd~)+1)12(f(R1~,,Rd~)+γ)12i=1dδvi(βRi~+1)(Ri~+γ)
    f(0,,0)γi=1dδviγ.

    如果viS代表所有vi,那么由于上述论点,δviλγdi代表免费的vi,其中0dikvi的孩子数量;简单地 δvi=0 表示固定的 vi 因此,ϵvλγd+1i=1dλγdi+1λ32dγd+32max{1γk,1} 如果γ>1,那么很容易看出dγd+32以常数为界,因此成立ϵv=O(1);如果k是有界的,那么dk也是有界的,因此显然是ϵv=O(1)

  3. 3.

    然后我们分析ϵv的逐步收缩。 定义 yv=φ(Rv) 以及相应的 yvi=φ(Rvi)1id 然后是yv+ϵv=φ(Rv+δv)yvi+ϵvi=φ(Rvi+δvi)1id 我们有

    yv = φ(f(φ1(yv1+ϵv1),,φ1(yvd+ϵvd))),
    yv+ϵv = φ(f(φ1(yv1),,φ1(yvd))).

    应用中值定理。 存在yi~[yvi,yvi+ϵvi]和对应的Ri~[Rvi,Rvi+δvi]满足yi~=φ(Ri~),1id,使得

    ϵv =φ(f(φ1(yv1),,φ1(yvd)))φ(f(Φ1(yv1+ϵv1),,φ1(yvd+ϵvd)))
    =φ(f(φ1(y1~),,φ1(yd~)))(ϵv1,,ϵvd)
    =Φ(f(R1~,,Rd~))i=1dfRi1Φ(Ri~)ϵi
    =(1βγ)(λi=1dβRi~+1Ri~+γ)12(βλi=1dβRi~+1Ri~+γ+1)12(λi=1dβRi~+1Ri~+γ+γ)12i=1dϵviRi~12(βRi~+1)12(Ri~+γ)12
    max1id{ϵvi}(1βγ)(λi=1dβRi~+1Ri~+γ)12(βλi=1dβRi~+1Ri~+γ+1)12(λi=1dβRi~+1Ri~+γ+γ)12i=1dRi~12(βRi~+1)12(Ri~+γ)12
    =α(d;R~1,,R~d)max1id{ϵvi}.

附录C启发式地找到一个好的势函数

也许我们证明中最神秘的一步是势函数Φ(x)=1x(βx+1)(x+γ)的选择。 正如在许多应用势分析的情况下一样,没有用于搜索合适的势函数的标准例程。 另一方面,我们不太可能在没有任何提示的情况下猜出这样一个相当复杂的公式。 在这里,我们提出了一种启发式方法,引导我们发现良好的潜在函数。 这部分并不严格,对于我们结果的合理性来说,逻辑上也是不必要的。 然而,这种启发式方法足够通用且有趣,并且可能会在其他场景中找到应用,因​​此值得在这里进行阐述。

启发式方法包括三个步骤:

  1. 1.

    求势函数的必要条件,即与势函数在某一点相关的方程。

  2. 2.

    (启发式步骤)通过假设方程适用于变量的整个范围来增强条件,这给了我们一个微分方程。

  3. 3.

    求解微分方程并得到势函数。

我们首先假设系统处于唯一性的边界。 这意味着 d 是临界度,f(x^)=1 是正不动点,其中 f(x)=λ(βx+1x+γ)dx^=f(x^) 是正不动点。

我们有以下身份:

x^=λ(βx^+1x^+γ)d and 1=λd(1βγ)(βx^+1)d1(x^+γ)d+1,

这一起意味着另一个身份

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

我们的目标是找到一个势函数,使 α(x)=|f(x)|Φ(f(x))Φ(x)1 适用于所有 x 另一方面,我们有α(x^)=|f(x^)|Φ(f(x^))Φ(x^)=1Φ(x^)Φ(x^)=1 所以α(x)x=x^时达到最大值。 作为一个可微函数,它满足α(x^)=0,即

[f(x)Φ(f(x))Φ(x)]|x=x^=0.

我们有以下计算:

[f(x)Φ(f(x))Φ(x)]x=x^=0
[f(x)Φ(f(x))]Φ(x)|x=x^=f(x)Φ(f(x))Φ(x)|x=x^
[f′′(x^)Φ(f(x^))+f(x^)Φ(f(x^))f(x^)]Φ(x^)=f(x^)Φ(f(x^))Φ(x^)
f′′(x^)Φ(x^)+Φ(x^)=Φ(x^)
f′′(x^)2=Φ(x^)Φ(x^)=(ln(Φ(x^))),

我们使用 x^=f(x^)f(x^)=1 的事实。

f(x)求二阶导数,我们有

f′′(x)=λd(βγ1)(βx+1)d2(x+γ)d+2((d1)β(x+γ)(d+1)(βx+1)).

该方程已经给出了 Φ 的恒等式。 但它太复杂了。 我们可以在x=x^点进一步简化上面的公式。

f′′(x^) = λd(βγ1)(βx^+1)d2(x^+γ)d+2((d1)β(x^+γ)(d+1)(βx^+1))
= 1(βx^+1)(x^+γ)((d+1)(βx^+1)(d1)β(x^+γ))
= d+1x^+γ(d1)ββx^+1
= d(1βγ)(βx^+1)(x^+γ)+1x^+γ+ββx^+1
= 1x^+1x^+γ+ββx^+1.

所以我们有

(ln(Φ(x^)))=f′′(x^)2=12(1x^+1x^+γ+ββx^+1). (17)

我们应用启发式方法并假设 (17) 对于所有 x 均成立。 这给了我们一个微分方程

(ln(Φ(x)))=12(1x+1x+γ+ββx+1).

该微分方程的解为

ln(Φ(x))=12ln(x(x+γ)(βx+1))+C1,

这给出了

Φ(x)=C2x(βx+1)(x+γ),

其中 C1,C2 是一些常量。 这给出了我们在论文中使用的势函数。

对于固定点x^,还有其他等式成立。 选择x^的哪个方程启发式扩展到所有x可能会影响我们获得的势函数。 例如,(16) 意味着

1dx^=1βγ(βx^+1)(x^+γ)=1x^+γββx^+1,

因此我们可以将 f′′(x^) 重写为

f′′(x^)=1x^+1x^+γ+ββx^+1=d+1dx^+2ββx^+1.

这给出了

(ln(Φ(x^)))=d+12dx^ββx^+1.

我们还可以将此方程视为变量 x 的微分方程,求解它可以得到

Φ(x)=cxd+12d(βx+1),

这是[15]中使用的势函数。