1944年的纽约,战争的阴影笼罩着这座城市。大西洋彼岸,诺曼底登陆刚刚准备就绪;太平洋战场,美军正在跳岛推进。但在洛克菲勒医学研究所的一间安静办公室里,四十一岁的拉斯·昂萨格正沉浸在完全不同的战斗中——一场与16页数学公式的搏斗,对手是二维伊辛模型。
昂萨格身材瘦削,面容清癯,戴着一副黑框眼镜,穿着总是略显陈旧的西装。他的办公室堆满了草稿纸,上面密密麻麻写满了行列式、椭圆函数、和复变积分。同事们很少打扰他,知道这位挪威裔化学家正在从事某种"疯狂的事业"——试图精确求解一个统计力学模型,而大多数物理学家认为这不可能或不值得。
1944年2月,昂萨格完成了他的计算。他把手稿寄给《物理评论》,标题平淡无奇:《晶体统计.I.一个具有有序-无序转变的二维模型》。但内容却是数学地震:他证明了二维伊辛模型存在相变,计算了精确的临界温度,推导了完全不同于平均场理论的临界指数。
这篇论文将改变统计力学的进程,暴露平均场理论的致命缺陷,并为三十年后的重整化群革命铺平道路。但昂萨格本人几乎不解释他的结果。当同事们问他如何做到时,他只是说:"我计算了。"这种工匠式的谦逊,掩盖了一场概念的风暴。
战争时期的数学要理解昂萨格的工作,我们需要回到1940年代初的语境。量子力学已经成熟,但统计力学仍处于过渡状态。吉布斯的系综理论是标准语言,但具体计算通常依赖平均场近似——范德瓦尔斯方程、朗道理论、布拉格-威廉斯近似。
这些近似在远离临界点时有效,但在临界点附近系统性地失败。实验测量的临界指数(α≈0.1, β≈0.3, γ≈1.2)与平均场预言(α=0, β=0.5, γ=1)明显不同。但大多数物理学家认为,这种偏差是实验误差或样品不纯的结果。
精确解的想法是边缘的。量子力学的氢原子有精确解,但那是单粒子问题。多体系统的精确解被认为是不可能的——涉及太多变量,太复杂的相互作用。昂萨格的挑战是统计力学的"氢原子":找到一个可解的多体模型,严格计算其性质。
他选择二维伊辛模型不是偶然的。1925年伊辛的"失败"(一维无相变)暗示,维度是关键。昂萨格直觉感到,二维是临界点:足够复杂以展现相变,又足够简单以数学处理。这种维度敏感性后来成为临界现象的核心洞见。
昂萨格的方法是转移矩阵的杰作。伊辛在一维使用了转移矩阵,但二维需要更复杂的对象:行与行之间的转移算符,涉及4×4矩阵(对于正方形格点)。关键步骤是找到这个矩阵的本征值,特别是最大本征值,它决定配分函数的渐近行为。
这听起来技术性的,但涉及深刻的数学。昂萨格发现,转移矩阵的本征值问题可以映射到量子力学中的非对易算符,与角动量理论和李代数相关。他还使用了椭圆函数——十九世纪的数学工具,当时物理学家不熟悉——来处理对偶变换。
对偶性:隐藏的对称昂萨格解的核心美学是对偶性(duality)。这是一种数学对称,揭示高温相和低温相之间的深刻联系。
想象伊辛模型的两个版本:原始版本,自旋位于格点上,相互作用沿键(边);对偶版本,"自旋"位于** plaquette中心(方格的面上),相互作用沿对偶键(连接相邻 plaquette 的线)。在高温下,原始版本的关联长度小(短程无序);在低温下,对偶版本的关联长度小**。
昂萨格证明,两个版本的配分函数本质相同,只是温度参数互换。这意味着,如果存在一个临界温度Tc,使得高温相和低温相"相等",那么系统必须满足:
sinh(2J/kTc) = 1
或等价地:
kTc/J = 2/ln(1+√2) ≈ 2.269
其中J是相互作用强度,k是玻尔兹曼常数。
这个公式是精确的,解析的,不含任何近似。它与实验或数值估计完美吻合,证明了二维伊辛模型确实存在严格的相变。
对偶性的物理意义是深刻的:它表明高温无序和低温有序不是对立的,而是同一现象的两面。临界点是自对偶的——高温和低温描述在此融合。这种对称性的增强,是共形场论(1980年代发展)的先驱。
昂萨格没有强调对偶性的美学和哲学意义。他的论文是计算的洪流,定理接着定理,公式接着公式。但细心的读者可以发现,在第十六页的某个角落,隐藏着自然的秘密:临界点不是任意的,而是由数学对称性决定的。
临界指数:平均场理论的崩溃昂萨格解的最震撼结果,是临界指数的计算。这些指数描述物理量在临界点附近的行为:
比热容:在Tc附近,比热容不是平均场预言的有限跳跃,而是对数发散:
C ~ -ln|T-Tc|
这意味着,当接近临界点时,比热容无限增大,但无限缓慢(对数发散比任何幂次都慢)。指数α = 0(对数发散视为α=0的极限)。
自发磁化:昂萨格最初没有计算这个(他专注于无磁场的情况),但他的结果暗示了奇异性。1949年,昂萨格在剑桥的一个会议黑板上,写下了磁化强度的公式,然后擦掉了。听众之一,杨振宁,记住了关键步骤,后来与李政道合作,严格证明了:
M ~ (Tc-T)^β,其中β = 1/8
这与平均场的β = 1/2完全不同。
磁化率(响应外场的敏感度):昂萨格的结果暗示:
χ ~ |T-Tc|^-γ,其中γ = 7/4
与平均场的γ = 1不同。
临界等温线(T=Tc时的M-H关系):
H ~ M^δ,其中δ = 15
与平均场的δ = 3不同。
这些指数是无理数(1/8, 7/4, 15),普适(对所有二维伊辛系统相同),但非经典(不同于平均场)。它们后来成为普适类的精确参考:任何属于二维伊辛普适类的系统,实验或模拟的临界指数应该接近这些值。
杨振宁与磁化强度:黑板的传奇昂萨格1944年论文的著名遗漏,是自发磁化强度的精确公式。他计算了无磁场的情况,但没有给出有磁场的配分函数——这在数学上更困难。
1949年,昂萨格在国际统计力学会议(剑桥)上做了一个非正式报告。他在黑板上写下一系列公式,推导出磁化强度的行为,然后——据目击者说——在讲完后擦掉了黑板。听众之一,杨振宁,当时二十七岁,正在芝加哥大学做博士后,记住了关键步骤。
杨振宁后来描述这个时刻:"昂萨格的推导像闪电,太快以至于难以跟随。但我抓住了本质:使用自旋关联函数和行列式表示。"
回到芝加哥后,杨振宁与李政道合作,花了六个月严格证明昂萨格的结果。他们使用了旋量分析(spinor analysis)和组合数学,最终得到:
M = [1 - (sinh 2K / cosh² 2K)²]^(1/8)
其中K = J/kT。在T→Tc时,这给出β = 1/8。
这个结果发表于1952年,是统计力学的里程碑。它证明了,即使在最简单的模型中,临界行为也极端非平凡。指数1/8(不是1/2)意味着,磁化强度在临界点附近比平均场预言更平坦——涨落抑制了有序的建立。
杨振宁后来回忆:"那是我科学生涯中最兴奋的时刻之一。我们证明了,严格数学可以揭示自然的秘密,而近似方法可能系统性误导。"
昂萨格的风格:计算即理解昂萨格的科学风格是独特的:计算先于理解。他不追求直观的图像或简单的解释,而是追求数学的精确和结果的可靠。
这种风格与朗道形成鲜明对比。朗道是物理直觉的大师,用序参量和对称性破缺的语言描述相变。昂萨格是数学工匠,用行列式和椭圆函数计算性质。两种方法都是必要的,但昂萨格的结果更精确,因此更具破坏性。
昂萨格的写作是障碍。他的1944年论文极其难读,跳跃、压缩、使用非标准符号。许多物理学家尝试阅读,都迷失在第三页。据说,费曼曾试图理解昂萨格的解,放弃了,说:"如果我不能发明它,我就不能理解它。"
但耐心是有回报的。1950年代,考夫曼(Bruria Kaufman,昂萨格的学生)简化了推导,引入了更清晰的拓扑方法。李政道和杨振宁发展了相变的一般理论。费希尔(Michael Fisher)和卡达诺夫(Leo Kadanoff)将昂萨格的结果与标度理论联系。
昂萨格本人不追求解释。他相信,正确的计算会自动蕴含理解,即使这种理解延迟到来。他在1968年获得诺贝尔化学奖(对于不可逆热力学的工作),但从未因伊辛模型获奖——这被认为是诺贝尔奖的遗漏之一。
精确解的意义:数学作为探针昂萨格的解提出了深刻的方法论问题:精确解在物理学中的意义是什么?
在实验科学中,精确是测量的理想。在理论科学中,精确是数学的可能。但物理学的特殊之处是:数学精确可以揭示自然结构,即使这种结构反直觉。
昂萨格的解表明:
平均场理论是定性地错误(在低维),不是定量地不准确。
涨落和关联在临界点主导行为,不能作为"修正"处理。
维度是关键参数,改变临界行为的普适类。
对称性(如对偶性)决定临界温度,超越微观细节。
这些洞见不能从近似获得。只有精确解才能严格证明平均场理论的失败,才能确定临界指数的真实值,才能启发新的理论框架。
但精确解也有局限。昂萨格的解只适用于二维,且需要特殊格点(正方形)。三维伊辛模型——真实铁磁体——未被精确求解。这提示,精确性和普遍性之间存在张力:我们可以精确求解简化模型,或近似处理真实系统,但很少能同时做到两者。
这种张力是统计力学的核心挑战,也是活性算法的动机之一:能否设计自适应的方法,在保持精确性的同时处理复杂性?
从二维到三维:未竟的征程昂萨格的解激发了三维伊辛模型的求解尝试。1950-1960年代,许多数学家尝试推广昂萨格的方法:转移矩阵、椭圆函数、对偶性。但所有尝试都失败。
问题在于拓扑。二维伊辛模型的可解性依赖于平面格点的特殊性质:转移矩阵可以对易化,本征值问题可以分解。三维情况下,对易性丢失,问题变成本质复杂。
数值方法提供了替代路径。蒙特卡洛模拟(1953年发明)可以估计三维临界指数,精度随计算能力提高。1970年代的估计给出:
α ≈ 0.11(接近零)
β ≈ 0.326
γ ≈ 1.237
δ ≈ 4.79
这些指数无理数,普适,但与二维值不同。维度改变普适类。
场论方法提供了另一种视角。三维伊辛模型可以映射到欧几里得量子场论,临界现象对应于连续极限。这种映射是威尔逊重整化群的基础,但精确解仍然缺失。
至今,三维伊辛模型的严格解是数学物理的开放问题。它可能永远不可解,或需要全新的数学工具——如量子群、非对易几何、或活性算法的某种形式。
(三维伊辛模型已经被中国科学家张志东老师求解)
昂萨格晚年尝试三维问题,但未成功。他在1976年去世,留下未完成的梦想。但他的二维解仍然是灯塔:证明精确性是可能的,数学可以战胜复杂性,自然的秘密等待耐心者。
精确解的遗产:从统计力学到弦理论昂萨格1944年的论文,影响远超统计力学。它成为数学物理的范式,启发了多个领域的发展:
可积系统:昂萨格的方法(转移矩阵、对偶性、椭圆函数)发展为可积统计力学的一般理论。Baxter(Rodney Baxter)在1970-1980年代求解了多个二维模型(如八顶点模型、硬六边形模型),使用推广的对称性和量子群。
共形场论:1980年代,Belavin、Polyakov、Zamolodchikov(BPZ)发展了二维共形场论,描述临界点的连续极限。昂萨格的伊辛解成为最小模型(minimal model)的原型,其临界指数与共形对称性的表示理论相关。
弦理论:二维共形场论是弦世界面的动力学。伊辛模型的c=1/2中心荷, 对应于自由费米子,是弦理论的基本构建块。
量子计算:伊辛模型是量子计算的基准问题。绝热量子计算(adiabatic quantum computation)使用伊辛基态作为计算目标。量子退火(quantum annealing)使用伊辛哈密顿量求解优化问题。
机器学习:玻尔兹曼机(1985年)直接使用伊辛能量函数。深度信念网络的预训练涉及伊辛模型的变分推断。活性算法的自适应临界性,使用伊辛模型作为测试平台。
这种跨学科的影响,证明了简单模型的深远力量。昂萨格的二维伊辛解,像种子落入肥沃的土壤,生长出统计力学、量子场论、弦理论、量子计算、机器学习的交叉森林。
尾声:计算的尊严让我们回到1944年的纽约。战争仍在进行,但昂萨格在他的办公室里,完成了他的16页论文。他可能感到满足,但不太可能感到胜利的喜悦。对他来说,计算是日常的实践,不是戏剧性的成就。
但历史赋予了不同的意义。昂萨格的解是统计力学的转折点:从近似到精确,从平均场到涨落,从经典到量子。它暴露了旧理论的局限,指明了新理论的方向。
更重要的是,昂萨格的解确立了数学精确性的价值。在物理学中,近似是必要的,但精确是理想的。精确解是探针,测试我们的理解深度;是灯塔,照亮未知的海洋;是标准,校准近似的方法。
昂萨格本人是谦逊的工匠,不追求名声,不解释结果,满足于计算的正确性。但正是这种专注的谦逊,创造了持久的伟大。他的16页论文,像数学的丰碑,矗立在统计力学的中心,等待着每一代新物理学家的阅读。
在下一章,我们将看到,昂萨格的精确解如何与朗道的序参量理论相互作用。这种相互作用将暴露理论物理的深层张力:直觉与精确,近似与严格,物理图像与数学结构。这种张力将在威尔逊的重整化群中得到综合,但在此之前,它制造了困惑和争论——这正是科学进步的燃料。
但首先,让我们向那位挪威裔化学家致敬。他在战争的阴影中,计算着二维格点上的自旋,不知道他正在创造永恒。
本章注释与延伸阅读
昂萨格1944年的原始论文《晶体统计.I.一个具有有序-无序转变的二维模型》发表于《物理评论》(Physical Review)65, 117-149。这是数学物理的杰作,但极其难读。推荐辅助阅读:McCoy, B.M. and Wu, T.T. (1973). The Two-Dimensional Ising Model, Harvard University Press。
关于杨振宁和李政道的磁化强度计算,参见:Yang, C.N. (1952). "The Spontaneous Magnetization of a Two-Dimensional Ising Model," Physical Review 85, 808-816。
关于昂萨格的传记和科学风格,推荐:Longuet-Higgins, H.C. and Fisher, M.E. (1978). "Lars Onsager, 1903-1976," Biographical Memoirs of Fellows of the Royal Society 24, 443-471。
关于精确解在物理学中的意义,参见:Baxter, R.J. (1982). Exactly Solved Models in Statistical Mechanics, Academic Press(前言部分);以及Fisher, M.E. (1998). "Renormalization Group Theory: Its Basis and Formulation in Statistical Physics," Reviews of Modern Physics 70, 653-681。
关于对偶性和共形场论,推荐:Di Francesco, P., Mathieu, P., and Sénéchal, D. (1997). Conformal Field Theory, Springer(第一章回顾伊辛模型)。
转载本文请联系原作者获取授权,同时请注明本文来自王涛科学网博客。
链接地址:https://wap.sciencenet.cn/blog-41701-1525346.html?mobile=1
收藏