（5）多体量子态与统计力学基础

1.多体系统量子态与量子算符

1.1量子态系数

例：两个自旋构成的基矢为四个4维向量, 可定义为
$$|0⟩|0⟩,|0⟩|1⟩,|1⟩|0⟩,|1⟩|1⟩$$
其中, $|i⟩|j⟩=|ij⟩=|i⟩\otimes |j⟩(\otimes$ 称为直积、张量积、外积或克伦内克积,$\otimes$符号可省略 ), 例如:
$|1⟩={\left[\begin{array}{ll}1& 0\end{array}\right]}^T$
$|11⟩={\left[\begin{array}{ll}1& 0\end{array}\right]}^T⨂{\left[\begin{array}{ll}1& 0\end{array}\right]}^T={\left[\begin{array}{llcc}1& 0& 0& 0\end{array}\right]}^T$
( 等号代表左边态的系数等于右边的张量)
任意的二自旋量子态可写成基矢的线性叠加:
$$|\phi ⟩={\phi }_{00}|00⟩+{\phi }_{01}|01⟩+{\phi }_{10}|10⟩+{\phi }_{11}|11⟩=\sum _{ij=0}^1{\phi }_{ij}|ij⟩$$
二自旋量子态 $|\phi ⟩$ 的系数可看作是 $4\times 1$ 的向量[ ${{\phi }_{00}{\phi }_{01}{\phi }_{10}{\phi }_{11}]}^T,$ 或
$2\times 2$ 的矩阵 $\left[\begin{array}{cc}{\phi }_{00}& {\phi }_{01}\\ {\phi }_{10}& {\phi }_{11}\end{array}\right],$ 二者相差一个reshape操作

1.2单体算符的运算

对于$N$自旋体系，对应希尔伯特空间维数为$2^N$ ，即量子态的系数为$2^N$维张量（$N$阶张量，有$N$个指标，每个指标的维数是2，总维数是$2^N$），算符的系数为$2^N$ × $2^N$ 维张量（2$N$阶张量,有2$N$个指标，每个指标的维数是2）
图形表示:

量子态

量子算符

---

定义单体算符：作用到某一个自旋上的算符，例如泡利算符，系数维数为2 × 2

---

单体算符作用到多体量子态的规则（以三自旋系统为例）：定义在第1个自旋空间中的算子${\hat{O}}^{(1)}$（即该算子仅作用在第1个自旋上），其对应的系数维数为2 × 2 ，三自旋量子态$|\phi ⟩$对应的系数维数为2 × 2 × 2，将${\hat{O}}^{(1)}$作用到$|\phi ⟩$上的公式可写为：
$$|{\phi }^{{}^{\prime }}⟩={\hat{O}}^{(1)}|\phi ⟩={\hat{O}}^{(1)}\otimes {\hat{I}}^{(2)}\otimes {\hat{I}}^{(3)}|\phi ⟩$$
其中，${\hat{I}}^{(n)}$为定义在第n个自旋空间的单位算符（单位算符的系数矩阵为单位阵）

---

注：对于多自旋态，严格而言，无法定义对某一个自旋的单独操作，相关算符也需定义在多自旋希尔伯特空间中；${\hat{O}}^{(1)}\otimes {\hat{I}}^{(2)}\otimes {\hat{I}}^{(3)}$类似的与单位阵的直积可看作是单体算符需满足的形式。

---

在公式$|{\phi }^{{}^{\prime }}⟩={\hat{O}}^{(1)}|\phi ⟩={\hat{O}}^{(1)}\otimes {\hat{I}}^{(2)}\otimes {\hat{I}}^{(3)}|\phi ⟩$中,${\hat{O}}^{(0)}\otimes {\hat{I}}^{(1)}\otimes {\hat{I}}^{(2)}$的维数为$2^N$ × $2^N$,可以以指标收缩的形式作用到维数为$2^N$的量子态$|\phi ⟩$上(矩阵与一个向量做一个矩阵乘)
但是,我们实际上不用按上述方式进行$2^N$ × $2^N$维矩阵与$2^N$维向量的矩阵积计算,而是作如下计算:
设$|\phi ⟩$与$|{\phi }^{{}^{\prime }}⟩$的系数分别为三阶张量${\phi }_{ijk}$和${\phi }_{i^{{}^{\prime }}jk}^{{}^{\prime }}$
设${\hat{O}}^{(1)}$的系数为二阶矩阵$O_{i{i}^{{}^{\prime }}}^{(1)}$,则有如下公式:
$${\phi }_{i^{{}^{\prime }}jk}^{{}^{\prime }}=\sum _i{O}_{i^{{}^{\prime }}i}^{(1)}{\phi }_{ijk}$$
将定义在第n个自旋的算符${\hat{O}}^{(n)}$作用到自旋多体态上,仅需将算符与第n个指标进行收缩,对应的图形表示如图:

注：虽然仅进行第n个指标的收缩，但实际上，所有张量元可能被改变，并非仅有第n个指标对应的张量元发生改变，无法定义第n个指标对应的张量元，这与“无法定义对某一个自旋的单独操作”这一事实是一致的

1.3多体算符的运算

对于多体算符，当该算符可以写成多个定义在不同空间的单体算符的直积时，计算算符作用到多体态上时，仅需进行多次单体算符的作用即可

---

由于单体算符定义在不同空间，算符之间相互对易（即可以交换作用的顺序,${\hat{O}}^{(m)}{\hat{O}}^{(n)}={\hat{O}}^{(n)}{\hat{O}}^{(m)}\iff [{\hat{O}}^{(m)},{\hat{O}}^{(n)}]=0$），故作用的顺序不影响结果。

---

例：将定义在第1个和第2个自旋空间中的算符$\hat{O}={\hat{O}}^{(1)}\otimes {\hat{O}}^{(2)}$作用到三自旋态$|\phi ⟩$上,得到的量子态$|{\phi }^{{}^{\prime }}⟩$，相应的系数满足：
$${\phi }_{i^{{}^{\prime }}jk}^{{}^{\prime }}=\sum _{ij}{O}_{i^{{}^{\prime }}i}^{(1)}{O}_{j^{{}^{\prime }}j}^{(1)}{\phi }_{ijk}$$

---

一般情况下的图形表示如图：

---

如果算符不能分解成多个单体算符直积的形式，则根据分解的情况进行收缩；如果存在不同算符作用在相同自旋上，则重复上述规则，由下至上依次将各个算符所用到量子态上

---

下图给出一个图形表示示例：

1.4拓展：多体算符

如果量子算符为幺正算符（$U{U}^{†}=I$或者说算符的系数为幺正矩阵），则这些算符构成一个作用在多体态上的大的幺正操作，称之为量子线路（注：特殊情况下可不满足幺正性）
量子线路是可运行于量子计算机的模型（类似于逻辑门线路与经典计算机间的关系），张量网络为量子线路提供一个给定基底下的数学表示.

---

作用在8个自旋上，每个自旋初态为0

2.经典热力学基础

对于经典平衡态，系综理论（描述经典热力学的经典理论）的核心是：对于一个全同粒子构成的系统，该系统处于某一种状态（或构型，记为($s_1,s_2,...$)）的概率𝑷 ，由该状态的能量𝑬决定（设玻尔兹曼常数与普朗克常数为1），满足：
$$P(s_1,s_2,...;\beta )=\frac{e^{-\beta E(s_1,s_2,...)}}{Z}\dots \dots （1）$$
其中:
$\beta =1/T$（物理中为倒温度，在机器学习中为超参数）
$Z$被称为配分函数（partition function），等于所有可能构型概率之和，满足
$$Z=\sum _{s_1,s_2,...}{e}^{-\beta E(s_1,s_2,...)}$$$Z$可理解为概率的归一化因子。
式（1）在机器学习中被称为玻尔兹曼机（Boltzmann machine）。

---

热力学量即对应物理量的概率平均值:
$$O(\beta )=\sum _{s_1,s_2,...}P(s_1,s_2,...;\beta )O(s_1,s_2,...)$$其中：$O(s_1,s_2,...)$为在当前构型下对应的物理观测量的取值。
可见，建立描述给定物理系统热力学性质的关键，在于建立能量𝑬与状态之间的函数关系（$E(s_1,s_2,...)$）

---

例：
定义Ising模型（最简单的概率统计模型）：由N个Ising自旋构成一个图（graph）,如：

