复杂网络重整化

重整化群（renormalization group, RG）理论是研究多自由度系统的强有力框架，通过不同尺度下系统状态、模型参数和耦合常数的描述，实现分辨率层次转换，可识别相变临界点并分析临界行为。传统物理应用中，重整化群基于同质性、对称性、几何结构和局部性构建尺度变换方案。近年来，该理论向复杂网络拓展，但网络系统存在显式几何坐标缺失、节点与子图统计异质性及类格点对称性缺乏三大特征，导致统一重整化流程构建面临挑战。本文梳理网络重整化主流方法，评述理论突破，探讨关键科学问题。

本文要点

重整化群方法在理解和分类多自由度物理系统整体行为中不可或缺； 复杂网络缺乏度量距离、几何规律性或拓扑规则性，直接应用传统物理重整化群技术存在困难； 尽管存在挑战，重整化群在若干场景仍取得显著进展，本综述系统梳理这些成果； 存在诸多开放性问题，论文基于第一性原理总结构建普适网络重整化框架的理论依据与实现路径。

背景知识

重整化和可重整化是现代物理学理解物质现象的基础概念，其核心假设是在足够大的时空尺度上，小尺度细节可被平均掉。重整化方法的形式化历经数十年，促成了对物理现象普遍性的深刻洞察。知识框1总结了传统物理系统中重整化和重整化群变换的基本概念。

除传统物理系统外，自然界和社会中的大型复杂系统包含多尺度相互作用，无先验最佳分辨率。例如，企业间商品流动建模中，将企业聚集成经济部门时关系的一致性表示，以及网络流行病传播中聚合分析与仓室模型结果的差异，均需重整化方法介入。

将物理系统研究经验移植到复杂网络时，需重新审视关键假设，如物质相互作用的局部性和空间规则性。此外，网络动态过程的粗粒化需保持动力学性质。因此，重新表述复杂系统中的重整化方法至关重要。

本文系统综述重整化方法在复杂网络中的应用，分析真实网络缺乏显式几何结构、多尺度异质性和跨尺度交互连接的特征，揭示传统重整化方法的改造需求，回顾粗粒化方案，探讨将重整化群核心原理推广至异质性网络的路径，厘清开放问题。

|知识框1|统计物理和场论中的重整化
物理系统中相互作用的局部性与因果性限制小尺度涨落影响，使宏观尺度仅需少量参数刻画。该思想规范化始于20世纪60年代，Leo Kadanoff将粗粒化概念应用于统计力学系统，理解二阶相变临界点邻域行为[118]。
Ising模型是经典范例，晶格自旋经历无序态与有序态相变，通过引入“块自旋”实现粗粒化，重整化群理论形式化块自旋相互作用的标度行为[111]。早期基于晶格的数值计算后发展为变分方法[169–172]。
重整化群理论通过实空间平均微观细节，定义任意尺度的有效自由度与相互作用，成为统计力学突破。Kenneth Wilson将其发展为连续场论方法[113,173]，不仅为临界现象标度行为奠定基础，还系统化处理场论中的发散量。此后，重整化群方法揭示微观细节迥异的系统在临界点附近的普适行为[175–179]。
过去五十年，重整化群技术扩展至从宇宙学到粒子物理的广泛领域。统计力学中，推广至同质分形格点的非平衡系统[180–186]，其局部平移不变性简化计算。重要范例包括流行病动力学[187]、生物群体行为[188–190]、神经网络[191–194]、动态与有向逾渗[195,196]、观点动力学的选民模型[197]，以及一维链[198]、格点[199]和层级树[200]上的非线性振子同步模型。此外，在张量网络[192,201–203]及信息论框架[204–208]中亦获应用。
非平衡系统中决定普适类的要素往往未知，高阶重整化群计算通常不可行。非微扰重整化群作为替代方案[209–212]，已应用于格点流行病扩散动力学[213]。| |---|

复杂网络重整化的挑战

重整化群方法已成功发现一类系统的普适性特征，但真实世界问题中相互作用极少呈现类格点规则性。结构无序与空间非均匀相互作用导致拓扑几何异质的网络模式，其特征包括：度分布多样性、小世界特性、非零集聚系数、模块化或层级结构 。这些不规则特征与系统有限性引起的边界效应和有限尺度效应 共同构成推广标准重整化群方法的挑战。

重整化群分析通常包含三个步骤：定义粗粒化变量、平均掉局部涨落、重整化耦合和参数。第一步在平移不变度量空间（如规则格点）中直接可行，但在高度非均匀和小世界网络中，如何定义块节点或慢模态成为问题。此类网络常基于受数据限制的经验分辨率构建，直接在聚合水平建模可能无法捕捉细尺度动态，如2008年金融危机和COVID-19疫情中的冲击传播 。

第二步平均局部波动在同质格点中可借助坐标邻近性定义连通性，但在不规则网络化空间中，局部性概念难以直接推广。网络的稀疏性、度分布异质性及小世界特性中的长程相互作用，使复杂网络中局部性定义复杂化。

