铺砌、替换与准晶

规则只有两行字，图案却不是随便什么图案。从一个 A 出发，每轮把每个 A 换成 AB、每个 B 换成 A，第 12 轮得到一条长 377 的串。A 和 B 的个数比是 1.618056。这不是巧合，它由替换矩阵 M = [[1,1],[1,0]] 的最大特征值决定，而这个矩阵正是斐波那契数列的那个矩阵，特征值就是 φ = 1.618034。第三行读数在核对另一件事：串长每轮乘以 1.618026，也等于 λ。"增长率"和"两种符号的比例"是同一个特征值的两种表现。换成图厄–莫尔斯规则：比例变成 1.000000，λ = 2。这条串有个著名的性质——它不含任何"立方"（连着三次的重复段），被用在国际象棋的和棋规则争论、公平分配（Thue–Morse 顺序比"你一个我一个"更公平）甚至无线电跳频序列里。换成倍周期规则：比例 2.000733，λ = 2。这条串描述的是通往混沌的倍周期道路上的轨道次序。现在换成康托三分：A→ABA，B→BBB。B 是空隙，一段被挖空中间三分之一，再对剩下两段重复。第四行读数停在 0.630930，一级都不动。这正是 log2/log3——康托集的维数。换成康托五分（五段留三段）：0.682606 = log3/log5。这两个数是精确的，不是拟合出来的：盒子数是 2^L（或 3^L），盒子尺度是 3^(−L)（或 5^(−L)），相除的对数当然就是常数。对照之下，斐波那契和图厄–莫尔斯的这个读数慢慢爬向 1（第 12 轮 0.9189，第 16 轮 0.9388）——因为它们的 A 占了正比例的位置，根本不是分形集，它们的有趣之处在"顺序"而不在"维数"。（这里有一个很值得记住的分类：同样是替换系统，有的生成分形集（康托类），有的生成非周期但满密度的序列（斐波那契类）。区分它们的是替换矩阵——看它有没有一个符号"只生成自己"，那就是空隙。）

这张图只用了两种形状：一种胖菱形（锐角 72°），一种瘦菱形（锐角 36°）。它铺满了整个平面，而且永远不周期——你没法把它平移一段距离后和自己重合。生成方法叫"充气"（deflation）：把每一块按黄金比例切开，换成更小的一批块，每细分一级，边长乘以 1/φ = 0.618034（第四行读数）。看第二行：胖瘦两种块的数量比是 1.617647（第 5 级），第 7 级是 1.618026。它在收敛到 φ。第三行更漂亮：面积比收敛到 φ² = 2.618034。这两个数决定了一件很深的事：胖瘦菱形的数量比是无理数，而一个周期性图案里任何两种基元的数量比必然是有理数。所以这张图不可能是周期的——不是"我们还没找到周期"，是数学上不存在。把放大倍数拖到 4 倍，在图中间慢慢看。你会发现一件很奇怪的事：任何一小块图案，都会在别处无限多次重复出现（这叫"重复性"，repetitivity），但整体永远不重复。局部处处似曾相识，全局永不闭合——这正是准晶最难被接受的那一点。还有一个很实际的追问：这样的铺法能不能"排"出来？彭罗斯给了两块瓷砖和一套边界匹配规则，按规则拼就只能拼出非周期图案。但局部规则不足以保证能一直拼下去——贪心地拼很容易走进死胡同。这说明"非周期铺砌存在"和"能用局部规则一步步生成"是两回事，而这正是非周期铺砌与不可判定问题（王浩铺砌问题）相连的地方。（2023 年数学界找到了"帽子"和"幽灵"——单独一块就能非周期铺满平面的形状，这个悬了六十年的问题才算真正闭合。）

