本文介绍了自动微分技术在张量网络中的应用，详细阐述了如何利用自动微分来优化张量网络收缩算法，包括张量重整化群（TRG）和角转移矩阵重整化群（CTMRG）方法。通过反向自动微分计算高阶导数，实现了张量网络状态的高效优化，并在自由能高阶微分和iPEPS基于梯度的优化中展示了其能力。同时，文章讨论了稳定后向传播的策略，如通过洛伦兹展宽处理数值稳定性问题，以及使用checkpointing和不动点迭代来节省内存。
本文为文献 Differentiable Programming Tensor Networks 的原创译文，转载请注明出处。
Author：小岚岚
目录
I.引言
张量网络是研究经典统计物理和量子多体物理问题的重要方法[1-3]。近年来，它的应用迅速扩展到不同的领域，包括量子电路的仿真与设计[4-7]、量子误差校正[8,9]、机器学习[10-14]、语言建模[15,16]、量子场论[17-20]、nolographic二象性[21,22]。
张量网络在一般情况下的优化是众多研究方向的核心问题之一。尽管一维矩阵积状态的优化方案非常成功[23-28]，但二维或高维张量网络的优化一直是一个具有挑战性的课题。这个困难部分原因是张量收缩的计算成本高，还要部分原因是在高维情况下缺乏有效的优化方案。对于无限平动不变量子系统，优化张量网络状态的挑战尤其紧迫。在这些情况下，同一个张量以多种方式影响变分能量，从而导致高度非线性的优化问题。基于近似虚时算符的优化方案[29-33]努力解决目标函数中的非局部依赖性。文献[34,35]采用基于梯度的优化，结果明显改善。
然而，即使对于一个简单的物理哈密顿量，从解析的角度推导张量网络状态梯度也很麻烦，而且容易出错。这一事实限制了基于梯度的张量网络状态优化在更复杂系统中的广泛应用。替代方法如使用数值导数计算梯度，精度和效率有限;因此，它们只适用于变化参数较少的情况[36,37]。使用链式法则手动推导梯度只适用于精心设计的简单张量网络结构[38]。
可微编程通过以完全可微的方式组合整个张量网络程序并使用自动微分来优化它们，为这些问题提供了优雅而有效的解决方案。在本文中，我们提出基本的自动微分技术，计算张量网络程序的(高阶)导数有效地数值精度。这一进展为一般情况下基于梯度(甚至基于hessian)的张量网络状态优化打开了大门。此外，计算张量网络算法输出的(高阶)梯度提供了一个简单的方法来计算物理量，如比热和磁化率。可微规划方法不知道详细的格点几何、哈密顿量和张量网络收缩方案。因此，该方法的通用性足以支持广泛的张量网络应用。我们的重点是涉及二维无限张量网络的应用，与传统方法相比，可微规划技术提供了显著的好处。我们证明，在通过奇异值分解(SVD)解决数值稳定微分和不动点迭代的内存高效实现等重大技术挑战后，可以获得张量网络状态变分优化的最先进结果。
本文的组织结构如下。在第二节中，我们在张量网络算法的背景下引入自动微分，并将张量网络收缩表示为计算图。在第三节中，我们给出了张量网络算法的稳定和可扩展可微规划的关键技术。在第四节中，我们证明了该方法的能力，并将其应用于无限正方晶格上的经典Ising模型和量子海森堡模型。最后，我们在第五节对未来的研究方向进行了展望。我们的代码实现是公开的(参见参考[39])。
II．一般理论
A．自动微分
通过计算图自动微分是一个统一的框架，涵盖了机器学习训练神经网络，优化量子物理张量网络等等。我们首先回顾自动微分的核心思想，然后解释它在以计算图形式表述的各种张量网络收缩算法中的应用。
自动微分机械地计算用计算机程序[40]表示的计算过程的导数。与数值微分不同，自动微分计算的导数值精确到机器精度。其性能有、普适的理论保证，不超过原程序的算法复杂度[41,42]。
图1所示。计算图的反向自动微分。黑色箭头表示从输入到输出的正向评价函数。红色箭头表示伴随反向传播的后退步骤。(a)链图。(b)一个更一般的计算图。在向后传递过程中，给定节点的伴随点根据公式(2)计算。
自动微分是现代深度学习应用的计算引擎[43,44]。在量子最优控制[45]和量子化学计算中也有应用，如计算力[46]，优化基参数[47]等。自动微分的中心是计算图的概念。计算图是由基本计算步骤组成的有向无环图。图的节点表示数据，数据可以是标量、向量，矩阵，或者张量[48]。图连通性表示计算过程中数据流的依赖性。最简单的计算图是图1(a)所示的链。例如，从输入参数向量theta开始，可以计算一系列中间结果，直到最终的输出L，我们假设它是一个标量。所谓的正向求值就是按顺序遍历链图，即 θ → T 1 → . . . → T n → L \theta \rightarrow T^1\rightarrow...\rightarrow T^n \rightarrow \mathcal{L} θ→T1→...→Tn→L。
为了计算目标函数相对于输入参数的梯度，我们可以利用链式法则
∂ L ∂ θ = ∂ L ∂ T n ∂ T n ∂ T n − 1 . . . ∂ T 2 ∂ T 1 ∂ T 1 ∂ θ . ( 1 ) \frac{\partial \mathcal{L} }{\partial \theta}= \frac{\partial \mathcal{L} }{\partial T^n} \frac{\partial T^n }{\partial T^{n-1}}... \frac{\partial T^2 }{\partial T^1} \frac{\partial T^1}{\partial \theta}. (1) ∂θ∂L=∂Tn∂L∂Tn−1∂Tn...∂T1∂T2∂θ∂T1.(1)
由于我们考虑的情况是，输入维数大于输出维数，在Eq.(1)中通过使用一系列雅可比向量乘积从左到右相乘来计算梯度是更有效的。根据图1(a)所示的计算图，将图向后遍历，并将梯度信号从输出端传播回输入端。用这种方法计算导数叫做反向自动微分。这种方法通常称为反向传播算法[43]，被认为是训练深度神经网络最成功的方法。引入伴随变量 T ˉ = ∂ L / ∂ T \bar T=\partial \mathcal{L} /\partial T Tˉ=∂L/∂T 来表示最终输出L相对于变量T的梯度。我们可以看到，反向自动微分将伴随变量从 T n ‾ = L ˉ ∂ L / ∂ T n \overline {T^n} =\bar \mathcal{L} \partial \mathcal{L}/\partial T^n Tn=Lˉ∂L/∂Tn， L ˉ = 1 \bar \mathcal{L}=1 Lˉ=1
一直传播到 T i ‾ =
最低0.47元/天 解锁文章
864
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
