中科大蒋彬课题组开发 FIREANN,分析原子对外界场的响应

内容一览: 使用传统方法分析化学系统与外场的相互作用,具有效率低、成本高等劣势。中国科学技术大学的蒋彬课题组,在原子环境的描述中引入了场相关特征,开发了 FIREANN,借助机器学习对系统的场相关性进行了很好的描述。

关键词: 化学物理学 分子动力学 外界场

作者 | 雪菜

编辑 | 李宝珠

化学系统与外场的相互作用在物理、化学及生物过程中至关重要。以电场为主的外场可以与原子、分子和凝聚态物质相互作用,导致电子或自旋极化,或是改变系统的空间取向。

密度泛函理论 (DFT) 和从头算分子动力学 (AIMD) 已被用于研究外加电场下的复杂周期和非周期系统。然而,AIMD 的应用要求很高,尤其是在核量子效应 (NQEs) 比较重要的系统中很难使用。

经验力场分析效率很高,但准确率有限,而精确的场相关量子散射计算只适用于非常小的系统

与此同时,机器学习 (ML) 在解决高维化学问题中取得了不俗的成果。然而,大多数机器学习模型将势能和系统对电场的响应分割看待,忽略了系统的场相关性。

为此,中国科学技术大学的蒋彬课题组,在原子环境的描述中引入了场相关特征,开发了场诱导递归嵌入原子神经网络 (FIREANN) 。FIREANN 不仅可以准确描述外场强度和方向的变化时,系统能量的变化趋势,还能对任意阶数的系统响应进行预测。这一成果已发表于「Nature Communication」。

在这里插入图片描述

这一成果已发表于「Nature Communication」

论文链接:

https://www.nature.com/articles/s41467-023-42148-y

FIREANN 模型链接:

https://github.com/zhangylch/FIREANN

关注公众号,后台回复「相互作用」获取论文完整 PDF

FIREANN:REANN + 伪原子场向量

FIREANN 的基础是 REANN 模型,后者通过嵌入原子密度 (EADs) 对原子环境进行描述。外加场时,电子密度会重新分布,系统的旋转不变性也会被破坏。系统和外场的相互作用显然会受到电场强度和方向的影响。

在这里插入图片描述

FIREANN 的架构

FIREANN 会为每个原子赋一个模拟真实原子行为的伪原子场向量,随后二者结合得到场相关嵌入原子密度,作为神经网络的输入,最终输出原子力、偶极矩、极化率等物理量。

每个原子的伪原子场向量可以表示为:

在这里插入图片描述

随后,将场相关轨道和高斯轨道 (GTOs) 组合成场诱导 EADs (FI-EAD) 矢量

在这里插入图片描述

这里,每个原子受到的外加场由伪原子相对于该原子的位置矢量表示。FI-EAD 便可以由原子间距离和封闭角重写:

在这里插入图片描述

实验验证

玩具模型:水分子的精准预测和外推

研究人员首先以水分子作为玩具模型 (Toy System) ,验证 FIREANN 对系统-外场相互作用的预测。在 yz 面上有一个水分子,x 方向上有强度为 0.1 V/Å 的电场。

由于外场和分子平面始终正交,分子的势能不会发生变化。FIREANN 准确预测到了这一结果。

同时,FIREANN 对分子沿 y 轴旋转时的偶极-电场相互作用进行了精准预测

在这里插入图片描述
FIREANN 对 yz 面水分子受外界电场影响的预测

a:水分子沿 x 轴旋转的情况;

b:水分子沿 y 轴旋转的情况;

c:电场强度变化时 DFT、FIREANN 和 FieldSchNet 的预测结果。

FIREANN 还有很强的外推能力,仅用单一的训练数据,便推理得到了电场强度在 -0.2-0.2 V/Å 时,分子势能的变化。这是传统的 FieldSchNet 模型做不到的。

NMA:IR 光谱的准确预测

FIREANN 的一个典型特征就是,它可以一步预测化学系统在有无外场时的能量和响应特性

研究人员在 N-甲基乙酰胺 (NMA) 上进行了测试。当外界电场在 0.0-0.4 V/Å 变化时,FIREANN 能够对 NMA 分子的能量、偶极矩和极化率进行有效预测,均方根误差 (RMSEs) 分别为 0.0053 eV、0.028 Debye 和 0.51 a.u

在这里插入图片描述

FIREANN 和 DFT 对 NMA 的能量 (a)、偶极矩 (b) 和极化率 © 的预测相关性图

