张量网络的构成与MPO模型的转化

问题：如何进行下图中的转换：

在 “TEBD算法计算一维格点模型基态”以及 “MPO与一维热力学计算”中我们都提出了以上的问题，现在让我们来看看这个问题是如何解决的。

---

通过研读论文“Linearized Tensor Renormalization Group Algorithm for Thermodynamics of Quantum LatticeModels”（量子点阵模型热力学的线性张量重整化群算法（LTRG））我们将对整个LTRG算法会更加的熟悉。
LTRG算法:采用了无限时间演化的分组抽取技术，使我们能直接处理转移矩阵，使其更具有可扩展性。

---

step1：得到张量网络

①首先，我们给出一维量子多体模型中哈密顿量的形式：
$$\begin{array}{rl}H& =\sum _{i=1}^N{h}_{i,i+1}=H_1+H_2\\ H_1& =\sum _{i=1}^{N/2}{h}_{2i-1,2i},H_2=\sum _{i=1}^{N/2}{h}_{2i,2i+1}\end{array}$$我们将$H_1,H_2$展开可知：
$$\begin{gathered} H_1=h_{1,2}+h_{3,4}+h_{5,6}+...+h_{N-1,N} \\ {H}_2=h_{2,3}+h_{4,5}+h_{6,7}+...+h_{N,N+1} \end{gathered}$$即：$H_1$为作用在奇偶上的阿密顿量的和，$H_2$为作用在偶奇上的阿密顿量的和（即将哈密顿量按格点编号分为两部分）
其中，$N$为格点的个数且为偶数
我们添加$2K$（$K$很大）个完备基：
$$\left\{|{\sigma }_i^j⟩\right\}\left({\sigma }_i^j=1,\cdots ,D\right)$$其中，$i$为格点指标， $j$为$Trotter$指标。
②于是我们可以将该模型的配分函数表示为：
Z N ≃ Tr ⁡ [ e − β H 1 / K e − β H 2 / K ] K = ∑ { σ i j } ∏ j = 1 K ⟨ σ 1 2 j − 1 … σ N 2 j − 1 ∣ e − β H 1 / K ∣ σ 1 2 j … σ N 2 j ⟩ × ⟨ σ 1 2 j … σ N 2 j ∣ e − β H 2 / K ∣ σ 1 2 j + 1 … σ N 2 j + 1 ⟩ \begin{aligned} Z_{N} & \simeq \operatorname{Tr}\left[e^{-\beta H_{1} / K} e^{-\beta H_{2} / K}\right]^{K} \\ &=\sum_{\left\{\sigma_{i}^{j}\right\}} \prod_{j=1}^{K}\left\langle\sigma_{1}^{2 j-1} \ldots \sigma_{N}^{2 j-1}\left|e^{-\beta H_{1} / K}\right| \sigma_{1}^{2 j} \ldots \sigma_{N}^{2 j}\right\rangle \\ & \times\left\langle\sigma_{1}^{2 j} \ldots \sigma_{N}^{2 j}\left|e^{-\beta H_{2} / K}\right| \sigma_{1}^{2 j+1} \ldots \sigma_{N}^{2 j+1}\right\rangle \end{aligned} ZN​​≃Tr[e−βH1​/Ke−βH2​/K]K={ σij​}∑​j=1∏K​⟨σ12j−1​…σN2j−1​∣∣∣​e−βH1​/K∣∣∣​σ12j​…σN2j​⟩×⟨σ12j​…σN2j​∣∣∣​e−βH2​/K∣∣∣​σ12j+1​…σN2j+1​⟩​该式还可以继续展开为：
Z N ≃ ∑ σ i j ⟨ σ 1 1 ⋯ σ N 1 ∣ e − β H 1 / K