第二十六年的整个秋天,江临都在做一件看起来很奇怪的事。
看sweet那条孤独的电流片,像一条横亘在物理学史上的冰冷伤痕。
看petschek那个被x標记的几何,优雅却带著理想化的脆弱。
看plasmoid那串绿色的磁岛標记,如同沸腾断裂的珠串。
看hall那个被打上不可触及的高点的位置。
看湍流重联旁边那行主动选择不补的红字,透著某种无奈的妥协。
看西南角那个最小最暗的新標记——【?-multi-loop-coupled】。
看那行更小的字——【变电站工程规格异常。铭牌残字:tm-,网,配电,长城。含义未明。】
每一次看完,他的目光最终都会停在那个新標记上。
那些从变电站废墟里带回来的剩磁数据,像是一根根细小的倒刺,扎在他过去几十年建立起来的理论框架上。
第二十六年十月的一个下午,他终於做了一件事。
他在桌前坐下,铺开一张大白纸,在纸的正中央写了一个词。
【磁重联】
然后围绕著开始画。
工作站上当然有更好的绘图软体,有更精密的拓扑建模工具,但他这一次不想要软体帮他整理思路。
软体的底层逻辑是確定的算法,而他此刻需要的,是思路自己长出来。
第一条线从中心向上,线的末端写sweet-parker。
第二条线向右上,末端写petschek。
第三条线向右,plasmoid磁岛串。
第四条线向右下,hall重联。
第五条线向下,无碰撞重联。
第六条线向左下,湍流重联。
第七条线向左,撕裂模/锯齿崩塌。
第八条线向左上,?-multi-loop-coupled。
八条线,八个流派。
每一个都是磁重联在某种边界条件下的一种具体实现。
如果在现实世界的大学讲堂里,教材会把这八个流派排成一条清晰的歷史时间线。
从1957年的sweet-parker开始,提出长电流片模型。
到1964年的petschek,引入慢衝击波加速重联率。
至2007年的plasmoid,发现长电流片的次级失稳。
於2016年的mms探测器在电子扩散区拿到实测数据。
……
按时间,按发现顺序排序。
在这种敘事里,一个流派的出现,往往被理解为对另一个流派的超越或修正。
后人总是比前人更接近真理。
但江临在过去十年里读完八十二篇a类核心论文之后,他知道这种排序方式太单线了。
它適合教学,不適合描述一个真实灾难现场里多尺度机制的同时打开。
变电站的剩磁数据告诉他,真实灾难现场里,八个流派不一定同时完整出现。
但真实的物理系统,绝不可能像做单选题一样,只选择其中一个机制。
它会在不同尺度,不同边界,不同耗散条件下,让不同机制先后打开,重叠,竞爭。
相互嵌套,相互干涉,谁也不超越谁,谁也不取代谁。
如果把变电站这个一万二千平方米的局部残骸系统,外推到整个废土星球的导电网络,甚至外推到太阳风与地球磁层交匯的宏观空间。
多数机制都有可能在不同尺度上先后打开,重叠,竞爭。
它们未必同时完整出现,但真实系统绝不只选择其中一个。
sweet-parker的缓慢电流片可能存在於某个宏观的背景尺度中,默默积蓄能量。
petschek的x型几何可能出现在另一个次级尺度的能量释放点。
plasmoid磁岛串可能在第三个更小的尺度上,像连珠炮一样裂开。
hall效应在离子惯性长度附近开始让电子和离子的运动分道扬鑣,再往下,电子扩散区才是真正不可触及的高点。
湍流可能像一场永不停歇的风暴,贯穿从宏观到微观的所有尺度。
撕裂模描述的是某些受边界和共振面约束的局部崩塌行为,在托卡马克里,它尤其清楚。
它们都在,同时在,组成一个多通道,多尺度,相互嵌套的系统。
而江临,即便已经把物理学史翻了个底朝天,所有这些流派加起来,显然都还是不够。
因为还有一个流派,在他北墙的研究史地图上,连一个完整的定义都没有。
那个最小最暗的【?-multi-loop-coupled】。
它的意思是,当系统的边界条件不再是理想化的无限大等离子体,而是一整张如同人类电网般错综复杂的闭合导体网络时,会发生什么?
这就像是问,当雷霆劈中一张巨大的铁丝网,这张网会如何反噬天空?
教材没有討论这个问题。
至少在他翻烂了的那些a类核心论文里,没有一篇把这个问题作为主轴。
学术界更喜欢乾净的等离子体腔室,而不是生锈的变电站残骸。
但变电站的数据,用一种冰冷的物证形式,逼问了这个问题。
江临盯著这张八角图,第一次模糊却又无比强烈地意识到,自己需要的不是在某一条线上再往前走一步,去搞出一个稍微精確一点的新模型。
他需要的,是一种新的语言。
一种能够描述这八个流派如何同时活动,如何相互耦合的新语言。
当脑海中浮现出这个念头並写下这句话时,他自己都有点惊讶。
他不知道这种语言长什么样,是张量微积分的某种变体?
还是某种全新的拓扑代数?
他只知道这种语言应该存在,因为废土上发生的真实物理,就在那里。
如果某种语言不存在,那绝不是因为它不可能,仅仅是因为,还没有人去发明它。
於是,江临开始动手寻找这种语言。
寻找物理学的语言,往往不能从纯理论开始。
他决定停掉所有论文阅读和实验室分析,先回到工程,回到那些看得见摸得著,能转动,会发热的金属造物中去。
每天花六到八个小时,坐在工作檯前,对著窗外的风机二號和远处的观测点a,画一张又一张的结构图。
这些图虽然形態各异,但標题都是同一个。
【前哨站电力系统的多通道结构】
风机二號,本身就是一个绝佳的多通道耦合系统。
叶轮捕获废土上肆虐的风能,第一个通道——气动通道。
叶轮通过轮轂和主轴直接拖动低速永磁发电机转子,第二个通道——机电转换通道。
发电机输出不稳定的三相交流电,经过整流桥转换成带纹波的直流,第三个通道——电力变换通道。
直流电经过控制器的智能分流,一部分充进电池(储能通道),一部分流向负载(消耗通道),一部分流向卸荷电阻(泄能通道)。
每一个通道都有自己的特徵时间尺度。
气动通道,风速变化的时间尺度,从秒到分钟,缓慢而宏大。
机电转换,电机电磁转矩的响应时间,从毫秒到秒。
电力变换,半导体开关的开断频率,微秒到毫秒,极快。
储能,电池化学反应的充放电时间,从分钟到小时,甚至几天。
消耗,负载的隨时切换,秒到分钟。
泄能,卸荷电阻的热平衡和散热时间,秒到分钟。
六个通道,六个截然不同的时间尺度,同时影响电力系统的整体行为。
教科书上写风力发电系统的时候,从来不会把所有六个通道清清楚楚地画在一张图上,並討论它们的耦合。
要么只讲气动,那是流体力学课的范畴,里面全是纳维-斯托克斯方程。
要么只讲机电,那是电机学,討论的是磁链和安培力。
要么只讲电力电子,那是功率器件课,全是功率开关器件的导通损耗和开关损耗。
要么只讲控制逻辑,那是系统工程课,玩的是pid算法。
每一门课,都像盲人摸象一样,只讲一个通道。
整个系统在真实环境下的多通道非线性行为,很少有一本入门教材会把这些通道放在同一张图里,认真討论它们在故障边界附近的耦合。
但江临在前哨站建立,维护,修復风机二號的过程中,亲眼看见过,亲手处理过这个系统的多通道耦合。
风速突变时,气动通道首先响应,叶轮转速剧烈波动。
转速波动通过机电通道传到发电机端,三相电压波动。
电压波动通过电力变换通道,导致直流母线电压跳动。
跳动通过储能通道,让电池充电电流变化。
变化通过控制器逻辑,触发或不触发卸荷。
整个链条,从风到卸荷电阻发热,所有通道都在同时活动。
而决定整个系统最终是安然无恙地稳定运行,还是在过载中轰然烧毁的,从来不是任何一个单独的通道。
而是它们之间的耦合关係。
哪个通道响应快,哪个通道响应慢。
哪个通道的输出溢出,会触发另一个通道的崩溃閾值?
哪个通道在过载时,会被另一个通道的冗余容量吃掉?
哪个通道有硬体保护,哪个通道是脆弱裸露的?
这就是真实的工程系统。
江临把电力系统的多通道结构图掛在研究史地图旁边。
於是他有了两张图。
一张是磁重联领域六十年的研究史。
一张是前哨站电力系统的工程结构。
它们一边是等离子体天体物理,一边是低端电气工程,看起来风马牛不相及,但在多通道耦合这件事上,两者说的是同一件事。
第二十六年深冬。
风雪肆虐,江临在石屋这个坚固的堡垒里面,开始尝试把前哨站电力系统的多通道思维,应用到磁重联上。
第一个尝试是,如果磁重联也是一个多通道系统,那它有几个通道?
每个通道对应著什么具体的等离子体物理过程?
每个通道的特徵尺度是什么?
通道之间到底是如何耦合的?
他在工作站里建了一个新文件夹。
【multi-channel-framework(多通道框架)】
文件夹里,目前只有一个可怜的文本。
【draft_v1.txt】
文件的第一行是——【磁重联的多通道分解尝试 v1】
往下,他庖丁解牛一般,一一列出候选的通道。
【通道一:大尺度mhd供能通道】
特徵:庞大的磁能从外部大尺度系统(如太阳日冕、地球磁尾),像水库放水一样,源源不断地输入到重联区域。
特徵尺度:系统宏观尺度l。
对应理论模型:sweet-parker的大尺度部分。
【通道二:中尺度plasmoid层裂通道】
特徵:当mhd通道输送的能量让长电流片不堪重负时,电流片自身发生撕裂失稳,碎裂成一长串被称为磁岛的结构。
特徵尺度:电流片长度的某个分数。
对应理论模型:loureiro (2007) 的链式不稳定性、bhattacharjee (2009) 的快重联机制。
【通道三:小尺度hall转捩通道】
特徵:在极微小的尺度下,沉重的离子跟不上轻盈电子的步伐,两者在磁场中退耦,產生了强烈的霍尔电场,改变了重联区的几何拓扑。
特徵尺度:离子惯性长度d_i。
对应理论模型:gem挑战赛確认的机制、burch (2016) mms探测器的实测。
【通道四:跨尺度湍流通道】
特徵:磁场本身的无序湍流,像一锅沸腾的热水,让大量微小尺度的重联事件在整个空间內隨机且同时发生,极大增加了耗散。
特徵尺度:从宏观能量注入尺度,一直到微观能量耗散尺度的整个惯性级联区间。
对应理论模型:lazarian-vishniac(1999)的隨机湍流重联模型。
【通道五:边界封闭通道】
特徵:在人造或特定的封闭几何(如托卡马克装置的环形室)中,磁场被容器强制保持某种拓扑约束,重联只能在特定的共振面上发生。
特徵尺度:边界容器的特徵尺度。
对应理论模型:托卡马克装置內的撕裂模。
【通道六:多迴路集体响应通道】
特徵:多个闭合的人造或天然导体迴路相互嵌套,在极端电磁扰动下,像一个庞大的lc振盪网络一样,產生集体响应和磁化烙印。
特徵尺度:整个导电网络的宏观尺度(可能覆盖整个行星地表)。
对应理论模型:尚无。
最后一行,带著一个刺眼的问號。
江临在【draft_v1】的下方,又补充了一行字。
“六个通道,前五个,都有现存的成熟理论可以支撑。第六个,是变电站的废墟用剩磁给我的残酷提示,但在人类的理论物理上,这依然是一片令人窒息的空白。”
写完保存,关掉编辑器,江临揉了揉发酸的眉心。
框架是搭起来了,但也仅仅是万里长征的第一步。
第二十七年,春。
冰雪消融,苔蘚焕发生机。
江临却陷入了试图把这六个通道之间的耦合关係说清楚的泥潭。
每两个通道之间的耦合,他都尝试在草纸上画一张小图。
这些图不是简单的连线,蕴含著严格逻辑因果的物理流程图。
他分析通道一(mhd)和通道二(plasmoid)的耦合,大尺度mhd把磁能输入,拉长电流片。
当电流片长宽比超过临界值,失稳发生,plasmoid產生。
这是一对一,非常经典的串行耦合。
他分析通道一(mhd)和通道三(hall)的耦合,
大尺度的mhd流体运动,绝不会直接驱动微观的hall效应。
hall效应只在离子惯性长度这种极小尺度上才会自发点燃。
但大尺度mhd却握有生杀大权。
它决定了hall效应有没有机会出现。
如果大尺度挤压出的电流片太短太厚,hall尺度还没来得及形成,整个系统就已经在电阻耗散中退化为了缓慢的sweet-parker重联。
这是阴险的间接的条件性耦合。
……
画到最后,只剩下通道六——多迴路集体响应。
江临拿著笔,悬在半空中很久,迟迟无法落笔。
这一个通道,他画不出任何耦合图。
他不知道庞大的人造地表网络,该怎么和等离子体內部的微观hall效应耦合?
他甚至无法在数学上证明它是否真的作为重联的一个主导变量而存在。
嘆了口气,他在图纸上把通道六这一栏画了一个巨大的虚线框,空著,留给未来那个可能更聪明的自己,或者留给下一次更危险的远征。
第二十七年的整个春天,江临都在画这些耦合图。
隨著图越画越多,他意识到,文字的描述是苍白的。
如果不能用数学语言把这些耦合关係钉死,那他所做的就只是一堆玄学。
每一组耦合关係,他都开始尝试用一组无量纲数来表徵其物理閾值。
通道一和通道二的耦合閾值,由电流片的长宽比决定。
通道一和通道三的耦合閾值,由电流片长度与离子惯性长度的比值决定。
通道二和通道三的耦合閾值,由plasmoid尺寸与离子惯性长度的比值决定。
通道四的贯穿强度,由等离子体的湍流强度决定。
通道五的封闭程度,由边界条件的开放/封闭比决定。
通道六,依然是令人挫败的未知。
每一个无量纲数,都对应著一个清晰的物理直觉。
每一个直觉,都是江临从过去十年的文献和工程经验里提炼出来的。
他深知,自己並没有发明这里面的任何一个无量纲数。
l/d_i在bhattacharjee的论文里写得明明白白,湍流马赫数在lazarian的公式里隨处可见。
它们都在浩如烟海的文献里。
但是,把它们像搭乐高积木一样,精准地咬合在一起,组合成一个完整的,能够描述整个等离子体系统多通道状態的网络,这件事,在人类已知的学术文献里,没有看到任何一篇论文做过。
最接近这个伟大构想的,是cassak在2006年发表的一篇著名综述。
在那篇综述里,cassak提出了一个令物理学界著迷的0.1重联率问题。
即为什么在太阳耀斑,地球磁尾,甚至是实验室里,那些被定义为快速重联的系统,无论底层的物理机制是hall效应还是plasmoid,其无量纲重联率数值,总是神秘地落在0.1左右这个量级?
cassak討论了这种神奇的收敛,但他那篇文章並没有给江临想要的多通道网络框架。
不过,它给了江临一个极大的启发。
不同尺度不同机制的背后,必定存在某种更高维度的,共同的速率选择约束。
就像万川归海。
江临站起身,在北墙上那张密密麻麻的研究史地图上,郑重地把cassak这篇2006年的论文位置用红笔圈了出来。
第二十七年的夏天,气温飆升,废土进入了一年中最酷热难熬的季节。
而江临,则在犹如蒸笼般的石屋里,开始尝试做一件在理论物理学界极具野心也极其危险的事。
他想用一个数字,仅仅是一个无量纲数。
去刻画整个复杂的多通道系统,从缓慢的电阻主导,瞬间转向狂暴的快速重联主导的临界几何条件。
这个数,必须是六个通道(或者是已知的五个通道)耦合关係的精妙组合。
它不应该是任何单一已有的旧无量纲数,否则就是多此一举。
在物理学史上,坐在黑板前拍脑门发明一个新的无量纲数很容易,隨便找几个参数乘乘除除就行了。
但要发明一个有深刻物理意义的,能够真正预测未知真实物理现象的新无量纲数,难如登天。
纵观歷史,绝大多数由野心勃勃的学者自创的无量纲数,最后在经过同行的无情审查后都会被发现,要么只是已有某个数的简单代数变形,毫无新意。
要么就是个数学游戏,无法在实验中对应任何可测量的物理现象。
江临提醒自己,先尝试,再审讯。
如果尝试失败,毫不犹豫地承认失败,绝不迷恋自己的造物。
如果尝试看起来成功了,立刻去找最刁钻的反例去攻击它。
他开始动手。
第一次尝试,他把lundquist数s、电流片长宽比l/a、离子惯性长度比l/d_i、湍流强度σ、边界开放比η,五个量加权相加。
每个量取对数,得到一个组合数,命名为g_1。
他打开工作站,调出已知样本库,用g_1去擬合那些经过无数次观测的经典事件.
太阳耀斑,地球磁层亚暴,托卡马克装置里的撕裂模破裂,以及大学实验室里小型重联装置的实验数据。
结果出来后,g_1在所有样本上都给出了似乎很接近的数值。
但这恰恰是最大的问题:接近,不等於物理上的一致。
太阳耀斑g_1≈2.3。
磁层亚暴g_1≈2.7。
托卡马克g_1≈1.9。
实验室装置g_1≈1.4。
如果g_1这个参数真的准確描述了系统快慢转变的本质,那么所有快速重联繫统的g_1应该都大於某个极其明確的临界值g_c,而慢速系统的g_1应该都小於g_c。
从分类上看,太阳耀斑和磁层亚暴绝对属於快速重联,托卡马克撕裂模勉强算中等,实验室装置绝对是慢速。
它们的g_1数值確实呈现出快>中>慢的趋势。
但这有什么用呢?
临界值g_c又在哪里?
如果g_c=2.0,那么实验室装置(1.4)和托卡马克(1.9)都被分到慢速,磁层亚暴(2.7)和太阳耀斑(2.3)被分到快速。
这种划分在统计学上勉强成立,但在物理学上,江临看著这堆数字,感到很不满意。
不是因为数值分布不够漂亮,而是g_1的定义本身,没有任何坚实的物理根基。
为什么是五个量加权相加,权重的分配凭什么定成这个比例,为什么是取对数而不是取指数或者其他非线性函数?
这些他都答不上来。
g_1看起来更像是一个为了凑出结果而事后诸葛亮般擬合出来的数学工具,而不是一个基於第一性原理先验定义出来的物理量。
事后擬合的工具是很危险的。
就像一个量身定製的谎言,能完美擬合你给它的任何已知旧样本,但一旦遇到从未见过的新样本,它就会立刻崩溃,无法做出任何有效的预测。
江临毫不犹豫地把写满g_1推导过程的十几页手稿全部扔到一边。
他决定摒弃这种已有无量纲数加权拼接的方法,从物理过程本身,从最底层的因果关係开始推导。
风机的每一个通道,都有一个特徵驱动时间和一个特徵响应时间。
驱动时间——指的是这个通道接收到外部施加的扰动,所需的时间尺度。
响应时间——指的是这个通道內部的机制,对这个扰动做出反馈並產生动作所需的时间。
如果驱动时间远大於响应时间(外界变化慢,自己动作快),这个通道就能游刃有余地主导局面。
如果驱动时间远小於响应时间(外界变化极快,自己动作跟不上),这个通道就会严重滯后,甚至被拋弃。
这正是风机稳定运行的判断標准。
风速变化(驱动,几秒钟)远大於偏航电机的动作时间(响应,几百毫秒),所以风机总能稳稳地迎著风。
但如果废土上颳起诡异的微秒级阵风(驱动时间极小),偏航机构根本来不及响应(响应时间太大),风机就会被动挨打,甚至断轴。
江临意识到,如果把磁重联的通道也这样拆解,每一对通道的驱动/响应时间比,都是一个有著清晰物理意义的无量纲数。
他开始计算。
如果用之前列出的六个通道两两组合比对,根据排列组合,会得到整整十五对驱动/响应比。
理论上,这十五个比值构成的高维矩阵,才是真正完整的耦合网络描述。
它不是某种粗暴的加权,而是十五个相互独立的,有血有肉的物理比值。
但问题隨之而来。
十五个数,太多了。
人脑根本无法直接处理十五维参数空间里的相变过程。
这在数学上也许很美,但在工程和预测上,毫无实用价值。
江临想,能不能在这个庞大的十五维矩阵里,做一次主成分分析?
能不能提炼出一两个起著决定性作用的关键比值,作为整个系统快慢转变的咽喉代表?
他重新审视这十五个比值的理论推导。
然后发现其中有些比值,在无论是恆星还是实验室的所有样本里,都无限接近於1。
这意味著这两种机制总是同步的,它们之间的耦合併不构成阻碍系统演化的瓶颈,可以忽略。
有些比值在所有样本里都极大或极小,这意味著某个通道始终处於绝对压制或绝对休眠状態,远离快慢相变的临界点,也可以忽略。
经过近一个月的残酷筛查和计算剥离,只有少数几个比值,在不同的样本环境之间,表现出了剧烈而显著的数量级差异。
也就是说,只有这几个比值,才是真正卡住系统脖子,决定其是在缓慢耗散中死亡,还是在狂暴重联中爆发的关键枢纽。
他小心翼翼地把这三个关键比值挑了出来。
第一个:电流片长度l与plasmoid次级不稳定性临界激髮长度l_c的比值Π_1 = l/l_c。
它决定了系统会不会碎裂出磁岛。
第二个:电流片长度l与离子惯性长度d_i的比值Π_2 = l/d_i。
它决定了系统能不能触碰到无碰撞的霍尔物理.
第三个:等离子体的湍流强度σ与某个临界湍流强度σc的比值Π_3 = σ/σ_c。
它决定了背景噪声是推波助澜,还是喧宾夺主。
三个比值,每一项都有坚如磐石的物理意义,每一项都对应著一个关键通道的驱动/响应匹配生存条件。
江临开始把这三个比值,用一种高度凝练的代数拓扑形式组合起来。
gate index,g。
门控指標。
为了处理湍流那种非单调影响,他引入了一个归一化响应函数f(σ/σc)。
这个函数在湍流强度適中时取峰值,在过弱和过强两端同时下降。
g的v0.1形式,被他暂时写成——
【plasmoid项+hall可达项+湍流响应项】
具体代数形式暂时不进入正文,只进入工作日誌。
第一个log在物理上相当於一个数学开关,对应著plasmoid通道是否被强行开启。
第二个log对应著微观的hall通道在当前尺度下是否可达。
第三个带负號的log,则精確地描述了湍流背景是否已经强到了足以压制其他所有单一重联机制的地步。
公式中出现的负號,以及括號里那个看似奇怪的分式结构x/(1-x)並非江临凭空捏造。
那是因为湍流的影响在物理上是非单调的,適度的湍流会增加耗散,帮助快速重联发生。
这种精妙的非单调代数结构,江临是在回忆lazarian-vishniac在1999年发表的那篇经典论文某个不起眼的附录里,见过类似的推导形式。
他从不认为自己是凭空创造的神,自己只是在借用巨人的砖块,並修改了水泥的配方。
g参数的定义终於出炉。
大为兴奋的江临强迫自己冷静下来,进入残酷的交叉检验阶段。
他用新的g参数,重新去擬合那些已知样本。
太阳耀斑:g≈4.1。
磁层亚暴:g≈3.8。
托卡马克撕裂模:g≈1.2。
锯齿崩塌:g≈1.5。
实验室小型重联装置:g≈-0.3。
江临看著这些数字,心跳开始加速。
如果设临界值g_c=2.0。
g > 2.0的,毫无爭议地全都是宇宙级別的快速重联(太阳耀斑、磁层亚暴,数值高达4左右)。
g落入[1,2]区间的,全都是处於临界状態、具有阵发性特徵的事件(托卡马克的破裂现象)。
g < 1甚至是负数的,则是那些死气沉沉的慢速事件(实验室装置)。
这个分类结果,比之前那个胡乱拼凑的g_1清晰锐利了无数倍。
更重要的是,g的定义有深厚的物理根基。
那三个log项,每一个都对应著具体的通道耦合物理方程的渐近展开,绝不是事后为了好看而凑出来的常数。
看著这组漂亮的擬合数据,江临並没有兴奋地跳起来欢呼。
因为他知道这仅仅是第一步,他还没有去主动寻找反例来攻击它。
如果g这个门控指標真的描述了宇宙中等离子体快慢转变的真理,那它必须且一定要满足以下三个苛刻的条件,缺一不可。
第一,必须能准確预测新样本。
把g应用到他之前定义公式时完全没有见过的新实验或观测数据上,它必须仍然能给出正確的物理分类。
不能见光死。
第二,必须对参数变化具备物理敏感性。
如果某个系统的外部参数被人为地从慢速区间调整到了快速区间,g的数值必须在数学上相应地从小变大,且变化趋势必须符合微积分的单调性预期。
第三,必须与已有的经典观测现象兼容。
比如cassak提出的那个著名的0.1重联率现象,应该在g≈g_c(临界閾值)附近,作为一种数学推论自然而然地涌现出来。
接下来的第二十七年,整整一个夏天和一个秋天,江临像一个不知疲倦的机器,全部耗在了这三项残酷的检验上。
第一项检验:盲测新样本。
他极其严谨地从资料库里,剔除了所有之前用来构建g参数的训练数据。
然后从c类索引库里,隨机挑出了十几篇他之前几乎没细看过的冷门观测论文。
这里面包括mms卫星在地球磁层顶探测到的几个异常高解析度微观重联事件,几篇统计太阳耀斑不同爆发功率层级的论文,以及几篇关於受控核聚变装置在特殊磁场位形下发生重联的偏门统计。
他把这些全新样本的物理参数,一个个硬生生地代入g参数公式中进行盲算。
mms观测的某特定重联事件:g≈3.5。
根据论文描述,该事件爆发极快,符合g>2.0的快速分类。
某次超大型 x 级太阳耀斑:g≈4.4。
显然极快,符合。
某次规模较小的m级耀斑:g≈2.8。
依然快速,符合。
某种被故意压制不稳定性的特殊位形托卡马克实验:g ≈ 0.7。
论文描述其演化极慢,符合g<1.0的慢速分类。
某实验室提出的一种极具爭议的新型重联宣称实验:g ≈ 2.1。
这恰好卡在临界值边缘,这解释了为什么该实验在学术界引发了到底是快还是慢的巨大爭议。
十几个完全盲测的新样本里,绝大多数精准命中了g的预测区间。
只有一两个边界极端情况因为测量误差需要更精细的重新判定。
第二项检验:参数空间的单调性扫描。
江临利用石屋工作站那少得可怜的算力,写了一个脚本,做了一系列暴力的参数空间扫描。
在代码里,他固定住其他所有参数,然后逐个连续地改变lundquist数,电流片长度,和湍流强度,观察g值隨之变化的三维曲面趋势。
结果显示:
g隨lundquist数s的增加呈现单调递增。
s越大,流体越理想,发生次级不稳定性导致快速重联的概率越高,g越大,物理直觉吻合。
g隨电流片宏观长度l的变化同样是单调递增的,符合理论预期。
g隨湍流强度σ的变化是非单调的,在曲面上呈现出了那个他刻意设计的非单调拋物线穹顶结构。
σ处於適中值时g达到最大,辅助重联;σ极度过大或过小时g迅速减小。
契合复杂等离子体物理事实。
第三项检验:歷史悬案的兼容性。
他把资料库里所有计算出g≈g_c(在临界值2.0附近游荡)的实验和观测样本全部单独提取出来,去查看这些论文里实际报告的无量纲重联率数值。
奇蹟出现了。
这些处在临界状態的样本,其报告的重联率数值,密密麻麻地集中在了0.05到0.15之间这个狭窄的带状区域里。
兼容了cassak综述中那个困扰了学界多年的0.1 重联率谜团。
而反观那些g远大於g_c的样本,比如耀斑,重联率如同脱韁野马显著高於0.15。
g远小於g_c的样本,重联率则像死水一样显著低於0.05。
三项检验都没有给出致命反例。
第二天清晨,阳光透过窗欞打在屏幕上。
一夜未眠的江临搓了搓僵硬的脸颊,在数字日誌上以克制到甚至带有警告意味的笔触写下了一段话。
【第二十七年九月十七日】
【g参数(门控指標)定义出炉,三项苛刻的物理交叉检验全部通过。】
【但必须牢记:通过检验,绝不等於这就是真理。】
【目前的g,极有可能只是在我目前已知、能够获取的有限样本范围內成立的一种高级现象学规律。】
【一旦进入黑洞吸积盘或者中子星表面那种更极端的相对论参数空间下,g很有可能瞬间失效。】
【在面对不同的、尤其是开放性的宇宙边界条件下,g的方程形式可能需要极其痛苦的拓扑修正。】
【现在的g,只能被称为v1版。】
【它仅仅是人类在理解复杂系统相变道路上的一个落脚点起点,绝不是最终的答案。】
写完,保存日誌,他走到东墙边,在那份改变了他认知轨道的【recon-vc-01】(变电站剩磁观测)下面,新写了一行字。
【theory-g-01:门控参数(g-index)定义v1,三项交叉检验通过。】
【適用范围:具有多通道耦合特徵的常规磁重联等离子体系统。】
【未解决的阿喀琉斯之踵 1:通道六(宏观多迴路网络集体响应)的数学形式过於复杂,尚未能成功纳入g的对数框架。】
【未解决的阿喀琉斯之踵 2:临界閾值g_c仅仅是基於统计擬合得出,极度缺乏第一性原理的严密数学物理推导。为什么是 2.0,而不是π或者e?】
【未解决的阿喀琉斯之踵 3:g在趋於无穷大或无穷小的极端参数渐近行为未知。】
四行,第一行是耀眼的胜利结果,而后面三行,全是对g参数本身挖出的新坑。
这大概就是科研的常態,任何一个真正有价值的新理论工具被发明出来,在解决一个旧问题的同时,往往会像潘多拉魔盒一样,一口气打开更多新问题。
不过在江临看来,一个好的科学工具,本来就应该在斩断旧有迷局的同时,劈开新的战场。
如果一个物理理论只负责完美解决当下的问题,却不再衍生出任何新的未知,那它要么是存在隱蔽的逻辑谬误,要么就是一门行將就木的死学问,已经被时代拋弃。
只有那种能不断催生,且同时打开崭新物理图景的工具,才是真正能够推动人类认知边界向前碾压的利器。
g参数也许並不完美,甚至可以说是千疮百孔。
也许在百年之后,它会被后人发明的一个更加优雅更加高维的新张量工具无情地扫进歷史的垃圾堆。
但至少在这一刻,在这个荒凉的废土前哨站里,它悍然打开了三个全新的物理学方向。
单凭这一点,它就值得继续追。
第二十七年十月,秋风渐起。
江临进入了漫长而琐碎的整理期。
他把过去整整一年里,在无数个不眠之夜画过的所有通道耦合图,所有那些被他证明是错误从而揉成一团的g_1手稿,所有针对g参数的每一项交叉检验的底层数据表格,以及那些针对反例绞尽脑汁的审讯记录日誌,全部按时间顺序进行了数位化归档。
他建立了一个庞大的主档案文件夹:【theory-g-01】
子文件夹结构——
【01_六通道现象分解】
【02_耦合响应矩阵推演】
【03_失败尝试_g_1的废墟】
【04_g参数v1版代数定义】
【05_三项交叉盲测检验数据集】
【06_异常反例审讯与辩护记录】
【07_未解决致命问题清单】
每一个子文件夹里,都附带了详尽的推导过程和思考中间態记录。
整理完成后,他在工作站桌面上把【theory-g-01】这个文件夹放在了桌面正中央这个最显眼的位置。
在这个中央文件夹的下方,犹如眾星拱月般排列著六个旧有流派的理论文件夹。
【mr-sp-01】、【mr-pk-02】、【mr-num-03】、【mr-plm-04】、【mr-edr-05】、【mr-turb-06】。
而在屏幕的右上角,孤高地悬掛著那个一切疯狂想法的物理起源——【recon-vc-01】(变电站物证数据)。
新的【theory-g-01】端坐在正中央。
这不仅仅是图標的摆放,更是江临大脑中物理学认知架构的外化投影。
中央是诞生於废土的全新理论尝试,下方是它所依託吸收並试图统合的所有人类科学遗產。
右上角则是无言的宇宙大自然,用灾难废墟给予他的真实启示。
江临双手交叉垫在脑后,静静地看著屏幕上这个充满內在逻辑张力的布局结构。
他並没有狂妄地觉得自己已经解决困扰人类半个多世纪的磁重联世纪难题。
但他有一种真切的踏实感,他觉得自己第一次在这片研究地图上,留下了一道属於自己的铅笔痕。
时光荏苒。
第二十八年到第三十二年。
整整五年的时间在废土的风沙中悄然流逝。
在这五年里,江临展现出了令人髮指的科研定力。他没有急功近利地立刻去死磕g参数留下的那三个极度诱人的未解决数学难题。
相反,他选择沉下心来,做了一件在聪明人看来极度笨拙,但在严谨科学体系构建中却最为基础的苦力活。
他要把g门控指標,强制应用到他硬碟文献库里能搜刮到的,人类歷史上发表过的所有磁重联观测和实验文献样本上,进行最疯狂的普適性压力测试。
不是十几个。
也不是几十个。
是全部,不留死角的穷举。
a类核心论文八十二篇,b类重点论文四百二十一篇,再加上c类海量索引库里,只要能够通过字里行间拼凑追溯出物理参数的全部案例记录。
最初粗略筛查下来,他的统计表格里赫然堆积了一千八百多个真实的物理发生场景。
太阳色球层爆发的各种级別耀斑。
地球磁尾在太阳风挤压下引发的复杂磁层亚暴。
mms编队卫星在地球磁层顶探测到的几十个微小的电子扩散区重联事件。
世界各地数十台托卡马克装置里爆发的撕裂模。
聚变等离子体核心发生的锯齿崩塌现象。
各大高校实验室里,用脉衝放电產生的二维重联装置实验报告。
超算中心耗费数百万机时跑出来的plasmoid高解析度数值模擬网格数据。
人为施加外部扰动的湍流背景下的重联物理实验。
所有这些来自天南海北,跨越天文尺度到实验台尺度的现象,全部被江临无情地丟进了一个由python驱动的巨大索引资料库表里。
面对这海量的数据,江临並没有急切地让电脑直接去画一张大一统的散点图。
那样画出来的图是用垃圾数据污染好数据。
他所做的第一步,是数据分级清洗。
他在那个巨大的索引表的最后一列,加上了一个决定样本生死存亡的栏位名称。
【data_quality(数据置信度等级)】
a级样本,文献必须同时清晰地给出了明確的lundquist数推算依据,无可爭议的电流片长度,基於实测的离子惯性长度,有谱分析支撑的湍流强度估计,毫不含糊的系统边界条件说明,以及最关键的,该论文对重联率的定义必须在物理公式上可严密追溯,不能只是正文里缺乏定义的经验数值。
这种条件严苛到变態,能符合的样本犹如凤毛麟角。
一千八百多个场景里,在经过江临如同拿放大镜般的逐篇审查后,最后竟然只剩下了可怜的三百多个。
b级样本,允许论文在某些次要参数上存在缺失或含糊,但缺失的参数必须能够通过江临利用等离子体基本流体力学方程,从论文给出的其他物理量里反向推导出来,並且这种反推的误差限带,必须能在表格里写得清清楚楚。
这类样本是学术界的主力,数量最多,经过整理,大约有七百多个。
c级样本,数据严重缺失,只能做数量级上的粗略估算。
比如很多观测论文只豪言壮语地给出了观测到的重联率,却根本没有给出电流片的宏观长度。
或者有些实验报告只给了一张彩色几何示意图,却没有提供任何可供同行覆核归一化的特徵尺度標尺。
这些样本被江临保留在索引库里作为参考,但被禁止进入核心参数的擬合计算。
d级样本,直接剔除。
原因包括:有的是重联率定义与其他通用论文完全不在一个维度上根本不具备可比性,有的是实验边界条件交代得一塌糊涂,有的是关键参数无法恢復。
还有一些早期数值论文受时代限制,人工耗散和网格处理交代得不够细,无法进入同一套统计口径。
第三十年的严冬,鹅毛大雪封堵了石屋的窗户。
在温暖的炉火旁,江临终於开始了数据可视化的时刻。
他决定第一次,仅仅只使用a级样本来绘製图表。
代码运行,图像生成。
x轴:他定义的骄傲——g指標数值。
y轴:那些顶尖论文中报告的实测无量纲重联率。
当蓝色的散点图在黑色背景的屏幕上跳跃出来时,江临的手指悬在滑鼠上,盯著屏幕看了很久很久。
出现在屏幕上的,並不是他在梦中期待过的那种,能够作为绝对真理去威慑整个学术界的,如同教科书般完美的狭窄曲线。
更不是一条能用简单线性方程描述的漂亮直线。
那是一片云。
一片非常厚,分布散乱,充满著现实世界观测误差和不同系统稟赋带来巨大方差的散点星云。
但江临的眼睛,透过这片看似混沌的云层,敏锐地捕捉到了其深处的灵魂。
这片云,有极其明確的宏观物理方向。
g大的样本,散点云整体剧烈抬升,重联率的数值绝对偏高,它们在狂野地燃烧。
g小的样本,散点云像一滩死水般紧贴著底部,重联率整体低迷,它们在缓慢地耗散。
g接近临界区间的样本,它们的重联率如同受到了某种神秘力量的约束,密密麻麻地挤压在y轴0.05到0.15这个狭窄的水平带状通道之间。
这绝对不是一个严丝合缝的数学证明。
甚至算不上是具有排他性的强证据,方差太大了。
但这確实是在三百多个最纯净的a级人类物理样本库里,呈现出的一个统计学相变趋势。
江临深吸了一口气,没有急著按保存键留下这张最终图。
他还要继续施压。
切出代码,换上包含更多方差的b级样本重新运行。
由於b级样本数量庞大且带有反推误差,图表上的云层瞬间变得更厚,离散度急剧扩大,有些不守规矩的异常点甚至飘得极远,像是在嘲笑理论的脆弱。
但是,如果你眯起眼睛,退后两步看,那条快中慢的整体斜向物理演化方向,宛如一条脊椎骨,依然撑在图表的中央没有断裂。
紧接著,他生猛地將a级和b级样本重叠合併。
一千多个数据点砸在图上。
趋势,依然坚如磐石。
最后,他深吸一口气,把那些仅供参考的c级样本,用黯淡的灰色,如叠印在最下层。
整张宏大的图表瞬间变得更加脏了。
密密麻麻的数据点交织在一起,就像是深空中一片被充满星际尘埃的混沌星云,让人看一眼都觉得眼晕。
可即便在如此极端的暴力干扰下,那片云的整体轮廓仍然没有被完全打散。
依然固执地大致沿著同一个从小到大,从平缓到跃升的方向,横亘在坐標系中。
时间来到第三十二年的初春。
废土迎来了新一轮的復甦。
江临在工作站前,为这张耗费了五年时间,匯聚了人类半个世纪智慧与谬误的终极星云图表,敲下了最终的定名標题。
【g-index v0.1:trend only(门控指標 v0.1:仅展示全局统计趋势)】
隨后,他又在判词下方,以冷静客观口吻,写下了四条针对图表巨大方差的自我解释。
第一,这种顽强趋势的出现,极有可能意味著g指標在底层的代数拓扑上,確实幸运地捕捉到了等离子体宇宙中跨尺度的多通道相变门控效应。
第二,图表中巨大的数据离散度,不可避免地来源於b级样本复杂的参数反推理论误差。
第三,离散度同样也可能来自不同年代不同作者在撰写论文时,对无量纲重联率在流体力学定义上的內在使用不一致。
第四,最深层的可能,这种离散说明目前的g指標v0.1版存在严重的理论缺陷,它遗漏了某个足以改变局部权重的关键物理通道,尤其是那个尚未找到数学途径纳入公式的【通道六:多迴路集体响应网络】。
这四条解释,让它不至於被巨大的数据误差直接证偽。
但这四条解释,同样也在宣布,它现在的形態,充其量只是一张大航海时代早期,航海家凭著经验和残破的指南针画出来的远洋地图。
模糊不清,做工粗糙,充满了令人心惊肉跳的统计学噪声,甚至有些岛屿的形状都被画歪了。
但幸运的是,它指示的大方向,那条通往物理学新大陆的航线,並没有完全错。
第三十二年,九月十七日下午。
他定义出g参数整整五年后的纪念日。
远处,那几个承载著废土环境数据的观测点a和b的指示灯,依然在顽强地记录著时间的流逝。
而在石屋內部,工作站的高清屏幕上,那张匯集了一千八百个样本,由蓝色和灰色斑点组成的混沌散点图,正安静地停留在那里。
江临靠在有些掉皮的电脑椅上,目光深邃地看著屏幕上那片被框在x轴和y轴之间的星云。
在长久的凝视中,他仿佛突然顿悟了某种更深层次的科研宿命。
在接下来漫长的岁月里,他不能再像这五年一样,继续去扩大索引库的样本量了。
样本数量哪怕再增加一万个,十万个,把整个图表填得密不透风,也永远无法发生从量变到质变的飞跃,无法把一片建立在统计学概率上的混沌云,萃取升华成一条可以用严格等式连接的纯粹数学定理。
基於宏观数据的经验趋势擬合,到这一步,已经走到了人类认知的尽头。
如果他內心深处真的渴望知道这个亲手捏造出来的g指標里,到底有没有触碰到大自然不可撼动的物理法则,他就必须对它动刀子。
把它从宏大敘事里剥离出来,把它拆到一个没有任何实验噪声干扰的几何模型里去。
目前的g指標妄图同时包含六个狂暴的物理通道,適用范围太宏大,数据来源太脏,著太多前哨站修风机式的野路子工程气息。
而在物理学那座名为真理的圣殿里,真正能够被毫无感情的数学逻辑死死咬住,能经得起几百年推敲的,从来不是什么宏大的综合模型。
江临的目光渐渐锁定了那六个通道中的一个。
plasmoid演化通道。
二维不可压缩电阻mhd。
有限长宽比的长薄电流片。
能量方法。
sobolev空间里的先验估计。
严酷的能量泛函变分方法。
以及那些隱藏在高维拓扑中的sobolev空间解析先验估计。
只有用这些纯粹的数学武器,在这个最狭窄的点上凿穿它,才能赋予g参数真正的灵魂。
江临移动滑鼠,毫不留恋地关掉了那张耗费他五年心血的散点星云图。
打开研究日誌。
【第三十二年归档总结】
【g-index v0.1在全量级分层样本资料库中,確实显示出不可忽视的强统计趋势。】
【阶段性结论:其底层现象学有效性,目前未被证偽。】
【限制警告:这种高可信样本中存在的趋势,仅仅是统计学上的妥协,绝非第一性原理的严密理论证明。它目前仍是一个非普適性的粗糙工具,严禁將其盲目外推到未知的,尤其是具有相对论效应的极端参数空间。】
【悬而未决的已知地雷:通道六的拓扑响应未被纳入,临界閾值缺乏解析推导,极端参数区间行为未知,处理非单调性的湍流项代数形式並不具备数学上的唯一性。】
【战略下一步动作,也是最终动作,放弃全局擬合。从庞杂的g参数体系中,只抽出那条最狭窄的门缝(电阻mhd条件下的理想电流片撕裂相变)。】
【我將不再用统计数据去支持它。】
【我將用纯数学,去证明它。】
本站所有小说均来源于会员自主上传,如侵犯你的权益请联系我们,我们会尽快删除。