FIREANN 还对场内的分子光谱进行了预测。当电场强度以 0.1 V/Å 的步长逐渐从 0.0 增加至 0.4 V/Å 时,C-O 伸缩带的变化最为明显。随着电场强度增加,C-O 伸缩带的 P/R 分支逐渐消失,吸收峰变得更加尖锐。

此外,FIREANN 还预测,外加电场会降低化学键的强度,导致 CO 的拉伸振动红移,距离与电场强度成正比。

在这里插入图片描述

300 K 下 NMA 在不同强度电场下的 FIREANN 预测结果

液态水:周期系统的高度吻合

为验证 FIREANN 模型预测周期系统对外加电场响应的能力,研究人员在液态水中进行了测试。与分子系统不同,周期系统的极化强度 (单位体积的偶极矩) 是一个多值量,导致其存在多个平行分支,引起偶极矩的突变。

AIMD 预测中我们可以清楚地看到偶极矩的突变带来的结果的不连续,在外加电场之后,这种突变会更加频繁,这为传统的机器学习算法带来了挑战。

在这里插入图片描述

AIMD、修正后的 AIMD 和 FIREANN 对无场 (a) 和有场 (b) 时液态水的偶极矩分析

而在 FIREANN 框架中,由于模型只在有电场的情况下进行了原子力的训练,系统的能量梯度其实是不受影响的,因此轻松绕过了这一问题。

为此,研究人员构建了一个包含 64 个水分子,x 方向上电场强度 0.6 V/Å 的模型,以原子力作为唯一的预测对象,称为 FIREANN-wF。模型对原子力的预测与实验高度一致,均方根误差仅 39.4 meV/Å

FIREANN-wF 对液态水的无场径向分布函数 (RDFs) 的预测与 DFT 和实验结果也是吻合的。

在这里插入图片描述

FIREANN 对 O-O (a)、O-H (b) 和 H-H © 的径向分布函数的预测,及与实验结果的对比

偶极矩对于 IR 谱具有重要影响。由于 FIREANN-wF 模型会对核量子效应进行分析,它能够对势能面 (PES) 和偶极矩面进行正确的预测,与 DFT 的结果一致

随后,利用 FIREANN-wF 对外加 0.4 V/Å 后的 IR 谱进行预测。由于电场降低了 O-H 键的强度,同时诱导水分子重定向与电场平行,光谱中 O-H 伸缩带出现了明显的红移。

在这里插入图片描述

FIREANN 对无场 (a) 和有场 (b) 情况下液态水的光谱预测及与实验结果的对比

对比 REANN:外推能力与高速训练

虽然此前已有与 FIREANN-wF 训练方式类似的模型,然而它们对外场的处理方式完全不同,导致这些模型无法对高阶的相互作用进行预测

而在 FIREANN 中,引入场相关原子轨道后,模型可以通过轨道间的相互作用捕获电子密度对外场的响应

前文已经对比过 FIREANN 和 FieldSchNet 在水分子中的差异,这一差异在周期系统中依然存在。

研究人员利用水分子和 x 方向的电场建立了测试体系。FIREANN 和 FieldSchNet 的预测均方根误差分别为 54.5 meV/Å 和 245.4 meV/Å。与之前的结果类似,FIREANN 可以将预测外推到 ±2 V/Å,而 FieldSchNet 不具备这一能力

在这里插入图片描述

DFT、FIREANN 和 FieldSchNet 对液态水体系能量在电场变化时的外推结果

训练时间上看,在单张显存为 80 GB 的 A100 上, FieldSchNet 的单个 epoch 为7.6 分钟,而 FIREANN 仅用 2.4 分钟

分子-场相互作用:微观系统的遥控器

化学系统和外场的相互作用为人们研究微观系统提供了窗口,更为微观系统的操纵提供了有力的工具。通过调控外加电场,可以改变物质的化学结构、促进电子转移、控制物质相变和生物分子的构象变化、调整催化剂的选择性,甚至影响冷化学反应的量子动力学。

在扫描隧道显微镜的尖端和金属表面之间施加电场,可以将金属表面的偶氮苯衍生物可逆的反式-顺式异构化。

在这里插入图片描述
金 (111) 表面偶氮苯的反式-顺式异构

同样的,改变电场的取向,可以在纳米尺度上改变分子的混合模式。

在这里插入图片描述

