海归学者发起的公益学术平台

分享信息，整合资源

交流学术，偶尔风月

霍尔

世间谁晓士为谁

穿越横行自有规

可叹百年君未老

芸芸一岁赋芳菲

---谨以此文致敬2018

---祝量子材料新年快乐

1.引子

2018年12月17日，英国《自然》(Nature)杂志刊登了美国麻省理工学院(MIT)两个实验小组：Pablo实验组和Gedik实验组的一项实验工作[1]。这一工作是关于在双层WTe2中观察到具有时间反演对称的非线性霍尔效应之最新实验结果，引起轰动！

乍一看，看君也许会有些狐疑：霍尔效应只是一个经典电磁学的现象。对一些相对复杂的半导体体系，霍尔效应也许不是那么好的线性关系，会展现出非线性霍尔效应。反常霍尔效应可能也算其中一种。这里所谓非线性是什么节奏？事实上，我们这里要接触的话题是：具有时间反演对称的非线性霍尔效应。几个关键词：时间反演、非线性、霍尔，叠加在一起，估计就不再是线性叠加了，而是某种崭新的现象与物理。

笔者经过一番纠结迷糊，好像有点明白什么是非线性霍尔效应了！可以夸张一点说：这是一个全新的研究方向！有了这个宣言，笔者就敢于撰写本文，评述一番什么是非线性霍尔效应，并在科普层次上触及最新理论进展[2]。

2.霍尔效应

我们在高中就学习过霍尔效应：在磁场(磁感应强调B)作用下，电流I感受到垂直于磁场和运动方向的洛伦兹力，会发生偏转，并在样品两端形成电势差V。这一现象称为霍尔效应(如图1所示)。

图1.霍尔效应原理示意图。三个字：特简单！

霍尔效应一直是凝聚态物理研究的一个主流方向。迄今为止，已经有四次诺贝尔物理学奖和霍尔效应有关(图2)，分别是：1985年量子霍尔效应[3]，1998年分数量子霍尔效应[4,5]，2016年拓扑物相和拓扑相变[6]，2010年获奖的石墨烯[7]之重要实验证据也是半整数量子霍尔效应[8]。当然，这些获奖成果之深度与广度不再是“霍尔效应”这四个字所能容纳与表达的，其中物理之妙与美远超霍尔效应本身，或者说霍尔效应只是这些新的物理最简单的一种表达而已。这么说的底气，很容易从几个方面来展示。

图2.和霍尔效应有关的诺贝尔物理学奖。

3.几何相位

0有关，大小为360o(1-θ0)。可以证明，这个角度正好等于“傅科摆每天随地球自转走过的路径所覆盖的面积”除以“地球半径的平方”。注意到，圆球球面面积除以半径的平方，就是空间角。因此，我们可以说，傅科摆一天变化的角度竟然等于单摆随地球自转所覆盖的空间角。

这一等价关系颇为神奇，是几何相位在经典系统中最直观的例子。

图3.上：北京天文馆的傅科摆。下左：球面上的圆环代表傅科摆随着地球自转一圈走过的路径，θ0是傅科摆所在的维度。圆环所覆盖的面积除以地球半径的平方，等于傅科摆每天摆动方向改变的角度。下右：量子力学经典例题“电子自旋磁矩在绝热转动的外磁场B中积累的几何相位”[参考,IntroductiontoQuantumMechanics，第10章]。

4.量子几何相位

其次，笔者假定自己已经理解了量子力学，虽然这个假定可能是一个几率(^_^)。我们来解释一下什么是量子中的几何相位。1984年，英国物理学家MichaelBerry撰写了一篇文章[9]，阐明了量子力学系统中的几何相位问题。这篇文章很像扔在水中的一颗小石头，激起的涟漪深远而美妙。

