24小时热门版块排行榜    

查看: 1047  |  回复: 1

styrenes

铁杆木虫 (著名写手)


[资源] AI求解薛定谔方程,兼具准确度和计算效率,登上《自然-化学》 已有1人参与

AI求解薛定谔方程,兼具准确度和计算效率,登上《自然-化学》
2020-12-23 12:39
编辑:杜伟、 魔王 、小舟

作为量子力学的基础方程之一,薛定谔方程一直广受关注。去年,DeepMind 科学家开发一种 新的神经网络来近似计算薛定谔方程 ,为深度学习在量子化学领域的发展奠定了基础。今年九月份,柏林自由大学的几位科学家提出了一种新的深度学习波函数拟设方法,它可以获得电子薛定谔方程的近乎精确解。相关研究发表在 Nature Chemistry 上。



即使并非物理学界人士,我们也对薛定谔这个名字并不陌生,比如「薛定谔的猫」。著名物理学家埃尔温 · 薛定谔是量子力学奠基人之一,他在 1926 年提出的薛定谔方程(Schrödinger equation)为量子力学奠定了坚实的基础。薛定谔方程是描述物理系统的量子态怎样随时间演化的偏微分方程,是量子力学的基础方程之一。

在经典力学里,人们使用牛顿第二定律描述物体运动。而在量子力学里,类似的运动方程为薛定谔方程。薛定谔方程的解完备地描述物理系统里微观尺寸粒子的量子行为,包括分子系统、原子系统、亚原子系统。微观系统的状态由波函数来描写,薛定谔方程即是波函数的微分方程。若给定了初始条件和边界的条件,就可由此方程解出波函数。另外,薛定谔方程的解还可完备地描述宏观系统,可能乃至整个宇宙。

求解薛定谔方程可以为化学反应提供线索。化学反应的结果基本上与电子以及它们环绕原子和分子的方式有关。而控制事物反应的能量以及电子在分子中的轨道的差异决定了化学物质的形状,也由此决定了其性质。计算这一能量的方式就是求解薛定谔方程。换句话说,求解出薛定谔方程,就可以知道化学反应的结果。

然而,这并非易事。此前,我们可以精确求解的原子只有氢原子——仅具备一个质子和一个电子。

最近,来自柏林自由大学的科学家提出利用人工智能计算薛定谔方程的基态解,相关研究发表在 Nature Chemistry 上。



论文地址:https://www.nature.com/articles/s41557-020-0544-y

用 AI 求解薛定谔方程

量子化学旨在预测分子的化学和物理性质,它仅利用分子在三维空间中的原子排列来完成。这可以减少对资源的需求,并加快实验速度。理论上,这可以通过求解薛定谔方程来完成,但在实践中这往往非常困难。目前,人们仍无法高效求得任意分子的精确解。

最近,来自柏林自由大学的科学家提出一种深度学习方法,达到了前所未有的计算效率和准确度权衡。



该研究作者之一 Frank Noé 教授表示:「我们认为这一方法或将极大地影响量子化学的未来。」

无需在准确度和计算成本之间做出取舍

波函数是量子化学和薛定谔方程的关键所在,是一种描述分子内电子行为的函数。它是一种高维实体,这使得捕获编码特定电子之间相互影响方式的频谱变得极度困难。

量子化学领域中的许多方法不再只是尝试以数学方式获得特定分子的能量,但这需要近似值,限制了预测的质量。还有一些利用大量简单数学构造块表示波函数的方法,但这些方法过于复杂,难以针对较多原子计算波函数。

该研究一作 Jan Hermann 设计了新方法的关键特征,他表示:「避免在准确度和计算成本之间进行权衡是量子化学的最高成就。」

将物理属性引入 AI 神经网络

Hermann 表示:「到目前为止,最流行的方法是极具成本效益的密度泛函理论。我们认为我们提出的深度『Quantum Monte Carlo』方法至少可以达到同样好的效果。该方法以可接受的计算成本提供了前所未有的准确度。」

该研究设计了一个深度神经网络来表示电子的波函数,这是一种全新的方法。Noé 解释说:「我们没有使用用相对简单的数学成分组成波函数的标准方法,而是设计了一种人工神经网络,它能够学习电子围绕原子核运动的复杂模式。」

Hermann 表示:「电子波函数的独特性在于反对称性。在交换两个电子时,波函数需要改变符号,我们必须将这种特性引入到神经网络架构中才能使之奏效。」

