针对这个难题，中国科学院理论物理研究所张潘课题组系统地发展了任意图的张量网络方法[1-4]，涵盖了从严格到近似、从固定线路到通用的张量网络网络的众多情形。进一步，张潘研究员及其硕士生王一佳、伦敦大学学院博士生张钰雯与新加坡国立大学博士后潘峰（此前为中国科学院理论物理研究所博士生）于近期提出了张量网络消息传递（TNMP）算法。TNMP方法结合了张量网络和消息传递算法的优点，同时解决了短圈和长圈结构所带来的问题。具体来说，TNMP算法在处理局部的多短圈子图时采用张量网络进行缩并，对于全局的多长圈环境则采用消息传递算法进行近似，既利用了张量网络在处理短圈结构时的优势，又保留了消息传递在处理长圈时的高效性，大幅提升了现有算法在宏观物理量计算上的精确度上限与规模上限。TNMP算法也可以被视为终极的圈展开平均场方法，利用张量网络方法严格地处理系统中圈的效应，圈的大小远远超过了例如Plefka展开，Kikuchi变分自由能方法等传统的圈展开平均场理论。
TNMP方法具体展示如下：如图1所示，在计算变量i的单点磁化率（即边缘概率）时，TNMP将相互作用系统对应的张量网络划分成两个部分：虚线圈中彩色的部分所对应的中心点i的“邻域”，以及圈外灰色部分所对应的“环境”。其中，邻域中的点与中心点有着更强的关联，于是在边缘概率的计算中扮演着更重要的角色，同时，其作为网络的一个局部有着足够小的规模，可以在可接受的计算开销下被精确处理。因此，TNMP使用张量网络缩并来对邻域进行严格地计算。而对于邻域之外的环境，虽然其所包含的信息由于关联的衰减而在边缘概率的计算中产生更小的影响，但获得严格的环境却具有和全局宏观量计算几乎一致的计算复杂度，故而具有很大的近似空间。基于此，TNMP沿用了消息传递的架构对环境进行近似——用邻域边界上所有“空腔”张量（图中的紫色圆圈）的直积来近似环境所对应的指数规模的张量，并通过迭代过程来得到满足邻域边界条件的空腔张量。如图2所示，通过数值实验可以看到，TNMP在人工合成网络与现实世界网络上均展现出了显著优于MCMC算法、传统消息传递算法与近期提出的含短圈修正的消息传递算法[5-6]的计算能力。
图1：消息传递算法(a)、张量网络方法(b)、张量网络消息传递算法计算单变量边缘概率(c)与计算空腔张量(d)的图示。
图2：TNMP在人工合成网络[7]上的铁磁Ising模型（左图）与现实世界网络[8]上的自旋玻璃模型（右图）中均表现出了优越的计算精确度，并有能力在可接受的计算成本下达到机器精度。
除了在多种图结构下的统计力学问题中展现出显著的优势之外，TNMP更是可以直接应用于贝叶斯推断问题、包括渗流问题与社区发现问题在内的网络科学问题、经典或量子的冗余码纠错问题与组合优化问题等广泛的领域与场景之中。同时，通过将张量网络上的代数延拓到更丰富的代数类型，TNMP可以在量子多体计算、复本对称性破缺下的近似计算等更大的舞台上发挥自己的潜力。论文于近期发表于Physical Review Letters并被选为编辑推荐。
https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.132.117401
论文代码：
https://github.com/YijiawWang/TNMP
参考文献
[1] F. Pan,P. Zhou,S. Li,and P. Zhang,Contracting arbitrary tensor networks: General approximate algorithm and applications in graphical models and quantum circuit simulations,Phys. Rev. Lett. 125,060503 (2020).
[2] F. Pan and P. Zhang,Simulation of quantum circuits using the big-batch tensor network method,Phys. Rev. Lett. 128,030501 (2022).
[3] F. Pan,K. Chen,and P. Zhang,Solving the sampling problem of the sycamore quantum circuits,Phys. Rev. Lett. 129,090502 (2022).
[4] X. Xu,S. Benjamin,J. Sun,X. Yuan,and P. Zhang,A herculean task: Classical simulation of quantum computers,arXiv:2302.08880.
[5] G. T. Cantwell and M. E. Newman,Message passing on networks with loops,Proc.Natl. Acad.Sci.U.S.A. 116,23398 (2019).
[6] A. Kirkley,G. T. Cantwell,and M. E. J. Newman,Belief propagation for networks with loops,Sci.Adv. 7,eabf1211 (2021).
[7] S. A. Williamson and M. Tec,Random clique covers for graphs with local density and global sparsity,in Uncertainty in Artificial Intelligence (PMLR,Cambridge,2020),pp. 228–238.
[8] T. A. Davis and Y. Hu,The university of florida sparse matrix collection,ACM Trans. Math. Softw. 38,1 (2011).
