而这个黑黢黢的小东西……用了不到三秒。
第304章 顺手的事
第二天清晨。
李东合上电脑时，窗外的天色已经微微亮了。
草稿纸散落一地，他揉了揉有些发胀打的太阳穴，然后看向了屏幕右下角。
小黑头上的进度条已经来到了【0.18】。
【主人，一百九十六万条，全吃完啦。嗝……】
李东失笑，敲下键盘。
【真信号多少？】
【八万六千三百九十二条，其余全是幻峰。】
李东点点头。
【小黑，挂载昨天的脚本，叠第二层筛子。】
【在十几个TeV能段，找形状符合AD 774/5外推预测的事件，χ2/d.o.f.卡在1.1以内。】
小黑没有废话，又开始蹦蹦跳跳起来。
到二十秒，一行字弹了出来。
【主人，挑好啦，一共十一笔。】
李东嘴角抽了一下。
十一笔。
跟昨晚他用那三十七笔“待复核”事件包外推出来的数字，分毫不差！
如果说昨晚是还是弱证据的话，那现在，用悟空号八年全量数据跑出来的结果，就是强证据了。
深吸一口气，李东知道，这第一关算是过了。
接下来才是真正的难搞的事。
他要确认十八个月、卡灵顿级的窗口期到底是不是正确的，甚至他还想将三宅级的时间也给推出来。
他抽出桌底那沓空白草稿纸。
李东压根没去等金陵大学超算的批文。
他要走的这条路，跟算力没什么关系。
十一笔真信号的时间节点被他重新铺开。
但是样本其实还是太少了点，传统的傅里叶变换拉不开频率分辨率。
于是他在纸上飞快写下了一行字。
【Lomb-Scargle Periodogram】
Lomb-Scargle周期图。
这是天文学界专门拿来对付“不均匀采样时间序列”的老把式，五十年前由两位天文学家先后独立提出。
它不要求数据点之间间隔均匀，而是在频率空间里直接做最小二乘拟合。
低样本量、不均匀间隔、信号叠加在噪声里，这三件事，是天文学家最难受的事。
而Lomb-Scargle就是为这三件事量身定制的。
李东拿起笔，开始倒推频率。
T0≈ 73.4天。
也就是说，如果这十一笔事件背后藏着一只钟，那这只钟，每七十三天敲一下。
对应的就是太阳巨大的磁场势阱底下，单位钟摆周期内，外面的高能粒子撞击势垒一次的频率。
它没法直接告诉你主峰何时爆发。
但只要再叠上一层东西，主峰自身的“激发概率”就能反推出来。
这层东西就是……
伽莫夫的量子隧穿衰变框架。
不过需要改上一点……
核内势阱改“太阳磁场势阱”、α粒子改“高能耀斑外溢通量”……李东在草稿纸上魔改着微分方程。
十一笔事件的振幅递增曲线被套入方程，MCMC（马尔可夫链蒙特卡洛）算法瞬间跑完五万步采样。
后验分布的均值，最后钉死在了一个坐标上。
【主峰抵达窗口期：T0+ 18.3个月（± 2.4个月）】
【量级：卡灵顿级】
【后验置信度：78.6%】
和昨晚的结果几乎一样，但这一次，它的可信度更高了。
就在这时，桌上的手机突然震动了起来，是个陌生号码。
“李东老师吗？我是金陵大学超算中心的，您申请的大规模回测节点批下来了，A100满规格，请问您什么时候来开账号？”
李东愣了一下，看了一眼屏幕右下角打饱嗝的小黑。
“不用了，老师。”李东笑了笑。
“计划有变，这节点让给别人排吧，我已经跑完了。”
电话那头一时间安静的有些诡异。
两百万条全量数据，你一天不到跑完了？拿什么跑的？算盘吗？
没等对方回过神，李东直接挂断了电话。
草稿纸开始散落在房间的各个地方，书桌、地板……
李东要把±2.4个月的窗口期继续往里压！
他开始引入羊八井观测站三年来的“待复核”异常邮件。
三笔同向事件叠入MCMC！
置信度飙升至83.4%！窗口期压至±1.9个月！
然后再引入LHAASO（拉索）观测站的TeV级宇宙线月报。
置信度推到86.7%！窗口期死锁在±1.5个月！
李东所有基础属性直接拉满，彻底进入了那种“洞穿真理”的专注状态。
天黑天亮，门外的风吹草动，全被他屏蔽了。
他的世界里只剩下那只在地球上空敲了八年的“钟”。
他正一笔一笔地，把死神敲门的倒计时，从宇宙的噪声里给硬生生抠出来！
……
两天后。
金陵某高层公寓。
“咔哒”一声，备用钥匙拧开了防盗门。
沈澈和张文平眉头紧锁地走了进来。
屋子里的气味有些发霉。
书桌前，一个人影正埋在几乎要将他淹没的草稿纸堆里。
“李东？”
李东没有回应。
两人这才走到他的身前，沈澈拍了拍他的肩膀。
李东猛地一抖，然后缓缓抬起头。
“沈老师……张老师……”
他的声音很沙哑。
“你这两天连电话都不接，就为了算这个？”
张文平看着满地的草稿纸，又好气又好笑。
“咕噜……”
李东的肚子突然叫了一声，他有些不好意思的说道。
“算入神了，忘了吃东西。”
沈澈看着他这幅亢奋又憔悴的模样，心里暗暗叹息。
国际上喊“狼来了”七八次，都没兑现。
李东把这推论看得太重，万一算空了，这孩子不得崩溃？
“李东啊，那个预测的事你先放一放把。”
沈澈语气温和，带着劝慰。
“山上还有一堆正篇论文的待复核数据等着你把关，你先好好休息下吧。”
可李东却摇了摇头。
“沈老师，不用了。”
沈澈一愣，不用了？破罐子破摔了？
“咱们预研组的事。”
李东指了指桌面右侧一叠码放整齐的草稿纸。
“我在推预测时间的时候顺手给推出来了。”
沈澈下意识的接话。
“就算你推出来了，可……”
他话还没说完，突然就傻住了。
因为他看到了那叠草稿纸最上面的一行字。
沈澈猛地扑过去，一把抓起那叠纸。
第一页：李判据沿能段往上爬的完整轨迹图。
25 GeV到4.6 TeV，锚定边界清清楚楚。
第二页：一百九十六万条全量数据复核。
【真信号：四万一千七百九十二条。】
【主物理结构集中区间：0.871 TeV— 1.000 TeV。】
【转折能量Eb=914± 98 GeV。拟合优度χ2/d.o.f.=23.3/18。】
【相对单一幂律模型偏好显著性：6.6σ。】
最后一页：一张系统误差修正后的最终能谱曲线。
那一段从25 GeV一路平滑下来的能谱到了0.9 TeV那一处，“咔嚓”一声，斜率突然变陡了。
不是幻峰。
是真的。
是一个真真切切的、卡在TeV那一档的转折结构。
可能外人根本无法理解，这个在0.9 TeV处斜率突然从γ1≈3.1剧烈转折到γ2≈3.9的“断点（Break）”，到底意味着什么。
毫不夸张地说，这是近十几年高能天体物理学界所有人都在眼红的“圣杯”！
上一篇：我家艺人太没上进心了
下一篇：返回列表