Berry在文章中举了一个例子：电子自旋磁矩在绝热转动的外磁场中会积累几何相位(图3下右)。自旋是电子携带的角动量，而角动量是描述转动的物理量。简单地说，一个物体旋转的快慢可用角动量大小来描述。带电粒子，比如电子，自转的时候会产生一个磁矩，就像环形电流产生磁场一样。如果再施加一个外磁场，电子的自旋磁矩会沿着磁场方向极化。我们熟知的电动机就是类似的原理。简单地说，如果外加一个磁场，磁场转动足够慢的话，电子自旋会一直跟着磁场方向转动，并且会积累一个量子力学“相位”。注意，这是量子力学相位！

图4.文中几个物理量的数学定义。

众所周知，相位是描述波的一个概念。比如湖面上的水波，我们可以用振动频率、振动幅度和振动相位来描述。说得粗鲁一点，量子力学的波动学说就是用波的那一套来描述微观粒子。Berry的文章指出：电子的自旋磁矩方向被磁场拉着转一圈后，会多出一个相位，其大小为1/2个空间角，和傅科摆类似。特别是，这个相位γ可用一个优美的面积分表达出来(前方高能，数学较多)，如图4之方程(1)所示，其中S就是转动的磁场在单位圆上划出的面积、ds是面积微分元。注意，被积分函数Ω被称为Berry曲率(Berrycurvature)，其定义由图4中的方程(2)给出。这里|ψ是自旋沿着磁场方向的某个瞬时本征量子态(根据其平行或反平行于磁场方向，称作自旋上或自旋下)，▽是关于磁场的三维矢量微分算符。

5.霍尔效应的几何来源

讲了这么多，数学已经看到两眼发黑，而我们才开始接触霍尔效应的几何本源。首先要回答的是：这个几何相位和霍尔效应有什么关系？

当不考虑无序的时候，霍尔电导σH的主要贡献可以写成类似的对Berry曲率的积分，也就是图4中的方程(3)，其中e是电子电荷、h是普朗克常数，对N求和表示有多少能带被电子填充。

Thouless等人1982年发现了这件事情，了不起！[6]。他们发现，当图4之(3)式中的积分包裹整个球面的时候(可借助图3去想象一下，能不能想得通是能不能学物理的一个判据^_^)，给出的结果刚好是1/2个圆球的空间角,也就是2π。由此，所谓的霍尔电导就量子化为图4中的方程(4)，其中e2/h被定义为量子电导、n必定是整数。

注意到，这里包裹球面的积分数学上实际上是一个拓扑数，称作拓扑陈数。“陈”代表著名数学家陈省(xing三声)身。这是拓扑这门数学和凝聚态物理学的一次优美的合作。

我们通常了解的霍尔效应是需要磁场激发的。最近十五年，Berry曲率相关的这套数学描述被推广到没有磁场的本征材料体系[10,11]，导致量子自旋霍尔效应的发现[12-14]、三维拓扑绝缘体的发现[15-18]、量子反常霍尔效应的发现[19,20]、拓扑半金属的发现[21,22]、及至整个拓扑物理学的广阔领域。所以，我们说“霍尔效应”这个名词现在看起来太沉重了，要负载这么多新物理。

事实上，现在开始有人用“拓扑物理学”来代替这个名词，以还原霍尔效应本来的面目。拓扑物理学是凝聚态物理的一个主要前沿方向，已无异议。当今，研究凝聚态物理或多或少都会涉及到拓扑的问题，这些概念也被推广到光学和原子分子物理。未来之路能走多远、还有多少意外发现，实未可知。

6.Berry曲率和弯曲空间

有了Berry曲率(图4之式(2))这一概念和物理图像，我们就有了强大的手段，可以深入拓扑量子物理之领地了。这也表明Berry曲率这个概念之重要，而Berry老先生耄耋之年依然未能等到斯德哥尔摩的电话是不大公平的。