第三步重整化耦合和参数的复杂性与前两步相关。在无序格点系统中，相互作用拓扑规则但耦合常数异质，可通过随机变量建模。而异质网络中，相互作用拓扑本身无规则，重整化需将细粒度图的概率分布变换到粗粒度图，理论上定义重整化流及其不动点。

受此启发，研究者致力于建立复杂网络重整化的通用框架，下文总结不同策略，先概述粗粒化方法，再深入探讨异质网络中重整化步骤的推广方法。

网络粗粒化方法

真实世界复杂网络缺乏基本格点结构，催生了多样化的粗粒化方案。

基于最短路的方法

基于最短路长度定义节点间距，借助分形几何的盒覆盖思想 实现迭代重整化群变换。该方法以最短路距离为度量，将指定距离内的节点聚合成粗粒化节点，生成重整化群序列，可揭示网络的自相似性与分形特征 ，并用于新生长机制[22]和网络普适类分类。但受小世界特性限制，该方法难以有效解耦和揭示标度行为，度分布自相似性不伴随关联函数标度律，尤其在网络集聚方面。

研究者还提出在规则格点网络添加长程连接实现实空间重整化群变换，用于研究Watts-Strogatz模型临界点行为和格点无标度网络模型，但需预知节点坐标或区分短程与长程连接，而真实网络中此类信息常未知。

谱方法

基于网络谱性质的方法通过保持特定过程动力学行为进行粗粒化。谱粗粒化最初用于研究网络随机游走，保留随机矩阵最大特征值及对应特征向量以保持大尺度行为，剔除快模态同时保留慢模态[26]，已推广至二部网络[27]和耦合振子同步动力学研究[28–30]，也用于有向网络同步[31]和可控性分析[32]。

利用图拉普拉斯矩阵的特征向量与特征值，研究者探索了场论重整化群方法处理有序-无序现象的可能性[33,34]，为建立基于拉普拉斯谱的图重整化理论框架[35]奠定基础，后文将详细讨论。

拓扑方法

基于拓扑性质的粗粒化方法包括节点中心性度量（如度[36]和广义度[37]）、基于邻居重叠的节点相似性度量[38]，以及网络模体分析[41]等。

值得关注的方法包括基于度阈值的重整化[42]和k核分解[43]。前者通过移除度低于阈值的节点定义嵌套子图层级结构，后者将网络分解为满足最小度约束的极大子图，有效揭示真实系统中超越度分布标度不变性的层级自相似性，如聚类谱的标度律。

此外，基于多尺度社区检测中模块化概念的策略[44]，类似于计算机科学中的多层图划分算法，包含基于节点或边属性的粗粒化启发式方法，如重边启发式算法[45]和基于社区的粗粒化方案[46]。拓扑方法依赖特定结构模式，如层级社区或结构等价节点，特征缺失时会失效。

对称方法

利用对称性定义的等价类构造粗粒化方案，在生物科学到计算机科学研究中，表征网络动态过程局部对称性的方法为刻画系统基本模块提供普适原则[47]。网络在节点变换下保持输入树不变时涌现均衡划分[48]，用于定义常微分方程的成块性[49,50]、布尔后向等价[51]和纤维染色性质[52]等概念。

对称性诱导的均衡划分可将系统约化为保留整体动力学与同步特性的动力学变量集合[53]，但通常仅能确定单一划分，无法迭代应用且在对称性缺失体系中缺乏普适性。

工程技术和计算机科学中的粗粒化

开发互联网及其他通信与交通网络的缩尺工程模型是通信技术领域的长期目标[54–56]。近年来，图嵌入技术和图神经网络等深度学习工具用于图结构数据压缩，粗粒化方案旨在降低计算成本与复杂度[57–63]，在设计高效通信与数据存储纠错编码中具有价值[64]。这些方法受网络重整化思想驱动，但聚焦特定过程，对图结构数据与复杂系统多尺度本质的底层原理揭示有限。

信息论方法

信息论方法探索网络不同尺度间的相互作用。复杂网络中更大尺度通常提供更丰富的连接结构信息，即因果涌现现象[65]，可降低节点间关系观测的不确定性。针对图上动力学过程，信息论方法有助于评估状态空间划分对约化后动力学行为的影响[66]。

重整化群方法的基本原理

真实世界网络异质性可能导致不同方法结果矛盾，解决这一问题需基于第一性原理建立普适框架。构建框架需重新审视规则格点和均匀度量系统上的重整化方案，注意前述挑战。

粗粒化变量定义中，引入$$N_\ell \times N_\ell$$邻接矩阵$$A^{(\ell)}$$表示给定分辨率尺度的网络，$$\ell=0$$为最精细层级，$$\ell>0$$为粗粒化迭代层级。加权网络需权重矩阵$$W^{(\ell)}$$。