每个Ising自旋为图中一个节点(node)，其可取状态$s_i$为1或-1；对于给定状态，其能量满足:
$$E(s_1,s_2,...)=\sum _{⟨i,j⟩}{J}_{ij}{s}_i{s}_j$$
其中， $⟨i,j⟩$代表图中任意一对相连的Ising自旋，$J_{ij}$称为对应连接的耦合系数（图中求和的项数等于图中边的个数，每条边上对应一个$J_{ij}$）
故，Ising模型由图（节点与边）定义；当每个节点可取的状态S数大于2时（比如：状态数取1，0，-1等），模型推广为S态Potts模型；后面我们会介绍如何使用张量网络计算Ising或Potts模型热力学

注：这里并不建议将Ising模型当作一种物理上的存在，它应是用来近似描述一大类物理现象的一个数学模型，因而这里使用了概率图理论的术语来描述Ising模型；物理上的概念并不是物理实在自身，它们都是用来解释已知物理实在、预言未知物理实在的数学模型而已。

3.量子格点模型

将经典热力学推广到量子热力学

---

3.1热力学基础

描述量子系统的热力学理论，应与经典热力学理论相容。
量子系统的热力学由有限温密度算子给出，定义为：
$$\hat{\rho }(\beta )=e^{-\beta \hat{H}}/Z$$
其中， $\hat{H}$ 为系统哈密顿量（在经典力学里面哈密顿量起源于哈密顿力学）, $Z$ 为量子配分函数。
对于量子系统，给定状态（量子态）下的能量满足（能量与$\hat{H}$的关系）：
$$E\left(s_1,s_2,\dots \right)=⟨s_1{s}_2\dots |\hat{H}|s_1{s}_2\dots ⟩$$（$\hat{H}$在构型下的能量）
与经典热力学理论相同， 定义处于 $|s_1{s}_2\dots ⟩$ 的概率为$$P\left(s_1,s_2,\dots ;\beta \right)=\frac{e^{-\beta E\left(s_1,s_2,\dots \right)}}{Z}$$且配分函数满足 $Z={\sum }_{s_1,s_2,\dots }{e}^{-\beta E\left(s_1,s_2,\dots \right)}$
可见定义量子模型即定义哈密顿量

---

将能量表达式代入得量子配分函数
$$Z=\sum _{s_1{s}_2\dots }{e}^{-\beta ⟨s_1{s}_2\dots |\hat{H}|s_1{s}_2\cdots ⟩}$$根据基矢的正交完备性 ${\Sigma }_{s_1{s}_2\dots }|s_1{s}_2\dots ⟩⟨s_1{s}_2\dots |=I,$ 得:
$$Z=\sum _{s_1{s}_2\dots }⟨s_1{s}_2\dots \left|e^{-\beta \hat{H}}\right|s_1{s}_2\dots ⟩=\mathrm{Tr}\left(e^{-\beta \hat{H}}\right)$$如果证明需要运用泰勒展开加上基底完备性。
证明如下：

多个热力学量可由配分函数关于温度的导数求得，例如自由能、能量、嫡等, 因此，求解配分函数是求解热力学问题的关键一步
算符平均值可由密度矩阵计算获得。

3.2基态问题

基态：能量最低的那个态

---

当系统温度极低时 $(\beta \to \infty ),$ 系统密度算符由哈密顿量最低的本征态 (记为 $|g⟩)$ 给出, 称为系统的基态 （ground state），对应的本征值 $E_g$ 称 为基态能
$$\begin{array}{c}{\lim }_{\beta \to \infty }{e}^{-\beta \hat{H}}/Z=|g⟩⟨g|\\ \hat{H}|g⟩=E_g|g⟩\end{array}$$(注：考虑基态非简并情况; 证明可参考最大本征值的幕级数解法)
基态观测量满足： $O(\beta )=\mathrm{Tr}\left(\hat{O}{e}^{-\beta \hat{H}}/Z\right)=⟨g|\hat{O}|g⟩,$ 与量子态观测量公式一致。
基态求解即求解哈密顿量对应矩阵的最低本征态及本征值，对应于如下最优化问题（回顾：最大本征问题对应的最优化问题）：
$$E_g=\underset{⟨g\mid g⟩=1}{\min }⟨g|\hat{H}|g⟩$$