毋庸置疑，诺贝尔化学奖，是我们的了

　　许铭把笔记本电脑打开，调出了肖宿写的代码。

　　两个人凑近，一左一右站在许铭身后。

　　屏幕上，肖宿写了四行注释安安静静地立在文件头部。

　　方清的目光从第一行扫到第四行，然后停住了。

　　“电子密度分布。”

　　他念出这几个字，声音有点干涩，“我们第一步算的是原子核坐标。”

　　许铭点头。

　　“三个多月，我们从第一步就错了。”

　　万汇杨把肖宿的代码从头到尾看了一遍，看到辛流形构造部分的时候，手指在屏幕前虚虚地点了一下。

　　“findSaddlePoints。”

　　他念出这个函数名，然后转过头看着许铭和方清，“肖宿在论文里写过：电子密度的鞍点构成了拉格朗日子流形的骨架。

　　但我怎么也没意识到这句话是实现的钥匙。”

　　“把VASP的测试数据准备好，用同一组分子体系，同样的计算精度要求，跑一遍对比测试。”

　　“好。”

　　方清很快把测试环境搭好了。

　　三台工作站并排运行，左边跑VASP，右边跑他们之前的旧框架，中间准备用来跑肖宿的新算法。

　　测试体系选的是十种常见分子：水、甲烷、苯、氨、甲醛、乙烯、乙炔、二氧化碳、甲醇、乙醇。

　　计算目标统一设定为基态电子结构，精度要求取化学精度。

　　“开始吧。”

　　许铭把肖宿的代码导入中间那台工作站，点击运行。

　　屏幕上的进度条开始快速移动起来。

　　第一个水分子实验，VASP用时四十七秒，旧框架用时三十二秒但结果偏差较大，而肖宿的新算法用时九秒，结果和实验值偏差低于零点零三电子伏特。

　　三个人互相看看，都能从对方的眼中看到震惊。

　　接下来是甲烷，VASP用时一分零三秒，肖宿的新算法用时十四秒，偏差低于零点零二电子伏特。

　　再接着是苯，这是最有挑战性的一个。

　　苯环的共轭π电子体系是DFT的传统软肋，VASP跑苯的基态电子结构需要两分四十秒，而且结果和实验值有大约零点二电子伏特的系统偏差。

　　而用肖宿的新算法跑苯，只用了三十一秒，结果和实验值的偏差还低于零点零一五电子伏特。

　　方清的手都有些颤抖了。

　　剩下的七个分子一个接一个跑完，每跑完一个，方清就在旁边的白板上记录下用时、偏差、和VASP的对比。

　　写完最后一个乙醇分子，尽管三人都有心理准备，可还是被惊到了。

　　十种分子，肖宿新算法的平均用时是VASP的四分之一不到，平均精度偏差是VASP的五分之一以下，算力消耗平均不到VASP的三分之一。

　　而且这还不是终点。

　　肖宿在代码里留了大量可以进一步优化的空间。

　　他给的是一个框架，一个方向，沿着这个方向往下走，性能还能再提升。

　　“这已经不是比VASP强了。”

　　万汇杨的声音有些沙哑，“这是重新定义了计算化学的底层逻辑啊。”

　　三个人站在白板前，看着上面那组数据，很久都没有说话。

　　谁都没法想象，这个实验数据公布之后，化学界将会迎来多么可怕的震动。

　　“我想，”万汇杨转过身，看着许铭和方清，脸上的疲惫一扫而空，“毋庸置疑，诺贝尔化学奖，是我们的了