第269章 老师,它很物理!(三合一)
玄浑道章小说推荐阅读:Apop之我在首尔当外教、赘婿出山、泥泞、股神传奇、乡村中医李茂阳、开局账号被盗,反手充值一百万、极品全能学生、推拿、最强狂兵、带着超市重返年代、都市:我的器官能自己强化、不装了,我有五个天仙姐姐!、国宝神鉴
组装体了,铖那一边的Tb(pcam-Br)s配合物,苏老师已经在排合成线了。
Au(111)单晶基底已经定了俄铅冶那家。
超高真空腔区已经开了。
锁相放大器和X射线斩波系统,在苏老师他们手底下走“接口对线”。
光束线那一头,等陆明远教授从京城回来再说。
所以现在卡在中间这一摊的,全是软的……
反演算法、配位场张量、谱拟合、相位标定。
这一摊要不动,前面那一摊样品做出来了也没意义。
後面陆明远的光束线调的再好,也只能拿到一组没人能看懂的隧穿谱。
“所以你的活儿。”
“今天上午,去跑一下数学那边。”
李东微微一愣。
“现在吗?”
“嗯。”
张燕点了点头。
“咱们组里数学这一块,本来就是让你来挑大梁的。”
“现在一开工,你的压力可不小哟。”
“合城那边,列旺教授现在做的就是数学反演这一块。”
“咱们这边再不下场,等他们做出来,咱们就只剩下“参考别人结论’的份儿了。”
张燕继续说。
“所以今天上午,你别管别的。”
“先去隔壁那间小办公室。”
“那间屋子隔音比较好,外头吵不到。”
“你先把那两个数学瓶颈,自己再过一遍。”
“先按现在公认的“反演’这条路推。”
“别管你心里那个“方向不对’。”
张燕看了他一眼,眼神很认真。
“我知道你昨天心里有想法。”
“但你想法没成型之前,咱们这边的进度不能停。”
“哪怕你的想法是对的,你也得先证实“反演’这条路真的不通……”
“否则你自己都说服不了你自己。”
“明白吗?”
李东看了她一下。
心里那一句“原来你不光直觉好,思路也清”差点从嘴里溜出来。
他点了点头。
“行。”
“我去过一遍。”
“争取就是尽快跑出来,不卡在前面那一两步。”
张燕点了点头。
“嗯,前面那一两步,你能跑多远跑多远。”
“现在国内国外这麽多组,全都卡在那。”
“你能往前挪半步,咱们组就比别人快半步。”
李东“嗯”了一声,就朝着办公室走去。
来到房间後,李东也有自己的计划。
他要做的可不只是跑一下数学那麽简单。
他还要顺手验一件事。
他清晨想的那一套物理切口,必须要有数据支撑。
不是说他物理那一头脑补得再漂亮,就一定是对的。
他得拿数学这边能跑出多远来做参照。
万一数学方向真的只是小卡一下,後面一路畅通呢?
那他这一晚上的物理切口,就成了无用功了。
他不能拿一个直觉去推翻整个组的工程进度。
所以他得先去把数学这边推到尽头。
李东把包放在桌上,也不开电脑。
他先把草稿纸和笔拿了出来。
李东是搞纯数的。
反问题,严格的说他不算专业。
可是这些天,鄂伟南院士给的那一份书单,他已经啃下了一大半。
Hadamard的良态性条件、Tikhonov正则化、谱截断、Backus-Gilbert方法…这些名字加上对应的核心定理,全都紮紮实实地挂在他记忆宫殿的【应用数学/反演】那一格。再加上他理解东西的速度,本来就在常人之上。
因此他也有信心不比一般走应用数学的学者差。
他坐下来,开始推导。
第一步:测量算子K:T→S的形式定义,T是配位场张量空间,S是隧穿谱空间。
第二步:核函数展开,把K拆成本征函数序列。
第三步:上吉洪诺夫正则化,调整参数α。
第四步:把K的奇异值谱分布,对应到T的低秩近似。
第五步……
他一步一步地推。
每一步,他都对照着昨晚已经搭好的那一套物理参数,看数学这一头能不能合上。
基本上前面五六步,都没有出大问题。
Tikhonov那一段他卡了不到十分锺,就用一个加权范数把它处理掉了。
基函数的选择那一段,他换了一组Whittaker型的局部基底,规避了一部分相位混淆。直到第七步。
李东看着自己写出来那一行,皱起了眉头。
【在(L2(R,du),〈·,·〉_K)上对相位算子做谱分解,要求三组基底的过渡矩阵收敛於一个无相位歧义的西变换……】
“是这样吗?”
“不对呀。”
这一段,本来是要把不同基底的展开“硬合”在一起。
按列正常思路,是用一个隐性约束量去反卡这个过渡矩阵的上界。
可是李东推到这里,就是过不去。
他换了三种正则化,结果一种比一种糟糕。
他试着把鄂伟南院士点出来的那一根“隐性约束”硬塞进去。
卡得住,又卡不死。
第三峰从糊里头能拎出来一点点边,但拎到一半,又掉回去了。
李东盯着那一行字,叹了一口气。
“……果然啊。”
“全国十几个组,全都卡在这。”
就在他第三次放下笔的时候,背後突然冒出一个声音。
“卡住了?”
李东猛地回过头。
吴开教授不知什麽时候已经站在了他的身後。
李东赶紧站起来。
“吴老师。”
“您来了。”
吴开摆摆手,走过来在他对面那张椅子上坐下。
“我刚从化院开完会过来。”
他看了一眼李东桌上的草稿纸。
看见最上头那一张推到一半的过渡矩阵谱分解。
吴开点了点头。
这一块,他也是看得懂的。
不是因为他数学多好。
是因为这一块,全圈儿的人都卡在这儿。
卡得太多了,他不熟也熟了。
吴开笑了笑。
“李东,你别急啊。”
“这个东西,已经卡好几年了。”
“国内国外大家都卡在你这一步呢。”
“你今天上午就推到这儿,已经比绝大多数人快了。”
“再说,你也别给自己压力。”
他温和地看着李东。
“我也没想过你一来就能解决问题啥的。”
“咱们组的进度计划,是按“数学这一块磨一两年’排的。”
“前面那一堆样品和针尖,咱们可以先把後续的衔接做了。”
“等你这边反推回来,咱们再把它接上,也是可以的。”
吴开这话说得是真心的。
他怕李东受挫。
李东毕竞还年轻。
基础数学那一头,他跑得是真的太快了。
这种年少成名的孩子,最怕的就是在另一个领域里被一面墙撞得头破血流。
吴开过来的不光是问候。
他是来给李东托个底的。
李东点了点头。
他听出了吴开的好意。
心里挺感激。
他正要开口说“老师我会好好磨的”。
结果话还没出口,他想起来一件事。
他刚才推到第七步卡住的时候,没再继续往第七步那一堵墙上撞………
他用了另一种“绕”的方法,从纯数那一边借了一招过来,把第七步先放到一边,跳到了第八步。第八步往下,他一口气推了好几页。
也不知道这样行不行,毕竟大家都没用纯数的方法来,万一在这里不行呢?
所以他问道。
“老师,这後头那一段我用的自己的方法,您给看一下这样行不行。”
李东说完,就低头去翻草稿纸最下面那几张。
然而听到这话的吴开整个人就在旁边僵住了。
啥意思啊?你推过去了?
吴开数学不算精通。
他这辈子混的是化学合成。
但数学这一头,因为这个项目卡了好几年,他这一两年也开始“陪着读”。
反演这里头,每一步在哪儿、为什麽会卡、卡在哪个具体的点上…
他门儿清。
刚才李东桌面上那一张过渡矩阵谱分解。
他一进门就看出来卡在哪了。
可是李东刚才嘴里冒出来的那一段……
“後头那一段”?
哪一段“後头”?
这一步根本就没“後头”啊。
所有人都是卡死在这一步。
吴开开口道。
“等一下。”
“李东,你刚才说的,是什麽意思?”
李东抬起头。
“啊?”
“哦。”
他指了指自己手里翻到一半的那几张纸。
“就是这一步过去以後,会遇到什麽样的一个问题。”
“我本来是想,您给看一下我这边推得对不对。”
李东说得理所当然。
吴开整个人愣了一下。
“这一步过去以後?”
他有点不太确定。
李东没听出他语气里的反应。
他的注意力还在自己那几页推导上。
“嗯。”
“您看这一段……这里我用的是局部Whittaker模型那一边的过渡函子,把三组基底的相位差先冻在一个固定的代数结构里……”
吴开听到“局部Whittaker模型”五个字,他整个人都是懵的。
这是什麽??
他从应用数学那一头听过基函数、算符、反演,可从来没听过局部Whittaker模型啊。这不是反演里的语言。
这是……
朗兰兹纲领的语言。
吴开的手不自觉地往李东那堆草稿纸抬去。
“你这里………”
“不止一张纸啊?”
他这一抬手,才发现李东桌面上那一遝草稿纸,最上头那一张底下,还有厚厚一摞。
吴开下意识地把最上面那一张过渡矩阵谱分解抽走,露出底下第二张。
他低头一看。
第二张上头,密密麻麻的全是符号。
最顶上那一行,他认得……
【设过渡矩阵U_ij之间的相位歧义为中∈R/2mZ,由偶极宇称选择定则有sgn(M_dipole) ∈ {+1,-1}固定】
吴开看到“偶极宇称选择定则”这八个字的时候,眉毛挑了一下。
这一路是物理那一边的玩法。
可下头紧接着又是一行……
【在GLs的局部基变换函子下,将三组基底的相位差归约为一个有限维代数表示的同构类,由L因子的局部根数决定】
这一行下面,李东用了一个李东自己常用的对关联函数符号F_r(α),把中的可能取值死死卡在一组离散的值上。
最後一步………
他把F_n(α)的零点和宇称选择定则给的sgn合起来,得到一个干干净净的结论:
【中=0】
吴开盯着那一行【o=0】。
看了三秒。
看了五秒。
他自己都没意识到,他端保温杯的那只手,已经把杯子里的水洒到了袖口上。
第三张。
第四张。
每一张推到後面,李东开始展开第三配位壳层那一套。
奇怪的是,李东这边推第三壳层,没有再用任何一种正则化。
他用的是另一组语言。
吴开看不太懂………
但他认得最关键的那几个名词。
“WKB”。
“势垒衰减”。
“伽莫夫因子”。
吴开“嗖”地一下抬起头,盯着李东。
他的喉结动了一下,吸了一口气。
“李东。”
“你……”
吴开下意识想问……
你这是把後面的反演也做了吗?
可这一句又咽回去了,因为他感觉这不是反演。
所以他换了一句最接近他直觉的话。
“你把这一步……推过去了?”
李东抬起头。
他眨了眨眼,明显没明白吴开问的是哪一步。
“什麽推过去了?”
“没有啊。”
李东很认真地摇头。
吴开:……
吴开把手里那一张草稿纸往桌上一摁,指着第二张【@=0】那一行。
“那这是啥?”
李东低头瞄了一眼。
“哦,这个啊。”
“我刚才在第七步这儿,把基函数相位歧义拆掉了。”
“再後头第八步、第九步,是接着第三配位壳层那一段往下推的。”
“再後面,就过不下去了。”
吴开盯着李东那一张坦然的脸,足足愣了三秒。
他终於反应过来了。
李东刚才嘴里说的那一句“过不下去”
不是他、不、不是国内外都卡死的那一堵墙!
李东说的“过不下去”
是这一整摊反演问题作为一个整体,过不下去!
而眼前这一面所有人头破血流的墙,他……
吴开干脆把【中=0】那一行拍出来,直接指着最上面那一张推导。
“这一步……”
“你看一下,就这一步啊!”
“你推过去了吗?”
李东低头看了一眼。
他抬起头,对着吴开说了一句。
“嗯。”
“这一步不是很难。”
吴开:???
吴开心里立刻就开始吐槽了。
不是很难?
不是很难,全国十几个组会卡在这一步整整三年?
李东压根没意识到吴开心里那一阵翻江倒海。
他指了指自己那一行【o=0】下面的推导。
“这里我是用Ls的局部基变换函子做的。”
“就是把三组基底,先看成同一个全局自守表示在不同局部位置上的Whittaker模型。”“再用L因子的局部根数把过渡矩阵的相位差锁死在有限维代数表示里。”
“最後用偶极宇称选择定则给的实符号收尾。”
“基本上不算太绕。”
李东说得很自然。
吴开听得整个人都僵了。
头一句“L.”,他还能跟。
第四个字开始,“局部基变换函子”……他直接断片了。
纯数那一头,朗兰兹纲领、自守表示、Whittaker模型………
这些东西离他十万八千里。
吴开无语了好半天。
最後他憋出一句。
“李东。”
“你说人话。”
李东:……
李东这才算是反应过来。
原来这一面墙,全国十几个组卡在这儿。
不是卡在“反演”这一整块出不来。
而是卡在“反演”里头具体的一步。
就是基函数相位差那一步。
他刚才一直以为,大家是卡在反演这一摊的总瓶颈上。
所以吴开教授刚才说“卡好几年了”,他下意识地是把“几年”往整个反演那个大问题上想。结果不是。
大家的几年,是卡在他半小时绕过去的那一步上。
李东也不能笑话人家。
他自己很清楚。
他能半小时绕过去,是因为他手里的工具,本来就不是从反演里借的。
他手里的工具,是朗兰兹纲领。
对一个搞应数的学者来说,他根本不会想到要从GL.的局部基变换函子那一头,往一个谱反演问题里头钻。
这是两套语言。
这是两个圈子。
如果不是李东这个本来就是给朗兰兹纲领“夯地基”的人,又恰好进了这个项目组。
这两条路永远都不会装上。
因为有一个搞纯数的学者会跑来这里高应用数学。
所以这一步在外人看来是天堑。
对他来说,是顺手。
李东想通了这一点。
他抬起头,看着吴开。
深吸了一口气。
“吴老师。”
他认认真真地说。
“这一步,其实不是真正的问题。”
吴开还在震惊李东的进度,下意识的问道。
“你说什麽?”
李东语气平静。
“咱们项目,乃至全世界这十几个组,真正卡的·……”
“不是这一步。”
“用数学反演这条路,本身就走不通。”
屋子里静了下来。
换做是十分锺以前,李东在他面前讲这一句话一
吴开顶多礼貌地笑一笑,表示尊重,然後心里默默给这小子记一笔“还得磨”。
毕竟,那是一句没有任何依据的话。
可现在不一样了。
现在李东已经把全球十几个课题组卡了三年的那一步,在一个上午里推了过去。
不光推过去,还顺手在後头推了好几页。
李东这一句“反演这条路本身就走不通”。
就不是一个少年人凭直觉甩出来的话了。
这是一个已经站在墙顶上的人,回过头来,对底下还在撞墙的人说的话。
吴开让自己先冷静下来,然後才开口问道。
“那你跟我说说。”
“为什麽你觉得不行?”
李东也没绕。
他组织了一下语言说道。
“老师。”
“反演这一块,本质是一个病态算子的逆问题。”
“病在什麽地方?”
“病在它的核是不完整的,从隧穿谱反推回配位场张量的过程里,输入端任何一个噪声,都会被那个算子的小奇异值放大几个数量级。”
“我们现在的所有数学工具。”
“吉洪诺夫、谱截断、Backus-Gilbert本质都是在压那个奇异值谱。”
“压的方法不一样,但是压的逻辑是一样的……”
“我们硬塞一个先验进去,告诉算子:解应该长什麽样。”
吴开点了点头。
这一段,是他这一两年陪着读应数最熟的一段。
李东接着说。
“问题在於,所有的先验,都是数学先验。”
“什麽叫数学先验?”
“就是,我希望解是光滑的、我希望解是稀疏的、我希望解的范数最小。”
“这些希望,跟单原子上那个铁原子周围的电子云长什麽样,没有任何物理上的关系。”
“我们用一个数学的偏好,去压一个物理的解。”
“两边的逻辑根本对不上。”
“所以第三配位壳层的那一根峰,永远糊。”
吴开微微皱起了眉。
李东又说。
“还有就是………”
“反演问题本身,是把过程的结果倒过来推过程。”
“你想清楚没有,这台机器测的是什麽?”
“测的是一个X射线激发、电子跃迁、空轨道伸出、隧穿到针尖这样一个完整的物理过程产生的电流。”“它本来就是过程的结果。”
“我们却在反过来,从结果倒推过程。”
“那中间任何一段过程的物理细节,都会被压成一个有效算子,然後我们再去硬解这个算子。”“这不就等於我们把一段活生生的物理,先杀成一组冷冰的算符。”
“再去问这组算符背後那段物理长什麽样?”
“这中间天然有一道墙。”
“墙不是物理给的,是我们自己造的。”
吴开听到这儿,眉头皱得更紧了。
他没说话。
李东深吸了一口气。
“老师,您问我有没有解决办法。”
“我有一个方向。”
“我有点不确定。”
“但我是想试着说一下。”
吴开看着他。
“你说说看。”
李东把今天清晨在草稿纸上搭好的那一套切口,简单地说了出来。
“我们不要做反演。”
“我们做正向。”
“把那一组测量算子的逆问题,重新写成一个用物理参数化的正向模型。”
“参数不再是配位场张量的抽象分量。”
“参数是……”
“偶极跃迁矩阵元、未占据轨道的对称性、轨道相对於针尖的取向、以及真空段的WKB衰减常数K。”“再把X射线偏振矢量正大光明地搬进来,作为模型的输入边界条件。”
“偏振矢量在源头上选择哪条空轨道被激发……”
“这件事本来就是物理在做的,不需要让基函数去吵。”
“再有就是,第三配位壳层那一根糊掉的峰。”
“我们不去拎它。”
“我们让它自己长出来。”
“每一个配位壳层,对应在隧穿这一关,有自己一根独立的KB衰减常数K。”
“越外的壳层《越小,从针尖底下伸出去的尾巴越长。”
“换而言之,第三壳层在隧穿电流的衰减曲线上,本来就有一段独属於自己的尾巴。”
“这段尾巴不需要去拎。”
“它在距离一档一档拉远的扫描里,会自己显出来。”
“靠的就是那一根伽莫夫因子。”
吴开听到“伽莫夫因子”四个字的时候,神色微微动了一下。
李东最後说。
“老师,简单总结一下。”
“两堵墙我们不推。”
“因为两堵墙之间,本来就有一道门。”
实验室外头实验台上的“嘀嘀”声远远地传过来。
办公室里没人说话。
吴开放在桌上的那只手,慢慢的握了起来。
他第一反应是不可能。
怎麽可能?
这这个方向已经一百年麽没人提出过新的思路了。
因为走不通,所以最後大家会一起转向数学反演
这是一条所有人筛了又筛、筛剩下来的路。
怎麽可能突然冒出来一句“两堵墙中间本来就有一道门”?
怎麽可能这道门是物理?
吴开忍不住开口。
“李东。”
“我先不说你那一套数学物理上能不能合得上。”
“我说一个最朴素的事。”
“你说要正向。”
“那意味着,你要在源头上,把那个X-ERT过程从头到尾,都用物理量参数化。”
“可这里头,有一段东西”
“针尖上方那一小段真空隧穿,再加上配体场和Au(111)衬底的相互作用-……”
“它的有效势垒形状,谁都没真正测过。”
“你那个伽莫夫因庙里的K,到底丫多宁?”
“你不知道,我不知道,没人知道。”
“我们之所以走反演,就丫因为这一段我们灭不过去。”
“我们没法直接去测那一段势垒。”
“所以我们才愿意承认,反演虽然病抗……”
“但起码,每一项里都是“实验值’,再病态它也有数据撑着。”
吴开顿了一万。
“你那一头要的“乘理参数化’。”
“它从一开始,源头那一段就丫空的。”
“你拿什麽去填?”
吴开越说越觉得自己说在了点上。
这一块的核心矛盾就在这儿。
大家不丫不想做正向,丫源头那一段没法直接测。
所以才退一步,做反演。
所以他最後摇了摇头。
他的声4轻,但很笃定。
“李东,你这套思路,整体上看丫漂亮。”
“但丫从化学这一行的逻辑……”
“它太悬了。”
“我跟你直说。”
“这一套很不化学!”
屋里又安静了丐来。
李东低头看着自己那一堆草稿纸。
最上面那一张,丫他拿红笔涂的那一团。
红色那一团旁边,他写的那五个字。
“让原庙自己说”。
他抬起头来,看着吴开。
眼神平静,没有半点不服。
“老师,我懂您的意思。”
“在化学这条路上,看上去丫没有路的。”
“但………”
李东也很笃定的说道。
“老师。”
“它很乘理!”
http://www.xuanhundaozhang.com/yt127592/48690861.html
请记住本书首发域名:www.xuanhundaozhang.com。玄浑道章手机版阅读网址:www.xuanhundaozhang.com
Au(111)单晶基底已经定了俄铅冶那家。
超高真空腔区已经开了。
锁相放大器和X射线斩波系统,在苏老师他们手底下走“接口对线”。
光束线那一头,等陆明远教授从京城回来再说。
所以现在卡在中间这一摊的,全是软的……
反演算法、配位场张量、谱拟合、相位标定。
这一摊要不动,前面那一摊样品做出来了也没意义。
後面陆明远的光束线调的再好,也只能拿到一组没人能看懂的隧穿谱。
“所以你的活儿。”
“今天上午,去跑一下数学那边。”
李东微微一愣。
“现在吗?”
“嗯。”
张燕点了点头。
“咱们组里数学这一块,本来就是让你来挑大梁的。”
“现在一开工,你的压力可不小哟。”
“合城那边,列旺教授现在做的就是数学反演这一块。”
“咱们这边再不下场,等他们做出来,咱们就只剩下“参考别人结论’的份儿了。”
张燕继续说。
“所以今天上午,你别管别的。”
“先去隔壁那间小办公室。”
“那间屋子隔音比较好,外头吵不到。”
“你先把那两个数学瓶颈,自己再过一遍。”
“先按现在公认的“反演’这条路推。”
“别管你心里那个“方向不对’。”
张燕看了他一眼,眼神很认真。
“我知道你昨天心里有想法。”
“但你想法没成型之前,咱们这边的进度不能停。”
“哪怕你的想法是对的,你也得先证实“反演’这条路真的不通……”
“否则你自己都说服不了你自己。”
“明白吗?”
李东看了她一下。
心里那一句“原来你不光直觉好,思路也清”差点从嘴里溜出来。
他点了点头。
“行。”
“我去过一遍。”
“争取就是尽快跑出来,不卡在前面那一两步。”
张燕点了点头。
“嗯,前面那一两步,你能跑多远跑多远。”
“现在国内国外这麽多组,全都卡在那。”
“你能往前挪半步,咱们组就比别人快半步。”
李东“嗯”了一声,就朝着办公室走去。
来到房间後,李东也有自己的计划。
他要做的可不只是跑一下数学那麽简单。
他还要顺手验一件事。
他清晨想的那一套物理切口,必须要有数据支撑。
不是说他物理那一头脑补得再漂亮,就一定是对的。
他得拿数学这边能跑出多远来做参照。
万一数学方向真的只是小卡一下,後面一路畅通呢?
那他这一晚上的物理切口,就成了无用功了。
他不能拿一个直觉去推翻整个组的工程进度。
所以他得先去把数学这边推到尽头。
李东把包放在桌上,也不开电脑。
他先把草稿纸和笔拿了出来。
李东是搞纯数的。
反问题,严格的说他不算专业。
可是这些天,鄂伟南院士给的那一份书单,他已经啃下了一大半。
Hadamard的良态性条件、Tikhonov正则化、谱截断、Backus-Gilbert方法…这些名字加上对应的核心定理,全都紮紮实实地挂在他记忆宫殿的【应用数学/反演】那一格。再加上他理解东西的速度,本来就在常人之上。
因此他也有信心不比一般走应用数学的学者差。
他坐下来,开始推导。
第一步:测量算子K:T→S的形式定义,T是配位场张量空间,S是隧穿谱空间。
第二步:核函数展开,把K拆成本征函数序列。
第三步:上吉洪诺夫正则化,调整参数α。
第四步:把K的奇异值谱分布,对应到T的低秩近似。
第五步……
他一步一步地推。
每一步,他都对照着昨晚已经搭好的那一套物理参数,看数学这一头能不能合上。
基本上前面五六步,都没有出大问题。
Tikhonov那一段他卡了不到十分锺,就用一个加权范数把它处理掉了。
基函数的选择那一段,他换了一组Whittaker型的局部基底,规避了一部分相位混淆。直到第七步。
李东看着自己写出来那一行,皱起了眉头。
【在(L2(R,du),〈·,·〉_K)上对相位算子做谱分解,要求三组基底的过渡矩阵收敛於一个无相位歧义的西变换……】
“是这样吗?”
“不对呀。”
这一段,本来是要把不同基底的展开“硬合”在一起。
按列正常思路,是用一个隐性约束量去反卡这个过渡矩阵的上界。
可是李东推到这里,就是过不去。
他换了三种正则化,结果一种比一种糟糕。
他试着把鄂伟南院士点出来的那一根“隐性约束”硬塞进去。
卡得住,又卡不死。
第三峰从糊里头能拎出来一点点边,但拎到一半,又掉回去了。
李东盯着那一行字,叹了一口气。
“……果然啊。”
“全国十几个组,全都卡在这。”
就在他第三次放下笔的时候,背後突然冒出一个声音。
“卡住了?”
李东猛地回过头。
吴开教授不知什麽时候已经站在了他的身後。
李东赶紧站起来。
“吴老师。”
“您来了。”
吴开摆摆手,走过来在他对面那张椅子上坐下。
“我刚从化院开完会过来。”
他看了一眼李东桌上的草稿纸。
看见最上头那一张推到一半的过渡矩阵谱分解。
吴开点了点头。
这一块,他也是看得懂的。
不是因为他数学多好。
是因为这一块,全圈儿的人都卡在这儿。
卡得太多了,他不熟也熟了。
吴开笑了笑。
“李东,你别急啊。”
“这个东西,已经卡好几年了。”
“国内国外大家都卡在你这一步呢。”
“你今天上午就推到这儿,已经比绝大多数人快了。”
“再说,你也别给自己压力。”
他温和地看着李东。
“我也没想过你一来就能解决问题啥的。”
“咱们组的进度计划,是按“数学这一块磨一两年’排的。”
“前面那一堆样品和针尖,咱们可以先把後续的衔接做了。”
“等你这边反推回来,咱们再把它接上,也是可以的。”
吴开这话说得是真心的。
他怕李东受挫。
李东毕竞还年轻。
基础数学那一头,他跑得是真的太快了。
这种年少成名的孩子,最怕的就是在另一个领域里被一面墙撞得头破血流。
吴开过来的不光是问候。
他是来给李东托个底的。
李东点了点头。
他听出了吴开的好意。
心里挺感激。
他正要开口说“老师我会好好磨的”。
结果话还没出口,他想起来一件事。
他刚才推到第七步卡住的时候,没再继续往第七步那一堵墙上撞………
他用了另一种“绕”的方法,从纯数那一边借了一招过来,把第七步先放到一边,跳到了第八步。第八步往下,他一口气推了好几页。
也不知道这样行不行,毕竟大家都没用纯数的方法来,万一在这里不行呢?
所以他问道。
“老师,这後头那一段我用的自己的方法,您给看一下这样行不行。”
李东说完,就低头去翻草稿纸最下面那几张。
然而听到这话的吴开整个人就在旁边僵住了。
啥意思啊?你推过去了?
吴开数学不算精通。
他这辈子混的是化学合成。
但数学这一头,因为这个项目卡了好几年,他这一两年也开始“陪着读”。
反演这里头,每一步在哪儿、为什麽会卡、卡在哪个具体的点上…
他门儿清。
刚才李东桌面上那一张过渡矩阵谱分解。
他一进门就看出来卡在哪了。
可是李东刚才嘴里冒出来的那一段……
“後头那一段”?
哪一段“後头”?
这一步根本就没“後头”啊。
所有人都是卡死在这一步。
吴开开口道。
“等一下。”
“李东,你刚才说的,是什麽意思?”
李东抬起头。
“啊?”
“哦。”
他指了指自己手里翻到一半的那几张纸。
“就是这一步过去以後,会遇到什麽样的一个问题。”
“我本来是想,您给看一下我这边推得对不对。”
李东说得理所当然。
吴开整个人愣了一下。
“这一步过去以後?”
他有点不太确定。
李东没听出他语气里的反应。
他的注意力还在自己那几页推导上。
“嗯。”
“您看这一段……这里我用的是局部Whittaker模型那一边的过渡函子,把三组基底的相位差先冻在一个固定的代数结构里……”
吴开听到“局部Whittaker模型”五个字,他整个人都是懵的。
这是什麽??
他从应用数学那一头听过基函数、算符、反演,可从来没听过局部Whittaker模型啊。这不是反演里的语言。
这是……
朗兰兹纲领的语言。
吴开的手不自觉地往李东那堆草稿纸抬去。
“你这里………”
“不止一张纸啊?”
他这一抬手,才发现李东桌面上那一遝草稿纸,最上头那一张底下,还有厚厚一摞。
吴开下意识地把最上面那一张过渡矩阵谱分解抽走,露出底下第二张。
他低头一看。
第二张上头,密密麻麻的全是符号。
最顶上那一行,他认得……
【设过渡矩阵U_ij之间的相位歧义为中∈R/2mZ,由偶极宇称选择定则有sgn(M_dipole) ∈ {+1,-1}固定】
吴开看到“偶极宇称选择定则”这八个字的时候,眉毛挑了一下。
这一路是物理那一边的玩法。
可下头紧接着又是一行……
【在GLs的局部基变换函子下,将三组基底的相位差归约为一个有限维代数表示的同构类,由L因子的局部根数决定】
这一行下面,李东用了一个李东自己常用的对关联函数符号F_r(α),把中的可能取值死死卡在一组离散的值上。
最後一步………
他把F_n(α)的零点和宇称选择定则给的sgn合起来,得到一个干干净净的结论:
【中=0】
吴开盯着那一行【o=0】。
看了三秒。
看了五秒。
他自己都没意识到,他端保温杯的那只手,已经把杯子里的水洒到了袖口上。
第三张。
第四张。
每一张推到後面,李东开始展开第三配位壳层那一套。
奇怪的是,李东这边推第三壳层,没有再用任何一种正则化。
他用的是另一组语言。
吴开看不太懂………
但他认得最关键的那几个名词。
“WKB”。
“势垒衰减”。
“伽莫夫因子”。
吴开“嗖”地一下抬起头,盯着李东。
他的喉结动了一下,吸了一口气。
“李东。”
“你……”
吴开下意识想问……
你这是把後面的反演也做了吗?
可这一句又咽回去了,因为他感觉这不是反演。
所以他换了一句最接近他直觉的话。
“你把这一步……推过去了?”
李东抬起头。
他眨了眨眼,明显没明白吴开问的是哪一步。
“什麽推过去了?”
“没有啊。”
李东很认真地摇头。
吴开:……
吴开把手里那一张草稿纸往桌上一摁,指着第二张【@=0】那一行。
“那这是啥?”
李东低头瞄了一眼。
“哦,这个啊。”
“我刚才在第七步这儿,把基函数相位歧义拆掉了。”
“再後头第八步、第九步,是接着第三配位壳层那一段往下推的。”
“再後面,就过不下去了。”
吴开盯着李东那一张坦然的脸,足足愣了三秒。
他终於反应过来了。
李东刚才嘴里说的那一句“过不下去”
不是他、不、不是国内外都卡死的那一堵墙!
李东说的“过不下去”
是这一整摊反演问题作为一个整体,过不下去!
而眼前这一面所有人头破血流的墙,他……
吴开干脆把【中=0】那一行拍出来,直接指着最上面那一张推导。
“这一步……”
“你看一下,就这一步啊!”
“你推过去了吗?”
李东低头看了一眼。
他抬起头,对着吴开说了一句。
“嗯。”
“这一步不是很难。”
吴开:???
吴开心里立刻就开始吐槽了。
不是很难?
不是很难,全国十几个组会卡在这一步整整三年?
李东压根没意识到吴开心里那一阵翻江倒海。
他指了指自己那一行【o=0】下面的推导。
“这里我是用Ls的局部基变换函子做的。”
“就是把三组基底,先看成同一个全局自守表示在不同局部位置上的Whittaker模型。”“再用L因子的局部根数把过渡矩阵的相位差锁死在有限维代数表示里。”
“最後用偶极宇称选择定则给的实符号收尾。”
“基本上不算太绕。”
李东说得很自然。
吴开听得整个人都僵了。
头一句“L.”,他还能跟。
第四个字开始,“局部基变换函子”……他直接断片了。
纯数那一头,朗兰兹纲领、自守表示、Whittaker模型………
这些东西离他十万八千里。
吴开无语了好半天。
最後他憋出一句。
“李东。”
“你说人话。”
李东:……
李东这才算是反应过来。
原来这一面墙,全国十几个组卡在这儿。
不是卡在“反演”这一整块出不来。
而是卡在“反演”里头具体的一步。
就是基函数相位差那一步。
他刚才一直以为,大家是卡在反演这一摊的总瓶颈上。
所以吴开教授刚才说“卡好几年了”,他下意识地是把“几年”往整个反演那个大问题上想。结果不是。
大家的几年,是卡在他半小时绕过去的那一步上。
李东也不能笑话人家。
他自己很清楚。
他能半小时绕过去,是因为他手里的工具,本来就不是从反演里借的。
他手里的工具,是朗兰兹纲领。
对一个搞应数的学者来说,他根本不会想到要从GL.的局部基变换函子那一头,往一个谱反演问题里头钻。
这是两套语言。
这是两个圈子。
如果不是李东这个本来就是给朗兰兹纲领“夯地基”的人,又恰好进了这个项目组。
这两条路永远都不会装上。
因为有一个搞纯数的学者会跑来这里高应用数学。
所以这一步在外人看来是天堑。
对他来说,是顺手。
李东想通了这一点。
他抬起头,看着吴开。
深吸了一口气。
“吴老师。”
他认认真真地说。
“这一步,其实不是真正的问题。”
吴开还在震惊李东的进度,下意识的问道。
“你说什麽?”
李东语气平静。
“咱们项目,乃至全世界这十几个组,真正卡的·……”
“不是这一步。”
“用数学反演这条路,本身就走不通。”
屋子里静了下来。
换做是十分锺以前,李东在他面前讲这一句话一
吴开顶多礼貌地笑一笑,表示尊重,然後心里默默给这小子记一笔“还得磨”。
毕竟,那是一句没有任何依据的话。
可现在不一样了。
现在李东已经把全球十几个课题组卡了三年的那一步,在一个上午里推了过去。
不光推过去,还顺手在後头推了好几页。
李东这一句“反演这条路本身就走不通”。
就不是一个少年人凭直觉甩出来的话了。
这是一个已经站在墙顶上的人,回过头来,对底下还在撞墙的人说的话。
吴开让自己先冷静下来,然後才开口问道。
“那你跟我说说。”
“为什麽你觉得不行?”
李东也没绕。
他组织了一下语言说道。
“老师。”
“反演这一块,本质是一个病态算子的逆问题。”
“病在什麽地方?”
“病在它的核是不完整的,从隧穿谱反推回配位场张量的过程里,输入端任何一个噪声,都会被那个算子的小奇异值放大几个数量级。”
“我们现在的所有数学工具。”
“吉洪诺夫、谱截断、Backus-Gilbert本质都是在压那个奇异值谱。”
“压的方法不一样,但是压的逻辑是一样的……”
“我们硬塞一个先验进去,告诉算子:解应该长什麽样。”
吴开点了点头。
这一段,是他这一两年陪着读应数最熟的一段。
李东接着说。
“问题在於,所有的先验,都是数学先验。”
“什麽叫数学先验?”
“就是,我希望解是光滑的、我希望解是稀疏的、我希望解的范数最小。”
“这些希望,跟单原子上那个铁原子周围的电子云长什麽样,没有任何物理上的关系。”
“我们用一个数学的偏好,去压一个物理的解。”
“两边的逻辑根本对不上。”
“所以第三配位壳层的那一根峰,永远糊。”
吴开微微皱起了眉。
李东又说。
“还有就是………”
“反演问题本身,是把过程的结果倒过来推过程。”
“你想清楚没有,这台机器测的是什麽?”
“测的是一个X射线激发、电子跃迁、空轨道伸出、隧穿到针尖这样一个完整的物理过程产生的电流。”“它本来就是过程的结果。”
“我们却在反过来,从结果倒推过程。”
“那中间任何一段过程的物理细节,都会被压成一个有效算子,然後我们再去硬解这个算子。”“这不就等於我们把一段活生生的物理,先杀成一组冷冰的算符。”
“再去问这组算符背後那段物理长什麽样?”
“这中间天然有一道墙。”
“墙不是物理给的,是我们自己造的。”
吴开听到这儿,眉头皱得更紧了。
他没说话。
李东深吸了一口气。
“老师,您问我有没有解决办法。”
“我有一个方向。”
“我有点不确定。”
“但我是想试着说一下。”
吴开看着他。
“你说说看。”
李东把今天清晨在草稿纸上搭好的那一套切口,简单地说了出来。
“我们不要做反演。”
“我们做正向。”
“把那一组测量算子的逆问题,重新写成一个用物理参数化的正向模型。”
“参数不再是配位场张量的抽象分量。”
“参数是……”
“偶极跃迁矩阵元、未占据轨道的对称性、轨道相对於针尖的取向、以及真空段的WKB衰减常数K。”“再把X射线偏振矢量正大光明地搬进来,作为模型的输入边界条件。”
“偏振矢量在源头上选择哪条空轨道被激发……”
“这件事本来就是物理在做的,不需要让基函数去吵。”
“再有就是,第三配位壳层那一根糊掉的峰。”
“我们不去拎它。”
“我们让它自己长出来。”
“每一个配位壳层,对应在隧穿这一关,有自己一根独立的KB衰减常数K。”
“越外的壳层《越小,从针尖底下伸出去的尾巴越长。”
“换而言之,第三壳层在隧穿电流的衰减曲线上,本来就有一段独属於自己的尾巴。”
“这段尾巴不需要去拎。”
“它在距离一档一档拉远的扫描里,会自己显出来。”
“靠的就是那一根伽莫夫因子。”
吴开听到“伽莫夫因子”四个字的时候,神色微微动了一下。
李东最後说。
“老师,简单总结一下。”
“两堵墙我们不推。”
“因为两堵墙之间,本来就有一道门。”
实验室外头实验台上的“嘀嘀”声远远地传过来。
办公室里没人说话。
吴开放在桌上的那只手,慢慢的握了起来。
他第一反应是不可能。
怎麽可能?
这这个方向已经一百年麽没人提出过新的思路了。
因为走不通,所以最後大家会一起转向数学反演
这是一条所有人筛了又筛、筛剩下来的路。
怎麽可能突然冒出来一句“两堵墙中间本来就有一道门”?
怎麽可能这道门是物理?
吴开忍不住开口。
“李东。”
“我先不说你那一套数学物理上能不能合得上。”
“我说一个最朴素的事。”
“你说要正向。”
“那意味着,你要在源头上,把那个X-ERT过程从头到尾,都用物理量参数化。”
“可这里头,有一段东西”
“针尖上方那一小段真空隧穿,再加上配体场和Au(111)衬底的相互作用-……”
“它的有效势垒形状,谁都没真正测过。”
“你那个伽莫夫因庙里的K,到底丫多宁?”
“你不知道,我不知道,没人知道。”
“我们之所以走反演,就丫因为这一段我们灭不过去。”
“我们没法直接去测那一段势垒。”
“所以我们才愿意承认,反演虽然病抗……”
“但起码,每一项里都是“实验值’,再病态它也有数据撑着。”
吴开顿了一万。
“你那一头要的“乘理参数化’。”
“它从一开始,源头那一段就丫空的。”
“你拿什麽去填?”
吴开越说越觉得自己说在了点上。
这一块的核心矛盾就在这儿。
大家不丫不想做正向,丫源头那一段没法直接测。
所以才退一步,做反演。
所以他最後摇了摇头。
他的声4轻,但很笃定。
“李东,你这套思路,整体上看丫漂亮。”
“但丫从化学这一行的逻辑……”
“它太悬了。”
“我跟你直说。”
“这一套很不化学!”
屋里又安静了丐来。
李东低头看着自己那一堆草稿纸。
最上面那一张,丫他拿红笔涂的那一团。
红色那一团旁边,他写的那五个字。
“让原庙自己说”。
他抬起头来,看着吴开。
眼神平静,没有半点不服。
“老师,我懂您的意思。”
“在化学这条路上,看上去丫没有路的。”
“但………”
李东也很笃定的说道。
“老师。”
“它很乘理!”
http://www.xuanhundaozhang.com/yt127592/48690861.html
请记住本书首发域名:www.xuanhundaozhang.com。玄浑道章手机版阅读网址:www.xuanhundaozhang.com