实空间重整化通过“块节点”定义粗粒化变量，对应划分$$\Omega_\ell$$，将$$\ell$$层级节点映射到$$\ell+1$$层级的粗粒化节点。k空间重整化则需推广对偶空间定义，通过邻接矩阵的显式算符表示分解为基本模态。

实空间中对精细细节的积分消去定义从$$A^{(\ell)}$$到$$A^{(\ell+1)}$$的满射，粗粒化矩阵元素仅依赖细粒度矩阵中包含于粗粒化节点的微观节点元，如： [math]\displaystyle{ A_{i_{\ell+1} j_{\ell+1}}^{(\ell+1)} = \mathcal F\left( \left\{ A_{i_\ell j_\ell}^{(\ell)} \right\}_{i_\ell \in i_{\ell+1}, j_\ell \in j_{\ell+1}} \right) }[/math] 其中$$\mathcal{F}$$为函数，例如二值层级下，当且仅当$$\ell$$层级块间存在连边时，粗粒化块间连边存在： [math]\displaystyle{ A_{i_{\ell+1}j_{\ell+1}}^{(\ell+1)} = 1 - \prod_{i_\ell \in i_{\ell+1}} \prod_{j_\ell \in j_{\ell+1}} \left(1 - A_{i_\ell j_\ell}^{(\ell)}\right) }[/math] 权重矩阵粗粒化类似： [math]\displaystyle{ W_{i_{\ell+1} j_{\ell+1}}^{(\ell+1)} = \mathcal G\left( \left\{ W_{i_\ell j_\ell}^{(\ell)} \right\}_{i_\ell \in i_{\ell+1}, j_\ell \in j_{\ell+1}} \right) }[/math] $$\mathcal F$$和$$\mathcal G$$的选择导致结果依赖方案，且满射性质对应小尺度信息丢失，使重整化群成为半群。对偶空间方法中对快模态的有损积分具有相同性质。

重整化推广到耦合系数和模型参数时，若网络为确定性图，取决于参数性质；若为随机图模型，需包括随机图参数的重整化，如几何网络和多尺度网络重整化部分所述，需分析图和模型参数上诱导概率分布的函数流。在不动点处，重整化流可能统计“可逆”，逆向流可产生细粒化版本。

完成梳理后，本文介绍三种互补方法解决引言中的三个步骤。

几何网络重整化

复杂网络中的空间局部性可通过几何框架恢复[67]，将真实网络嵌入双曲空间，转化结构异质性与非局部性为局部性特征，恢复几何重整化群方法的适用性[68–71]。

1. 1. 1. 几何模型

几何软配置模型[42,72]基于节点流行度与相似性定义度量距离，双曲几何成为嵌入层级结构与小世界结构的自然空间，距离编码节点连接概率，邻近节点更易连接，与随机几何图[73,74]和社交网络隐空间模型[75]思路相似。

在牛顿版本$$\mathbb{S}^D$$模型[42]中，$$\ell=0$$层每个节点$$i_0$$被赋予流行度变量$$\kappa_{i_0}$$（从幂律分布$$\rho(\kappa)$$抽取）和相似性空间坐标。维度$$D=1$$时，相似性空间为半径$$R_0$$的一维球面，节点$$i_0$$赋予角坐标$$\theta_{i_0}$$，定义角距离$$\Delta\theta_{i_0 j_0}$$。节点连接概率$$p_{i_0 j_0}$$遵循引力定律形式[76]： [math]\displaystyle{ p_{i_0 j_0} = \frac{1}{1 + X_{i_0 j_0}} }[/math] 其中[math]\displaystyle{ X_{i_0 j_0} = \frac{(R_0 \Delta \theta_{i_0 j_0})^{\beta}}{(\mu_0 \kappa_{i_0} \kappa_{j_0})^{\max(1,\beta)}} }[/math]，参数$$\mu_0$$决定平均度，节点期望度与$$\kappa$$值成正比，$$\beta$$调控网络拓扑集聚性，作为逆温度参数[76,77]。

该模型在$$\beta=1$$处呈现反常拓扑相变[77]，几何相（$$\beta>1$$）中集聚系数保持有限[42]，非几何相（$$0\leq\beta\leq1$$）热力学极限下集聚系数趋于零，但准几何区域$$0.5\lessapprox\beta<1$$内衰减极慢[78]，表明有限规模网络保留显著集聚性，与某些真实复杂网络性质相当。

$$\mathbb{S}^D$$模型的纯几何形式为$$\mathbb{H}^{D+1}$$模型[72]，隐藏度$$\kappa$$对应径向坐标$$r$$，连接概率转换为仅依赖双曲空间节点间距的费米-狄拉克函数。Mercator工具[79,80]通过统计推断与机器学习技术，将真实网络嵌入此空间，生成的映射可用于可视化、分析、导航及重整化等任务。

