第112章 你是怎么发现的?
肖宿走上讲台,插上u盘,打开ppt。
第一页,標题很简单:对称性约束下的加权度量构造及其在计算化学中的应用。
他看向台下,开口说:“之前,万教授找我帮忙解决一个计算化学的问题。”
没有寒暄,没有客套,直接进入正题。
“问题本身是关於含过渡金属的有机配合物的电子结构计算。这类分子的对称性高,电子关联强,传统方法处理不好。”
他调出一张分子结构图,是一个八面体的金属配合物,中心原子是铁,周围连著六个配体。
“万教授尝试用我上一篇论文里的加权度量方法,但没有成功。原因在於,他们在奇点附近的数值离散化破坏了对称性约束。”
肖宿在图上画了几个圈,標註出奇点所在的位置。
“对称性约束在量子化学计算中很重要。分子轨道属於群的不可约表示,哈密顿量在不同表示之间是对角的。如果数值处理破坏了这种对角性,计算结果就会出现偏差。”
台下,姜明微微点头。
这孩子讲得清楚,不仅知道问题在哪儿,还知道为什么重要。
肖宿继续讲那个加权函数的构造。
他用了大概二十分钟,把整个数学框架梳理了一遍。
从群作用的轨道分解,到奇点附近的局部坐標,再到加权函数的显式表达式。
讲完之后,他停顿了一下,说:“这个构造本身不难,但有个副產品。”
ppt切换到下一页。
“atkins的《molecular quantum chemistry》,第五章关於分子轨道对称性与化学键稳定性的推导,可以进行简化处理。”
台下突然安静下来。
atkins的书?
那个全世界几代化学学生都在用的经典教材?
肖宿调出一张截图,正是那本书第187页的內容。
“这里的推导,假设了某个积分在对称操作下保持不变。这个假设在阿贝尔群的情形下成立,但在非阿贝尔群,比如四面体群、八面体群的情形下,是不成立的。”
他用红圈標出那个关键步骤。
“正確的处理应该是用群轨道平均的方法,引入一个修正因子。”
ppt上出现了一个公式:
e_stab = Σ_i c_i x_i(g) · w(r)
c_i是轨道係数,x_i(g)是特徵標,w(r)是由对称群决定的修正因子。
“用这个修正公式重新计算一些高对称性的分子,结果会有明显变化。”
他调出几张图表,是甲烷、苯、一个八面体配合物的计算结果对比。
“甲烷的稳定性指標,修正后提高了大约百分之三。苯的变化小一些,百分之一点几。这个八面体配合物变化最大,修正后的稳定性能比修正前高了將近百分之十。”
他顿了顿,说:“对应的实验数据支持修正后的结果。”
台下彻底安静了。
姜明盯著屏幕上那些图表,眼神里闪过一丝复杂。
atkins的书他也读了四十年了。
给学生讲这门课讲了三十多年。
从来没有想过,那个被无数人引用、验证、再验证的推导,居然有缺陷。
而且这个缺陷,一个十五岁的孩子来指出来了。
华清大学的符教授举起手。
肖宿点点头:“请讲。”
符宴民站起来,语气很认真:“肖宿同学,你刚才说的这个修正,我大概理解了。但我想问的是,这个修正只影响高对称性的分子吗?有没有可能在一些低对称性的体系里,也存在类似的问题?”
肖宿想了想,说:“理论上,只要分子的对称群是非阿贝尔的,这个修正项就不为零。但实际影响的大小,取决於分子轨道的表示类型和相互作用的强度。”
他在白板上画了一个简单的示意图。
“阿贝尔群的不可约表示都是一维的,特徵標就是±1这种简单数值。非阿贝尔群有高维表示,特徵標是矩阵。
积分在群作用下混合的,就是这些高维表示內部的各个分量。”
“如果分子轨道都属於一维表示,那修正项就很小。如果有轨道属於高维表示,修正项就不可忽略。”
符宴民听完,若有所思地点点头,坐下了。
又一个研究员举手,问的是加权函数的数值实现细节。
肖宿答得很简洁,但每一步都清晰。
姜明一直没说话,只是安静地听著。
他看著台上那个少年,看著他在白板上写下一串串公式,看著他不紧不慢地回答每一个问题,心里涌起一种复杂的感觉。
他自己年轻时刚接触量子化学的时候,花了整整一年,才把群论在分子轨道中的应用搞明白。
而现在这个孩子,接触化学才两周,就已经能指出经典教材里的错误。
而且他指出的,不是那种印刷错误或者表述不清的小问题,而是一个在特定条件下会导致明显偏差的数学缺陷。
这意味著什么?
意味著这孩子读那本书的时候,不是像普通学生那样“接受知识”,而是在用自己的数学框架“审视知识”。
看到推导,想的应该是:这个假设成立吗?
在什么条件下成立?
如果不成立,会怎样?
然后他会用自己的工具去验证,去修正。
这才是真正做学问的態度。
提问环节进行到一半,一个年轻的研究员站起来,问了个问题:
“肖宿同学,你刚才讲的这些,都还是针对单个分子的静態性质。你有没有想过,把这些方法用到化学反应动力学上?比如模擬化学键断裂和形成的过程?”
肖宿沉默了几秒。
他確实没想过这个问题。
他看的那几本书,讲的都是平衡態的性质,像分子轨道、电子结构、光谱性质等等。
化学反应涉及的是动態过程,是分子从一个势能面跃迁到另一个势能面,是原子核的运动和电子的重新排布。
那需要的时间尺度、数学工具,和静態性质完全不同。
“还没想。”肖宿说,“但应该可以。”
年轻研究员笑了:“那我们等著你的下一篇文章。”
台下响起一阵笑声。
姜明这时候举起了手。
万匯杨眼睛一亮:“姜老,您请讲。”
姜明站起来,看向台上的肖宿。
他的头髮已经全白,但腰板挺得很直,两眼有神。
“肖宿同学,”他的声音有些沙哑,但中气很足,“我问你一个不一样的问题。”
肖宿点点头。
“你刚才说,你看atkins那本书,发现了这个推导的问题。我想知道的是,你是怎么发现的?是带著找茬的心態去读的,还是读著读著就觉得不对劲?”
第一页,標题很简单:对称性约束下的加权度量构造及其在计算化学中的应用。
他看向台下,开口说:“之前,万教授找我帮忙解决一个计算化学的问题。”
没有寒暄,没有客套,直接进入正题。
“问题本身是关於含过渡金属的有机配合物的电子结构计算。这类分子的对称性高,电子关联强,传统方法处理不好。”
他调出一张分子结构图,是一个八面体的金属配合物,中心原子是铁,周围连著六个配体。
“万教授尝试用我上一篇论文里的加权度量方法,但没有成功。原因在於,他们在奇点附近的数值离散化破坏了对称性约束。”
肖宿在图上画了几个圈,標註出奇点所在的位置。
“对称性约束在量子化学计算中很重要。分子轨道属於群的不可约表示,哈密顿量在不同表示之间是对角的。如果数值处理破坏了这种对角性,计算结果就会出现偏差。”
台下,姜明微微点头。
这孩子讲得清楚,不仅知道问题在哪儿,还知道为什么重要。
肖宿继续讲那个加权函数的构造。
他用了大概二十分钟,把整个数学框架梳理了一遍。
从群作用的轨道分解,到奇点附近的局部坐標,再到加权函数的显式表达式。
讲完之后,他停顿了一下,说:“这个构造本身不难,但有个副產品。”
ppt切换到下一页。
“atkins的《molecular quantum chemistry》,第五章关於分子轨道对称性与化学键稳定性的推导,可以进行简化处理。”
台下突然安静下来。
atkins的书?
那个全世界几代化学学生都在用的经典教材?
肖宿调出一张截图,正是那本书第187页的內容。
“这里的推导,假设了某个积分在对称操作下保持不变。这个假设在阿贝尔群的情形下成立,但在非阿贝尔群,比如四面体群、八面体群的情形下,是不成立的。”
他用红圈標出那个关键步骤。
“正確的处理应该是用群轨道平均的方法,引入一个修正因子。”
ppt上出现了一个公式:
e_stab = Σ_i c_i x_i(g) · w(r)
c_i是轨道係数,x_i(g)是特徵標,w(r)是由对称群决定的修正因子。
“用这个修正公式重新计算一些高对称性的分子,结果会有明显变化。”
他调出几张图表,是甲烷、苯、一个八面体配合物的计算结果对比。
“甲烷的稳定性指標,修正后提高了大约百分之三。苯的变化小一些,百分之一点几。这个八面体配合物变化最大,修正后的稳定性能比修正前高了將近百分之十。”
他顿了顿,说:“对应的实验数据支持修正后的结果。”
台下彻底安静了。
姜明盯著屏幕上那些图表,眼神里闪过一丝复杂。
atkins的书他也读了四十年了。
给学生讲这门课讲了三十多年。
从来没有想过,那个被无数人引用、验证、再验证的推导,居然有缺陷。
而且这个缺陷,一个十五岁的孩子来指出来了。
华清大学的符教授举起手。
肖宿点点头:“请讲。”
符宴民站起来,语气很认真:“肖宿同学,你刚才说的这个修正,我大概理解了。但我想问的是,这个修正只影响高对称性的分子吗?有没有可能在一些低对称性的体系里,也存在类似的问题?”
肖宿想了想,说:“理论上,只要分子的对称群是非阿贝尔的,这个修正项就不为零。但实际影响的大小,取决於分子轨道的表示类型和相互作用的强度。”
他在白板上画了一个简单的示意图。
“阿贝尔群的不可约表示都是一维的,特徵標就是±1这种简单数值。非阿贝尔群有高维表示,特徵標是矩阵。
积分在群作用下混合的,就是这些高维表示內部的各个分量。”
“如果分子轨道都属於一维表示,那修正项就很小。如果有轨道属於高维表示,修正项就不可忽略。”
符宴民听完,若有所思地点点头,坐下了。
又一个研究员举手,问的是加权函数的数值实现细节。
肖宿答得很简洁,但每一步都清晰。
姜明一直没说话,只是安静地听著。
他看著台上那个少年,看著他在白板上写下一串串公式,看著他不紧不慢地回答每一个问题,心里涌起一种复杂的感觉。
他自己年轻时刚接触量子化学的时候,花了整整一年,才把群论在分子轨道中的应用搞明白。
而现在这个孩子,接触化学才两周,就已经能指出经典教材里的错误。
而且他指出的,不是那种印刷错误或者表述不清的小问题,而是一个在特定条件下会导致明显偏差的数学缺陷。
这意味著什么?
意味著这孩子读那本书的时候,不是像普通学生那样“接受知识”,而是在用自己的数学框架“审视知识”。
看到推导,想的应该是:这个假设成立吗?
在什么条件下成立?
如果不成立,会怎样?
然后他会用自己的工具去验证,去修正。
这才是真正做学问的態度。
提问环节进行到一半,一个年轻的研究员站起来,问了个问题:
“肖宿同学,你刚才讲的这些,都还是针对单个分子的静態性质。你有没有想过,把这些方法用到化学反应动力学上?比如模擬化学键断裂和形成的过程?”
肖宿沉默了几秒。
他確实没想过这个问题。
他看的那几本书,讲的都是平衡態的性质,像分子轨道、电子结构、光谱性质等等。
化学反应涉及的是动態过程,是分子从一个势能面跃迁到另一个势能面,是原子核的运动和电子的重新排布。
那需要的时间尺度、数学工具,和静態性质完全不同。
“还没想。”肖宿说,“但应该可以。”
年轻研究员笑了:“那我们等著你的下一篇文章。”
台下响起一阵笑声。
姜明这时候举起了手。
万匯杨眼睛一亮:“姜老,您请讲。”
姜明站起来,看向台上的肖宿。
他的头髮已经全白,但腰板挺得很直,两眼有神。
“肖宿同学,”他的声音有些沙哑,但中气很足,“我问你一个不一样的问题。”
肖宿点点头。
“你刚才说,你看atkins那本书,发现了这个推导的问题。我想知道的是,你是怎么发现的?是带著找茬的心態去读的,还是读著读著就觉得不对劲?”