怎么计算z=x+y的概率密度_混合概率密度与拓扑不变量

前排提示：由于本人对傅立叶变换部分正负号搞的很迷糊，本文中可能有些地方正负号是有问题的，一切与论文中的图表与公式为准，等我学好了再改正！！！

本文介绍HTOF(hybrid time-of-flight)图的概念，并给出利用混合概率密度计算HTOF图的方法，并且说明其与拓扑不变量的关系。

---

一、拓扑不变量

陈数一般怎么算呢？首先由哈密顿量

得到波函数

，然后算贝里联络

，然后算贝里曲率

，最后在布里渊区积分得到

，当然还有数值的方法

Chern Numbers in Discretized Brillouin Zone: Efficient Method of Computing (Spin) Hall Conductances​journals.jps.jp

但我们还有一种直接画图看陈数的神奇方法，在文章

http://www.iop.cas.cn/xwzx/kydt/201311/P020131113340018357225.pdf​www.iop.cas.cn

中阐述了这个利用混合密度函数(hybrid density)

做出HTOF图，然后根据当

在一个

周期变化中函数变化了多少直接得出陈数。下面我们简单给出计算的过程。

二、模型

我们首先给一个简单的正方格子的Hofstadter model，其哈密顿量如下：

其中

为hopping，当电子顺时针绕小正方形四个顶点上的原子运动一圈时相位会变化

，我们取

，其中

是互质的两个整数。特别的，这里取

这样模型如图所示：

绿色椭圆圈起来的是原胞，一个原胞中有7个原子

我们首先计算一下两边都是周期性边条情况下的能带与陈数。其哈密顿量的核为：

也就是一个

的矩阵，其能带如下：

显然有7个能带，然后我们用数值算一下陈数，按照能量从低往高分别为(1,1,1,-6,1,1,1)

(*算能带代码*)

三、混合概率密度

接下来我们介绍一下混合概率密度：首先将以上模型在y方向改为有限边条，于是整个体系就有一个混合边界条件，这样在做傅里叶变换之后解哈密顿量得到的波函数为

，其概率即为

，将一个能带中的所有波函数模方求和再取平均就可以得到总的混合概率密度函数：

于是我们就按部就班地构造混合概率密度函数吧。首先得到哈密顿量（以前讲过了怎么算了，具体看程序吧），计算得到能谱：

出现了许多边界态

然后我们检查一下最低能带一共有10个态，于是把他们都平方再求和得到概率密度函数，继续利用Densityplot做出图像如下：

C=1

这个图怎么看？首先固定一个n，然后画一条直线平行于

轴，看这条直线与带相交多少次陈数的绝对值就是多少，而如果带是从左下往右上的，那么陈数是正的，反之为负。这里已经画了一个白线，只有一个相交，于是陈数为1.然后我们看一下第四能带的图像是什么样：

C=-6

从左上往右下，且有六个交点，已经画在图中了。我们还可以画一下HWF中心位置的图像，其定义为

其周期性也能看出陈数。如下：

重复了一个周期，其陈数为1

计算代码如下：

J

四、自由边界条件

当两个方向都是有限的情况，也就是所谓自由边界条件的情况下，不能将哈密顿量中的产生湮灭算符做傅里叶变换，然后再得到

了。这时我们可以解实空间的哈密顿量，然后得到实空间的波函数

。但是我们还是可以通过以下式子定义

：

所以我们只需要计算出相应的实空间的波函数，再在一个方向通过离散傅里叶变换得到相应的混合概率密度函数。

我们选取x方向

个原子，y方向上

个原子(由于一个原胞里有q个原子，故取

个原胞可以得到

个原子)，特别的，取

，首先得到mode number:

然后选取前100个位于最低能带的态，作图得到：

C=1

如上所述，图像基本一致，陈数为1.

接着我们再计算一种情况，就是当

时，这时陈数为-2.通过同样的方法我们得到下图：

C=-2

代码如下：

J

综上，我们在混合边界与自由边界情况下，利用构造混合概率密度函数得到了HTOF图，并且从图中直接得到了拓扑不变量。

参考文献：

http://www.iop.cas.cn/xwzx/kydt/201311/P020131113340018357225.pdf​www.iop.cas.cn