频谱分析是信号处理的看家本领。林薇知道，如果非高斯残差中真的隐藏着某种周期性信号，频谱分析会把它暴露出来——就像棱镜把白光分解成彩虹一样。
她选择了多窗谱分析（Multitaper Method）而不是传统的快速傅里叶变换。原因很简单：多窗谱分析在处理短数据和有色噪声时表现更好，而地磁残差恰好是这种情况——数据长度有限，背景噪声的功率谱不是平坦的。她用了5个Slepian窗，带宽参数NW=4，这是地磁数据分析中的标准配置。
计算跑了整整一个下午。当结果终于出现在屏幕上时，林薇的手指在键盘上停住了。
43个非高斯台站的残差频谱中，有一个共同的峰值——7.83Hz。不是所有43个台站都有这个峰值，但38个台站的频谱在7.83Hz处有统计显著的功率峰（超过背景噪声3倍标准差）。7.83Hz——这是舒曼共振的基频。
舒曼共振是地球电磁空腔的自然共振频率。1952年，德国物理学家Winfried Otto Schumann从理论上预测了这一现象：地球表面和电离层之间的空腔像一个巨大的共振腔，全球雷暴活动产生的电磁波在这个空腔中反射、叠加，形成驻波。基频约为7.83Hz，谐波频率约为14.3Hz、20.8Hz、27.3Hz和33.8Hz。舒曼共振是地球电磁环境中最基本的周期性信号之一，它像地球的心跳一样，从未停止。
"所以残差中包含舒曼共振成分，"林薇自言自语。"这有什么奇怪的？地磁数据当然会包含舒曼共振——它是全球性的电磁信号，地磁台站当然能记录到。"
但奇怪的不是舒曼共振本身的存在，而是它的出现方式。正常情况下，舒曼共振在地磁数据中的表现应该是：去除主磁场和日变化后，残差中的舒曼共振成分应该非常微弱——因为INTERMAGNET台站的采样频率通常是1分钟一次（1/60 Hz），远低于7.83Hz的奈奎斯特频率。7.83Hz的信号不应该出现在分钟级数据中——除非它通过某种非线性过程被"混叠"到了低频。
林薇重新检查了数据采样率。137个台站中，有12个台站同时记录了1秒级和1分钟级数据。她下载了这12个台站的1秒级数据，重新做了频谱分析。结果让她更加困惑：在1秒级数据中，7.83Hz的峰值确实存在，但它的幅度与1分钟级数据中的"混叠"峰值不成比例——1分钟级数据中的峰值太强了，强到无法用简单的混叠效应解释。
这意味着什么？这意味着分钟级数据中的7.83Hz峰值不仅仅是1秒级数据的混叠——它还包含了额外的贡献。某种东西在分钟级时间尺度上产生了与舒曼共振频率相关的信号。
她想到了一种可能性：调幅。如果舒曼共振的幅度被某种低频过程调制，那么在分钟级数据中，这种调幅信号就会表现为7.83Hz的"影子"——不是7.83Hz本身，而是7.83Hz的包络。
她做了希尔伯特变换，提取了7.83Hz信号的瞬时幅度。果然，瞬时幅度不是常数——它在缓慢地波动，波动的周期大约是几个小时。她对这个幅度包络做了频谱分析，发现了一个新的峰值：3.92Hz。
3.92Hz。恰好是7.83Hz的一半。
林薇盯着这个数字看了很久。7.83Hz的一半——次谐波。在非线性动力学中，次谐波的出现通常意味着系统发生了某种对称性破缺。一个以频率f驱动的系统，如果以f/2的频率响应，这意味着系统的动力学方程在时间平移T→T+π/f下不再不变——时间平移对称性被打破了。
但她立刻告诫自己不要过度解读。次谐波响应可以出现在很多非线性系统中，并不意味着时间晶体。一个简单的Duffing振子就可以产生次谐波——这并不稀奇。稀奇的是，这个次谐波出现在全球38个地磁台站上，而且它的幅度与舒曼共振功率之间存在定量关系。
她计算了3.92Hz信号幅度与舒曼共振功率的相关性。结果是：相关系数r=0.73，p<0.001。3.92Hz信号的幅度与舒曼共振功率强相关——舒曼共振越强，3.92Hz信号越强。这暗示了因果关系：舒曼共振驱动了3.92Hz信号。
但这仍然不能排除一个平凡解释：3.92Hz信号可能是舒曼共振的非线性混频产物。如果地磁台站的传感器或数据采集系统存在非线性响应，两个7.83Hz的信号可以混频产生3.92Hz的差频信号。她需要排除这种可能性。
她检查了台站的传感器规格。INTERMAGNET台站使用的磁通门磁力仪，其线性度在0.01%以内——这意味着非线性混频产生的3.92Hz信号应该比7.83Hz信号弱至少10000倍，远低于观测到的水平。传感器非线性不是3.92Hz信号的来源。
她又检查了数据采集系统的模数转换器。24位ADC，动态范围142dB——同样，非线性度远低于观测信号所需的水平。
排除了传感器和采集系统的非线性后，3.92Hz信号的来源只能是地球本身。某种地球物理过程正在以舒曼共振频率的一半振荡，而这个过程的幅度与舒曼共振功率相关。
林薇在笔记本上写下了推论链：残差非高斯→包含周期性成分→频谱分析→7.83Hz峰值（舒曼共振）→幅度调制→3.92Hz次谐波→与舒曼共振功率相关→舒曼共振驱动→排除非线性混频→信号来自地球本身。
她画了一个方框，在方框里写下了那个她一直在回避的问题："3.92Hz次谐波是时间晶体相的标志吗？"
她在方框旁边画了一个问号，然后合上了笔记本。她还没有准备好回答这个问题——甚至还没有准备好认真考虑它。但她知道，这个问题不会因为她回避就消失。
那天晚上，她第一次搜索了"discrete time crystal"的文献。Google Scholar返回了大约2300篇论文。她从被引最多的几篇开始读——Norman Yao 2016年的PRL、Monroe组2017年的Nature、Lukin组2017年的Nature、Google Sycamore 2021年的论文。她读到凌晨三点，读完了十二篇论文，做了七页笔记。
第二天早上，她顶着黑眼圈走进办公室，泡了一杯浓茶，然后给一个她从未谋面的人写了一封邮件。收件人是周远川，中国科学技术大学物理学院教授，凝聚态理论与量子信息方向。邮件的主题是："一个关于宏观时间晶体相的实验线索"。
邮件只有三行字："周教授，我在地磁数据中发现了一个3.92Hz的次谐波信号，它与7.83Hz舒曼共振的功率强相关。这个信号可能满足离散时间晶体相的次谐波响应判据。如果您有兴趣，我可以提供详细数据。"
她犹豫了整整十分钟才按下发送键。她知道这封邮件意味着什么——她在向一个陌生人暴露自己尚未完成的工作，冒着被嘲笑或被抢发的风险。但她也知道，如果这个信号真的是时间晶体相的标志，她一个人无法证明——她需要一个理论物理学家的帮助。
周远川的回复几乎是即时的："来合肥。带上你的数据。"周远川的回复几乎是即时的：来合肥。带上你的数据。
林薇看着屏幕上的这八个字，心跳加速了。她知道这封邮件和这个回复意味着什么——一个理论物理学家对一个实验线索产生了兴趣，愿意投入时间和精力来分析它。在物理学的历史上，很多重大发现都始于这样的一拍即合——实验学家发现了异常，理论学家提供了解释框架，两者结合，产生了突破。
但她也知道，这种一拍即合也可能导致灾难——如果理论学家过早地介入，可能会引导实验学家走向错误的方向。她需要保持独立思考的能力，不能因为周远川的权威就盲目接受他的判断。
她买了第二天去合肥的高铁票。G27次，北京南站到合肥南站，4小时23分钟。她带上了笔记本电脑、打印好的分析报告和一整盒方便面——她预感这会是一个漫长的讨论。
在高铁上，她重新整理了所有数据，把关键图表做成了PPT。她知道周远川是一个视觉思考者——他需要看到图表才能思考，纯文字的描述对他来说不够直观。她把频谱图、自相关函数图、空间分布图和相关性散点图都做成了大尺寸的高清图，确保每一个细节都清晰可见。
列车穿过华北平原，窗外的景色从城市过渡到农田，又从农田过渡到丘陵。她看着窗外飞速后退的风景，心里想着7.83Hz——这个频率在地球电磁空腔中已经回荡了数十亿年，从地球拥有大气层和电离层的那一天起，舒曼共振就开始了。它是地球最古老的声音——比任何生命都古老，比任何文明都持久。而现在，这个古老的声音似乎在诉说着一个更深层的故事——一个关于时间晶体的故事。合肥的中国科学技术大学校园比林薇想象的更美丽——绿树成荫，湖水清澈，教学楼散落在山丘之间，有一种宁静而专注的学术氛围。周远川的办公室在物理学院的四楼，门口贴着一张手写的告示：请敲门后等待三秒再进入——我可能在思考。
林薇敲了门，等了三秒，然后推门进去。周远川正站在黑板前，手里拿着一支粉笔，黑板上写满了公式。他四十出头，瘦高个子，戴着一副金丝边眼镜，头发有点乱——像是刚用手抓过。
林薇？请坐。他指了指办公桌前的椅子，然后继续在黑板上写了一会儿，才放下粉笔转过身来。
抱歉，我需要把那个推导写完——不然会忘。他笑了笑，一种带着歉意的笑。你带了数据？
林薇打开笔记本电脑，把分析报告和图表一一展示。周远川看得很仔细，不时停下来问问题——不是泛泛的问题，而是非常具体的技术问题：Shapiro-Wilk检验的功效在样本量N=1000时是多少？Moran's I检验的零假设是什么？3.92Hz信号的带宽是多少？
林薇一一回答。她发现周远川的提问方式非常独特——他不是在质疑她的结论，而是在检验她的推理。每一个问题都指向推理链中的一个薄弱环节，如果她能回答所有问题，就意味着推理链是完整的。列车穿过华北平原，窗外的景色从城市过渡到农田，又从农田过渡到丘陵。她看着窗外飞速后退的风景，心里想着7.83Hz——这个频率在地球电磁空腔中已经回荡了数十亿年，从地球拥有大气层和电离层的那一天起，舒曼共振就开始了。它是地球最古老的声音——比任何生命都古老，比任何文明都持久。
