第304章 顺手的事
赘婿出山小说推荐阅读：Apop之我在首尔当外教、赘婿出山、泥泞、股神传奇、乡村中医李茂阳、开局账号被盗，反手充值一百万、极品全能学生、推拿、最强狂兵、带着超市重返年代、都市：我的器官能自己强化、不装了，我有五个天仙姐姐！、国宝神鉴
第二天清晨。
李东合上电脑时，窗外的天色已经微微亮了。
草稿纸散落一地，他揉了揉有些发胀打的太阳穴，然後看向了屏幕右下角。
小黑头上的进度条已经来到了【0.18】。
【主人，一百九十六万条，全吃完啦。嗝隔……】
李东失笑，敲下键盘。
【真信号多少？】
【八万六千三百九十二条，其余全是幻峰。】
李东点点头。
【小黑，挂载昨天的脚本，叠第二层筛子。】
【在十几个TeV能段，找形状符合AD 774/5外推预测的事件，x/.f.卡在1.1以内。】小黑没有废话，又开始蹦蹦跳跳起来。
到二十秒，一行字弹了出来。
【主人，挑好啦，一共十一笔。】
李东嘴角抽了一下。
十一笔。
跟昨晚他用那三十七笔“待复核”事件包外推出来的数字，分毫不差！
如果说昨晚是还是弱证据的话，那现在，用悟空号八年全量数据跑出来的结果，就是强证据了。深吸一口气，李东知道，这第一关算是过了。
接下来才是真正的难搞的事。
他要确认十八个月、卡灵顿级的窗口期到底是不是正确的，甚至他还想将三宅级的时间也给推出来。他抽出桌底那遝空白草稿纸。
李东压根没去等金陵大学超算的批文。
他要走的这条路，跟算力没什麽关系。
十一笔真信号的时间节点被他重新铺开。
但是样本其实还是太少了点，传统的傅里叶变换拉不开频率分辨率。
於是他在纸上飞快写下了一行字。
【Lomb-ScarglePer iodogram】
Lomb-Scargle周期图。
这是天文学界专门拿来对付“不均匀采样时间序列”的老把式，五十年前由两位天文学家先後独立提出。
它不要求数据点之间间隔均匀，而是在频率空间里直接做最小二乘拟合。
低样本量、不均匀间隔、信号叠加在噪声里，这三件事，是天文学家最难受的事。
而Lomb-Scargle就是为这三件事量身定制的。
李东拿起笔，开始倒推频率。
To～ 73.4天。
也就是说，如果这十一笔事件背後藏着一只锺，那这只锺，每七十三天敲一下。
对应的就是太阳巨大的磁场势阱底下，单位锺摆周期内，外面的高能粒子撞击势垒一次的频率。它没法直接告诉你主峰何时爆发。
但只要再叠上一层东西，主峰自身的“激发概率”就能反推出来。
这层东西就是……
伽莫夫的量子隧穿衰变框架。
不过需要改上一点……
核内势阱改“太阳磁场势阱”、α粒子改“高能耀斑外溢通量”……李东在草稿纸上魔改着微分方程。十一笔事件的振幅递增曲线被套入方程，MCMC（马尔可夫链蒙特卡洛）算法瞬间跑完五万步采样。後验分布的均值，最後钉死在了一个坐标上。
【主峰抵达窗口期：T。+18.3个月（±2.4个月）】
【量级：卡灵顿级】
【後验置信度：78.6%】
和昨晚的结果几乎一样，但这一次，它的可信度更高了。
就在这时，桌上的手机突然震动了起来，是个陌生号码。
“李东老师吗？我是金陵大学超算中心的，您申请的大规模回测节点批下来了，A100满规格，请问您什麽时候来开账号？”
李东愣了一下，看了一眼屏幕右下角打饱嗝的小黑。
“不用了，老师。”李东笑了笑。
“计划有变，这节点让给别人排吧，我已经跑完了。”
电话那头一时间安静的有些诡异。
两百万条全量数据，你一天不到跑完了？拿什麽跑的？算盘吗？
没等对方回过神，李东直接挂断了电话。
草稿纸开始散落在房间的各个地方，书桌、地板……
李东要把±2.4个月的窗口期继续往里压！
他开始引入羊八井观测站三年来的“待复核”异常邮件。
三笔同向事件叠入MCMC！
置信度飙升至83. 4%！窗口期压至±1.9个月！
然後再引入LHAAS0（拉索）观测站的TeV级宇宙线月报。
置信度推到86.7%！窗口期死锁在土1.5个月！
李东所有基础属性直接拉满，彻底进入了那种“洞穿真理”的专注状态。
天黑天亮，门外的风吹草动，全被他屏蔽了。
他的世界里只剩下那只在地球上空敲了八年的“锺”。
他正一笔一笔地，把死神敲门的倒计时，从宇宙的噪声里给硬生生抠出来！
两天後。
金陵某高层公寓。
“哢哒”一声，备用钥匙拧开了防盗门。
沈澈和张文平眉头紧锁地走了进来。
屋子里的气味有些发霉。
书桌前，一个人影正埋在几乎要将他淹没的草稿纸堆里。
“李东？”
李东没有回应。
两人这才走到他的身前，沈澈拍了拍他的肩膀。
李东猛地一抖，然後缓缓抬起头。
“沈老师……张老师……”
他的声音很沙哑。
“你这两天连电话都不接，就为了算这个？”
张文平看着满地的草稿纸，又好气又好笑。
“咕噜………”
李东的肚子突然叫了一声，他有些不好意思的说道。
“算入神了，忘了吃东西。”
沈澈看着他这幅亢奋又憔悴的模样，心里暗暗叹息。
国际上喊“狼来了”七八次，都没兑现。
李东把这推论看得太重，万一算空了，这孩子不得崩溃？
“李东啊，那个预测的事你先放一放把。”
沈澈语气温和，带着劝慰。
“山上还有一堆正篇论文的待复核数据等着你把关，你先好好休息下吧。”
可李东却摇了摇头。
“沈老师，不用了。”
沈澈一愣，不用了？破罐子破摔了？
“咱们预研组的事。”
李东指了指桌面右侧一叠码放整齐的草稿纸。
“我在推预测时间的时候顺手给推出来了。”
沈澈下意识的接话。
“就算你推出来了，可……”
他话还没说完，突然就傻住了。
因为他看到了那叠草稿纸最上面的一行字。
沈澈猛地扑过去，一把抓起那叠纸。
第一页：李判据沿能段往上爬的完整轨迹图。
25 GeV到4. 6 TeV，锚定边界清清楚楚。
第二页：一百九十六万条全量数据复核。
【真信号：四万一千七百九十二条。】
【主物理结构集中区间：0.871 TeV-- 1.000 TeV。】
【转折能量Eb=914±98GeV。拟合优度 x/d. o.f.=23.3/18。】【相对单一幂律模型偏好显着性：6.60。】
最後一页：一张系统误差修正後的最终能谱曲线。
那一段从25GeV一路平滑下来的能谱到了0.9TeV那一处，“哢嚓”一声，斜率突然变陡了。不是幻峰。
是真的。
是一个真真切切的、卡在TeV那一档的转折结构。
可能外人根本无法理解，这个在0.9 TeV处斜率突然从y.～3.1剧烈转折到v.～3.9的“断点（Break）”，到底意味着什麽。
毫不夸张地说，这是近十几年高能天体物理学界所有人都在眼红的“圣杯”！
为什麽？
因为在这个TeV级转折的背後，极有可能藏着整个物理学界苦苦追寻了半个世纪的终极答案一一暗物质粒子湮灭或衰变的直接线索，亦或是某颗未知的邻近脉冲星！它能直接终结学术界关於“正电子异常”的百年争论。
为了找到这个拐点，全世界的科研机构早就杀红了眼。
欧洲的地面切伦科夫望远镜阵列（H. E. .）曾经“间接”瞥见过一点模糊的影子，但系统误差太大了美国的费米卫星（Fermi-LAT）在天上盯了那麽多年，硬是因为粒子背景污染太高、仪器能量分辨率不够，活生生错过了这个TeV级别的转折！
所有人都知道那里可能有东西，但在此之前，全世界没有任何一个团队，能拿出一张铁证如山的“直接探测”能谱！
直到今天。
直到悟空号（DAMPE）带着史无前例的高能量分辨率和极低本底，用这张被李东硬生生从底噪里洗出来的图谱，一锤定音！
这不仅是人类历史上第一次“直接探测”到TeV宇宙线电子能谱的断折。
这更是一张注定要轰动全球、预定《Nature》封面的王炸！
沈澈拿着纸的手都在微微颤抖，他倒退了一步，差点没站稳。
“老沈！怎麽了？”
张文平赶紧扶住他。
沈澈有些颤抖的说道。
“老张……我先回山上，你盯着他，别让他饿死了！”
张文平愣了半秒，无奈地点了点头。
看这个样子，李东好像是真的把老沈他们组那个憋了五六年的事，给一锅端了。
张文平虽然也是这个组里的一员，可他其实是来帮忙的，跟悟空号没什麽深的瓜葛。
他不像沈澈，是真的在这上面耗费了几年的心血。
所以他理解沈澈现在的心情。
也不在意自己一个堂堂院士，居然成了李东的临时保姆。
沈澈把那一摞草稿纸朝胸口一抱，几乎是用跑的冲出了房门。
门“砰”地一声关上。
张文平转过头，想问李东要吃点什麽。
可一回头，他就发现……
李东又埋头在算什麽东西去了。
张文平：……
他也没去打扰李东。
只是默默地掏出手机，给自己一个还在金陵大学读博後的学生拨了个电话。
“喂，小赵啊，去川香府打包两份套餐送过来……”
陈越：川香府？？？
日子过得飞快。
李东这一阵几乎就没出过那间公寓的门。
张文平的博後小赵每天定时送饭，送完就把外卖盒收走。
李东在公寓里头算。
沈澈在山上推进。
那一摞从李东桌子上抱走的草稿纸被沈澈连夜带回山上以後。
老周、老吴、老郑、几个助研、博後，那是连轴转。
把李东草稿纸上的每一个步骤、每一组数字，重新在悟空号的全量数据上独立跑了一遍。
跑到第十二天的时候。
每一组复核结果都跟李东草稿纸上的数字对上了。
跑到第十五天的时候……
十月一日国庆节到了。
紫金山，暗物质粒子探测协同创新中心。
主控大厅。
大厅最里头那块巨大的屏幕上，一条从25 GeV一路平滑下来的能谱曲线，正缓缓地往4.6 TeV那一头延伸。
每往前推一格，曲线下面那条系统误差带就在跟着收紧。
屋子里站满了人。
老周、老吴、老郑、几个助研、博後……
所有人都没说话。
整个大厅里安静得只能听见风扇转动的声音。
曲线一直往前推。
推到0.87TeV的位置斜率开始往下沉。
推到0.91 TeV的时候斜率“哢嚓”一声，明显地往下折了一档！
那一段转折结构，没有幻峰的毛刺。
就是一条清清楚楚的一转折。
最後一个像素点落定的瞬间。
屏幕底下又跳出来一行小字。
【单一幂律vs平滑断点幂律：6.6g。】
沈澈慢慢地转过身来。
朝大厅里头那一圈人，轻轻地点了一下头。
随着他的点头，整个大厅炸了。
欢呼声、拍桌子的声音、有人激动得直接给身边人来了一个熊抱、还有几个年纪大些的研究员当场抹起了眼泪……
老周抓着旁边那位博後的胳膊，连声音都在抖。
“立住了！”
“五年了！”
“五年了啊！”
“咱们这一群人，可算把这条曲线立住了！”
就在这时候。
大厅最角落里，常津朝着沈澈招了招手。
小会客间里。
常津望着窗外那一片刚开始转黄的梧桐，过了好半天，才开口。
“老沈啊。”
“你们这个项目前前後後……搞了多久？”
沈澈被这一问，愣了一下。
他下意识地朝身後那扇门看了一眼。
门里那些欢呼声还隐隐地透出来。
“常台长。”
“您是指悟空号正篇？”
“还是指………”
常津挥了挥手。
“我指的就是从你们把李判据这把尺子接进来，到今天这条曲线立住。”
“一共多久？”
沈澈眨了眨眼睛。
他在脑子里默默地算了一下。
“愿……”
“一个月零三天。”
常津半天只憋出两个字：牛逼！
http://www.zhuixuchushan.com/yt124292/48927900.html
请记住本书首发域名：www.zhuixuchushan.com。赘婿出山手机版阅读网址：www.zhuixuchushan.com
李东合上电脑时，窗外的天色已经微微亮了。
草稿纸散落一地，他揉了揉有些发胀打的太阳穴，然後看向了屏幕右下角。
小黑头上的进度条已经来到了【0.18】。
【主人，一百九十六万条，全吃完啦。嗝隔……】
李东失笑，敲下键盘。
【真信号多少？】
【八万六千三百九十二条，其余全是幻峰。】
李东点点头。
【小黑，挂载昨天的脚本，叠第二层筛子。】
【在十几个TeV能段，找形状符合AD 774/5外推预测的事件，x/.f.卡在1.1以内。】小黑没有废话，又开始蹦蹦跳跳起来。
到二十秒，一行字弹了出来。
【主人，挑好啦，一共十一笔。】
李东嘴角抽了一下。
十一笔。
跟昨晚他用那三十七笔“待复核”事件包外推出来的数字，分毫不差！
如果说昨晚是还是弱证据的话，那现在，用悟空号八年全量数据跑出来的结果，就是强证据了。深吸一口气，李东知道，这第一关算是过了。
接下来才是真正的难搞的事。
他要确认十八个月、卡灵顿级的窗口期到底是不是正确的，甚至他还想将三宅级的时间也给推出来。他抽出桌底那遝空白草稿纸。
李东压根没去等金陵大学超算的批文。
他要走的这条路，跟算力没什麽关系。
十一笔真信号的时间节点被他重新铺开。
但是样本其实还是太少了点，传统的傅里叶变换拉不开频率分辨率。
於是他在纸上飞快写下了一行字。
【Lomb-ScarglePer iodogram】
Lomb-Scargle周期图。
这是天文学界专门拿来对付“不均匀采样时间序列”的老把式，五十年前由两位天文学家先後独立提出。
它不要求数据点之间间隔均匀，而是在频率空间里直接做最小二乘拟合。
低样本量、不均匀间隔、信号叠加在噪声里，这三件事，是天文学家最难受的事。
而Lomb-Scargle就是为这三件事量身定制的。
李东拿起笔，开始倒推频率。
To～ 73.4天。
也就是说，如果这十一笔事件背後藏着一只锺，那这只锺，每七十三天敲一下。
对应的就是太阳巨大的磁场势阱底下，单位锺摆周期内，外面的高能粒子撞击势垒一次的频率。它没法直接告诉你主峰何时爆发。
但只要再叠上一层东西，主峰自身的“激发概率”就能反推出来。
这层东西就是……
伽莫夫的量子隧穿衰变框架。
不过需要改上一点……
核内势阱改“太阳磁场势阱”、α粒子改“高能耀斑外溢通量”……李东在草稿纸上魔改着微分方程。十一笔事件的振幅递增曲线被套入方程，MCMC（马尔可夫链蒙特卡洛）算法瞬间跑完五万步采样。後验分布的均值，最後钉死在了一个坐标上。
【主峰抵达窗口期：T。+18.3个月（±2.4个月）】
【量级：卡灵顿级】
【後验置信度：78.6%】
和昨晚的结果几乎一样，但这一次，它的可信度更高了。
就在这时，桌上的手机突然震动了起来，是个陌生号码。
“李东老师吗？我是金陵大学超算中心的，您申请的大规模回测节点批下来了，A100满规格，请问您什麽时候来开账号？”
李东愣了一下，看了一眼屏幕右下角打饱嗝的小黑。
“不用了，老师。”李东笑了笑。
“计划有变，这节点让给别人排吧，我已经跑完了。”
电话那头一时间安静的有些诡异。
两百万条全量数据，你一天不到跑完了？拿什麽跑的？算盘吗？
没等对方回过神，李东直接挂断了电话。
草稿纸开始散落在房间的各个地方，书桌、地板……
李东要把±2.4个月的窗口期继续往里压！
他开始引入羊八井观测站三年来的“待复核”异常邮件。
三笔同向事件叠入MCMC！
置信度飙升至83. 4%！窗口期压至±1.9个月！
然後再引入LHAAS0（拉索）观测站的TeV级宇宙线月报。
置信度推到86.7%！窗口期死锁在土1.5个月！
李东所有基础属性直接拉满，彻底进入了那种“洞穿真理”的专注状态。
天黑天亮，门外的风吹草动，全被他屏蔽了。
他的世界里只剩下那只在地球上空敲了八年的“锺”。
他正一笔一笔地，把死神敲门的倒计时，从宇宙的噪声里给硬生生抠出来！
两天後。
金陵某高层公寓。
“哢哒”一声，备用钥匙拧开了防盗门。
沈澈和张文平眉头紧锁地走了进来。
屋子里的气味有些发霉。
书桌前，一个人影正埋在几乎要将他淹没的草稿纸堆里。
“李东？”
李东没有回应。
两人这才走到他的身前，沈澈拍了拍他的肩膀。
李东猛地一抖，然後缓缓抬起头。
“沈老师……张老师……”
他的声音很沙哑。
“你这两天连电话都不接，就为了算这个？”
张文平看着满地的草稿纸，又好气又好笑。
“咕噜………”
李东的肚子突然叫了一声，他有些不好意思的说道。
“算入神了，忘了吃东西。”
沈澈看着他这幅亢奋又憔悴的模样，心里暗暗叹息。
国际上喊“狼来了”七八次，都没兑现。
李东把这推论看得太重，万一算空了，这孩子不得崩溃？
“李东啊，那个预测的事你先放一放把。”
沈澈语气温和，带着劝慰。
“山上还有一堆正篇论文的待复核数据等着你把关，你先好好休息下吧。”
可李东却摇了摇头。
“沈老师，不用了。”
沈澈一愣，不用了？破罐子破摔了？
“咱们预研组的事。”
李东指了指桌面右侧一叠码放整齐的草稿纸。
“我在推预测时间的时候顺手给推出来了。”
沈澈下意识的接话。
“就算你推出来了，可……”
他话还没说完，突然就傻住了。
因为他看到了那叠草稿纸最上面的一行字。
沈澈猛地扑过去，一把抓起那叠纸。
第一页：李判据沿能段往上爬的完整轨迹图。
25 GeV到4. 6 TeV，锚定边界清清楚楚。
第二页：一百九十六万条全量数据复核。
【真信号：四万一千七百九十二条。】
【主物理结构集中区间：0.871 TeV-- 1.000 TeV。】
【转折能量Eb=914±98GeV。拟合优度 x/d. o.f.=23.3/18。】【相对单一幂律模型偏好显着性：6.60。】
最後一页：一张系统误差修正後的最终能谱曲线。
那一段从25GeV一路平滑下来的能谱到了0.9TeV那一处，“哢嚓”一声，斜率突然变陡了。不是幻峰。
是真的。
是一个真真切切的、卡在TeV那一档的转折结构。
可能外人根本无法理解，这个在0.9 TeV处斜率突然从y.～3.1剧烈转折到v.～3.9的“断点（Break）”，到底意味着什麽。
毫不夸张地说，这是近十几年高能天体物理学界所有人都在眼红的“圣杯”！
为什麽？
因为在这个TeV级转折的背後，极有可能藏着整个物理学界苦苦追寻了半个世纪的终极答案一一暗物质粒子湮灭或衰变的直接线索，亦或是某颗未知的邻近脉冲星！它能直接终结学术界关於“正电子异常”的百年争论。
为了找到这个拐点，全世界的科研机构早就杀红了眼。
欧洲的地面切伦科夫望远镜阵列（H. E. .）曾经“间接”瞥见过一点模糊的影子，但系统误差太大了美国的费米卫星（Fermi-LAT）在天上盯了那麽多年，硬是因为粒子背景污染太高、仪器能量分辨率不够，活生生错过了这个TeV级别的转折！
所有人都知道那里可能有东西，但在此之前，全世界没有任何一个团队，能拿出一张铁证如山的“直接探测”能谱！
直到今天。
直到悟空号（DAMPE）带着史无前例的高能量分辨率和极低本底，用这张被李东硬生生从底噪里洗出来的图谱，一锤定音！
这不仅是人类历史上第一次“直接探测”到TeV宇宙线电子能谱的断折。
这更是一张注定要轰动全球、预定《Nature》封面的王炸！
沈澈拿着纸的手都在微微颤抖，他倒退了一步，差点没站稳。
“老沈！怎麽了？”
张文平赶紧扶住他。
沈澈有些颤抖的说道。
“老张……我先回山上，你盯着他，别让他饿死了！”
张文平愣了半秒，无奈地点了点头。
看这个样子，李东好像是真的把老沈他们组那个憋了五六年的事，给一锅端了。
张文平虽然也是这个组里的一员，可他其实是来帮忙的，跟悟空号没什麽深的瓜葛。
他不像沈澈，是真的在这上面耗费了几年的心血。
所以他理解沈澈现在的心情。
也不在意自己一个堂堂院士，居然成了李东的临时保姆。
沈澈把那一摞草稿纸朝胸口一抱，几乎是用跑的冲出了房门。
门“砰”地一声关上。
张文平转过头，想问李东要吃点什麽。
可一回头，他就发现……
李东又埋头在算什麽东西去了。
张文平：……
他也没去打扰李东。
只是默默地掏出手机，给自己一个还在金陵大学读博後的学生拨了个电话。
“喂，小赵啊，去川香府打包两份套餐送过来……”
陈越：川香府？？？
日子过得飞快。
李东这一阵几乎就没出过那间公寓的门。
张文平的博後小赵每天定时送饭，送完就把外卖盒收走。
李东在公寓里头算。
沈澈在山上推进。
那一摞从李东桌子上抱走的草稿纸被沈澈连夜带回山上以後。
老周、老吴、老郑、几个助研、博後，那是连轴转。
把李东草稿纸上的每一个步骤、每一组数字，重新在悟空号的全量数据上独立跑了一遍。
跑到第十二天的时候。
每一组复核结果都跟李东草稿纸上的数字对上了。
跑到第十五天的时候……
十月一日国庆节到了。
紫金山，暗物质粒子探测协同创新中心。
主控大厅。
大厅最里头那块巨大的屏幕上，一条从25 GeV一路平滑下来的能谱曲线，正缓缓地往4.6 TeV那一头延伸。
每往前推一格，曲线下面那条系统误差带就在跟着收紧。
屋子里站满了人。
老周、老吴、老郑、几个助研、博後……
所有人都没说话。
整个大厅里安静得只能听见风扇转动的声音。
曲线一直往前推。
推到0.87TeV的位置斜率开始往下沉。
推到0.91 TeV的时候斜率“哢嚓”一声，明显地往下折了一档！
那一段转折结构，没有幻峰的毛刺。
就是一条清清楚楚的一转折。
最後一个像素点落定的瞬间。
屏幕底下又跳出来一行小字。
【单一幂律vs平滑断点幂律：6.6g。】
沈澈慢慢地转过身来。
朝大厅里头那一圈人，轻轻地点了一下头。
随着他的点头，整个大厅炸了。
欢呼声、拍桌子的声音、有人激动得直接给身边人来了一个熊抱、还有几个年纪大些的研究员当场抹起了眼泪……
老周抓着旁边那位博後的胳膊，连声音都在抖。
“立住了！”
“五年了！”
“五年了啊！”
“咱们这一群人，可算把这条曲线立住了！”
就在这时候。
大厅最角落里，常津朝着沈澈招了招手。
小会客间里。
常津望着窗外那一片刚开始转黄的梧桐，过了好半天，才开口。
“老沈啊。”
“你们这个项目前前後後……搞了多久？”
沈澈被这一问，愣了一下。
他下意识地朝身後那扇门看了一眼。
门里那些欢呼声还隐隐地透出来。
“常台长。”
“您是指悟空号正篇？”
“还是指………”
常津挥了挥手。
“我指的就是从你们把李判据这把尺子接进来，到今天这条曲线立住。”
“一共多久？”
沈澈眨了眨眼睛。
他在脑子里默默地算了一下。
“愿……”
“一个月零三天。”
常津半天只憋出两个字：牛逼！
http://www.zhuixuchushan.com/yt124292/48927900.html
请记住本书首发域名：www.zhuixuchushan.com。赘婿出山手机版阅读网址：www.zhuixuchushan.com