该模型族能解释真实网络连接性的关键特征，包括异质度分布[42,72,82]、显著集聚性[42,72,82–84]、小世界现象[42,85–87]、渗流性质[88–90]、谱特征[91]及自相似性[42]，并可扩展至生长网络[92]、加权网络[93]、二部网络[94,95]等多种网络类型。

1. 1. 1. 几何重整化技术

$$\mathbb{S}^1$$或$$\mathbb{H}^1$$模型通过对网络的双曲映射迭代操作，为几何重整化群奠定基础[68-71]。在$$\ell=0$$层配置嵌入空间坐标后，通过粗粒化和重标度得到$$\ell=1$$层，涉及将相似性子空间划分为$$N_1$$个扇区，同一扇区节点成为超节点，若原始节点间存在连接则超节点间建立连边，迭代生成跨尺度重整化流。

几何重整化表明多尺度自相似性是真实世界网络的普遍对称性[68,71]，能准确预测人类大脑连接组多尺度重构中的自相似性[104]，可由$$\mathbb{S}^1$$或$$\mathbb{H}^2$$模型在几何（$$\beta>1$$）[68]和非几何（$$\beta\leq1$$）[71]机制下的可重整化性解释。隐藏度的尺度变换通过要求连接概率不变性获得，角度坐标流由保持原始模型旋转不变性的变换定义。记$$z_{i_\ell} = \kappa_{i_\ell}^{\max(1, \beta)}$$，保持连接概率和网络拓扑性质的变换为： [math]\displaystyle{ z_{i_{\ell+1}} = \sum_{i_\ell \in i_{\ell+1}} z_{i_\ell}, \quad \theta_{i_{\ell+1}} = \frac{\sum_{i_\ell \in i_{\ell+1}} z_{i_\ell} \theta_{i_\ell}}{\sum_{i_\ell \in i_{\ell+1}} z_{i_\ell}} }[/math] 全局参数演化如下：$$R_{\ell+1} = R_\ell / r_\ell$$，$$\mu_{\ell+1} = \mu_\ell / r^{\min(1, \beta)}$$，其中$$r$$为等大小超节点中的节点数，$$\beta$$尺度不变。该变换满足半群性质，概率$$p_{i_{\ell+1} j_{\ell+1}}$$保持原始形式，$$\beta>1$$时平均度$$\langle k \rangle$$按流$$\langle k \rangle_{\ell+1} = r^\nu \langle k \rangle_\ell$$演化，真实网络通常$$\nu>0$$，流趋向高连接度图；$$\nu<0$$趋向一维环；转变处$$\nu=0$$，平均度保持不变。$$\beta\leq1$$时，链接数在重整化下保持不变[71]。

1. 1. 1. 推广到带权网络

通过保持权重与度关系的变换，网络权重结构在重整化群流中得以保留，产生几何权重重整化（GRW）技术[70]，可对加权网络结构进行多尺度展开[105,106]，形成自相似层级，数学基础在于加权$$W\mathbb{S}^D$$模型的可重整化性[93]。实际操作中常选择重整化链接最大权重作为近似方法，在代谢网络、航班网络、科学合作网络等异质度分布加权网络中，使用sup-GRW策略或基于其理论的变换[70]进行粗粒化和重标度时，可观察到几何标度行为。

1. 1. 1. 自相似细粒化

几何分支生长（GBG）模型[69]作为几何重整化的逆向技术，能准确预测世界贸易网和期刊引文网络等增长网络的自相似行为，通过生成自相似度量推广复制网络原始结构，可扩展至任意大网络规模。

从层级$$\ell$$生成细粒度层级$$\ell-1$$时，将原始层级节点以概率$$p$$划分为$$r$$个后代，节点数量增加为$$N_{\ell-1} = N_\ell(1 + p(r - 1)) = bN_\ell$$，其中$$b$$为分支率。不分裂节点保留坐标，后代节点赋予新坐标，隐藏度遵循几何重整化规则，分布$$\rho(z)$$与祖先层级相同，要求度分布为稳定分布[107]，角坐标固定在祖先左右两侧，略有偏移以保持旋转顺序。

设$$\mu_{\ell-1} = b\mu_\ell$$，迭代应用GBG模型生成细粒化自相似层级，保留原始经验连接概率、度分布、集聚系数和社区结构，平均度减小。添加额外链接并设置$$\mu_{\ell-1} = a\mu_{\ell-1}$$（$$a>1$$）得到膨胀版本，避免平均度快速减小。GBG是积分重整化的统计逆过程，膨胀GBG是收缩几何重整化的统计逆过程，加权和非加权$$\mathbb{S}^1$$模型的重整化复现真实网络多尺度自相似性，表明相同原则决定不同尺度网络连接性。

1. 1. 1. 应用

几何重整化与GBG模型结合实现全尺度网络自上而下的多尺度展开，应用包括设计提升单尺度方法性能的多尺度导航[68]，生成保持统计特性的真实网络缩放副本[68,69]，用于研究尺寸依赖现象，如尺寸诱导的随机共振效应、通过有限尺度标度技术分析临界行为[69]，以及模型降阶以进行高保真模拟。

