通过收缩快速混合格劳伯动力学直至达到独特性
摘要
对于一般的反铁磁自旋系统,包括加权独立集上的硬核模型和反铁磁伊辛模型,当无限正则树位于李等人(2013)的唯一性区域时,最大阶数图上的分割函数有一个。 此外,在树非唯一性区域中,Sly(2010)表明,除非,否则不存在来估计配分函数。 算法结果遵循 Weitz (2006) 提出的相关衰减方法或 Barvinok (2016) 开发的多项式插值方法。 然而,运行时间只是常数的多项式。 对于硬核模型,Anari 等人(2020)最近的工作在树唯一性区域建立了简单单点马尔可夫链的快速混合,称为格劳伯动力学。 我们的工作通过考虑固定顶点 对其他顶点的总成对影响(而不是其他顶点对 的总影响)来简化他们对格劳伯动力学的分析,从而扩展他们对所有 2 旋转模型的工作并改进了混合时间。
更重要的是,我们的证明将三种不同的算法方法联系在一起:我们证明了所谓的树递归与合适的势函数的收缩,这是建立 Weitz 相关衰减方法和 Barvinok 多项式插值方法效率的主要技术,也建立了格劳伯动力学的快速混合。 我们强调,这种联系适用于所有 2 自旋模型(反铁磁和铁磁),并且相关衰减或多项式插值方法的现有证明立即意味着格劳伯动力学的快速混合。 我们的证明利用图配分函数是韦茨自回避游走树配分函数的除数。 这一事实催生了用于分析顶点影响的新工具,并且可能对复零的研究具有独立意义。
1简介
在最大度数 的一般图中的近似计数问题的计算复杂性与度数 的无限规则树(或向上)上的统计物理相变之间建立了显着的联系。在更一般的情况下为 )。 这种联系适用于二态反铁磁自旋系统——独立集上的硬核模型和伊辛模型是此类系统最有趣的例子。
给定一个 顶点图 ,2 个自旋模型的配置是对顶点的自旋 的 分配。 2-自旋系统由三个参数定义:边权重和顶点权重。 边缘参数控制相邻-自旋之间相互作用的(相对)强度,对应于相邻-自旋,并且 是应用于具有 -自旋的顶点的外部场。
每个自旋配置 都分配有一个权重
其中,对于自旋,是具有自旋的单色边的数量,是具有自旋的顶点的数量(按照标准,参数已标准化,因此我们可以避免两个额外的参数)。 自旋构型的吉布斯分布由 给出,其中 是配分函数。
有两个特别有趣的例子:硬核模型和伊辛模型。 当和时,唯一具有非零权重的配置是的独立集和独立集的权重是;这个例子被称为硬核模型,其中参数对应于逸度。
在的情况下,重要的量是单色边的总数,配置的权重是;这是经典的伊辛模型,其中参数对应于逆温度,是外部场(表示没有外场)。 请注意,当时,模型是铁磁,因为相邻顶点更喜欢具有相同的自旋,而是反铁磁伊辛模型。 在一般的自旋系统中,当时模型是铁磁的,当时模型是反铁磁的。 (当 时,我们得到一个简单的产品分布。)
基本算法任务是从吉布斯分布中采样并估计配分函数。 对于近似采样问题,我们给出一个图 和一个 ,我们的目标是从总变异距离内的分布 生成样本 的吉布斯分布 在时间 上的分布。 高效的近似采样算法意味着近似计数问题[JVV86,ŠVV09]的(完全多项式随机近似方案)。 回想一下,给定一个 顶点图 和 , 输出一个 近似值 的概率 在时间 中,而 是确定性模拟(即 )。
近似采样问题的标准方法是马尔可夫链蒙特卡罗 (MCMC) 方法;事实上,有一个简单的马尔可夫链,称为格劳伯动力学。 格劳伯动力学的工作原理如下:从时间 的配置 中,选择一个随机顶点 ,然后设置 为所有,最后我们从的条件分布中选择。 对于硬核模型的情况,如果当前没有邻居被占用,则被设置为被占用(即自旋),概率为,否则它被设置为空闲。
很容易验证格劳伯动力学是遍历的,吉布斯分布是唯一的平稳分布。 混合时间是保证从最差初始状态开始,的分布在总变化距离 的吉布斯分布。 目标是证明混合时间是中的多项式,在这种情况下,链被称为快速混合。
对于铁磁 Ising 模型(有或没有外场)的情况,Jerrum 和 Sinclair [JS93] 的经典结果通过 MCMC 方法为所有图给出了 。 这是通用图的有效算法的唯一情况。 对于反铁磁二自旋模型,该图与常规树上的统计物理相变密切相关。
硬核模型的情况很好地说明了唯一性/非唯一性相变。 考虑以 为根的无限 -正则树 ,并让 表示在第一个 此相变捕获 叶子处的配置是否在限制 内“影响”根。 对于硬核模型,我们可以考虑偶数高度树(对应于全偶数边界条件)与奇数高度树。 令表示吉布斯分布中根被占据的边际概率。 让 和 。 如果成立,我们就说树唯一性成立;如果不相等,我们就说树非唯一性成立。 对于所有,都存在临界逸度[Kel85],其中树唯一性保持当且仅当。
显着的联系是最大度 的一般图的算法相变发生在同一树临界点。 对于所有常量、所有、所有、所有最大度图、[Wei06] 提出了一个用于近似配分函数的。 另一方面,对于所有,所有,[Sly10,SS14,GŠV16]证明,除非,没有用于估计配分函数的。
一个重要的警告是,Weitz 算法的运行时间为 ,其中近似因子为 ,常数 多项式取决于间隙 (回想一下,)。 Sinclair 等人[SST14]将Weitz的相关衰减算法扩展到树唯一性区域中的反铁磁Ising模型,以及相应树唯一性区域中的所有反铁磁2-自旋系统(如下详述) ) 作者:Li、Lu 和 Yin [LLY13]。
Barvinok [Bar16] 提出了一种有趣的新算法方法,并由 Patel 和 Regts [PR17] 改进,利用复平面中配分函数零点的缺失来使用泰勒近似有效地近似配分函数对数的合适变换。 这种多项式插值方法被证明在与 Peters 和 Regts [PR19] 的 Weitz 结果相同的树唯一性区域中是有效的,尽管运行时间的指数取决于 。
长期以来,人们推测简单的格劳伯动力学正在树的独特性区域中迅速混合。 Anari、Liu 和 Oveis Gharan 最近证明了这一点[ALO20];他们证明,对于所有,只要,混合时间就是。 我们改进了这个结果。 首先,我们将混合时间从 改进到 ,如以下定理所述。
Theorem 1 (硬核模型)。
令 为整数,。 对于最大度的每个-顶点图和每个,硬核模型的格劳伯动力学的混合时间在 上,逸度 为 。
除非进一步改进 [AL20] 中的局部到全局参数,此界限是最佳的。 我们改进的结果源于更简单、更清晰的证明方法,该方法使我们能够将结果扩展到各种 2-spin 模型,与相关衰减算法的关键结果相匹配,并大大缩短运行时间。
我们的证明方法统一了近似计数的三种主要算法工具:相关衰减、多项式插值和 MCMC。 相关衰减和多项式插值方法的大多数已知结果都是通过在所谓的树递归上显示适当定义的势函数的收缩来证明的;树递归是韦茨的自回避游走树的结果,我们将在本文后面更详细地描述。 Shao 和 Sun 最近的一项工作[SS20] 统一了这两种方法,表明通常用于证明相关衰减算法效率的收缩也意味着(在一些额外的分析条件下)多项式插值方法是有效的。
在这里,我们证明了势函数的相同收缩也意味着格劳伯动力学的快速混合,以及我们改进的独立于 的运行时间;详细说明请参见4和5。 我们的证明利用了有关韦茨自回避游走树的几个新工具,这些工具在部分 3中有详细介绍。 特别地,我们证明了图的配分函数除以Weitz自回避游走树的配分函数;请参阅8。 这一结果对于确定具有复杂参数的配分函数不存在零点具有潜在的独立意义,因为它使人们能够考虑自回避游走树。 这一结果还产生了一个新的、有用的等价物,用于根据自回避树来限制图中的影响,这加强了 Weitz [Wei06] 先前已知的联系;详情请参见8。
作为一个简单的结果,我们在树唯一性区域中获得了反铁磁 Ising 模型的 Glauber 动力学的快速混合。 就边活动而言,-正则树上 Ising 模型的两个临界点位于 和 ;第一个位于反铁磁区域,而第二个位于铁磁区域。 如果,则唯一性适用于-常规树上的所有外部字段。
如前所述,对于铁磁 Ising 模型, 被称为通用图 [JS93]。 此外,Mossel和Sly[MS13]证明了铁磁Ising模型的Glauber动力学在时的混合时间。 然而,反铁磁伊辛模型在树独特性区域的快速混合尚不清楚。
我们为情况 提供以下混合结果。 请注意,当 时,外部字段 的某些值有一个额外的唯一性区域;该区域由 3 覆盖。
Theorem 2 (反铁磁伊辛模型)。
令 为整数,。 假设和。 然后对于最大度的每个-顶点图,伊辛模型的格劳伯动力学在上的混合时间边权重为,外场为。
我们的硬核模型和伊辛模型的结果适合一般反铁磁二自旋系统的更大框架。 回想一下,反铁磁情况是在 时。
对于一般的 2-自旋系统,适当的树相变更加复杂,因为存在树唯一性阈值在 中不是单调的模型。 因此,正如 [LLY13] 所考虑的那样,适当的概念是“最高唯一性”。 粗略地说,如果对于每个整数 ,距离 处的所有顶点,我们说间隙 的唯一性在 正则树上成立> 从根开始对根的边缘有总的“影响”。 如果带有间隙 的唯一性在所有 常规树上成立,则我们说最多 带有间隙 的唯一性成立;准确定义请参见部分 2。
Theorem 3 (一般反铁磁-自旋系统)。
令 为整数,。 令为实数,使得、、和。 假设参数 最多 是唯一的,间隙为 。 然后对于最大度的每个-顶点图,反铁磁-自旋的格劳伯动力学的混合时间 上带有参数 的系统是 。
1.1电位法混合
树递归在近似计数的研究中非常有用。 考虑一棵以 为根的树。假设 有 子级,用 表示。 对于 ,我们将 定义为以 为根的 的子树,其中包含 的所有后代。 让 表示根的边际比率,并且 对于每个子树。 由于的独立性,树递归是一个在给定的情况下计算的公式。 更具体地说,我们可以编写 ,其中 是一个多元函数,对于 ,
然而,在本文中,我们特别关注边际比率的对数。 原因是我们要仔细研究[ALO20]中介绍的吉布斯分布的成对影响矩阵 并定义为每个
在[ALO20]中,作者表明,如果的最大特征值有适当的界限,那么格劳伯动力学就是快速混合。 我们在本文中做出的一个重要观察是,影响 上的 可以被视为 相对于 处的日志外部字段的导数(参见 12 因此,我们使用对数比率更方便。 为此,我们将树递归重写为 ,其中 是一个函数,对于 ,
观察。 此外,我们定义
对于 ,以便 对于每个 。
为了证明我们的主要结果,我们使用了已广泛用于建立相关性衰减的势方法。 通过为对数比率选择合适的势函数,我们表明给定顶点的总影响随着距离呈指数衰减,从而建立了格劳伯动力学的快速混合。 让我们首先明确我们对潜力的要求。 对于每个整数 ,我们定义一个有界区间 ,其中包含 度顶点处的所有对数比率。更具体地说,当 时,我们让 ;当 时,我们让 。 此外,将定义为包含次数小于的所有对数比率的区间。
Definition 4 (-势函数).
在-势的定义中,应该将视为顶点处的对数边际比,势函数为。 以下定理在给定 -势函数的情况下建立了格劳伯动力学的快速混合。
Theorem 5。
令 为整数。 令为实数,使得、和。 假设对于某些 和 ,存在相对于 和 的 势。 然后对于最大度的每个-顶点图,-自旋系统的格劳伯动力学的混合时间带有参数 的 是 。
我们在部分3中概述了我们的证明。 请注意,在 4 和 5 中,常数 允许依赖于最大度数 和参数 总的来说。 例如,对 [LLY13] 中的势进行简单的黑盒应用将为 有界性 条件提供 ,从而产生 混合。 然而,这对于具有潜在无限度数的图来说是不可取的。 我们的贡献之一在于,我们证明了 Boundedness 条件在一个通用常数 与 和 无关的情况下成立。 因此,我们的混合时间是 ,除了 之外,指数中没有任何参数。
我们注意到,在潜在方法的所有先前工作中,结果和证明总是以 、 的树递归和 的形式呈现, 的潜在函数。事实上,我们的结果也可以翻译成的语言。 要看到这一点,从开始,如果我们选择,则可以直接检查,从而检查。 这意味着当且仅当相应的收缩条件对于成立时,4中的收缩条件对于成立。 有界性条件也可以等效地表述为。 尽管如此,在本文中我们选择使用 出于以下两个原因。 首先,正如前面提到的,事实是 是 的导数,因此可以很自然地考虑对数比率的树递归。 事实上,直接使用 呈现我们的结果和证明比切换到 更容易、更清晰。 其次,我们将使用的势函数 是通过转换 从 [LLY13] 中的精确势 获得的。111更准确的说,我们还乘了一个常数因子,它只是简化了我们的计算,并没有多大关系;另请注意,[LLY13] 表示 的势函数及其 的导数。 有趣的是,这个势能的导数就是。 那么收缩条件有一个很好的形式:; Boundedness 条件仅涉及 的上限。 这似乎揭示了 [LLY13] 中神秘的势函数 ,并且还表明 是树递归的一个有意义的变体考虑。 再添加一个证据,对于很多情况(例如,),选择了潜在的,这仅仅意味着我们可以选择恒等函数和 本身正在收缩,没有任何不平凡的潜力。
2021 年 7 月修订。
这篇论文在 FOCS 2020 上发表后,在[LLY13]中发现了关于反铁磁二自旋系统唯一性区域描述的一个小错误。 该错误已在最新版本的[LLY13]中修复。 在本次修订中,我们更新了部分7和附录E中相应的结果和证明 依赖于 [LLY13] 中的更改;特别地,36根据当前唯一性区域的描述进行调整。 我们指出,这些变化纯粹是技术性的,不会影响我们主要结果(例如5)的有效性。
致谢。
我们要感谢 Shayan Oveis Gharan 和 Nima Anari 激发了讨论。 我们也感谢匿名审稿人提出的有益的意见和建议。 我们感谢 Yitong Yin 与我们沟通[LLY13]的最新更新,并提供有关修改附录中的陈述和结果证明的有用说明E,特别是36。
2 预赛
混合时间和光谱间隙
令 为具有平稳分布 的有限状态空间 上遍历(即不可约且非周期)马尔可夫链的转移矩阵。 让表示从开始的步之后链的分布。 的混合时间定义为
如果对于所有,我们说是可逆的。 如果是可逆的,则只有实特征值,可用表示。 的光谱间隙定义为,的绝对光谱间隙为定义为。 如果 相对于内积 也是半正定的,则 的所有特征值都是非负的,因此 也是非负的。 最后,混合时间和绝对光谱间隙的关系为
(1) |
独特性
令 为整数或 。 令为实数,使得、、和。 对于 ,定义
并用表示的唯一不动点。 对于,我们说参数是up-to-唯一的,间隙如果 对于所有 。
比例及影响力
考虑图 上的 -自旋系统。 让 和 。 对于所有,我们将处的边际比率定义为
对于所有,我们定义对的(成对)影响:
为(成对)影响矩阵编写,其条目由下式给出 .
韦茨的自回避行走树
令 为连通图, 为 的顶点。 自回避行走(SAW)树定义如下。 假设顶点集存在全序。从 出发的自回避步行是一条路径 ,使得所有 都具有 。 SAW树是一棵以为根的树,由所有带有的自回避行走以及附加的关闭循环的又一个顶点(即 对于某些 ,例如 )。 请注意, 的顶点可能在 SAW 树中具有多个副本,并且除了叶子之外,顶点的度数将被保留。 请参阅图1了解示例。
我们可以使用相同的参数在上定义一个-自旋系统,其中一些叶子被固定到特定的自旋。 更具体地说,对于附加了 的自回避行走 ,我们将 修复为旋转 如果 相对于 上的总排序,如果 则旋转 。 对于每个 ,我们用 表示 中 的所有自由(未固定)副本的集合。 对于 和部分配置 ,我们通过将自旋 分配给 的每个副本 从 开始,并删除 的所有后代,为每个 定义具有条件 的 SAW 树。 请注意,一般来说,来自 的 的不同副本可以接收不同的自旋分配。 最后,如果每个顶点 都有不同的字段 ,则来自 的 的所有副本都将具有相同的字段SAW 树中的 。
3 主要结果的证明大纲
步骤 1 ([ALO20]):光谱独立意味着快速混合。
我们的证明建立在[ALO20]的基础上,他证明了从最大程度图上的硬核分布中采样的格劳伯动力学最多混合在步中每当 时。 他们证明的关键要素之一是他们称之为谱独立性的概念。 [ALO20]表明光谱独立性意味着快速混合。 请注意, 的对角线条目是 ,而不是 [ALO20] 中原始定义中的 。
Definition 6 (谱独立性[ALO20])。
我们说 顶点图 上的吉布斯分布 与 谱无关,如果对于大小为 和 的每个 、,都有 。
Theorem 7 ([ALO20]).
如果 是一个 谱独立分布,则从 采样的格劳伯动力学至少具有谱间隙
我们现在的主要目标是限制 的最大特征值。
步骤 2:自我回避行走树保留影响。
从标准线性代数,我们知道 的最大特征值的上限是 -范数 ,对应于顶点 的总影响,以及无穷范数 ,对应于的总影响。在[ALO20]中,作者使用作为的上限。 粗略地说,他们证明了在所有边界条件下,固定顶点 上的绝对影响总和的上限是以 为根的自避让行走树中 上的最大绝对影响。 在本文中,我们将使用 来代替上限 。 事实上,如果我们看看自回避树中 的影响,则事实更多。 我们证明对于每个顶点,影响 中的内容以影响力总和的形式保存在以为根的自回避游走树中 的所有副本 。
我们建立这一事实的方法是将配分函数视为 中的多项式。 事实上,考虑每个 具有任意外部字段 的更一般情况将很有用。 让 表示字段。 对于和,以为条件的的权重定义为,其中 是 的 中至少有一个端点的 - 边的数量。 此外,是以为条件的配分函数。 我们将和视为一些固定常量,并将视为变量。 从这个意义上说,我们将权重视为中的单项式,将配分函数视为中的多项式。 此外,边际比率及其影响 for 都是 中的函数。 我们的主要结果是 的配分函数除以每个 的 的配分函数。 由此,我们表明 SAW 树保留了根的影响,并重新建立了 Weitz 的著名结果[Wei06],请参阅13。
Lemma 8.
令 为连通图, 为顶点, 使得 连通。 令为以为根的的自回避游走树。然后对于每个,除以。 更准确地说,存在一个独立于 的多项式 ,使得
(2) |
作为推论,对于每个顶点 ,
(3) |
其中 是 中 的所有自由(未固定)副本的集合。
第三步:给予良好潜力的影响力的衰减。
树递归为我们提供了一个很好的工具来递归计算树的顶点(对数)比率。 正如我们在 12 中所示,影响 实际上是 处对数边际比的导数版本。因此,可以自然地使用树递归来关联这些影响。 然后,我们应用势方法(该方法已在文献中广泛使用)来建立相关性的衰减(强空间混合)。 下面的引理表明,对距离 的绝对影响之和随 呈指数衰减,可以将其视为成对影响的衰减。
Lemma 9.
如果存在关于和的-潜在函数,其中和,然后对于每个 、 和所有整数 ,
其中表示距距离为的所有自由顶点的集合。
第四步:找到好的潜力。
作为最后一步,我们需要找到 4 中定义的 势函数。 我们选择的潜在正是来自[LLY13]的那个,适应对数边际比率和树递归(参见部分6了解更多详情)。 我们证明,如果参数 最多 是唯一的,且有间隙 以及 或 ,则 是 势。
Lemma 10.
令 为整数。 令为实数,使得、、和。 假设 最多 是唯一的,间隙为 。 隐式定义函数
(4) |
如果,则是具有和的势函数。 如果 和 ,那么 与 和 是一个 电位;如果 ,我们可以进一步取 。
4 保留自回避行走树的影响
在本节中,我们将展示[Wei06]中介绍的自回避行走(SAW)树(另见[SS05]),保留了根的所有影响,从而建立8。 为此,我们将 的配分函数(视为外部字段 的多项式)除以 SAW 树的配分函数。 由此我们证明根顶点 对 中另一个顶点 的影响与对 使用我们的证明方法,我们表明根的边缘保留在 SAW 树中,重新建立了 Weitz 著名的结果 [Wei06],以及与 相关的所有成对协方差> 被保留。
Theorem 11.
令 为连通图, 为顶点, 使得 连通。 令为以为根的的自回避游走树。然后对于每个,除以。 更准确地说,存在一个多项式 使得
此外,多项式独立于。
Remark 2。
我们注意到,[Ben18] 证明了硬核模型的 11 的单变量版本,而 [LSS19] 对于零模型显示了类似的结果。具有均匀边缘权重的伊辛场模型。 我们的结果适用于所有 自旋系统和每个顶点的任意字段。 我们还可以以简单的方式将其推广到每条边的任意边权重。 至关重要的是,商多项式与根处的场无关,由此我们可以立即推导出根的边际保留和影响。
在证明11之前,我们首先给出它的一些结果。 对于所有 ,我们将 处的 marginal 定义为 (为了方便起见,我们将事件 写为 ),以及 和 的 协方差 > 作为
以下引理将我们感兴趣的量与(对数)配分函数的适当导数联系起来。 引理的第一部分和第二部分是民间传说。
Lemma 12.
对于每个图 、 和 ,以下内容成立:
-
1.
对于所有,
-
2.
对于所有,
-
3.
对于所有,
证明。
前两部分是标准的。 为了完整性,我们将这两个事实的证明包含在附录B中。 对于第 3 部分,我们从第 2 部分推断出
仍有待证明的是
这实际上适用于任何两个二元随机变量。 为了看到这一点,我们首先计算 根据定义:
同时,协方差可以写为
这表明 从而确立了第 3 部分。 ∎
Lemma 13.
令 为连通图, 为顶点, 使得 连通。 令为以为根的的自回避游走树。那么对于每个 我们有:
-
1.
([Wei06,定理3.1])根边缘的保留 :
-
2.
保留协方差和影响:对于每个,
其中,是在.中的所有自由(非固定)副本的集合。
证明。
4.1 11的证明
在提出证明之前,让我们首先回顾一下之前介绍的符号和定义。 用 表示所有顶点处的字段集。 对于和,以为条件的的权重由下式给出
其中对于 , 表示两个端点都接收自旋 且至少其中一个位于 中的边数>。 以为条件的配分函数定义为。 对于 SAW 树,我们以相同的方式定义条件权重和配分函数。 特别是,回想一下,当我们在 SAW 树上修复条件 时,我们还会删除每个 的 的所有后代。
对于每个 和 ,我们应编写 来表示配置集,使得 (即 ) 并让 是所有配置的权重之和 . 我们进一步扩展这个符号并写成 对于每个 和 。 对于 SAW 树,我们也采用相同的符号。
11 的证明。
我们将证明存在一个独立于 的多项式 ,使得
(5) |
方程 5的高级证明思路与[Wei06,定理3.1]中的相应结果类似>。 令 为 中至少有一个端点的边数。 我们对 使用归纳法。当 时,该语句自 起就很简单。 假设 Eq. 5 对于所有图形和所有小于 边的条件成立。 假设根 有 个邻居 。 定义为将顶点替换为顶点,然后连接得到的图对于。
首先考虑仍然连接的情况。 对于每个,让。 使用相同的参数 在 上定义 -自旋系统,加上顶点 固定为自旋的附加条件而固定为自旋;我们用 和 来表示这种条件。 然后,可以通过以下递归过程生成。 另请参阅图2了解说明。
算法:
-
1.
对于每个,让加上条件;
-
2.
通过将 连接为 ,使 成为 和 的并集;输出。
为了证明的目的,我们还考虑 上具有相同参数 的 -自旋系统,但我们让顶点 没有字段(即为 设置 而不是 )。 然后我们观察到
将自旋 替换为 也是如此。 对于,让表示条件和的并集,其中。 然后对于每个 我们有
请注意,两侧都独立于字段:对于左侧,所有在上都没有用于旋转系统的字段;对于右侧,请记住,我们不计算每个 的条件分区函数的固定顶点的权重。 现在定义 通过
它独立于 。 然后我们得到
使用类似的论证,我们也有
由于我们假设 是连通的,因此对于每个 来说,图 也是连通的。 然后,根据归纳假设,对于每个 都存在一个多项式 使得
这些多项式独立于 ,因为 的条件分配函数不涉及 。 现在如果我们让
那么从树递归可以得出
另一个平等 以同样的方式建立。 这样就完成了已连接的情况的证明。
如果有两个或多个连通分量,那么我们可以通过每个分量的SAW树构造。 回想一下, 是通过将顶点 拆分为图中 中的 个副本来定义的。假设 具有整数 的 连接分量。 令 为每个组件导出的子图,以及来自与其相邻的 的顶点。 对于每个 ,让 是通过将 的所有副本收缩为一个顶点 从 得到的图形,并让 . 观察一下,一旦我们收缩 的根 ,生成的树就是 。
我们在每个 上使用相同的参数 定义 -spin 系统,只是顶点 没有字段(即 而不是 )。 对于,令和为限制在上的配置。 然后 与每个 连接,并且由于 ,每个具有条件 的 的数量少于 边缘。 因此,我们可以应用归纳假设;即,对于,存在一个独立于的多项式,使得
我们将多项式 定义为
然后很容易检查
和类似地 . 定理如下。 ∎
5 对树木的影响范围
在本节中,我们研究根对树中其他顶点的影响。 我们给出根对固定距离处所有顶点的总影响的上限。 为此,我们应用了势方法,该方法已用于建立相关衰减属性(例如,参见[LLY12,LLY13,GL18])。 给定任意势函数 ,我们的上限是根据 的属性,涉及 和 上的界限,其中 。 然后,在 为 势的情况下,我们推导出 9。
假设是一棵以为根、最大度数至多的树。 令 和 是任意且固定的。 考虑 上的 -自旋系统,参数为 ,以 为条件。 我们需要限制影响力 从根到另一个顶点。 请注意,如果在删除 时 与 断开连接,则 由自旋系统的马尔可夫性质。 因此,我们可以假设,通过删除所有此类顶点, 只包含 的叶子。
对于顶点 ,令 为以 为根的 的子树,其中包含 的所有后代>;请注意。 我们将 写为 中距 距离为 的所有自由顶点的集合。 我们特别关注子树中处的边际比率,为了简单起见,写成。 通过树递归 关联。如果顶点 有 子节点,用 表示,则树递归由下式给出
其中对于 和 ,
另请记住,对于 ,我们定义
和所有 的 。
下面的引理允许我们使用任意势函数将所有影响的总和从根限制到距离 。
Lemma 14.
令 为具有图像 和导数 的可微且递增(势)函数。 用表示根的次数。 然后对于每个整数 ,
在哪里
在证明14之前,我们首先提出对树的影响的两个有用的属性。 首先,在[ALO20]中表明,影响满足以下形式的树链式法则。
Lemma 15 ([ALO20,引理B.2])。
假设 是三个不同的顶点,使得 位于从 到 的唯一路径上。然后
其次,对于树上的两个相邻顶点,一个对另一个的影响由函数给出。
Lemma 16.
令 和 为子树 中 的子级。 然后
证明。
我们现在准备证明14。
14 的证明。
然后我们推导出 9 作为推论。
6 验证良好的潜力:收缩
在本节中,我们第一步来证明10。 令 为整数。 令为实数,使得、、和。 回想一下,我们通过导数定义了我们的势函数
(1) |
我们在附录C中包含一个简短的证明,以表明是明确定义的。 If is up-to- unique with gap , then we show that satisfies the Contraction condition for . 这适用于唯一性区域中的所有参数,无需要求。 在后面的 附录 E中,我们建立了 当 时的 有界性 条件,完成了 10 的证明。 的情况比较复杂,留给Section7。
在给出证明之前,我们首先指出势函数 本质上与 [LLY13] 中使用的势函数 相同(请注意,[LLY13] 使用 作为势函数的符号,使用 作为其导数)。 回想一下,边际比率的树递归由函数 给出,其中 使得对于所有 ,
来自 [LLY13] 的势函数 通过其导数隐式定义为
以下引理解释了我们如何从 获得潜在的 。
Lemma 17.
我们有;即, 代表所有 。
证明。
检查起来很简单
所以,
结合[LLY13]引理12、13和14的结果,我们得到当在时,势函数满足以下梯度界限独特性区域。 请注意,这可以被视为收缩条件,但对于和而言。
Theorem 18 ([LLY13]).
令 为 的图像。 如果参数 最多为 唯一,且间隔 ,则对于每个整数 使得 和每个 ,
其中。
Lemma 19.
令 为 的图像。 如果参数 最多为 唯一,且间隔 ,则对于每个整数 使得 和每个 ,
其中。
7 剩余的反铁磁情况:和
在本节中,我们讨论 和 的情况。 正如[LLY13]中所研究的,在这种情况下,唯一性区域更加复杂。 例如,存在一个临界,使得具有的自旋系统处于任意图的唯一性区域中;也就是说, 最多是 是唯一的。 为了处理大度数,我们需要放宽4中的有界性条件,并定义更通用版本的势。 我们将看到 5 对于这个一般的 潜力仍然成立。 其背后的原因是,为了限制影响力矩阵的最大特征值,只需考虑具有大度数的顶点的绝对影响力的顶点加权和。
Remark 20。
回想一下,我们的目标是限制矩阵 的最大特征值。 我们可以通过限制绝对行总和的上限来做到这一点 固定的 ,从而为我们提供了 的有效上限。 然而,当和时,这种方法不起作用。 在这种情况下,势 无法成为独立于 的通用常数 的 势。 事实上,不存在这样的 势作为绝对行和 可以与一样大。 特别是,如果参数 最多为 唯一,这意味着自旋系统对于任意图都具有唯一性,则绝对行和 可以是 ,其中 。 我们举一个具体的例子来说明这种情况。
Example 21.
考虑由参数 、 和 在以 为中心、有 叶的星形图上指定的反铁磁性 2 自旋系统。 一个简单的计算表明 对于任何叶顶点。 因此, . 现在,从 开始,我们有
强迫 即使 位于唯一性区域。 然而,我们仍然有 因为 .
为了解决这个问题,人们可能需要考虑绝对列总和,涉及对固定顶点的绝对影响的总和。 然而,这不允许我们使用图和 SAW 树之间的美丽连接,如 8 所示。 相反,我们在这里考虑 绝对行和的顶点加权版本,它也是最大特征值的上限。
Lemma 22.
令 为顶点的正权重函数。 如果存在一个常量 ,对于每个 我们有
(8) |
然后。
证明。
让。 那么假设就等价于。 由此可见。 ∎
然后,我们修改 4 中 -势的定义,以允许较弱的 有界性 条件。 我们注意到 23 和 4 之间唯一的两个区别是:我们允许 ;并且有界性条件被放宽到我们所说的一般有界性。 回想一下,对于每个 ,我们在 时设置 ,在 时设置 。
Definition 23 (一般-势函数)。
请注意,一般有界性 是比有界性 更弱的条件。 要看到这一点,如果潜在函数 满足带有参数 的 有界性,则对于每个 和每个 其中 我们有
以下定理概括了 5 并表明一般的 势函数足以建立格劳伯动力学的快速混合。
Theorem 24.
令 为整数或 。 令为实数,使得、和。 假设对于某些 和 存在相对于 和 的一般 势。 然后对于最大度的每个-顶点图,-自旋系统的格劳伯动力学的混合时间带有参数 的 是 。
Lemma 25.
令 为整数。 令为实数,使得和。 假设 最多 是唯一的,间隙为 。 那么由Eq.4隐式定义的函数就是一个通用的-势函数与和;如果,我们可以进一步取。
24的证明可以在附录D中找到。对于25,的收缩条件来自19,并且一般有界性1>是与所有其他案例一起在附录3>E4>2>中得到证明。
8 铁磁案例
在铁磁情况下,[GL18, SS20] 中给出了最著名的相关衰减结果。 使用[GL18]和[SS20]中的势函数,我们显示以下两个结果,与已知的相关衰减结果相匹配。 事实上,[SS20] 的势函数是常量 和 的 势函数。
Theorem 26.
固定一个整数、正实数和,并假设满足以下三个条件之一:
-
1.
,是任意的;
-
2.
和;
-
3.
和。
那么恒等函数(基于[SS20]中给出的势)是-的势函数,并且。 此外,对于每个最大度数至多的-顶点图, 与参数 是 ,对于通用常量 。
Remark 3。
[GL18] 的势函数确实是一个 势函数,但不幸的是, 必须依赖于 。 我们得到以下结果,该结果弱于无界度图[GL18]中的相关衰减算法。
Theorem 27.
固定一个整数和满足、和的非负实数。 然后对于每个-顶点图,其最大度数至多为,铁磁2-自旋系统的格劳伯动力学在 与参数 是 ,对于常量 仅取决于 ,但不取决于 .
这些定理的证明在附录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 -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主要结果证明
5 的证明。
请注意,由于格劳伯动力学的转移矩阵 具有所有非负特征值,因此我们有 ,因此为了推导出混合,它足以满足下界 。 我们通过使用7来做到这一点。 对于足够小的,足以显示的光谱独立性。
要绑定,只需绑定 对于具有 顶点的所有图 以及 顶点的子集 上的所有边界条件 。 我们声明如下:
(9) |
其中 是仅取决于 的常数。 第一个上限 推导如下
(8; ) | ||||
(split the sum by levels) | ||||
(9) | ||||
第二个上限 更简单。 直观上,这意味着每个绝对成对影响 至多是某个常数,因此绝对影响力的总和上限为。 下面的两个主张给出了更精确的陈述,其证明在部分A.2中提供。
Claim 28 (反铁磁案例)。
固定整数和实数,并假设、、和。 那么,对于每个最大阶数至多为 的 顶点图 ,参数为 的 上的反铁磁 自旋系统是 光谱独立的,其常数 仅取决于 。 此外,如果 最多是 唯一,那么我们可以放弃对 的依赖。
Claim 29 (铁磁外壳)。
固定一个整数和正实数,并假设和。 那么,对于每个最大阶数至多为 的 顶点图 ,参数为 的 上的铁磁 自旋系统是 光谱独立的,常数 仅取决于 。
3 的证明。
2 的证明。
A.1 参数 paps 方面的唯一性差距
Claim 30 (硬核模型;来自 [ALO20] 的引理 C.1)。
固定整数 、 和 。 如果,则最多是唯一的,间隔。
Claim 31 (大)。
固定整数和。 如果 ,则 最多 是唯一的,对于所有 来说,间隔为 。 注意如果,这正是条件。
证明。
A.2 恒定大小图的谱独立界限
在本节中,我们证明具有少于 多个顶点的图的谱独立边界,因为对于具有如此少顶点的图,基于树递归收缩的边界变得微不足道。
28 的证明。
如果表示顶点的边际比率,则。 在的情况下,我们也有;如果,我们就有。 由此可见,我们立即就有了界限
对于所有。 请注意,这些常量小于 ,并且仅取决于 ,从而产生第一个声明。
29 的证明。
证明与反铁磁情况相同,我们在此省略。 ∎
附录 B 12 的证明(第 1 部分和第 2 部分)
12 的证明(第 1 部分和第 2 部分)。
为了看到第一个等式,我们直接计算并得到
对于第 2 部分,使用上面的结果,我们还可以得到
附录 C 的技术引理
以下引理意味着由 Eq. 4 给出的势 是明确定义的。
Lemma 32.
对于所有 这样的 ,我们有
证明。
对于 一侧,我们有
类似地,对于 一侧,我们有
附录D势法混合:24的证明
在本节中,我们以与 5 相同的方式证明 24,如 Section 3 中所述。 这里的主要区别是我们考虑绝对影响的加权和 其中 是权重函数。 这足以让我们限制影响矩阵的特征值,如 22 所示。 我们将选择顶点的权重为,度为。给定一般的势,以下引理为我们提供了对距离的绝对影响的加权和的上限。 特别是,它概括了9。
Lemma 33.
如果存在关于和的一般-潜在函数,其中和,然后对于每个 、 和所有整数 ,
其中表示距距离为的所有自由顶点的集合。
Lemma 34.
令 为具有图像 和导数 的可微且递增(势)函数。 用表示根的次数。 然后对于每个整数 ,
在哪里
附录E验证良好的潜力:有界性
在本小节中,我们将展示由 Eq. 定义的势函数 的有界性 或一般有界性 条件。 > 4 在不同的参数范围内。 结合19,我们完成了10和25的证明。
E.1 唯一性区域的预备知识
在本节中,我们简要描述参数的唯一性区域。 这里的所有结果及其证明可以在最新版本的[LLY13]的引理21中找到。
令 为整数, 为实数。 我们假设 、、 和 。 对于 定义
并用表示的唯一不动点。 回想一下,如果 对于所有 都是唯一的,那么参数 直至 都是唯一的,间隙为 。
当时,自旋系统称为硬约束模型。 在这种情况下,存在一个外部场的临界阈值,定义为
当且仅当 时,参数 最多为 唯一。 特别是,当 临界字段由下式给出时
当时,自旋系统称为软约束模型。 如果 ,则 对于所有 来说最多是唯一的 。 如果唯一性区域更加复杂,我们现在描述。 让
这样对于每个 我们有 ,对于每个 我们有 。 对于每个 ,我们将 定义为二次方程的两个正根
更具体地说, 和 由下式给出
在哪里
请注意, 适用于所有 。 对于 我们让
然后,当且仅当 属于以下状态时,参数 最多是唯一的
(10) |
特别是,当 存在两个关键阈值 时,参数 最多可达 唯一,当且仅当如果 或 (即 ),其中
以下关键字段的界限对我们稍后的证明很有帮助。
Lemma 35.
-
1.
如果 ,则对于每个整数 使得 我们有
-
2.
如果 和 ,则对于每个整数 这样我们有
其中。
35的证明被推迟到SectionE.3。
E.2 有界性证明
令 为整数。 令为实数,使得、、和。 回想一下,势函数 定义为
(1) |
令人惊讶的是发现,因为潜在的正是[LLY13]中的那个,如17所示>。 这似乎不是巧合,它提供了一些直觉,说明为什么 [LLY13] 的潜力有效。 更重要的是,这一事实有助于我们证明有界性和一般有界性。 回想一下,对于 和 ,我们让 为具有 子节点的顶点的对数边际比范围。 然后对于每个 和 其中 ,我们有
(11) |
下面的引理给出了 的上限,从中和 Eq. 11 我们推导出 有界性 和一般有界性 立即。 引理中的括号指示边界应用于哪个引理。
Lemma 36.
下面的引理很有帮助,其证明可以在 Section E.3 中找到。
Lemma 37.
对于每个 我们有
我们在这里给出36的证明。
E.3 技术引理的证明
35 的证明。
1. 对于每个 我们有
对于所有整数 ,最后一个不等式是从 得出的。
2. 对于每个 我们有
观察函数 在 时在 中单调递增,因此我们推断
所以,
对于所有整数 ,最后一个不等式是从 得出的。
第二部分可类似证明。 对于每个 我们有
因此,
然后我们得出结论:
对于所有整数 ,最后一个不等式再次从 得出。 ∎
37 的证明。
我们从 AM-GM 不等式推断出
附录F铁磁情况的证明
F.1 26的证明
F.2 27的证明
在本小节中,我们使用 [GL18] 的结果来证明 27。 它们的潜在函数由其边际比率的导数隐式定义为
对于仅依赖于 的常量 (有关精确定义,请参阅 [GL18])。 在我们的背景下,对数比率的相应潜力是
并且受常数限制,具体取决于 的 。
[GL18] 中的主要技术成果之一是表明树递归与势函数 收缩,导数 的边界为存在仅依赖于 的正常数 的感觉,使得 对于所有 而言。 [GL18]将此类函数称为通用势函数。
Theorem 38.
假设是满足、和的非负实数。 那么对于取决于 的常数 和取决于 的常数 来说,函数 是一个 势函数。
与5结合,这使得与仅取决于的常量混合。 我们注意到这比[GL18]中的相关性衰减结果弱,因为不依赖于,因此对于任意图表。
附录 G 混合速度稍快
在本节中,我们通过更仔细地考虑我们基于收缩证明的(非平凡)谱独立性界限与我们获得的(平凡)谱独立性界限之间的权衡,稍微优化某些反铁磁2自旋系统的混合时间结果在 Section A.2 中用于处理恒定大小的图形。
Proposition 39。
假设 子集上的分布 对于 来说是 谱独立的,对于某些 和。 那么从采样的格劳伯动力学具有至少的谱间隙
证明。
假设我们已经以 为条件,将元素的一部分“输入/输出”。 所得分布既是-光谱独立的,又是-光谱独立的。 边界优于的确切阈值由下式给出
我们注意到这样的 仅在 或等效的 时才有意义。 现在,我们基于固定至多 顶点分数,对所有条件分布应用 谱独立界限。 否则我们应用-谱独立性。 我们获得最终的光谱间隙下界
观察一下
我们还有
将它们放在一起,我们就得到了所需的下界。 ∎
有了这个结果,我们可以将其应用于具有 和 的反铁磁模型,因为查看 28 的证明,我们有这样的系统 -大致与光谱无关。
Corollary 40 (软约束)。
修复整数、。 令为满足和的非负实数。 进一步假设 最多可达 且具有间隙 。 然后,对于每个 顶点图 ,其最大度数至多为 ,从具有参数 以 步骤混合。
Corollary 41 (硬约束)。
修复整数 ,修复 ,并让 最多为 唯一,间隙为 。 然后,对于每个 顶点图 ,其最大度数至多为 ,从具有参数 -以步骤混合。