1,3,5-三(4-羧基苯基)苯和三聚硅酸在不同电压下的混合模式

可以说,分子与外场的相互作用便是微观系统的遥控器。理解这一相互作用,对于微观尺度的科学研究具有重要意义。FIREANN 能够准确分析周期系统和非周期系统与外场的相互作用,并对任意阶数的系统响应进行预测,为微观研究提供了新方法

参考链接:

[1]https://pubs.acs.org/doi/full/10.1021/ja065449s

[2]https://pubs.acs.org/doi/full/10.1021/acsnano.7b04610

—— 完 ——

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/177252.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

js执行异常处理 箭头函数 正则表达式

JS代码执行过程中,可能会出现异常 try-catch-finally try块包含可能抛出异常的代码,catch块包含处理异常的代码。finally 块始终执行,无论是否抛出异常。 try代码块出现异常 执行catch代码块 finally不管try代码有无异常 都会执行fina…

NX二次开发UF_CURVE_convert_conic_to_gen 函数介绍

文章作者:里海 来源网站:https://blog.csdn.net/WangPaiFeiXingYuan UF_CURVE_convert_conic_to_gen Defined in: uf_curve.h int UF_CURVE_convert_conic_to_gen(UF_CURVE_conic_p_t conic_data, UF_CURVE_genconic_t * gen_conic_data ) overview 概…

给定序列a,选k个数排成一排,从左往右扫,如果当前数小于上一个数,那么当前数变成上一个数,然后形成最后的序列,问形成的序列有多少种

题目 #include<bits/stdc.h> using namespace std; #define int long long const int maxn 3005, mod 1e9 7; int a[maxn], b[maxn], cnt_low[maxn];//cnt_low[i]表示小于i的数字个数 int sum[maxn][maxn];//前缀和 int f[maxn][maxn];//f[i][j]表示选了i个数&#x…

第二十章——多线程

Windows操作系统是多任务操作系统&#xff0c;它以进程为单位。一个进程是一个包含有自身地址的程序&#xff0c;每个独立执行的程序都称为进程。也就是说每个正在执行的程序都是一个进程。系统可以分配给每一个进程有一段有限的使用CPU的时间&#xff08;也可以称为CPU时间片&…

中非市场开发攻略,走进非洲唯一废除死刑的国家

中非共和国是位于非洲中部的一个国家&#xff0c;是世界最不发达国家之一&#xff0c;经济以农业为主&#xff0c;工业基础薄弱&#xff0c;80%以上的工业品靠进口&#xff0c;中国是其重要的贸易伙伴之一。今天就来跟大家聊聊中非市场的开发攻略。文章略长&#xff0c;大家点赞…

2311d,d运行时勾挂

原文 上周,我开始转换_d_arrayliteralTX为模板.勾挂有点类似_d_newarrayT,因为用它来分配数组字面,如[1,2,3]. 注意,与也会初化新创建数组的_d_newarrayT不同,_d_arrayliteralTX只是为它分配内存.来限制传递给勾挂的参数数. 如果它也要初化数组,它必须按额外变参接收内容.我坚持…

E云管家微信群聊机器人开发

请求URL&#xff1a; http://域名地址/modifyGroupRemark 请求方式&#xff1a; POST 请求头Headers&#xff1a; Content-Type&#xff1a;application/jsonAuthorization&#xff1a;login接口返回 参数&#xff1a; 参数名必选类型说明wId是String登录实例标识chatRo…

Linux学习笔记6-串口应用

到现在为止都是在开发板上运行的裸机程序&#xff0c;相当于之前学习STM32单片机时走过的路&#xff0c;还没有真正进入到核心的驱动开发部分&#xff0c;但这都是基础&#xff0c;所以慢慢来不着急。 接下来进入串口通信的学习&#xff0c;和GPIO一样&#xff0c;也是和单片机…

Docker Swarm总结+service创建和部署、overlay网络以及Raft算法(2/5)

博主介绍&#xff1a;Java领域优质创作者,博客之星城市赛道TOP20、专注于前端流行技术框架、Java后端技术领域、项目实战运维以及GIS地理信息领域。 &#x1f345;文末获取源码下载地址&#x1f345; &#x1f447;&#x1f3fb; 精彩专栏推荐订阅&#x1f447;&#x1f3fb;…

全新爱蜗影视优码双端影视源码v9.1/影视视频APP源码+支持代理/在线支付+支持对应苹果CMS