拉普拉斯网络重整化

复杂网络重整化的另一方法[35,108]利用节点间拉普拉斯扩散概念[109,110]，通过实时检测多尺度结构实现顶点和边的粗粒化，最终重整化网络上的动力学模型。该框架可建立直观的实空间版本，类似于统计物理中的Migdal-Kadanoff重整化群方法[111,112]，同时遵循Wilson重整化群思想[113]建立严格的对偶空间或k空间表述，两种表述均通过迭代粗粒化网络，在大尺度范围内保持扩散动力学性质。

1. 1. 1. 将拉普拉斯扩散视为近邻检测

构建异构无向网络拉普拉斯重整化群（LRG）方法的第一步是利用对称图拉普拉斯算子[110,114]$$L^{(0)}$$定义节点在固定尺度下的“等效”邻域。该算子在初始层级$$\ell=0$$由邻接矩阵$$A^{(0)}$$构造，矩阵元素为[math]\displaystyle{ L_{i_0j_0}^{(0)} = k_{i_0}\delta_{i_0j_0} - A_{i_0j_0}^{(0)} }[/math]，其中$$\delta_{i_0j_0}$$为Kronecker delta，$$k_{i_0}$$为节点度。

设$$X^{(0)}(\tau)$$为网络上扩散场在时刻$$\tau$$的值，遵循热扩散方程[math]\displaystyle{ \dot{X}^{(0)} = -L^{(0)}X^{(0)}(\tau) }[/math]，形式化为： [math]\displaystyle{ X^{(0)}(\tau) = K^{(0)}(\tau)X^{(0)}(0) = \mathrm{e}^{-\tau L^{(0)}}X^{(0)}(0) }[/math] 其中扩散演化算子$$K^{(0)}(\tau) \equiv \mathrm{e}^{-\tau L^{(0)}}$$，矩阵元素$$K_{i_0j_0}^{(0)}$$给出场在时间$$\tau$$内从节点$$j_0$$到$$i_0$$通过所有可能路径的扩散比例[115,116]。

1. 1. 1. 实空间中的拉普拉斯重整化群

基于物理诠释，利用$$K^{(0)}(\tau)$$将原始网络$$\ell=0$$的$$N_0$$个节点在任意时刻$$\tau$$划分为$$N_1$$个“扩散等效”子图，粗粒化为超节点定义下一层级$$\ell=1$$网络，通过剪裁微观连边关系确定超节点间重整化后的边。

网络拉普拉斯粗粒化方法概括如下：定义二值化拉普拉斯矩阵$$\zeta^{(0)}(\tau)$$，元素为： [math]\displaystyle{ \zeta_{i_0j_0}^{(0)}(\tau) = \begin{cases} 1, & \text{若 } K_{i_0j_0}^{(0)}(\tau) \geq \min\left[K_{i_0i_0}^{(0)}(\tau), K_{j_0j_0}^{(0)}(\tau)\right]; \\ 0, & \text{其他情况}. \end{cases} }[/math] 需注意$$\zeta^{(0)}(\tau=0)$$为单位矩阵，$$\tau\to\infty$$时所有元素为1。对固定$$\tau^*>0$$，矩阵$$\zeta^{(0)}(\tau^*)$$视为原网络超图的邻接矩阵，包含$$N_1$$个在尺度$$\tau^*$$上通过扩散连接的连通分支，将$$N_0$$个节点划分为$$N_1<N_0$$个超节点。$$\tau^*$$越大，连通分支数量越少且规模越大，$$\tau^*\to\infty$$时仅存单一分支。基于粗粒化加权邻接矩阵$$W^{(1)}$$可通过定义超节点间加权超边构造，权重等于连接对应细尺度簇的微边权重归一化总和，亦可采用统计物理粗粒化中的一致性裁剪方法获取二值邻接矩阵$$A^{(1)}$$，整个过程可迭代执行定义更高层级。

1. 1. 1. k-空间上的拉普拉斯重整化群

该公式化过程与Wilson在连续场论中的重整化群方法直接类比[119]，体现半群性质，是实空间LRG在傅里叶空间的对应，傅里叶基由拉普拉斯特征向量给出，波模态由拉普拉斯特征值替代，在规则格点上特征向量与平面波傅里叶基一致，特征值简化为$$k^2$$。

