自旋系统中相关性衰减至唯一性
摘要
我们给出了在一般图上具有强空间混合的二态反铁磁自旋系统的完整表征。 我们证明,当且仅当系统在无限正则树上具有唯一的吉布斯测度时,二态反铁磁自旋系统在所有最大度数最多 的图上具有强空间混合到 ,其中 可以是有界的,也可以是无界的。 因此,当在度数高达 的无限正则树上满足唯一性条件时,在最大度数最多 的图上,存在一个二态反铁磁自旋系统的配分函数 FPTAS 。 特别是,如果所有无限正则树都满足唯一性,则任意图都存在 FPTAS。 作为特殊情况,这涵盖了一般结构图上两态反铁磁系统的所有先前算法结果。
结合 Jerrum-Sinclair 和 Goldberg-Jerrum-Paterson 的二态铁磁自旋系统的 FPRAS,以及 Sly-Sun 和独立于 Galanis-S̆tefankovic̆-Vigoda 的最新硬度结果,这给出了完整的分类,除了相变边界,所有二态自旋系统的近似性,在有度限制的图族或所有图族上。
1简介
自旋系统在统计物理、应用概率和计算机科学领域得到了深入研究,作为捕捉局部相互作用和约束如何影响粒子系统宏观特性的本质的通用框架。 系统通常用图来描述,每个顶点处于固定数量的状态之一(称为自旋),边指定系统的邻域关系。
令 为图表, 为自旋状态数。 系统的配置是可能的分配之一。 每个配置都有一个能量作为所有边和顶点的总和,使得边的贡献由对称函数确定自旋状态和,顶点的贡献由其自旋状态的函数确定。 配置的权重是,其中是温度。 我们专注于两种状态的自旋系统。 在归一化之前,二态自旋系统完全由三个参数 捕获,其中 和 确定边缘贡献的对称函数,,也称为外部场,决定顶点贡献的函数。 吉布斯测度是所有配置的自然概率分布,配置 的概率为 ,其中归一化因子 称为配分函数。 配分函数编码了有关自旋系统宏观行为的丰富信息。 然而,对于几乎所有重要的设置来说,计算配分函数的精确值是非常困难的。
自旋系统最重要的特性之一是相关性衰减,它表示两个顶点的边缘吉布斯分布之间的相关性随着两个顶点之间的距离而快速衰减。 此属性也称为弱空间混合 [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。
对于任何有限或,二态反铁磁自旋系统在最大度图上最多具有强空间混合当且仅如果系统在所有的无限-正则树上表现出唯一性。
由于Weitz的自回避游走树构造[26],有度有界图上的二态自旋系统的强空间混合立即意味着配分函数的FPTAS。 事实上,我们展示了之前的工作[15]中引入的更强大的相关性衰减概念,即计算有效的相关性衰减,它不仅为有度有界图提供 FPTAS,还为任意图提供 FPTAS当满足相应的唯一性条件时,具有无界度。
Theorem 2.
对于任意有限的或,在最多的图上,二态反铁磁自旋系统的配分函数存在一个FPTAS > 如果对于所有,系统参数位于无限-正则树的唯一性区域的内部。
在上述两个定理中,情况代表无界度图。
由于 Sly 和 Sun [24] 针对一般二态自旋系统的最新硬度结果以及 Galanis、S̆tefankovic̆ 和 Vigoda [7] 针对一般二态自旋系统的独立结果不太通用的设置,违反唯一性条件意味着配分函数的不可近似性。
Theorem 3 (由于[24]和[7])。
对于任何有限的或,除非,不存在用于二态反铁磁自旋系统在图上的配分函数的FPRAS如果对于某些,系统参数位于无限-正则树的非唯一性区域的内部,则最大度至多。
[24] 中的原始定理适用于 - 具有固定 的正则图,这立即意味着有度有界图或任意图的硬度。 而 [24] 和 [7] 中的硬度条件要求在具有 的 不规则树上的非唯一性。 但无限2-正则树(即无限路径)上的唯一性对于任何二态反铁磁自旋系统始终成立,因此定理3中的条件就足够了。
对于最大次数为 2 或更小的图,可以在多项式时间内精确计算配分函数。 因此,定理2和定理3与两态铁磁自旋系统[12, 10]的FPRAS一起给出了几乎完整的(除了在相变边界)所有二态自旋系统的配分函数的近似性分类,在有度有界图族或所有图族上。
1.1.1 规律性和单调性
在统计物理学中,相关性衰减通常在正则图甚至结构对称图(例如 Bethe 晶格)上研究。 尽管对于仅考虑正则图的硬度结果将加强下界,但从算法的角度来看,考虑通用图上的自旋系统更为普遍。 [26]中硬核模型的近似算法和反铁磁Ising模型的[21]中的近似算法均适用于具有有限最大度的一般图。 正如[21, 24]中所讨论的,对于具有固定的规则图上的一般双态自旋系统,在参数平移之前,硬核模型和伊辛模型是完整的。然而,这并不包括最一般的情况,即一般图上有界或无界度的一般双态自旋系统。 造成这种差距的根本原因是一般自旋系统的非单调性。
充分研究的硬核和反铁磁 Ising 模型(以及所有具有 的双态自旋系统)都是单调自旋系统,从某种意义上说,无限 上的唯一性 -正则树意味着所有具有较小度数的无限正则树的唯一性。 这种单调性在一般的二态自旋系统中并不一定成立。
在[26]中,Weitz 在硬核模型中建立了以下含义:
Claim ([26]中的定理2.3)。
-正则树上的强空间混合意味着最大程度至多的所有图上的强空间混合。
在[26]中,Weitz在没有证明的情况下也指出这一推论对于所有二态自旋系统都成立(事实上,这在[21]中的反铁磁伊辛模型中得到了严格证明>)。 有了这一点,为有度有界图上的二态自旋系统设计近似算法就可以简化为验证 正则树上的强空间混合。 这已被接受为关于二态自旋系统的事实,并已成为此类系统的近似算法的构建块(参见[21]中的定理2.8以及[22中的讨论) ,23])。 [22]中还提出了该主张是否适用于一般多态自旋系统的猜想。
作为我们分析的副产品(参见第 5 节),我们发现 正则树和最大度图上的强空间混合之间的这种可信的含义大多数 仅适用于单调自旋系统(包括硬核和反铁磁 Ising 模型)。 对于一般的二态自旋系统,唯一性以及所有有度有界图之间的强空间混合的最坏情况确实是常规树,但可能不再是最高度的树。 这种复杂性的一个好处是,较高的度数可能会产生更快的相关性衰减,从而使 FPTAS 能够用于具有无界度数的图。
这些新现象表明,一般的二态自旋系统比经过充分研究的单调自旋系统(例如硬核和反铁磁伊辛模型)具有更丰富的结构。 前一种方法通过正则树上的强空间混合,在一般图上的单调自旋系统和正则图上的一般自旋系统上取得了成功,但在处理一般图上的一般自旋系统时遇到了障碍。 我们通过任意树而不是-常规树上的强空间混合,给出了一般二态自旋系统中相关性衰减的统一方法。 这种方法是在我们之前处理无界度图的工作[15]中发起的。 在本文中,我们设计了一种统一的基于势的分析,它适应任意树的不规则性和一般二态自旋系统的非单调性,并给出所有二态反铁磁自旋系统的紧密相关衰减结果关于图的有度族和所有图的族。
1.1.2 主要结果的意义
通过解决唯一性条件,我们可以以各种阈值形式重申我们的主要结果。 定理2作为特殊情况涵盖了一般结构图上二态反铁磁自旋系统的所有先前算法结果,并清除了以前未发现的情况。
在互动方面。
我们可以固定外场并讨论的可处理区域。 由于和的作用是对称的,我们可以进一步固定其中一个并讨论另一个的可处理范围。 [10, 15]中使用了该公式。
我们的主要结果可以重述如下:对于任何,对于无限正则树上的唯一性有一个临界阈值,直到程度,这样如果存在最大度数最多为的图的FPTAS;特别地,是绝对常数,因此如果则任意图都存在FPTAS。
作为特殊情况,这涵盖了[10]中关于反铁磁系统的所有算法结果以及[15]中的所有结果,扩展了这些先前作品中的可处理区域,并且也考虑有度有界图。
在外场方面。
在硬核伊辛模型和反铁磁伊辛模型研究的推动下,我们可以固定并讨论外场的可处理范围。
由于和的对称作用,我们可以假设。 我们的主要结果可以在特定设置下重述如下:
-
•
硬约束(当 时):对于任何 , 是无限正则树上唯一性的关键阈值,最高可达 度> 这样,如果 存在最大度数最多 的 FPTAS。
对于,临界阈值等于。 在无界度的无限正则树上不存在满足唯一性的外场。 这与没有度数限制的硬核模型的硬度结果[4, 23]一致。 一个特别有趣的情况是当时,在这种情况下,模型是具有逸度的硬核模型,并且是临界阈值。 这涵盖了[26]的结果。
对于 ,除了度受限图的结果外,还存在一个绝对正常数 ,它可以降低所有 的 的边界,这样,如果 ,就存在一个适用于任意图的 FPTAS。
-
•
软约束(当 时):如果 的唯一性条件对于任何外部域 都成立,并且对于最大度最多为 的图,总是存在一个 FPTAS。
现在假设。 令为满足的最小。 对于每个整数 ,对于无限 正则树上的唯一性,存在两个阈值 和 ,这样如果或对于所有整数,最大度数最多为的图存在一个FPTAS。
此外,如果,则上述唯一性条件可以简化为或,其中和。 特别是,对于反铁磁伊辛模型,我们有,并且当时唯一性成立。 这涵盖了[21]中反铁伊辛模型的结果。
对于无界最大度,如果,则在所有度的无限正则树上不存在满足唯一性的,这与硬度结果一致[10];如果,则当或对于所有时,所有无限正则树的唯一性成立,在这种情况下,存在用于任意图的 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 定义和基本知识
两态自旋系统由图描述。 系统的配置是顶点状态的可能分配之一。 我们还使用蓝色和绿色两种颜色来表示这两种状态。 配置的权重可以描述为各个边和顶点的贡献的乘积。 让 和 。 配置的权重由下式给出
吉布斯测度是定义的所有配置的概率分布。 归一化因子称为配分函数。
我们可以将 边和绿色顶点的贡献标准化为 1。 因此 代表某些 , 代表某些 。 由于蓝色和绿色的角色是对称的,我们可以假设而不失一般性。 三个参数与和完全指定了一个二态自旋系统。 具有的二态自旋系统是伊辛模型,具有或对称的二态自旋系统是硬核模型。
如果相邻顶点有利于不同的自旋,即如果,则二态自旋系统称为反铁磁。 不失一般性,我们重点关注的情况。
Definition 4。
如果、、和,则为反铁磁。
通过 和 的对称性以及情况 的平凡性,对于所有非平凡的二态反铁磁系统,该定义是完整的
2.1相关性衰减
吉布斯测度定义了每个顶点的状态边缘分布。 让 表示顶点 被着色为蓝色的概率。 令 为 的配置。 我们将顶点 固定顶点称为自由顶点。 我们使用表示在的配置固定为的条件下被着色为蓝色的边际概率。
Definition 5。
如果对于图族中的任意图 、任意 和 来说,图族上的自旋系统具有 强空间混合性、
其中 是 和 不同的子集, 是距 的最短距离到 中的任意顶点。
弱空间混合可以通过测量相对于而不是的衰减来定义。 空间混合特性在统计物理学中也称为相关衰减。
令 为一棵以 为根的树。将 定义为在施加条件 时根 为蓝色和绿色的概率之间的比率(约定为 当 时)。 假设 有 个子级,并且 是以第 个子级为根的子树。 由于子树的独立性,我们有一个简单的递归来计算:
(1) |
令 为图表。 类似地定义。 与树的情况相反,由于循环引起的依赖性,对于一般图来说,计算并不容易。 在[26]中,引入了一种称为自回避游走(SAW)树的结构,它将一般图中边缘分布的计算减少到树中的计算。 具体来说,给定一个图 和一个顶点 。 SAW 树 是一棵以 为根的树,具有新的顶点集 (它有效地枚举了源自 的所有路径) 并且可能包括固定叶子)。 此外,任何顶点集都映射到相应的,任何配置都映射到。 如果不会引起歧义,我们会滥用这种表示法并编写 和 。 给定图形 、 和 ,令 为 中距 到 中的任意顶点。
Theorem 6 (Weitz定理3.1[26])。
令 为图形,、 为 上的配置,以及 。 让。 其认为的最大阶数等于、和的最大阶数。 此外, 中 的任何邻域都可以在与邻域大小成正比的时间内构建。
2.2 唯一性条件
我们考虑无限 不规则树 上吉布斯度量的唯一性,其中由于 的对称结构,递归由 给出。
设为的正不动点,即。 已知[14, 18],上的二态反铁磁自旋系统在处发生相变,并且当。 这引发了以下定义。
Definition 7。
令 是反铁磁性的。 设为函数的正不动点。我们说 是 up-to- 唯一的,如果对于所有整数 ,
特别是,如果是最多唯一,那么它就是普遍唯一。
达到的唯一性相当于系统在无限正则树上的弱空间混合达到度。 唯一性条件可以用各种阈值形式来描述,在附录A中给出。
唯一性是通过要求 来定义的。 以下引理指出,只要唯一性条件成立, 就受到绝对常数的限制。
Lemma 8.
令 是反铁磁性的。 如果 最多为 唯一,则存在一个绝对常数 ,它仅取决于 、、 和 ,例如 适用于所有 。
3 一般图上的强空间混合
本节我们证明定理1。 唯一性条件的必要性是微不足道的,因为一般图上的强空间混合意味着常规树上的弱空间混合。 接下来就剩下证明下面的定理了。
Theorem 9。
令 是反铁磁性的。 对于任何有限的或,如果最多是唯一,则参数的二态自旋系统在最多的最大度图上具有强空间混合。
我们的方法是证明任意树上最大度至多的强空间混合,这通过自回避游走树构造隐含了定理。 请注意,与 [26] 和 [21] 不同,我们分析任意树而不是常规树的衰减。 这是因为对于一般的二态反铁磁自旋系统,最大程度的所有图之间强空间混合的最坏情况可能不再是-正则树。 我们将在5节中详细解释这一点。
给定最大度最多为 的任意图 、 上的任意配置 和任意 ,固定任意顶点 、根据定理 6,可以构造一棵自避走树 ,使得 的最大度受 、 和 的约束。 回想一下 因此 。 对于任何 和 ,
这引发了以下定义。
Definition 10。
令 为以顶点 为根的树, 为 上的配置, 为顶点集。 将 和 定义为:对于所有 , 只在 上与 不同。
通过为 构造这样的 和 并证明 就足以证明定理 9 >。
令 为一棵以 为根的树,它有 个子节点 , 为子树以 为根。 它认为
下限和上限和可以递归地构造如下。 基本情况是: (1) ,在这种情况下和; (2) ,即固定为所有中相同的颜色,此时和(或)如果固定为蓝色(或绿色)。 对于,由于对于反铁磁是单调递减的,
其中 和 是 、 相应的下限和上限。 特别是,当时,按照惯例,空乘积等于1,因此,这与是没有子节点的自由顶点的情况一致。
通过的单调性,很容易检验上面构造的和是否满足定义10的要求。 我们的目标是证明,当唯一性成立时, 的递归深度呈指数衰减。 一种简单的方法是尝试证明 在每个递归步骤中以恒定速率收缩。 但这不一定是真的才能保证指数衰减。 事实上,在某些情况下,误差不会单步衰减,而是长期衰减。 为了克服这个问题,我们使用潜力 来摊销收缩并表明以恒定速率收缩。 我们选择势函数为
我们正在分析任意树上具有不规则程度的腐烂。 为了适应这种不规则性,势函数不能以作为输入,而只能携带有关当前顶点的分布信息,但它必须能够为步骤提供正确的补偿 -在 的任何状态下以及对于满足足够唯一性的所有自旋系统而言,明智的衰减。 附录C中讨论了引导我们得到这个良好势函数的启发式过程。
令为满足的单调函数。 我们定义
相应地,、。
我们定义一个函数如下:
以下引理是通过应用中值定理获得的。 之前在 [15, 20] 中使用了类似的例程。 引理的证明在附录B中。
Lemma 11.
以下内容适用于。
-
1.
(与相关) 对于某些.
-
2.
(绝对界限) 假设 或 的最大度以常数为界,如果 则 ,如果 对于所有 则 。
-
3.
(逐步收缩)存在,,使得
To prove the strong spatial mixing, we first relate to by Item 1 of Lemma 11, and then apply induction on the depth in , with Item 2 of Lemma 11 as basis and Item 3 of Lemma 11 as induction step. 然后我们需要限制收缩率。 请注意, 对于 ,因此它无条件地适用于所有 、,即
(2) | ||||
(3) |
凭借唯一性,可以证明以下更紧的收缩界。
Lemma 12.
令 是反铁磁性的。 如果 最多为 唯一,则存在一个常量 ,使得对于任何整数 和任何 、,则成立。
定理证明9。
令 为 ,其最大度数至多为 。 那么 的最大度最多为 ,因此根 在 中最多有 个子顶点,而 中的每个其他顶点都有少于 个子顶点。 我们为 中的每个子树递归构造 、 和 。
请注意, 和 代表所有其他 。 假设 最多为 是唯一的。 如果 是有界的,那么根据定理 12,存在一个常数 ,使得 对于 ,并且 由于 (2);如果 ,那么根据定理 12,对于所有 。 在这两种情况下我们都有。
由于引理 11 的第 1 项, 对于某些 而言。 然后我们绑定和。 由于引理 21 的项 2, 最多 唯一这一事实意味着 是有界的或 。 请注意, 必须是自由的,否则该定理的证明很简单,并且 的子级都不在 中,因为 。 因此,根据引理11、和的项2,这意味着。
总之,如果最多是唯一,则存在一个常量,使得。 正如本节开头所讨论的,这证明了定理 9。
引理12的证明。
本节的其余部分致力于引理 12 的证明。 鉴于 最多 是唯一的,因此存在绝对常数 ,使得 对于任何 和 。
我们定义 的对称版本:
以下引理表明,通过使用 Cauchy-Schwarz 不等式以及算术和几何平均值,对称情况主导 的最大值。
Lemma 13.
令 是反铁磁性的。 那么对于任何整数和任何,都存在一个使得。
证明。
让。 然后是和。 用表达:
由于柯西-施瓦茨不等式,
由于算术平均值和几何平均值不等式,
让。 然后结合上面的计算,
假设 是 。 然后 并用 替换 ,我们有
∎
Lemma 14.
令 是反铁磁性的。 如果 最多为 唯一,则存在一个常量 ,这样对于任何整数 ,它都满足 对于所有 。
证明。
将 固定为任意正整数。 我们描述 达到最大值时 的值。 我们表示。对 对 求导,我们得到
其中,其导数为
由于 的范围超过 ,因此函数 在 中严格递减,范围从 到 ,函数在中严格递增,并且有一个有界范围。 因此,方程
(4) |
有唯一解,记为。 此外,它认为
(5) |
对于也是如此。 因此,对于任何固定的,在时达到最大值。
设为的正不动点,即。 然后我们声称,如果 最多是 唯一,则对于任何整数 , 是唯一的,要看到这个声明足以暗示引理,请注意,替换 后,我们有 。 并且由于引理 8,如果 最多是 唯一,则存在一个常量 使得 对于所有整数 。
然后我们证明这一主张。 假设 最多为 是唯一的且 。 那么足以证明,如果 ,则 正在减少,如果 ,则 正在增加。
情况 1:。 由于 (5),。 注意
由于唯一性,,因此。 与 结合,我们有 。 由于函数是单调递减的,并且是其不动点,因此。 由于满足(4),和必须同时为正或负,因此也成立。 那么和都是正数,并且在中单调递减。 因此,。
情况 2:。 根据与上面相同的论证,它成立、和。 因此,和都是负数,并且在中单调递减,因此它们的乘积是正数,并且在中递增。 因此,。 ∎
普通树上的强烈空间混合。
作为我们分析的副产品,我们证明了规则树的强空间混合定理。 当图 本身是一棵常规树时。 所有顶点(根除外)都具有相同的数量。 并且证明中出现的所有(摘录根的一个)都相同且等于该元数。 那么度数达到的所有无限正则树上唯一性成立的条件就可以用无限正则树上的唯一性来代替。
Theorem 15。
对于二态反铁磁自旋系统,在任何无限-正则树上,唯一性意味着强烈的空间混合。
4 算法含义
本节我们证明定理2。 也就是说,如果反铁磁性 到 是唯一的,那么对于最大度最多为 的任何图 ,存在分割函数 的 FPTAS,特别是普遍唯一性意味着任意图 的 FPTAS。
众所周知, 可以通过以下标准过程从 计算得出。 让 枚举 中的顶点。对于 ,令 为固定第一个 顶点 的配置,如下所示: for 和被固定,使得。 特别是, 是 的配置。对于 的吉布斯度量, 以及 都成立,因此 ,其中权重 可以针对任何特定的 以 的多项式时间精确计算。请注意, 等于 或 。 因此,如果可以在和中时间多项式的加性误差内近似,则配置 可以有效地构造,使得所有 都远离 0,因此乘积 可以在时间上在 因子内近似 和 中的多项式,这意味着
有界度图。
任意图形。
令 为任意图, 为任意顶点。 让。 我们使用[15]中引入的计算高效相关性衰减方法来处理无界度。 直观上,使用这种方法,我们观察到精细度量中的相关性衰减,而不是图距离,这样在这个新度量中,中等大小的邻域足以保证所需的相关性衰减。
类似地,我们使用3节中描述的递归过程来计算的上限和下限,但这次终止条件依赖于定义如下的新深度。
Definition 16。
令 为有根树, 为常量。 对于 中的任意顶点 ,定义 的基于 的深度,表示为 :如果 是根节点,则定义为 ;如果 是 的 子节点之一,则定义为 。
让 被修复。 用 表示具有基于 的深度 的所有顶点及其在 中的子代和孙子的集合。归纳起来可以验证。 当当前不在内时,递归将用于估计,直到不再在内,在这种情况下,将使用微分边界。
让 的定义如 3 节中所述。 重复应用定理 11 的 3 项,在不违背一般性的前提下,我们在 中有一条路径 从根 到 的 和 、使得为,其中为和、的子节点数。 克服无界度引起的爆炸的关键是观察收缩随着度的增加而急剧减小。
Lemma 17.
令 是反铁磁性的。 如果满足通用唯一性,则存在常量和,使得对于任何整数和任何,,它成立。
证明。
异质自旋系统。
我们在上一节和本节中的分析实际上适用于异构自旋系统,该系统允许每个顶点 具有不同的恒定外场 。
Theorem 18.
对于每个顶点处具有参数、和外场的二态反铁磁异质自旋系统,对于任何有限的或,如果对于所有,最多是唯一则自旋系统具有强空间混合性,并且在最大度图上具有至多的FPTAS。
5 一般二态自旋系统的非单调性
在本节中,我们证明以下定理。
Theorem 19.
存在二态反铁磁自旋系统,其在无限-正则树上表现出强空间混合,但在至多的所有最大度图上不表现出强空间混合。
在开创性的工作[26]中,Weitz证明了对于硬核模型,无限正则树上的强空间混合意味着最大度图上的强空间混合最([26]中的定理2.3)。 他进一步指出,同样的含义适用于所有两种状态的自旋系统。 那是,
对于任何二态自旋系统,无限正则树上的强空间混合意味着最大度数最多的图上的强空间混合。
这一主张对于当前对二态自旋系统相关衰减的理解以及为此类系统设计 FPTAS 发挥了重要作用。 [21] 中引用了该主张的算法形式作为所有二态自旋系统的定理([21] 中的定理 2.8),并被证明用于反-铁磁伊辛模型。 在[22]中提出了这个主张是否适用于多状态自旋系统的猜想。
在这里我们澄清这一说法仅适用于有限设置下的二态自旋系统,但并不适用于所有一般的二态自旋系统。 这反驳了[22]中的猜想,并表明了人们普遍认为-正则树代表了所有最大度图之间强空间混合的最坏情况,最多 不能推广到一般的二态自旋系统。 我们首先描述一个区域,证明该说法是正确的。
Lemma 20.
对于,无限-正则树上的强空间混合意味着最大度数最多为的树上的强空间混合。
证明。
给定最大度数为 的有根树,对于具有 的 子代的每个顶点,我们可以附加 虚拟子代具有固定的自旋态(分布)。 这是[26]和[21]中使用的方法:对于硬核模型,虚拟子项被固定为空闲,对于反铁磁Ising模型,虚拟子项被固定为空闲子级在自旋态上呈均匀分布。 在这两种情况下,虚拟孩子对其父母没有影响。 一般来说,我们将每个虚拟子节点的分布固定为 满足
当时,系统在中有解。 第 个子级的比率 、 和虚拟子级 的 的比率,父级的比率由递归给出
与原始数量相同。 ∎
由于自回避游走树构造(定理6),它认为对于,无限-正则树上的强空间混合意味着强空间混合和最大度至多图上的FPTAS。
对于,自旋系统表现出以下单调性质:无限正则树上的唯一性意味着所有更小度的无限正则树上的唯一性。 这可以通过以下推理来验证:根据定理15,在无限-正则树上,唯一性意味着强烈的空间混合,对于 ,意味着所有较小度的无限正则树上的强空间混合(包括唯一性)。
存在非单调的二态反铁磁自旋系统。 我们可以选择满足和的反铁磁,其中是Lemma 21的8项中给出的普遍唯一性的上临界值。 由于 Lemma 21的 3 项,对于 而言,在所有足够大的 中,唯一性在 不规则树上成立。另一方面,由于 Lemma 21 的 8 项,当 时, 不可能是普遍唯一的,这意味着存在一个有限的 ,使得系统在 不规则树上是非唯一的。
对于这样的非单调系统,由于定理15,唯一性意味着对于足够大的,正则树上存在强空间混合,但是强空间混合不适用于较小程度的规则树(因为较小树上的非唯一性)。 因此,正则树上的强空间混合与最大度至多图上的强空间混合之间的含义对于一般的二态自旋系统并不成立。
致谢。
我们要感谢 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.
最多为 2 个唯一。
-
2.
如果,那么唯一性不成立于无限-所有足够大的正则树。
-
3.
如果,则唯一性保持无限-所有足够大的正则树 。
-
4.
对于任何(包括),存在一个临界阈值,当且仅当时,直到是唯一的。特别是, 和 对于某个有限的.。
-
5.
如果,对于任何(包括),存在一个临界阈值,使得直到是唯一的,当且仅当.。
-
6.
如果,则对于任何外部字段、 最多为- 唯一。
-
7.
如果,对于任何(包括) ,的制度满足最多-0> 唯一性如下。
让。对于任意整数,令为方程的两个正根>,更具体地说,
让(其中)被定义为整数 .
那么 最多 - 唯一当且仅当 属于以下政权
(8) 如果此外,则可达- 唯一当且仅当 ,其中
特别是,当时,它成立,因此 最多 - 唯一当且仅当 。0>
-
8.
只有当 时, 才可能是普遍唯一的。如果,则存在有限正常数和,其中、使得在条件下是普遍唯一的,而只有在条件下是普遍唯一的。
Remark (SODA’13版本错误)。
引理21的证明。
令和为的正不动点。 然后
令 是反铁磁性的。 即 、 和 ,即 。
-
1.
很容易验证对于,对于任何正。 因此,当时,我们有,这意味着是最多2个唯一的。
-
2.
对于所有足够大的 ,它认为
根据矛盾,假设。 然后,
-
情况1:
。 然后。 因此,
这意味着。 然而,它认为
矛盾。
-
案例2:
。 然后。 因此,
这意味着。 然而,它认为
矛盾。
-
情况1:
-
3.
让。 固定点,以及
对于所有足够大的 来说,它严格小于 1。因此,对于所有足够大的,无限-正则树具有唯一性。
-
4.
我们首先表明,存在一个临界阈值 ,使得 最多 是唯一的,当且仅当 时。 足以证明,如果反铁磁 的唯一性可达 ,则 的唯一性可达 对于任何 和 都是唯一的。
回想一下, 是 的正不动点。还让 表示 的正解。
我们首先展示。 根据矛盾,假设。 由于对于反铁磁,是单调递减的,我们有
矛盾。
自从
我们有
对于,也成立
将以上两个不等式相乘,我们有
请注意,这些是当参数为和时在各自固定点处的绝对导数。 因此,如果 最多为 唯一,则 最多为 唯一。
- 5.
-
6.
我们注意到 在 中不是单调的。 它在处达到最大值。 因此,如果对于任何,对于,则最多是唯一对于任何。 当 对于所有 时,即当 时,此条件成立。
-
7.
让。 它对所有 都适用 。 因此, 是明确定义的,可以表示为:
回想一下 代表 。
下面的的单调性和无界性很容易验证。
Claim 22。
在 中单调递增,并随着 的增长而趋于无穷大。
证明。
观察到对于、和,它认为和在。 而且,
很容易验证 在 中不断增加,并且随着 的增长而无限制。 因此,我们可以得出结论, 在 中不断增加,并且随着 的增长而无界。 ∎
我们继续分析唯一性机制。 对于 ,它保持 ,在这种情况下,由于此引理的第 6 部分,它始终保持 。
回想一下。 仍然需要分析 中的那些 。 对于这样的 ,我们有 ,并且由于 和 是方程 的根>,当且仅当 或 时,才成立 。 请注意, 在 中单调递增。 因此:
这意味着对于任何,当且仅当。
然后我们证明,对于 ,唯一性机制 (8) 变为 。 我们需要以下断言,其证明推迟到本节末尾。
Claim 23。
对于,它总是成立
此外,如果 对于所有整数 ,则
Claim 24。
以下内容成立:
-
•
如果 ,则 对于所有整数 ;
-
•
如果,则存在有限的使得 对于所有。
-
•
-
8.
的情况在该引理的 5 部分中处理。 在这种情况下,由于引理的 5 部分, 最多是 唯一,当且仅当 其中。 特别是,如果 , 变为 0;如果 , 变为正且有限的 。 因此,仅当 时, 才是普遍唯一的,并且当 时, 是普遍唯一的,当且仅当 . 我们可以令 和引理如下。
Claim 25。
以下内容成立:
-
•
如果 ,则 接近 0,如 增长到无穷大;
-
•
如果 ,则 无界,如 增长到无穷大。
证明。
对于整数 ,让 ,并让
(9) 当时,很容易验证随着增长到无穷大,接近0。
此外,的上限以为界。 具体来说,对于,
由于当时在中增加,我们有
因此,对于所有 ,
(10) 因此,当时,随着增长到无穷大,趋近于0。
另一方面,当时,对于,
随着 的增长,这显然是无限的。 ∎
当时,由于上述说法,随着增长到无穷大,趋近于0,这意味着;根据声明22,随着增长到无穷大,也趋于无穷大,这意味着。 因此,只有当时,才可以是全局唯一的。
现在假设。 根据声明 25,随着 的增长, 是无界的。 因此, 在有限的 处达到其最小值。 显然,对于任何有限的 , 总是有限的且为正值。 因此,是有限且正的。 让 表示这个有限的正阈值。
另一方面,根据权利要求 24,对于 ,存在一个有限的 ,使得 适用于所有 。 让。 显然,。 并且 也是有限的。
-
•
∎
索赔证明23。
首先,很容易验证对于,它总是成立
因为对于所有,它总是成立
接下来,我们证明如果对于所有整数 都 ,则两种机制相等。 这是通过归纳法证明的。
当时,权利要求中的相等关系是微不足道的,因为在这种情况下,最多有一个整数 满足 。 那么等式两边就变成了
这就奠定了归纳基础。
对于归纳假设,假设 等式成立:
(11) |
我们在哪里表示
索赔证明24。
它还成立。 对于所有 ,我们有 。 由于在时在中单调递增,因此成立。 同时,对于所有,和成立。 所以,
(14) |
假设。 也不难观察到 在 中正在减少。 事实上,如果,那么在中单调递减;我们可以验证函数对于的单调性,如下所示:
这保证了 在 中单调递减,因此对于 , 在 中单调递减,因为 和 都是正数,并且 中是递减的。
附录 B 引理证明 11(中值定理方法)
-
1.
根据中值定理,存在 使得
-
2.
假设每个顶点最多有 个子节点。 然后由于假设 有界或 。
如果 ,则 ,其中 。 如果,则;如果是有限的,则。
根据中值定理,存在,,使得
如果代表所有,那么由于上述论点,代表免费的,其中是的孩子数量;简单地 表示固定的 。 因此,。 如果,那么很容易看出以常数为界,因此成立;如果是有界的,那么也是有界的,因此显然是。
-
3.
然后我们分析的逐步收缩。 定义 以及相应的 、。 然后是和、。 我们有
应用中值定理。 存在和对应的满足,使得
附录C启发式地找到一个好的势函数
也许我们证明中最神秘的一步是势函数的选择。 正如在许多应用势分析的情况下一样,没有用于搜索合适的势函数的标准例程。 另一方面,我们不太可能在没有任何提示的情况下猜出这样一个相当复杂的公式。 在这里,我们提出了一种启发式方法,引导我们发现良好的潜在函数。 这部分并不严格,对于我们结果的合理性来说,逻辑上也是不必要的。 然而,这种启发式方法足够通用且有趣,并且可能会在其他场景中找到应用,因此值得在这里进行阐述。
启发式方法包括三个步骤:
-
1.
求势函数的必要条件,即与势函数在某一点相关的方程。
-
2.
(启发式步骤)通过假设方程适用于变量的整个范围来增强条件,这给了我们一个微分方程。
-
3.
求解微分方程并得到势函数。
我们首先假设系统处于唯一性的边界。 这意味着 是临界度, 是正不动点,其中 和 是正不动点。
我们有以下身份:
这一起意味着另一个身份
(16) |
我们的目标是找到一个势函数,使 适用于所有 。 另一方面,我们有。 所以在时达到最大值。 作为一个可微函数,它满足,即
我们有以下计算:
我们使用 和 的事实。
对求二阶导数,我们有
该方程已经给出了 的恒等式。 但它太复杂了。 我们可以在点进一步简化上面的公式。
对于固定点,还有其他等式成立。 选择的哪个方程启发式扩展到所有可能会影响我们获得的势函数。 例如,(16) 意味着
因此我们可以将 重写为
这给出了
我们还可以将此方程视为变量 的微分方程,求解它可以得到
这是[15]中使用的势函数。