接下来，我们通过一个例子来展示Berry曲率的魅力。对处于电磁场中的固体而言，其中的带电粒子运动时，Berry曲率描述了动量空间的弯曲。这有点类似广义相对论在描述引力场时使用时空弯曲概念的味道(图5)。其意大概是说固体中的Berry曲率可以使动量空间弯曲。由此一来，电子看起来是直着行走，但实际上却走偏了，从而产生横向电荷积累和电压，也就是霍尔效应。很显然，这里的霍尔效应产生动力不再是外加磁场或内禀的固有磁矩，而是某种似乎可以等效为磁场的量子效应。

更重要的是，这种弯曲可以通过很多物理效应来测量[23-25]，而这些效应亦或对发展未来电子器件有指导意义。这就是量子霍尔效应博大精深之处！

图5.Berry曲率和广义相对论的时空弯曲的类比。

7.线性霍尔效应和非线性霍尔效应

行文至此，绕了很大一个圈圈，我们终于可以进入本文的主题，即非线性霍尔效应。

目前，我们所知道的霍尔效应都是线性的(图6)。线性霍尔效应测量方法是：在纵向施加一个电流，测量横向电压。为了滤波去除噪声，实验都会使用锁相放大器(Lock-inmeasurements)，意思是说纵向电流按照一个很慢的频率振荡，比如15赫兹振荡，测量横向电压的时候只把15赫兹的信号挑出来。这样一来，无关的噪音就被过滤掉。这里，我们用到了频率的概念！

所谓非线性霍尔效应，意思是说存在高频的霍尔信号。一种测量方法是挑选2倍频率测量横向电压，比如30赫兹(图7)，这样测得的霍尔信号称之为非线性霍尔效应。

好吧，从一般的物理动力学角度看，我们总是认为一个系统对激励的响应以线性响应为主。比如，很多动力学理论几乎无一例外都是去计算线性响应函数或者线性关联函数，从而获得动力学规律。非线性响应在物理上无非是将响应过程展开为激励量的泰勒级数，由此就有了二倍频、三倍频、四倍频等非线性响应。这一逻辑在非线性光学中被用到极致。现在，物理人将这个倍频概念用到了霍尔效应，美其名为“非线性霍尔效应”。

图6.各种线性霍尔效应。上：霍尔效应、自旋霍尔效应、反常霍尔效应[27]。中：量子霍尔效应[3]、量子自旋霍尔效应[12-14]、量子反常霍尔效应[19-20]。下:三维量子霍尔效应[26,27]。

图7.非线性霍尔效应的测量。

8.对称性破缺分析

那么什么情况下可以测量到非线性霍尔效应呢？或者说，什么物理使得倍频非线性霍尔效应很强？遗憾的是，要说清楚其中的子丑寅卯，我们需要一点数学和对称性分析，看君耐心一点，其实问题并不复杂。

前面提到，霍尔效应是Berry曲率的积分。Berry曲率有两个对称性性质：一是在时间反演下变号、二是在空间反演下不变号[28]。线性霍尔效应是Berry曲率的积分，如果有时间反演，积分的时候正动量和负动量部分会抵消。所以，要观测到霍尔电导，必须时间反演破缺。这就是为什么必须施加磁场或掺杂磁性杂质(看君一定明白，B和M都是破坏时间反演对称的)才能看到霍尔电导，而非磁性无磁场的自旋霍尔效应很难被观测到。

不同于线性霍尔效应，MIT的Sodemann和付亮[29]以及明尼苏达大学的Low[30]发现：非线性霍尔效应正比于一个被称作Berry曲率偶极子的量(简称Berrydipole)，其表达式为：

这一表达式物理上了不得，内涵非常丰富！例如，它有两种等价的表达方式：一是在费米面上的群速度乘以Berry曲率(第一个等号)，一是对费米面下所有Berry曲率的一阶导数求和(第二个等号)。第二个等号告诉我们，非线性霍尔效应是Berry曲率的高阶效应，所以是新物理，非常非常重要。

不过，第一个等号可以帮我们分析什么时候可以观测到非线性霍尔效应(图8)。