该方法为重整化过程提供更深刻洞见，允许网络粗粒化和重标度嵌入其中的统计动力学模型（如Ising模型或接触过程），以探测特征尺度及可能的相变。由于$$L^{(0)}$$与$$K^{(0)}(\tau)$$是自伴算子，采用量子力学符号书写： [math]\displaystyle{ L^{(0)} = \sum_{k=1}^{N_0} \lambda_k^{(0)} \left| \lambda_k^{(0)} \right\rangle \left\langle \lambda_k^{(0)} \right| }[/math] 特征向量完备性保证粗粒化变量定义的合理性，后续步骤按Wilson范式[113]实现重整化程序：设定时间$$\tau^*$$确定新分辨率尺度，将拉普拉斯谱划分为“快”特征值$$\lambda_i^{(0)} \geq \lambda^* = 1/\tau^*$$和“慢”特征值$$\lambda_i^{(0)} < \lambda^*$$，积掉“快”区域，重定义慢特征向量贡献的截断拉普拉斯算子，重新计算时间标度使$$\tau^*$$对应粗粒化拉普拉斯算子重整化至层级$$\ell=1$$，迭代生成所有层级的重整化拉普拉斯算子。

1. 1. 1. 扩散演化算子的统计物理学解释

算子$$K^{(0)}(\tau)$$可用于构建网络的“本征尺度扫描仪”，具备平衡态统计物理中熵敏感性或热容的数学特性，其显著极大值（无限尺寸极限下发散）表明网络结构组织中的相变。具体而言，$$K^{(0)}(\tau)$$的特征值$$\{e^{-\lambda_i^{(0)} \tau}\}$$用于定义实时概率密度函数$$p(\lambda; \tau)$$，对应的香农熵$$S(\tau)$$可形式化为冯·诺依曼熵，与量子正则密度算子相关，$$L^{(0)}$$扮演哈密顿量角色，$$\tau$$对应逆温度[130,131]，配分函数$$Z(\tau)$$与自由能$$F(\tau)$$相关。

该形式类比的直接推论是量[math]\displaystyle{ C(\tau)=-\frac{\mathrm{d} S(\tau)}{\mathrm{d} \log \tau} }[/math]具有熵敏感性，是热容[15,129]，节点数无穷大极限下，相变点解释为$$C(\tau)$$在特征尺度$$\tau^*$$处的奇异性，$$\tau^*$$是网络拓扑突变的内蕴尺度，类似Ising模型中关联长度尺度，$$C(\tau)$$可作为本征尺度探测器，通过识别网络结构拓扑转变的特征尺度诱导粗粒化过程[132]。

1. 1. 1. 拓扑标度不变网络

当拉普拉斯特征值密度$$\omega(\lambda)\sim\lambda^\nu$$时，网络具有拓扑标度不变性，其中$$\nu=d_S/2-1$$，$$d_S$$为拉普拉斯谱维度[133]，此时熵敏感性与$$\tau$$无关，保持恒定值$$C(\tau)=d_S/2$$[35]，这种$$\tau$$无关性在规则格点、随机树及分层模块化网络[134]中出现。根据拉普拉斯重整化群定义，所有拓扑尺度不变网络在拉普拉斯尺度变换和粗粒化下保持相同拓扑结构[135]，因此$$C(\tau)$$可作为实际网络尺度不变性的探测器，常数值直接用于度量谱维度。

多尺度网络的重整化

多尺度网络重整化方法[137–139]旨在引入在任意节点聚合下保持一致的随机图模型，动机源于两方面：给定网络的不同节点划分可能因不同原因相关，需在多种聚合方式下一致建模；真实世界网络数据中“节点”含义非同质，需唯一随机图模型在任意节点聚合或拆解后保持一致性。

该要求意味着多尺度重整化方法对粗粒化变量定义完全不可知，可任意“平铺”节点集，引入任何所需划分$$\Omega_0$$，不依赖度量或扩散邻近性概念。平均掉细微细节的过程允许不同选择，但考虑到构成节点间的连接对定义超节点间联系的重要性，选择特定方案[137]。多尺度方法的独特之处在于识别重整化群流不动点，在该点可对任意节点聚合方式进行耦合系数和参数值的精确重整化。

1. 1. 1. 节点聚合下的不变性

所求的聚合不变随机图概念类似于稳定随机变量[140-142]，即组合后仍服从相同分布的随机变量。考虑以概率$$P_\ell\left(A^{(\ell)} \mid \Theta_\ell\right)$$生成具体实现$$A^{(\ell)}$$的随机图系综，其中$$\Theta_\ell$$是模型参数集合。给定概率分布和任意划分$$\Omega_\ell$$，可在高一层级$$\ell+1$$导出结果图系综由特定概率分布描述，多尺度重整化方法的核心思想是识别在所有层次和所有可能划分下保持不变的图概率分布函数形式，仅允许（重整化）参数$$\Theta_{\ell+1}$$依赖层级$$\ell+1$$，通常通过$$\Omega_\ell$$与$$\Theta_\ell$$关联，该条件对应通过在图概率空间中识别重整化群流不动点实现参数的精确重整化。

1. 1. 1. 带独立边的多尺度模型

