有篇近作,7月2日投出,8月17日接受! European Physical Journal C ,算来也是物理学顶级一区刊物! 论文的基本框架都在这里: General covariant geometric momentum, gauge potential and a Dirac fermion on a two-dimensional sphere https://arxiv.org/abs/1905.07735 这篇论文不仅我们基本不懂,而且审稿人编辑肯定也不懂。这不妨碍我们一起莫名的激动并赞赏这件工作! —————————— Cover Letter Dear Editors of The European Physical Journal C, It is well-known that curved manifold can induce somegauge structure, but how this gauge structure appears, e.g., for a Dirac fermion on the curved graphene, and how it affects the motion of the particle,is not clear. Our approach starts with that the symmetry expressed by the Poisson or Dirac brackets in classical mechanics preserves in quantum mechanics ; and so the Hamiltonian is determined by the symmetry, and we establish a general framework to explore these problems. As two examples, we show that for the Dirac fermion on a two-dimensional sphere, the gauge field enters the momentum, but does not exhibit in the Hamiltonian. The details arein the following. In2011, we proposed a geometric momentum for non-relativistic particles constrained on the curved surface (PRA 84, 042101(2011)) as the proper form of momentum. This concept together with geometric potential is now widely accepted for it has stood both theoretical examinations and experimental testifications. Present paper develops such a geometric momentum to be general covariant, otherwise it is inapplicable to spin particles on curved surface. First of all, we need to establish a quite general formalism to accommodate the general covariant geometric momentum (GCGM). This is in fact the framework of Dirac canonical quantization, in which we make some close examinations, emphasizing the fundamental role of symmetry indicated by both the fundamental quantum conditions and dynamical ones. It has been recognized, through construction of quantum field on curved manifolds,that some gauge structures emerge, our GCGM makes the gauge structure transparent. Two important and interesting demonstrations of the GCGM are: 1, the generalized total angular momentum for a Dirac fermionon two-dimensional sphere, which reproduce the same result as that given by 2003 Nobel prize laureate Abrikosov via the purely consideration of dynamics of the particle ; 2, the GCGM together with the dynamical quantum conditions clears up a controversial issue as whether there is geometric potential for the Dirac Fermion on the two-dimensional sphere, and our answer to this question is that no geometric potential is acceptable. The problem we discuss in present paper is of broad interests, and present work is of fundamental contributions into physics. It relates to the topologic state on surface of the low-dimensional materials, and gauge theory of gravitation, and fundamental problem in quantum mechanics. So, it is appropriate for publicationin your journal. Sincerelyyours, Q.H. Liu ———————————— 一审 一审意见 Referee: 1 Comments to the Author The paper by Liu et al. deals with the definition of the general covariant geometric momentum, which is the first part of the paper, and the study of Dirac equation on a sphere, which cover the second part.. The paper is interesting, and deserves publication . I have just some minor points: 1) May be the authors wish to enlarge the presentation of Eq. (1). 2) From Eq. (1) it appears that the canonical commutation relation are modified. They have on the right handed side a quantity different from 1. This should have non trivial consequences, I suppose. Is it? 3) To better understand the model, the author should also comment the meaning of (1) in the case in which the index i is different from index j 4) The construction studied in the paper has any connection with the generalized uncertainty principle? Referee: 2 Comments to the Author The paper deals with the quantization of generally covariant systems when constrained on submanifolds of flat spaces. After a general introduction to this problem, and offering quite an extensive list of cogent references, it faces a specific problem, namely whether Dirac fermions on a sphere do carry the so-called 'geometric potential'. In other words, whether the 'standard' quantization in curved space techniques (leading there to their eq (28), or the Hamiltonian (29)), need be corrected when the two dimensional manifold is thought of as coming from constraints in a larger three-dimensional flat space. I found the paper sound, the calculations seem correct, the logic is convincing. I have two possibly minor concerns: i) could they improve a little the English, by checking some expressions (e.g., right in the Abstract 'Once the momentum is made general covariance'...), and by making their sentences shorter, sharper and clearer? ii) since they cite work where Dirac fermions are studied on 2d surfaces of negative curvature, could they address this problem here too? or at least, comment on that? Besides this, I believe that the paper deserves publication in EPJC. —————— 接下来,论文被接受! Dear Prof. Liu: It is a pleasure to accept your manuscript entitled Generally covariant geometric momentum, gauge potential and a Dirac fermion on a two-dimensional sphere in its current form for publication in the European Physical Journal C. The paper will also be published electronically under the journal's title. Thank you for your fine contribution. On behalf of the Editors of the European Physical Journal C, we look forward to your continued contributions to the Journal. Sincerely, Prof. Kostas Skenderis Editor in Chief, European Physical Journal C Associate Editor Comments to the Author: (There are no comments.) Referee: 1 Comments to the Author (There are no comments.)
# 编者信息 熊荣川 明湖实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz More often than not, when a phylogenetic dataset is divided into smaller partitions, each one gives rise to trees that have different topologies. One can draw two conclusions from this result: (1) one, some or all of the trees are wrong and the partitions share the same history, (2) or the trees are correct and the different partitions have experienced distinct evolutionary histories. Distinguishing between these options requires statistical testing to determine if the differences in topology are likely to have been observed simply by chance. There are many different tests of incongruence available in the field of systematic biology that all use comparable measurements and ideas. However, tests differ, sometimes subtly and sometimes drastically, in their assumptions, implementation, and interpretation. These details can be difficult to discern in the disjointed literature and controversies surrounding these tests. Incongruence tests may be broadly classified into tests that consider character information (character incongruence) and those that only consider tree shape or topology (topological incongruence). Character congruence analyses are particularly useful and powerful because they take both the tree topology and the underlying support for the tree topology into account. Topological congruence techniques have the advantage of being able to compare trees derived from data that may not be strictly comparable or easy to include in the same analysis. 通常情况下,当一个系统发育数据集被划分成更小的分区时,每个分区都会产生具有不同拓扑结构的树。从这个结果可以得出两个结论:( 1 )部分或全部的树是错误的,而各分区具有相同的历史;( 2 )或者树都是正确的,不同的分区经历了不同的进化历史。区分这些选项需要进行统计检验,以确定拓扑结构的差异是否可能只是偶然观察到的。 在系统生物学领域有许多不同的差异检验,都使用可比较的测量和想法。然而,检验在假设、实现和解释方面有时有微妙的,有时是显著的不同。在杂乱无章的文献和围绕这些检验的争论中,很难辨别出这些细节。 不一致性检验可以大致分为基础特征不一致的检验和只考虑树形状或拓扑(拓扑不一致)的检验。特征一致性分析特别有用和强大,因为它们同时考虑了树拓扑和对树拓扑的底层支持(数据来源)。拓扑同余技术的优点是能够比较来自数据的树,这些数据可能不具有严格的可比性,也不容易包含在同一分析中。 Planet P J . Tree disagreement: Measuring and testing incongruence in phylogenies . Journal of Biomedical Informatics, 2006, 39(1):86-102.
Creating Majorana fermions in topological insulators http://nsr.oxfordjournals.org/content/1/1/36.extract HOT ! 由上海交通大学物理与天文系刘灿华、贾金锋教授共同撰写的观点文章Creating Majorana fermions in topological insulators(译名:在拓扑绝缘体中创建Majorana费米子)已在 National Science Review (NSR ) 2014 年3 月份的创刊号 上发表。 这篇观点文章先简要地介绍了Majorana费米子在凝聚态物理学领域当中的概念、特性和存在形式,然后介绍了如何在拓扑绝缘体的表面创建Majorana费米子并通过实验加以探测的研究现状及科研前景。 Majorana 费米子这一概念最初是于1937年,在粒子物理学领域当中为了描述与自身的反粒子同体的一类费米子而提出来的。但至今,尚未有严格的实验证据表明某个实体粒子是Majorana费米子。在凝聚态物理学领域当中,存在着许多准粒子,当中的某些可视为粒子与反粒子。在特定的条件下,这些互为正反的准粒子由于具有完全一样的自旋以及零电荷和零能量,从而符合描述Majorana费米子的方程。而且,成为Majorana费米子的准粒子由于具有非阿贝尔的统计特性,而且对去相干具有极强的抵抗性,所以很适合用于量子计算,从而在近年引起了广泛的研究热潮。 比如,超导体的超导能隙当中的对称能级处的电子激发态与空穴激发态,就可视为一对正反准粒子,他们在费米面处的同权叠加态,就可构建出Majorana费米子。不过,一般的超导体由于其库珀对的轨道角动量多为零(s波超导体),所以是无法在费米能级处形成零能量的电子与空穴的同权叠加态。只有在一些特殊的超导体当中,库珀对的轨道角动量为1的时候(p波超导体),才能构建出符合Majorana费米子特性的准粒子。但p波超导体材料并不多见。有理论科学家发现,将s波超导体放到拓扑绝缘体表面,由于该表面态电子的自旋与角动量是藕合在一起的,所以通过近邻效应引入到该表面的库珀对就可形成类似于p波的配对形态,从而在特定条件下(比如外加磁场下的vortex中心部),构建出Majorana费米子。 近年来这方面的实验进展非常迅速。比如中科院物理所吕力研究员带领的团队,就率先在拓扑绝缘体的表面探测到了因近邻效应而引入的超导特性 。上海交大贾金锋教授的团队与合作者一起,则率先在s波超导体的表面生长出了高质量的拓扑绝缘体(TI/sSC异质结),并且从实验上证实了拓扑表面态与超导态的共存 。最近,该研究团队还成功的在该TI/sSC异质结表面观测到了磁场下的vortex,并在vortex中心部检测到了包括Majorana费米子在内的bound states和Majorana费米子存在的迹象 。而清华大学薛其坤教授的研究团队与周树云教授合作,将高质量的拓扑绝缘体薄膜生长到了高温超导体的表面,并观测到了拓扑表面态的超导能隙,从而为探测Majorana费米子创造了更为有利的实验平台 。Majorana费米子探测到工作方兴未艾,相信在不远的将来人们能够捕捉到这种神秘的粒子。 图1. 拓扑表面态与超导的共存,以及拓扑绝缘体薄膜表面的vortex和零偏压峰。(a) 6层Bi2Se3 薄膜上测得到拓扑表面态。(b) 扫描隧道谱测得的费米面附近的态密度表明存在一个超导能隙。(c) 生长在超导体NbSe2 衬底上的拓扑绝缘体Bi2Se3薄膜的STM形貌图。 (d) 在vortex中心测得的芯能级峰,Majorana费米子就隐藏在这些芯能级中。(e) 单个的vortex。 引文列表: F. M. Qu et al., Sci. Rep., 2, 339 (2012). M. X. Wang et al., Science 336, 52 (2012). J. P. Xu, et al., arXiv: 1312.3713 (2013). J. P. Xu, et al., arXiv: 1312.7110 (2013). E. Wang et al., Nat. Phys. 9, 621 (2013). 原文详见 National Science Review 网站: SPECIAL TOPIC: TOPOLOGICAL INSULATORS- PERSPECTIVES Canhua Liu and Jin-Feng Jia Creating Majorana fermions in topological insulators d oi :10.1093/nsr/nwt036 (2 pages) o Abstract o Full Text (HTML) o Full Text (PDF)
31.拓扑绝缘体(续)(系列完结篇) 上节中介绍的石墨烯,由于它独特的物理性质而引起了人们的兴趣。它的无质量的相对论性准粒子,被观察到的整数及分数量子霍尔效应,为基础物理研究的许多方面,提供了理论模型和实验依据。它优异的电子输运性质,又使其在自旋电子学等工程领域可能得到广泛的实际应用。 图 31.1列出了石墨烯及量子霍尔态等几种物态在费米能级附近的能带图。 从图 31.1中的(a)和(b),我们可以看到双层和单层碳原子结构能带形状的不同。前者是抛物线型接触,而后者是线性的。(必须提醒注意的是,我们所说的这两种石墨烯能带图都是指在二维空间中能无限延伸的理想晶体之能带图。) 那么,量子霍尔态的能带形状又如何呢? 图 31.1:两种石墨烯及量子霍尔态等能带图之比较 图 31.1c是量子霍尔态的能带示意图。它的导带及价带在费米能级附近的形状,接近抛物线,类似于普通绝缘体。但是,我们在上一节中也说过,量子霍尔态体内虽然是绝缘体,但它们由于边缘态的存在而导电。在图中,量子霍尔态的边缘态是一条连接导带和价带的直线。因此,量子霍尔态在低能态附近的行为,和石墨烯相仿,能量和动量的关系也是线性的,也存在无质量的相对论性准粒子。 因为量子霍尔态的实现需要强大的外磁场,由此人们将兴趣转向不需要磁场的量子自旋霍尔效应,并且在实验室里已经多次观察到了此种现象。对量子自旋霍尔态而言,不同的自旋有不同的边界态,因此,在图 31.1d所示的自旋霍尔态能带图中,有两条直线连接导带和价带,它们分别对应于自旋上和自旋下的边缘电流。这种情形下的能带图,看起来与理想石墨烯的能带图更为类似了。 普通的绝缘体,也可能产生边缘态而形成边缘导电,但却和前面两种情形下的边缘态有本质的区别。图 31.1e画出了普通绝缘体的能带。图中的边缘态曲线与费米能级相交,意味着在此绝缘体中可以存在边缘电流。 再仔细对照一下 c、d、e三个图边界态的异同点便不难发现,即使从这三个简单图中,也可以看出一点刚才所谓的“本质区别”来:普通绝缘体的那个边缘态的导电性是不稳定的:系统的缓慢连续变化可以使导电性增加或消失。比如说,在缓慢变化下,这个边缘态可以降低到与价带相交而增加导电性,但也可能渐渐升高而脱离费米能级线,最后被归类到导带中,而使得边缘失去导电性。但是,图c和图d所示两种量子效应下的边缘态,却是一条直线,直通通的从上到下,将导带和价带绑到一起,这个连接方式不会因为系统的缓慢连续变化而改变。或者说,图c或d,与e的不同之处,可以用一句话概括:边缘态的拓扑结构不同。图e所示边缘态的拓扑结构是平庸的,而图c或d的边缘态则非平庸,其导电性能受其拓扑性质所保护,这一类的量子物态,便被称为“拓扑绝缘体”,以区别于平庸的普通绝缘体。真空属于普通绝缘体。 前面的叙述中,为什么总是要加上一句“系统缓慢连续变化”呢?这句话的意思,在数学上是为了保证系统的拓扑性质不变,在物理上则是保证系统不发生量子相变。比如说:一坨类似球形的面团,如果被你缓慢连续地揉来揉去,仍然是类球形的一坨面。但如果你把它从中间挖了一个洞,那就不是保持拓扑性质不变的“缓慢连续”变化,而是“相变”了。 刚才是用简单的图像来说明拓扑绝缘体与普通绝缘体的基本不同点。现在让我们在这条路上走得更远一些。其实,图 c、d、e中边界态的拓扑性质只是表面现象,并不足以解释拓扑绝缘体的本质,边界态表现不同的更深层原因,是由于体材料能带拓扑的不同。 当两个拓扑特征不同的绝缘体放在一起,就会产生导电的边界态。界面变成金属性,才能实现两种拓扑特征的连续变化。 既然是用拓扑性质来区分量子态,那么,便需要找一个拓扑不变量来表征不同的态。这个拓扑不变量通常对应于参数空间中不可积的贝里相位,贝里积分是在体材料的动量空间中进行,与边缘态无关。由此再次证明,是体材料的能带拓扑结构决定了边缘态的拓扑结构,从而才又决定了拓扑绝缘体的那种“被拓扑保护、不受杂质和缺陷干扰”的边缘导电性。 对整数量子霍尔态而言,这个拓扑不变量就是在动量空间计算出来的“第一陈数”,它同时也等于与经典朗道能级有关的填充因子 n。朗道能级是由外磁场而产生的,所以,正如我们从描述整数量子霍尔效应的电阻平台示意图所见,实验中观察到的n与外磁场强度有关。但是,在量子自旋霍尔效应中,外磁场强度等于0。也就是说,量子自旋霍尔效应的n值只能为0,换言之,不能再用第一陈数来表征量子自旋霍尔态了。 那么,有什么其它的不变量,能用来表征量子自旋霍尔态呢? 量子自旋霍尔态的特点是不存在外加磁场,因而,在一定的条件下可以具有时间反演对称性。“时间反演”,什么意思?顾名思义嘛,那就是将时间的流逝方向反过来。当然,真实的世界中时间是不会倒流的,但是电影技术为我们提供了一个用想象来检验时间反演特性的最佳场所。如果将一个个的电影画面反过来放,就能模拟时间反演的过程。从倒放的电影中我们会发现:有些东西(物理量)是正放反放不变的,而有些是改变的。比如说,我们考虑电磁场中的运动电子所涉及的几个物理量:位置将不受时光倒流的影响,但速度要反向;电子的电荷是时间反演不变的,但因为速度反过来了,所以电流要反向;电场强度 E是时间反演不变的,而磁场B要反向。磁场反向的原因是因为磁场是由电流产生的,时间倒过来时,电流反向了,因而磁场也反向了。 由上可知,磁场不是时间反演不变的。量子自旋霍尔态没有磁场,因而便有可能保持系统的时间反演对称性。或许可以利用这点来找出表征量子自旋霍尔态的守恒量? 相关于时间反演不变性, Kane and Mele提出用Z2不变量来区别拓扑绝缘体和普通绝缘体 【 1】 (Z2是指有两个元素的循环群)。在他们的模型中,将自旋霍尔态看成两个(自旋上和自旋下)边缘电流方向相反的整数霍尔态的合成,见图31.2。 图 31.2:自旋下的IQHE加自旋上的IQHE等于QSH 两个整数量子霍尔态相加,外磁场互相抵消了,剩下两个方向相反的自旋流,表现为量子自旋霍尔态。这两个 IQHE,可以分别用自旋陈数n 上 (自旋上)和n 下 (自旋下)来表征。Kane等人证明,时间反演对称性要求:n 上 + n 下 = 0,所以,总陈数为零。但是, n S = (n 上 - n 下 )/2 不会等于 0。并且,他们还证明,可以用n S 的奇偶性来描述合成量子态的非平庸性:当n S 为奇数时,系统是非平庸的拓扑绝缘体;当n S 为偶数时,系统是平庸的普通绝缘体。 因此,类似于 IQHE中的陈数n,定义一个Z2拓扑不变量 n = n S mod(2)。不变量 n 便可以用来表征二维拓扑绝缘体。这个概念还可以扩展到三维的拓扑绝缘体,即用 4个Z2不变量来表征三维拓扑绝缘体 【 2】 。 与文小刚提出的属于长程整体纠缠的拓扑序概念不同 【 3】 ,拓扑绝缘体和量子自旋霍尔态是属于更局域的短程量子纠缠态。它们也可以看作是被某种对称性所保护的拓扑序的例子:拓扑绝缘体被电荷守恒和时间反演所保护;而量子自旋霍尔态则被电荷守恒和 z方向自旋守恒所保护 【 4】 。 前面讨论的量子自旋霍尔态,是假设材料中两种自旋的密度在费米能级附近是相等的。反之,如果某一个方向的自旋被抑制,比如说,在某种材料中掺入某种铁磁性杂质,这样,就将破坏时间反演对称性,并有可能得到另外一种也不需要强大外加磁场的量子物态:量子反常霍尔效应。 刚才几句话说起来容易,实现起来却是非常困难。中国科学院院士薛其坤带领的团队, 2013年在世界上首次发现了量子反常霍尔效应。对此我们不再作更多的介绍,请见参考资料 【 5】 。 拓扑绝缘体及各种量子物态拓扑分类的理论中,仍有许多尚待解决的问题。这其中涉及的概念,既关联到基础物理思想,也包括不同领域的数学理论。总之,大门已经敞开,理论还需完善,精度日益提高的实验技术,也将供给我们越来越精确的数据。随着越来越多的不同量子态被研究、被发现,物理学必将继续造福文明社会。 参考资料: 【 1】C. L. Kane and E. J. Mele,Z2Topological Order and the Quantum Spin Hall Effect,Phys.Rev. Lett. 95, 146802 (2005)。 【 2】 余睿、方忠、戴希, Z2拓扑不变量与拓扑绝缘体,《物理》2011年第7期 462-468页 【 3】Topological order from wiki: http://en.wikipedia.org/wiki/Topological_order 【 4】Quantum spin Hall state from wiki: http://en.wikipedia.org/wiki/Quantum_spin_Hall_effect 【 5】Chang C Z, Zhang J, Feng X, etal. Experimental Observation of the Quantum Anomalous Hall Effect in a MagneticTopological Insulator . Science, 2013. (全文完) 上一篇: 系列科普目录
前面几篇简要涉及了代数拓扑、微分拓扑、低维拓扑,如果大伙儿还不知道这些是啥,请复习拙著 :) 。其实,拓扑的概念和方法已经渗透到整个数学,而不仅限于拓扑学研究本身。 很多文献会提到,拓扑学起源于柯尼斯堡七桥问题,以及与此相关的一笔画问题。欧拉解答了这个问题。同样是欧拉,给出了第一个拓扑不变量,多面体表面的欧拉数(点数-边数+面数)。可以认为欧拉是第一个研究代数拓扑学的人(虽然莱布尼兹曾经臆想过)。 传说中高斯当年是德国国土局的领导,负责丈量土地。他由此观测到山脊附近的地面弯曲性质可以用内在的测量方法得到(所谓曲面的内蕴几何),他还定义了度量曲率大小的量(高斯曲率),并且将这个局部定义的、可以用微积分计算的量同整体定义的欧拉数联系起来(高斯-博内定理),所以他可以被追认为研究微分拓扑的第一人。当然,高斯的主业其实是政府工作、开创现代数论、以及电磁学。业余时间研究一下各种误差的分布啊,欧几里德的几何原本有什么错误啊之类的小问题作为消遣。 第二尊菩萨,黎曼同学,除了开创流形的几何学、发展了傅立叶分析和积分理论、研究了一下素数分布提出世界第一难题黎曼猜想之外,他不到40岁的短暂一生的其余时间基本上都在思考复变函数的问题。为什么有的复变函数是多值的?比如平方根和对数。能否把它们以某种方式变成单值函数?人们的思维定势是,既然函数是多值的,就像一条横着的抛物线,一个 x 对应到两个 y, 要变成单值很容易啊,砍掉抛物线的一半就行了。黎曼不这么看。他觉得,如果把函数图像本身作为定义域,每个点当然只对应到一个 y. 这样函数就变成单值的了,而且没有丢掉任何信息。如果是复变函数,其图像就是一个二维曲面,这就是黎曼曲面。黎曼曲面上有很多复杂的现象,这些现象催生了诸如连通性、单连通性、复叠空间这些拓扑概念,以及奇点、除子、函数域等等一些代数几何的概念。抛开代数几何不说,黎曼也许可以被屈尊为研究低维拓扑的第一人。 牛顿莱布尼兹发明了微积分之后,大家对无穷小无穷大这两个概念很不放心。为了让我们用得更安心,柯西和威尔斯特拉斯等人后来把无穷小解释得非常透彻,基本上就是说两个东西越来越近。离得近这个概念从而成为分析中最核心的概念。它正是所有拓扑学分支的共同基础点集拓扑的起源。拓扑学在现实生活中的应用多数跟点集拓扑有关,即,通过分析收敛性体现在应用数学中。所以,最牛的牛顿被屈尊为研究点集拓扑第一人(有好事者不以为然,一定要扯到古希腊的阿基米德,这就见仁见智了。) 总而言之,拓扑学有着高贵的血统。当然,好汉不提当年勇,拓扑学的现在和将来如何?20世纪中叶,代数拓扑学朝着高度抽象的方向发展,两位名不见经传的数学工作者(艾伦伯格-麦克雷恩)突发奇想,从中总结出一套抽象语言。为了体现这套语言之形而上,他们重载了先贤亚里士多德的概念范畴。继解析几何与微积分以来人类数学又一次在概念上经历了大变革。范畴论诞生了。20世纪数学的上帝格罗登迪克,用7000页的数学圣经将拓扑学和范畴论全面而深遂地渗透到代数几何和数论研究中,改变了整个数学的风貌。与此呼应,世纪之交的物理学也在经历变革,量子场论和弦论呈现出绚丽多姿的数学结构。20世纪物理学的耶稣(上帝被爱因斯坦附身了)爱德华.威顿,身负绝世的拓扑神功,一经施展,整个理论物理学为之色变。拓扑学三位一体,必须要像生物学一样将21世纪纳入自己的势力范围。所以,在这个系列的结尾,让我骄傲地替拓扑学宣称:21世纪是拓扑的世纪!! I am just kidding.
前面讲到有理面上的奇异性,提到:正因为这个薄层非常薄,给我们处理问题反而带来极大的方便:第一,我们可以把轮胎形状的有理面沿着赤道平面切一刀,再沿任一环向角切一刀,展开成一个平面;第二,我们可以利用边界层理论来处理这一问题。这两点是研究像磁重联这样在大尺度下存在奇异性的数学物理问题的关键。 其中第一点将复杂的三维问题变成了一维问题! 把有理面展开成一个平面,则在两个延伸方向(大环和小环方向)上都有周期条件。做相应的 Fourier 展开之后,在每个特定的有理面上只有一个 Fourier 分量。且对应的微分符号变成代数符号。这样三维的微分方程在两个方向上代数化了,只剩下小环径向的变化。在物理上,这是由于有理面上的奇异性导致垂直有理面的特征尺度远远小于其它两个周期条件方向的尺度。 而第二点则为解决这样的问题提供了常用的方法。 我们知道,对连续介质中没有外部驱动时物理量 F 随时间演化的典型数学物理方程基本上都是带有耗散项的抛物型方程,比如 dF/dt = DF/Dt + vDF/Dx = l D 2 F/Dx 2 这里 D 表示偏微分。方程的右边是所谓耗散项, l 是耗散系数。如果 F 是磁通(磁矢势的某一分量),则这个方程就是欧姆定律:左边是电场(包括 v x B 部分), l 是电阻, D 2 F/Dx 2 是电流。等等。如果耗散很弱,耗散系数 l 非常非常小,则我们可以做理想情况下的近似,令 l = 0 。物理上就是,如果系统的特征尺度是 L ,那么,对应耗散( dissipation )的特征时间显然就是 T D =L 2 / l 。 l 趋于 0 对应于物理量 F 被 dissipated 的时间趋于无穷大。所以近似有 dF/dt = 0 ,或者说, F 基本保持不变。 这样的近似在绝大多数区间是适用的。但是如果区间里存在着奇异性,问题就来了:这个奇异性存在于一个很小的区间,引进一个非常小的特征尺度 D 。则其对应的特征时间尺度 t D = D 2 / l 成为一个可以和系统特征运动时间 T 0 相比拟的有限值!相应的,在这个奇异面上,方程右边的耗散项就不能忽略。这时,理想近似下得到的所谓奇异面因为这个耗散效应的存在变成了奇异层。显然这个奇异层的厚度 D ~ ( l T 0 ) 1/2 (我们又看到了 Sweet-Parker 模型的 1/2 方关系)。 这样一类在 奇异层外部可以用理想近似求解,但是在其内部必须考虑耗散项的数学物理问题,我们称之为边界层问题 ( Boundary Layer Problems )。求解的方法称为边界层方法。
Sweet-Parker 模型虽然很成功,但是存在一个致命的问题:磁重联的速率太慢。事实上,太阳大气等离子体的电阻率大约在 1/10 10-12 ,由此得到的重联率 v x ~ h 1/2 ~ 10 -5 -10 -6 Alfvn 速度,远远不能解释像日耀斑这样的快过程。 人们注意到 Sweet-Parker 模型重联率所以相对比较慢,原因是其重联区的拓扑结构近似是一维的,即我们前面说的:等离子体携带磁力线进入扩散区的方向(入流方向,通常选作为 x- 方向)上的特征尺度远远小于磁力线重联以后携带等离子体离开扩散区的方向(出流方向,通常选作为 y- 方向)的特征尺度。这样,由于出口太小、进口太大,导致已经 merging 到扩散区附近的磁力线的排队等候,物理学家用的词汇叫 magnetic flux piled up 。(这样的过程会在重联区形成很薄的强电流片,其物理效应我们以后再谈。)因此,有人( Petschek , 1964 )提出一种快磁重联模型:认为重联区的拓扑是呈具有 X 分形线的二维结构,这样入流区(在 y- 方向上的)长度与出流区(在 x- 方向上的)宽度大约在同一个数量级。而出流的喇叭口形状会形成一个如钱江潮的慢激波(道理相似,但相对运动方向相反)。根据这个 slow shock (慢激波)上下游的连接条件,可以得到磁重联的速率 v x ~ - ln h h 1/2 h 这个几乎与电阻无关的重联率基本上可以很好解释日耀斑这样的快过程。 但是人们后来发现,在 Petschek 模型的物理讨论所依赖的电阻磁流体框架下,无法得到 X 型的磁场几何结构,除非电阻很大 h 10 -3 ( W. Park, et al, 1984: Phys. Fluids 27 , 137; D. Biskamp, 1986: Phys. Fluids 29 , 1520; Z. W. Ma et al, 1995: Phys. Plasmas 2 , 8 ) 。而对于这么大的电阻, - ln h 与 h 1/2 的重联率几乎没有可以明显区分的差别!而在对应实际物理世界的电阻很小、 Sweet-Parker 和 Petschek 这两种模型的结果有可以明显区分的差别的情况下, 1980 年代以后发展的高精度的数值模拟结果告诉我们:即使初始条件取 Petschek 模型的磁场分布,我们也总是得到 Sweet-Parker 的电流片几何位形和 h 1/2 的重联率! 后来人们才意识到:尽管 Petschek 模型的磁场拓扑结构是出于增大重联率的正确考虑,但是使用的电阻磁流体的物理模型是错误的!正确的快磁重联模型依赖于 1990 年代无碰撞磁重联理论的发展。