源码简介&#xff1a; 这个是最新爱蜗影视优码双端影视源码&#xff0c;作为实用方便的影视视频APP源码&#xff0c;它不仅支持代理/在线支付&#xff0c;而且也支持对应苹果CMS。 爱蜗影视优码双端影视支持对应苹果CMS支持代理在线支付 带图文教程&#xff0c;全新美化多功能…

统信UOS上使用localsend传输文件

原文链接&#xff1a;统信UOS上使用localsend传输文件及文件夹方法 hello&#xff0c;大家好&#xff0c;今天带给大家一个全新的主题&#xff1a;在统信UOS上如何高效使用 localsend 软件进行文件和文件夹的传输。这篇文章将为您提供详细的指导和操作步骤&#xff0c;让您能够…

Couchdb 命令执行漏洞复现 (CVE-2017-12636)

Couchdb 命令执行漏洞复现 &#xff08;CVE-2017-12636&#xff09; 1、下载couchdb.py 2、修改目标和反弹地址 3、Python3调用执行即可 couchdb.py文件下载地址: https://github.com/vulhub/vulhub/blob/master/couchdb/CVE-2017-12636/exp.py ‍ 在VULFocus上开启环境 …

Java之面向对象《ATM自动取款机》

一、前言&#xff1a; 关于上次我写的博客文章中"Java之《ATM自动取款机》(面向对象)"&#xff0c;里面还不够完善&#xff0c;因为在各个服务功能相互跳转时&#xff0c;会出现混乱问题。这次我对其进行了修改和改进&#xff0c;若还有其它在大家测试时出现的bug请及…

Linux网络

Linux网络 1、Linux网络配置文件 ​ 查看第一张网卡的网卡信息&#xff1a; [rootlocalhost yum.repos.d]# cat vi /etc/sysconfig/network-scripts/ifcfg-ens33 cat: vi: 没有那个文件或目录 TYPEEthernet PROXY_METHODnone BROWSER_ONLYno BOOTPROTOstatic DEFROUTEyes I…

HarmonyOS-Service服务开发(一)

文章目录 创建新项目启动Serviceets获取service的bundleName DataAbility开发指导开发Data步骤创建Data 创建新项目 ServiceAbility开发指导 在config.json中也有配置出现 启动Service ets获取service的bundleName 项目的bundleName service的bundleName 这里serviceAbil…

vs2019 - MFC对话框程序的工程名称不支持下划线命名法

文章目录 vs2019 - MFC对话框程序的工程名称不支持下划线命名法概述笔记备注END vs2019 - MFC对话框程序的工程名称不支持下划线命名法 概述 正在写账单分析程序, 用MFC 对话框. 因为比较习惯下划线命名法, 就在向导中给工程名称起了一个my_test这样的名称(下划线命名法, 小…

车载通信架构 —— 传统车内通信网络CAN(可靠性为王)

车载通信架构 —— 传统车内通信网络CAN(可靠性为王) 我是穿拖鞋的汉子,魔都中坚持长期主义的汽车电子工程师。 老规矩,分享一段喜欢的文字,避免自己成为高知识低文化的工程师: 屏蔽力是信息过载时代一个人的特殊竞争力,任何消耗你的人和事,多看一眼都是你的不对。非…

antv/g6的学习总结

新建一个简单实例 1、使用命令行在项目目录下执行以下命令 cnpm install --save antv/g6 2、创建容器 <div id"mountNode"></div> 3、在需要用的 G6 的 JS 文件中导入 import G6 from antv/g6; 4、 数据准备 引入 G6 的数据源为 JSON 格式的对象。…

mysql 计算同比、环比,查询结果拼接成列、行、

一、理论 国家统计局同比、环比计算公式 增长速度是反映经济社会某一领域发展变化情况的重要数据&#xff0c;而同比和环比是反映增长速度最基础、最核心的数据指标&#xff0c;也是国际上通用的指标。在统计中&#xff0c; 同比和环比通常是同比变化率和环比变化率的简称&…

手把手教会你--渗透实战--最开始的靶机演示--项目二:Lampiao

有什么问题&#xff0c;请尽情问博主&#xff0c;QQ群796141573 前言。。。知识1.1 开头介绍1.2 22/80端口1.3 nmap1.4 IP地址1.5 开启靶机--网络模式1.6 信息收集(1) 查看源代码(2) 查看常见漏洞(3) 随便点一点\改一改 1.7 robots.txt。。。实战1.1 确定目标(1) 查看自己的ip地…