首先，在时间反演下，速度和Berry曲率都反号，因此速度和Berry曲率的乘积不反号，这样我们在有时间反演的时候也有可能观测到非线性霍尔效应，这和线性霍尔效应需要破缺时间反演很不同。

其次，在空间反演下，速度反号，Berry曲率不反号，这样速度和Berry曲率的乘积反号，积分的时候正动量和负动量的部分会抵消。如果想观察到非线性霍尔效应，必须破缺空间反演。

总结一下：(1)线性霍尔效应需要破缺时间反演才能观测到，比如施加磁场或借助磁矩。(2)非线性霍尔效应需要破缺空间反演才能被观察到。

图8.对称性分析。非线性霍尔效应在时间反演破缺时也可能存在，非线性霍尔效应需要破缺空间反演。

9.实验与理论进展

最近有多个实验小组报道了具有时间反演的非线性霍尔效应的观测结果，包括MIT的Pablo组和Gedik组的双层WTe2实验(并列一作为马琼、徐苏扬、沈汇涛)[1]、Cornell的Mak和Shan实验组的多层WTe2实验[31]及Dzsaber等人的实验[32]。

等一下，这里又是WTe2？看君一定在很多场合都看到这个材料。它最近的确挺火：三维块材是第二类Weyl半金属[33,34]，单层体系可能会有高温量子自旋霍尔效应[35-39]或者超导[40,41]，多层体系有铁电性[42]。

在MIT的Pablo组和Gedik组有关双层WTe2

10.无序贡献的非线性霍尔效应

参考文献

1.(online,2018).

2.,,*,2.,266601(2018).

3.,494(1980).

4.,,,1559(1982).

5.:,1395(1983).

6.,,,,405(1982).

7.,197(2005).

8.',201(2005).

9.(1802),45(1984).

10.2,146802(2005).

11.,,,,136602(2005).

12.,226801(2005).

13.,,andS.-,1757(2006).

14.önig,,üne,,,,X.-,andS.-,766(2007).

15.,045302(2007).

16.,,,,,,,970(2008).

17.,C.-,X.-,,,andS.-2Se3,Bi2Te3,andSb2Te3,438(2009).

18.,,,,,,,,,,,398(2009).

19.,,H.-,S.-,,,61(2010).

20.C.-,,,,,,,,,L.-,Z.-,,,,,,,S.-,,,,X.-,andQ.-,167(2013).

21.[21],205101(2011).

22.3As2.,125427(2013).

23.,739(2012).

24.,,andH.–*.Negat,166601(2017).

25.H.––,127201(2017).

26.,H.-,H.–*,,136806(2017).

27.,,,,,,,,,,,,,,H.-,,3As2.Nature(online,2018)..

28.,,,1959(2010).

29.,216806(2015).

30.,,,235447(2015).

31.,,,,2.arXiv:1809.08744.

32.:1811.02819.

33.,,,,,,495(2015).

34.,,,,,,,,,,,H.-,,,,type-IIWeylsemimetalWTe2.NatureCommunications7,13142(2016).

35.,,,,1344(2014).

36.,,,,,,,,,,,,J.-,,’-WTe2.AdvancedMaterials28,4845(2016).

37.,,,,,,,,,2.NaturePhysics13,677(2017).

38.,,,,H.-,,,,,,,,,,,C.-,,,,,,,,,S.-,andZ.-’-WTe2.NaturePhysics13,683(2017).

39.,,,,,,,76-79(2018).

40.,,,,,,,,,,922(2018).

41.,,,,,,,,,926(2018).

42.,,,,,,,,,336(2018).

43.,,,087206(2009).

44.,,,,,,217203(2015).

45.,,H.–*,:1812.08377.

备注：

1.题头小诗乃Ising添加。文中得罪和挪喻之语及文法错误归于Ising，与笔者无关。

2.本文科普素材由卢海舟老师提供，本文也采用了若干卢海舟老师课题组的图表。