使用最大纠缠态来制备MPS（with periodic boundary conditions）

前提背景

Affleck, Kennedy, Lieb, and Tasaki（AKLT）模型：
AKLT模型：是一维的延伸量子（海森堡自旋模型）
通过构造，AKLT哈密顿量的基态是带有连接每个相邻位点对的单价键的价键固体。在图形上，这可以表示为：

在这里，实心点表示自旋1/2，并处于单重态。连接自旋1/2的线是表示单线态的价键。椭圆是投影算子，它们将两个自旋1/2一起“绑定”为一个自旋1，突出自旋0或单重态子空间，并仅保留自旋1或三重态子空间。符号“ +”，“ 0”和“-”标记了标准自旋1基态（$S^z$运算符）。

正文：

对于一个具有𝑛个自由度的量子多体系统：


问题：一共$d^N$个元素这是一个巨大无比的张量计算机根本无法存储
解决方案：将其化为一个近似的特殊的张量网络表达，下面以MPS态表示
矩阵乘积态有两种不同的MPS形式，在本文中，我们将运用最大纠缠态来制备周期边界条件的MPS（PBC-MPS）如图所示：

---

注：其中每个粒子各自的取值情况数为2，分别为0，1
$a_i{b}_i$这两个粒子是两个不相关的粒子
$$\frac{1}{\sqrt{2}}(|01⟩-|10⟩{)}_{b_i{a}_{i+1}}$$

第一步

我们对这个纠缠做一个变换：
$\frac{1}{\sqrt{2}}(|01⟩-|10⟩{)}_{b_i{a}_{i+1}}=\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}0& 1\\ -1& 0\end{array}\right]}_{b_i{a}_{i+1}}(|00⟩+|11⟩{)}_{b_i{a}_{i+1}}$
证明：

---

①方法一
因为：
$|0⟩=(1,0{)}^T,|1⟩=(0,1{)}^T$
所以$|01⟩=(0,1,0,0{)}^T,|10⟩=(0,0,1,0{)}^T,|00⟩=(1,0,0,0{)}^T,|11⟩=(0,0,0,1{)}^T$
$|01⟩-|10⟩=(0,1,-1,0{)}^T,|00⟩+|11⟩=(1,0,0,1{)}^T$
将$|01⟩-|10⟩=(0,1,-1,0{)}^T$reshape为矩阵
$$\left[\begin{array}{cc}0& 1\\ -1& 0\end{array}\right]$$将$|00⟩+|11⟩=(1,0,0,1{)}^T$reshape为矩阵
$$\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$$因为：
$$\left[\begin{array}{cc}0& 1\\ -1& 0\end{array}\right]\times \left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]=\left[\begin{array}{cc}0& 1\\ -1& 0\end{array}\right]$$将$$\left[\begin{array}{cc}0& 1\\ -1& 0\end{array}\right]$$reshape为向量即为$(0,1,-1,0{)}^T=|01⟩-|10⟩$
②方法二
由①可知$|00⟩+|11⟩=(1,0,0,1{)}^T$reshape出来的是一个单位阵，与任何矩阵乘积不改变原来矩阵形式

---

令$M_{b_i{a}_{i+1}}=\frac{1}{\sqrt{2}}{\left[\begin{array}{cc}0& 1\\ -1& 0\end{array}\right]}_{b_i{a}_{i+1}}$
则原式$=M_{b_i{a}_{i+1}}(|00⟩+|11⟩{)}_{b_i{a}_{i+1}}$
因为$|00⟩+|11⟩$为单位阵，对于它的作用相当与对$b_i{a}_{i+1}$的作用
则对于整个模型而言即为：${\sum }_{b_i{a}_{i+1}}{M}_{b_i{a}_{i+1}}|b_i{a}_{i+1}⟩$

第二步

注：再次声明——这两个粒子是两个不相关的粒子。各自的取值情况数为2，因此总的希尔伯特空间维度为4。但是我们可以将它们投影到一个更小的空间上去。

①：定义投影算符 ( σ i ∈ { + 1 , 0 , − 1 } ) \left(\sigma_{i} \in\{+1,0,-1\}\right) (σi​∈{+1,0,−1}) :
$$\begin{array}{l}{P}_i={\sum }_{a_i{b}_i{\sigma }_i}{N}_{a_i{b}_i}^{{\sigma }_i}|{\sigma }_i⟩⟨a_i{b}_i|\end{array}$$其中：$$\begin{array}{l}{N}^{(+1)}=\left[\begin{array}{ll}1& 0\\ 0& 0\end{array}\right],N^{(0)}=\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ \frac{1}{\sqrt{2}}& 0\end{array}\right],N^{(-1)}=\left[\begin{array}{ll}0& 0\\ 0& 1\end{array}\right]\end{array}$$

②：将投影算符$P_i$展开观察：

$P_i={\sum }_{a_i{b}_i{\sigma }_i}{N}_{a_i{b}_i}^{{\sigma }_i}|{\sigma }_i⟩⟨a_i{b}_i|$
$=N_{00}^{(+1)}|+1⟩⟨00|+N_{01}^{(0)}|0⟩⟨01|+N_{10}^{(0)}|0⟩⟨10|+N_{11}^{(-1)}|-1⟩⟨11|$
$$\begin{gathered} =\left[\begin{array}{cc}1& 0\\ 0& 0\end{array}\right]|+1⟩⟨00|+\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ \frac{1}{\sqrt{2}}& 0\end{array}\right]|0⟩⟨01|+\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ \frac{1}{\sqrt{2}}& 0\end{array}\right]|0⟩⟨10| \\ +\left[\begin{array}{cc}0& 0\\ 0& 1\end{array}\right]|-1⟩⟨11| \end{gathered}$$

