作者：胖胖的小橘
而這個黑黢黢的小東西……用了不到三秒。
第304章 順手的事
第二天清晨。
李東合上電腦時，窗外的天色已經微微亮了。
草稿紙散落一地，他揉了揉有些發脹打的太陽穴，然後看向了螢幕右下角。
小黑頭上的進度條已經來到了【0.18】。
【主人，一百九十六萬條，全吃完啦。嗝……】
李東失笑，敲下鍵盤。
【真訊號多少？】
【八萬六千三百九十二條，其餘全是幻峰。】
李東點點頭。
【小黑，掛載昨天的指令碼，疊第二層篩子。】
【在十幾個TeV能段，找形狀符合AD 774/5外推預測的事件，χ?/d.o.f.卡在1.1以內。】
小黑沒有廢話，又開始蹦蹦跳跳起來。
到二十秒，一行字彈了出來。
【主人，挑好啦，一共十一筆。】
李東嘴角抽了一下。
十一筆。
跟昨晚他用那三十七筆“待複核”事件包外推出來的數字，分毫不差！
如果說昨晚是還是弱證據的話，那現在，用悟空號八年全量資料跑出來的結果，就是強證據了。
深吸一口氣，李東知道，這第一關算是過了。
接下來才是真正的難搞的事。
他要確認十八個月、卡靈頓級的視窗期到底是不是正確的，甚至他還想將三宅級的時間也給推出來。
他抽出桌底那沓空白草稿紙。
李東壓根沒去等金陵大學超算的批文。
他要走的這條路，跟算力沒什麼關係。
十一筆真訊號的時間節點被他重新鋪開。
但是樣本其實還是太少了點，傳統的傅立葉變換拉不開頻率解析度。
於是他在紙上飛快寫下了一行字。
【Lomb-Scargle Periodogram】
Lomb-Scargle週期圖。
這是天文學界專門拿來對付“不均勻取樣時間序列”的老把式，五十年前由兩位天文學家先後獨立提出。
它不要求資料點之間間隔均勻，而是在頻率空間裡直接做最小二乘擬合。
低樣本量、不均勻間隔、訊號疊加在噪聲裡，這三件事，是天文學家最難受的事。
而Lomb-Scargle就是為這三件事量身定製的。
李東拿起筆，開始倒推頻率。
T?≈ 73.4天。
也就是說，如果這十一筆事件背後藏著一隻鍾，那這隻鍾，每七十三天敲一下。
對應的就是太陽巨大的磁場勢阱底下，單位鐘擺週期內，外面的高能粒子撞擊勢壘一次的頻率。
它沒法直接告訴你主峰何時爆發。
但只要再疊上一層東西，主峰自身的“激發機率”就能反推出來。
這層東西就是……
伽莫夫的量子隧穿衰變框架。
不過需要改上一點……
核內勢阱改“太陽磁場勢阱”、α粒子改“高能耀斑外溢通量”……李東在草稿紙上魔改著微分方程。
十一筆事件的振幅遞增曲線被套入方程，MCMC（馬爾可夫鏈蒙特卡洛）演算法瞬間跑完五萬步取樣。
後驗分佈的均值，最後釘死在了一個座標上。
【主峰抵達視窗期：T?+ 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。擬合優度χ?/d.o.f.=23.3/18。】
【相對單一冪律模型偏好顯著性：6.6σ。】
最後一頁：一張系統誤差修正後的最終能譜曲線。
那一段從25 GeV一路平滑下來的能譜到了0.9 TeV那一處，“咔嚓”一聲，斜率突然變陡了。
不是幻峰。
是真的。
是一個真真切切的、卡在TeV那一檔的轉折結構。
可能外人根本無法理解，這個在0.9 TeV處斜率突然從γ?≈3.1劇烈轉折到γ?≈3.9的“斷點（Break）”，到底意味著什麼。
毫不誇張地說，這是近十幾年高能天體物理學界所有人都在眼紅的“聖盃”！
上一篇：我家艺人太没上进心了
下一篇：返回列表