对于带独立边的图模型，可找到显式不动点$$P\left(A^{(\ell)} \mid \Theta_\ell\right)$$，概率分布分解为节点对$$i_\ell j_\ell$$的连接概率$$p_{i_{\ell} j_{\ell}}\left(\Theta_\ell\right)$$： [math]\displaystyle{ P\left(A^{(\ell)} \mid \Theta_\ell\right)=\prod_{i_\ell,j_\ell}^n\left[p_{i_{\ell} j_{\ell}}\left(\Theta_\ell\right)\right]^{A_{i_\ell,j_\ell}^{(\ell)}}\left[1-p_{i_{\ell} j_{\ell}}\left(\Theta_\ell\right)\right]^{1-A_{i_\ell,j_\ell}^{(\ell)}} }[/math] 多尺度模型（MSM）[137]定义为连接概率空间中重整化群流的唯一非平凡不动点，满足尺度不变性要求： [math]\displaystyle{ p_{i_{\ell} j_{\ell}}\left(\Theta_{\ell}\right)= \begin{cases} 1-\mathrm{e}^{-\delta x_{i_\ell} x_{j_\ell} f\left(d_{i_\ell j_{\ell}}\right)}, & i_{\ell} \neq j_{\ell}; \\ 1-\mathrm{e}^{-\delta x_{i_\ell}^2 f\left(d_{i_\ell j_{\ell}}\right)/2}, & i_{\ell}=j_{\ell}, \end{cases} }[/math] 其中参数$$\Theta_{ij}$$分解为全局参数$$\delta>0$$调节网络整体连接密度，$$N_\ell$$个正定节点特异性参数$$\{x_{i_\ell}\}$$（适应度值）决定节点建立连接的个体倾向，以及$$N_\ell^2$$个二元参数$$\{d_{i_\ell j_\ell}\}$$表示节点对效应，函数$$f$$是任意单调正值函数。

1. 1. 1. 参数重整化与聚合不变性

使用划分$$\Omega_\ell$$将网络粗粒化至$$\ell+1$$层级时，构型$$A^{(\ell+1)}$$的概率仍由上述方程描述，其中$$\ell$$替换为$$\ell+1$$且参数重整化为： [math]\displaystyle{ x_{i_{\ell+1}} = \sum_{i_{\ell} \in \Omega} x_{i_{\ell}}, \quad d_{i_{\ell+1},j_{\ell+1}} = f^{-1}\left(\dfrac{\sum_{\substack{i_{\ell} \in i_{\ell+1}, j_{\ell} \in j_{\ell+1}}} \delta x_{i_{\ell}} x_{j_{\ell}} f(d_{i_{\ell},j_{\ell}})}{\sum_{\substack{i_{\ell} \in i_{\ell+1}, j_{\ell} \in j_{\ell+1}}} \delta x_{i_{\ell}} x_{j_{\ell}} }\right) }[/math] 其中$$\delta$$尺度不变，保持不变。当$$\{d_{i_\ell j_\ell}\}$$为超度量时，重整化后仍保持超度量性，另一特例是令$$f(d) \equiv 1$$消去所有二元效应，使模型完全由适应度驱动，关键是适应度参数不可消去，代表模型中不可约简的特征。

值得注意的是，对应配分函数$$Z(\delta)$$可精确计算，证明在任意计算层级取值相同[137]，证实多尺度模型位于重整化流的不动点上。与几何或拉普拉斯方法不同，多尺度模型既不需要节点坐标的度量空间嵌入，也不需要扩散邻近性概念，允许任意划分，聚合不变性较邻近驱动的尺度不变性更显著。

1. 1. 1. 真实世界网络重整化的淬火变体

在多尺度模型的“淬火”变体中，节点适应度解释为确定性属性，如观测量或潜在节点特征[137,139,150,151]，可用于基于经验量（如适应度模型族[152]）或潜在变量（如节点嵌入算法族[153]）对真实世界网络建模。然而，这些算法族的通用模型在固定分辨率层级构建，而多尺度模型要求通过重整化参数在所有聚合层级一致描述同一系统，是唯一在聚合后保持适应度模型特征的适应度可加性模型[137,150]，也是在任意粗粒化下生成一致节点嵌入的方法[151]。

该模型在多个层级成功复现国际贸易网络[137,139]、企业间网络[150]和投入产出网络[151]的特性，当二元和/或适应度量从数据输入时，仅剩全局参数作为自由参数，可在特定分辨率层级拟合，提供所有其他层级的预测，多尺度预测与多层级经验网络特性高度吻合[137,139,150]。

1. 1. 1. 集群无标度网络与细粒度化的退火变体

在退火变体[137,138]中，节点适应度视为随机变量，类似于非均匀随机图家族[154]，作为“隐”变量不与特定真实世界特征关联，聚合不变性要求也适用于适应度，任何层次$$\ell$$的节点适应度应从相同概率密度函数抽取，适应度的加性表明其为$$\alpha$$-稳定随机变量[141]，密度函数在大$$x$$时衰减为$$x^{-1-\alpha}$$，为确保适应度正性，指数$$\alpha \in (0,1)$$，导致适应度均值及更高阶矩发散，此范围内唯一已知$$\alpha$$-稳定分布为Lévy分布（$$\alpha=1/2$$），但对任意$$\alpha$$值可推导严格结果，且基本与$$\alpha$$无关[138]。