③：将投影算符$P_i$作用在粒子上
$|\Psi ⟩=P_1{P}_2\cdots P_N{\sum }_{b_1{a}_2}{M}_{b_1{a}_2}|b_1{a}_2⟩{\sum }_{b_2{a}_3}{M}_{b_2{a}_3}|b_2{a}_3⟩\cdots {\sum }_{b_N{a}_1}{M}_{b_N{a}_1}|b_N{a}_1⟩$
$$\begin{gathered} ={\sum }_{a_1{b}_1{\sigma }_1}{N}_{a_1{b}_1}^{{\sigma }_1}|{\sigma }_1⟩⟨a_1{b}_1|{\sum }_{a_2{b}_2{\sigma }_2}{N}_{a_2{b}_2}^{{\sigma }_2}|{\sigma }_2⟩⟨a_2{b}_2|\cdots {\sum }_{a_N{b}_N{\sigma }_N}{N}_{a_N{b}_N}^{{\sigma }_N}|{\sigma }_N⟩⟨a_N{b}_N| \\ {\sum }_{b_1{a}_2}{M}_{b_1{a}_2}|b_1{a}_2⟩{\sum }_{b_2{a}_3}{M}_{b_2{a}_3}|b_2{a}_3⟩\cdots {\sum }_{b_N{a}_1}{M}_{b_N{a}_1}|b_N{a}_1⟩ \end{gathered}$$
注：$⟨a_1{b}_1|⟨a_2{b}_2|\cdots ⟨a_N{b}_N|与|b_1{a}_2⟩|b_2{a}_3⟩\cdots |b_N{a}_1⟩$计算为$⟨a_1{b}_1{a}_2{b}_2\cdots a_N{b}_N|b_1{a}_2{b}_2{a}_3\cdots b_N{a}_1⟩$（其中：$⟨a_N|a_N⟩$ = $⟨b_N|b_N⟩$ =1）
所以其最后的值为1
$={\sum }_{a_i{b}_i{\sigma }_i}\left(N_{a_1{b}_1}^{{\sigma }_1}{M}_{b_1{a}_2}\right)\cdots \left(N_{a_N{b}_N}^{{\sigma }_N}{M}_{b_N{a}_1}\right)|{\sigma }_1,\cdots ,{\sigma }_N⟩$

$$\begin{array}{c}|\Psi ⟩={\sum }_{a_i{b}_i{\sigma }_i}\left(N_{a_1{b}_1}^{{\sigma }_1}{M}_{b_1{a}_2}\right)\cdots \left(N_{a_N{b}_N}^{{\sigma }_N}{M}_{b_N{a}_1}\right)|{\sigma }_1,\cdots ,{\sigma }_N⟩\end{array}$$令 $\begin{array}{c}{N}_{a_1{b}_1}^{{\sigma }_1}{M}_{b_1{a}_2}=A_{a_1{a}_2}^{{\sigma }_1},\cdots ,N_{a_N{b}_N}^{{\sigma }_N}{M}_{b_N{a}_1}=A_{a_N{a}_1}^{{\sigma }_1}\end{array}$
$$\begin{array}{c}|\Psi ⟩={\sum }_{a_i{b}_i{\sigma }_i}\left(A_{a_1{a}_2}^{{\sigma }_1}\right)\cdots \left(A_{a_N{a}_1}^{{\sigma }_1}\right)|{\sigma }_1,\cdots ,{\sigma }_N⟩\end{array}$$

第三步

由第二步中$\begin{array}{c}{A}_{a_1{a}_2}^{{\sigma }_1}=N_{a_1{b}_1}^{{\sigma }_1}{M}_{b_1{a}_2},\cdots ,A_{a_N{a}_1}^{{\sigma }_1}=N_{a_N{b}_N}^{{\sigma }_N}{M}_{b_N{a}_1}\end{array}$我们来看看算符$A$的值：
$A^{(+1)}=N^{(+1)}M=\left[\begin{array}{cc}1& 0\\ 0& 0\end{array}\right]\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ -\frac{1}{\sqrt{2}}& 0\end{array}\right]=\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ 0& 0\end{array}\right]$
$A^{(0)}=N^{(0)}M=\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ \frac{1}{\sqrt{2}}& 0\end{array}\right]\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ -\frac{1}{\sqrt{2}}& 0\end{array}\right]=\left[\begin{array}{cc}-\frac{1}{2}& 0\\ 0& \frac{1}{2}\end{array}\right]$
$A^{(-1)}=N^{(-1)}M=\left[\begin{array}{cc}0& 0\\ 0& 1\end{array}\right]\left[\begin{array}{cc}0& \frac{1}{\sqrt{2}}\\ -\frac{1}{\sqrt{2}}& 0\end{array}\right]=\left[\begin{array}{cc}0& 0\\ -\frac{1}{\sqrt{2}}& 0\end{array}\right]$
矩阵$A_{a_i{a}_{i+1}}^{{\sigma }_i}$的值取决于$M$的值，即最终MPS中各张量的关联关系与每一个纠缠态的纠缠特性有关系，而这个特性就是由$M$矩阵决定的。

注：
${\sigma }_i$是开放的indices，在一些文献中称为“fixed collection of physical indices” 。
而$a_i$很明显，是粒子的编号，在这里叫做“connecting bond indices”。