受泡利不相容原理启发,PauliNet 方法诞生

受到「泡利不相容原理」(Pauli exclusion principle)的启发,研究者将他们的方法命名为「PauliNet」。它是一种深度学习波函数拟设,可以获得电子薛定谔方程的近乎精确解。PauliNet 具有一个作为基线的内置多参考哈特里-福克(Hartree–Fock)解,集成有效波函数的物理特性,并使用变分量子蒙特卡罗方法(variational quantum Monte Carlo, VMC)进行训练。

PauliNet 拟设架构的信息流如下图所示:



在实验部分,研究者采用了用于 DeepWF(Han et al., 2019)的相同系统,具体为氢分子(H_2)、氢化锂(LiH)、铍(Be)以及硼(B)和线性氢链 H_10。研究者将 PauliNet 与 SD-VMC(singledeterminant variational, 标准单行列式变分蒙特卡罗)、SD-DMC(singledeterminant diffusion, 标准单行列式扩散蒙特卡罗)和 DeepWF 进行了比较。

结果表明,PauliNet 的性能优于这三种用于原子、双原子分子和强相关氢链的 SOTA VMC 拟设方法,并且具有较高的计算效率。下表 1 为使用这四种不同方法时,H_2、LiH、Be、B 和 H_10 五种系统的基态能量对比:
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
czyzsu2楼
2020-12-25 09:08   回复  
送红花一朵
谢谢分享!
相关版块跳转 我要订阅楼主 styrenes 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26申博自荐 +8 烷基八氮- 2026-03-06 10/500 2026-03-08 01:44 by wszdzda
[考研] 材料325求调剂 30+6 mariusuki 2026-03-02 12/600 2026-03-07 20:12 by xindong
[考研] 286求调剂 +10 Faune 2026-03-06 11/550 2026-03-07 17:42 by fantasticxu
[考研] 274求调剂 +7 cgyzqwn 2026-03-01 14/700 2026-03-07 14:51 by Miko19
[考研] 085600材料与化工 298 调剂 +11 小西笑嘻嘻 2026-03-03 11/550 2026-03-07 10:30 by 一直走不要停
[考研] 085701环境工程295求调剂,一志愿A区211 +7 等春来, 2026-03-06 8/400 2026-03-06 22:04 by 星空星月
[考研] 材料学硕318求调剂 +21 February_Feb 2026-03-01 24/1200 2026-03-06 16:08 by @飒飒飒飒
[考研] 一志愿武理085601专硕347分求调剂 +4 啊欧欧欧 2026-03-04 5/250 2026-03-05 19:27 by Leding1356
[考研] 302材料工程求调剂 +7 Doleres 2026-03-01 8/400 2026-03-05 09:59 by oxidpl
[考研] 281电子信息求调剂 +5 jhtfeybgj 2026-03-02 9/450 2026-03-05 09:36 by 085400调剂党
[考研] 0856材料与化工,270求调剂 +17 YXCT 2026-03-01 20/1000 2026-03-05 09:22 by 一切OK
[考研] 材料与化工,291分,求调剂 +9 咕噜咕噜123123 2026-03-04 11/550 2026-03-05 09:03 by 斩魂滴兔子!
[考研] 一志愿武汉理工大学-085602-总分296分-求调剂 +7 紫川葡柚 2026-03-04 7/350 2026-03-04 21:04 by kakakapanpan
[考研] 286 +6 ksncj 2026-03-04 6/300 2026-03-04 20:49 by lature00
[考研] 322,求调剂 +3 菜菜爱玩 2026-03-04 3/150 2026-03-04 12:15 by xiongyaxuan
[基金申请] 请问大家,研究风险与应对措施那里, 大家都怎么写呢 ? +3 cauasen 2026-03-02 3/150 2026-03-03 19:43 by Tide man
[考研] 化学0703求调剂 学硕 理/工科均可 总分279 +3 1一11 2026-03-03 5/250 2026-03-03 12:37 by 1一11
[考研] 338求调剂 +5 18162027187 2026-03-02 6/300 2026-03-03 07:52 by njzyff
[考研] 求调剂 +11 yunziaaaaa 2026-03-01 13/650 2026-03-02 21:59 by sunny81
[考研] 一志愿中石油(华东)本科齐鲁工业大学 +3 石能伟 2026-03-02 3/150 2026-03-02 18:54 by caszguilin
信息提示
请填处理意见