若将节点聚集成大小相等的块，适应度保持相同密度，仅需重新缩放全局参数。退火变体中，多尺度模型生成任意大网络，其性质可在节点数量发散的渐近极限下研究[138]。适应度的无穷均值特性使退火多尺度模型的期望拓扑性质与有限均值适应度模型不同，可严格证明[137,138]，期望度分布$$P(k)$$具有普适的幂律尾衰减$$k^{-2}$$，与$$\alpha$$值无关。

由于$$\alpha$$-稳定随机变量具有无限可分性，退火方法中可将节点无限细分为子节点，每个子节点拥有独立同分布的适应度，意味着重整化流不仅定义自底向上的半群，还定义双向进行的群[137]。如何将该理论特性转化为实用工具，细化仅在粗粒化水平观察的现实网络，仍是公开挑战。

展望

尽管网络重整化方法在构建几何与系综异质性通用框架方面取得进展，仍存在若干挑战与开放性问题。

1. 1. 1. 分辨率水平

真实世界网络实现中，拓扑结构异质性使将几何和拓扑组织归于“功能单元”的操作具挑战性，比同质格点或规则树复杂。此外，节点定义的合理分辨率尺度通常非系统内蕴，源于观测限制或数据聚合粗糙程度，如金融关系网络和流行病网络中的信息不全。理想情况下，需确定并区分网络属性固有特征尺度与数据分辨率相关尺度。

1. 1. 1. 与动力学过程相容的重整化

对任意网络上动力学结构和过程的重整化，主要挑战在于如何以与底层图结构和系综随机性重整化兼容的方式重整化动力学。此前研究限于规则或分形格点上的动态过程，粗粒化结构部分可自然定义。

具体而言，精确重整化群计算已应用于分形网络上的高斯场[155]和随机游走[156]，扩展至研究Hanoi网络上Ising模型的不动点[157]，揭示不同区域的临界非普适行为。利用裁剪技术探索网络增长模型上渗流临界行为，发现与无关联网络偏差[158]，也用于分析神经元网络活动[159]。

一般而言，真实世界网络及其粗粒化版本的结构异质性可能导致动力学过程与拓扑结构重整化间的耦合，开发独立于具体方案的框架，兼容地同时重整化网络动态及底层几何或系综结构，仍是挑战。

1. 1. 1. 广义临界性？

真实世界网络连接的强异质性与非局部性可能影响动力学模型在临界点附近的局部激发与集体激发转变性质，引发更复杂行为，需要更广义的临界性定义。例如，真实接触网络上疫情传播动力学中，集体大流行态与局部地方病流行态间的转变可能比近似格点的理想化接触网络复杂，相互作用参数可能存在完整区域，使临界行为在宏观但局部的集体子网络中逐步显现，由强拓扑异质性及潜在层级结构决定，表明拓扑复杂性可能要求将临界点概念扩展至Griffith相范畴[160,161]。

1. 1. 1. 信息论视角下的参数相关性与无关性

粗粒化意味着信息丢失，这是重整化群技术的核心工具，用于识别最有价值时空尺度的参数。重整化方法最初设计用于研究大自由度物理系统，识别参数从微观到宏观尺度的流动，跨尺度信息流的关键在于随尺度增大，某些参数相关性增强，其他逐渐无关。

在复杂网络背景下，探索这种信息随参数的自下而上流动，可为理解哪些量和控制参数在最大时空尺度最重要提供第一性原理解释。通过信息论（Fisher信息与信息几何形式体系），可从纯经验视角获得互补的自上而下方法。

1. 1. 1. 信息几何与参数流

给定任意数据集，可构建由假设模型生成的似然函数，最大似然估计识别使似然函数最大化的参数集。对数似然相对于参数的二阶导数形成度量参数空间距离的结构，捕捉似然偏离最大值的变化方式，这种推断的几何化（信息几何[162]）使在参数流统一框架下重构统计推断与重整化成为可能[163]。

基于数据集构建的似然对应的Fisher信息度量具有显著的特征值层级结构，最大特征值对应参数微小变化对观测结果影响最大，直接映射重整化群意义上的参数相关性概念[164]，解释了为何信息有限时简单模型常能捕捉大部分相关特征。

当对网络数据集或模型嵌入模拟实现不断粗粒化后，粗粒化似然函数的部分特征值相关性增强，其他减弱，与重整化群相关性精确对应。基于此，仅通过研究数据粗粒化下的参数流，即可系统推断不同相与控制参数的存在，唯一假设是粗粒化方案选择、底层几何/系综随机性的参数化方式及其先验分布设定，该方法已用于识别异质性化学系统的不同相[166]，在复杂网络背景下亦具研究前景。