"自相似铺砌"是一个很朴素的想法：一块形状，能被切成几块和它自己形状相同、只是变小的块。任何三角形都可以——连四条中位线一切就是 4 块，边长减半；正方形可以切成 4 块，也可以切成 9 块。把"留下几块"拖到最大（全留），维数是 2，因为它就是原来那块形状，一点没少，实测盒维数 1.9401 / 2.0000 / 2.0000。现在把留下的块数减一。图案立刻变成分形。三角形留 3 块，理论维数 log3/log2 = 1.5850，这就是谢尔宾斯基三角形；正方形切 9 留 8（丢中间那块），理论 log8/log3 = 1.8928，这是谢尔宾斯基地毯。这一页想让人亲手体会的，是"丢一块"的后果有多大。全留时面积是 1；每级丢掉九分之一，第 4 级面积只剩 (8/9)⁴ = 0.624，而继续下去，面积趋于 0——图看上去还在，面积已经没了。一个维数 1.8928 的东西，既不是曲线（维数 1）也不是面（维数 2），它的"长度"是无穷，"面积"是零。这不是文字游戏：多孔材料、肺泡、催化剂载体正是靠这个性质工作——在有限的体积里塞进近乎无限的表面积。最后一件事，这一台最值得动手的地方：把「细分几级」从小拖到大，盯住实测维数。它不是稳定在理论值上，而是从上方一路降下来，中途穿过理论值，然后继续往下掉。谢尔宾斯基三角形（理论 1.5850）实测：第 3 级 1.7192、第 4 级 1.6330、第 5 级 1.5607、第 6 级 1.5014、第 7 级 1.4421。地毯（理论 1.8928）：第 2 级 1.9442、第 3 级 1.9055、第 4 级 1.8747、第 5 级 1.8501。两头错的原因不一样，这才是重点。级数太少时结构还没长出来，图形基本还是实心的，盒计数量到的接近 2，所以偏高；级数太多时每个小块已经比一个像素还小，整张图退化成一堆孤立的点，量到的接近 0，所以偏低。只有中间那一段——小块比像素大、而层次又够多——测出来的数才可信。（留 2 块那一档更极端：理论 1.0000，第 4 级实测 1.6443，第 6 级 1.2286，第 7 级才落到 1.0000。）这是所有盒计数测量的通病：报一个「实测分形维数」而不报它是在哪一段尺度上测的，等于没报。真实的研究会把双对数图上那段直线的范围一并画出来，让人自己判断它到底直不直。

衍射图是判定一个固体"有没有秩序"的黄金标准。上面是原子链，下面是它的衍射强度。先选"等间距"：一根根又细又高的尖峰，最强峰强度正好 1.0000，位置 k·a/2π = 2.0000——整数，这就是布拉格条件。再选"完全随机"：尖峰全没了，只剩一片起伏的背景，最强处只有 0.0423。而且把原子数拖大，这个值还会继续掉——非晶态没有布拉格峰，这是定义性的。现在选"斐波那契链"。这条链的间距只有两种：长的 L 和短的 S，排列顺序由上一个实验台那条斐波那契串给定，它绝对不是周期的。按当时所有人的共识，它应该给出弥散的背景。而读数是：最强峰强度 0.7739，位置 k·a/2π = 1.6184。尖锐的峰，清清楚楚。而那个位置正是 φ = 1.618034。把 k 的范围拖到 20，最强峰跑到 2.6185 —— φ² = 2.618034。这就是准晶最核心的事实：它有长程序（所以有尖锐衍射峰），但没有周期（所以峰的位置是无理数）。晶体的峰在整数处，准晶的峰在 m + nφ 这样的位置上，需要两个整数指标才能标定一个峰，这就是"晶体学要用六个指标"的由来。最后选"等间距 + 随机抖动"做个对照：峰还在原来的整数位置，但强度掉到 0.2847。无序削弱峰的高度，却不移动峰的位置；非周期改变峰的位置，却不削弱它的尖锐。这两件事完全不同，而 1982 年之前的晶体学把它们混为一谈——看到五重对称的尖锐斑点，第一反应是"这一定是孪晶假象"，而不是"周期性假设需要修改"。谢赫特曼花了两年才发表那篇论文，国际晶体学联合会在 1992 年重写了"晶体"的定义：从"原子周期排列的固体"改成"能给出清晰衍射图的固体"。定义跟着证据走，这件事本身比准晶更值得记住。