请选择 进入手机版 | 继续访问电脑版
本论坛在12月4日~12日期间暂停发新帖功能。给带来的不便敬请谅解!

计算化学公社

 找回密码
 现在注册!
查看: 1049|回复: 23

[综合交流] 计算分子内质子转移的轨迹,该用什么软件呢?大家来给我的想法提提意见吧

[复制链接]

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

发表于 2018-8-8 07:43:55 | 显示全部楼层 |阅读模式
苯环上有一个氨基(-NH3)、一个羧基(-COOH),质子在这两个基团之间转移,如果要计算这个质子的运动轨迹,有没有哪款免费软件能算全过程呢?
感觉已有的各种分子模拟软件都是针对周期性体系的,那我这个不是周期性的,就是单个分子(目前不考虑溶剂),不知能不能用那些软件,如果能用,用哪款最适合?并且,该选哪种算法,velocity verlet? Leapfrog? 还是别的?
如果真的没有一款软件能算这种东西,需要我手动一步一步的计算的话,大家看看我这个方案可行不可行:假设每一步的步长是1飞秒(这个是不是太长了呢?),我每一步用Gaussian里的DFT方法(DFT算这个是不是足够了?)算质子受力,然后根据这个力把质子移动到下一飞秒它该在的位置,然后再进行DFT计算,这样一步一步的算,直到质子最终和另一基团成键。
但是我对这个一步一步手算的方案有一个很大的疑问:每一步DFT计算时,Gaussian算出的都是能量最低的电子结构,但是在现实中,这个质子转移的过程中,不一定总是能量最低的结构(在几个飞秒内它不能迅速转换电子结构以达到能量最低),那我这样算出来的质子运动轨迹,是不是就和现实有很大偏差了?这种情况有解决办法吗?如果用分子模拟软件,会不会也出现这种情况?
谢谢大家花时间阅读我这么多问题!

84

帖子

0

威望

812

eV
积分
896

Level 4 (黑子)

无形之刃,最为致命

发表于 2018-8-8 08:30:19 | 显示全部楼层
个人感觉,既然说到分子受力,为何还要用DFT?分子力学不行吗?

评分

参与人数 1eV +1 收起 理由
算一算 + 1

查看全部评分

1万

帖子

25

威望

1万

eV
积分
36194

管理员

公社社长

发表于 2018-8-8 09:13:36 | 显示全部楼层
直接在Gaussian中通过ADMP或BOMD关键词做从头算动力学模拟即可

评分

参与人数 1eV +3 收起 理由
算一算 + 3

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

2330

帖子

23

威望

7813

eV
积分
10603

Level 6 (一方通行)

发表于 2018-8-8 15:56:09 | 显示全部楼层
你想通过“质子转移的轨迹”说明什么问题?

要用“分子模拟”的软件做,分子模拟即使用周期性盒子也可以做“孤立”的分子(稀溶液),盒子足够大就可以了,这不是问题。问题是大部分力场不会断键的。

简单的算受力、按受力走的方法显然有问题,因为如果从极小点出发、受力是0,哪儿都走不到。简单的从头算动力学的问题也不是“电子结构是能量最低的”,因为电子的弛豫快,往往认为电子结构处在基态造成的问题不太大。

真去做、你会发现你得到的不是质子转移的轨迹,大部分都是几何扭转之类的低能过程,气相的质子转移(其实应该是氢自由基转移)还不是什么特别优先的过程,所以还得等一堆低能的其他反应。不现实。

以及你确定你要的不是IRC而必须是动力学?

评分

参与人数 1eV +3 收起 理由
算一算 + 3

查看全部评分

34

帖子

0

威望

293

eV
积分
327

Level 3 能力者

发表于 2018-8-8 18:03:33 | 显示全部楼层
IRC分析就可以了

评分

参与人数 1eV +1 收起 理由
算一算 + 1

查看全部评分

166

帖子

2

威望

1437

eV
积分
1643

Level 5 (御坂)

发表于 2018-8-9 06:05:31 | 显示全部楼层
liyuanhe211 发表于 2018-8-8 15:56
你想通过“质子转移的轨迹”说明什么问题?

要用“分子模拟”的软件做,分子模拟即使用周期性盒子也可以 ...

想用MD研究这个过程的话可以从质子转移的过渡态开始吧
虽然不知道这样做能得到什么额外的信息

评分

参与人数 1eV +1 收起 理由
算一算 + 1

查看全部评分

1万

帖子

25

威望

1万

eV
积分
36194

管理员

公社社长

发表于 2018-8-9 06:54:33 | 显示全部楼层
alonewolfyang 发表于 2018-8-8 08:30
个人感觉,既然说到分子受力,为何还要用DFT?分子力学不行吗?

一般的分子力场(非反应力场)没法描述诸如质子转移这样有键的形成和断裂的过程
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

103

帖子

0

威望

899

eV
积分
1002

Level 4 (黑子)

发表于 2018-8-10 05:51:25 | 显示全部楼层
你的描述让我想起来14年的一篇Nature Chemistry用TeraChem这个软件来模拟原始大气产生氨基酸的。它的SI里面好像是有一些质子迁移的东西。他们用的方法好像是周期性地压缩体系的体积,得到短时间的高温,再解除对体积的限制,让温度降下来,然后再循环。但按照他们文中的描述的话,数据处理时有点麻烦,他们自己好像是写了一堆python脚本才完成的。
我的感觉是,如果想从过渡态开始研究质子迁移的轨迹,就得保证找到了质子迁移的决速步的过渡态,而且其能垒不能太低,然后才可以在过渡态附近采样,再往反应物、产物两个方向计算轨迹。这方面就得参考Texas Tech University的Hase、Texas A&M University的Singleton、大连化物所张东辉院士等人的工作。
如果你想从反应物开始研究的话,那么就得参考Emory University的Bowman的一些工作。只不过这样的话,最后要看的轨迹的数目有点小多。

评分

参与人数 2eV +6 收起 理由
算一算 + 4 很有帮助,谢谢
sobereva + 2

查看全部评分

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

 楼主| 发表于 2018-8-10 18:32:50 | 显示全部楼层
alonewolfyang 发表于 2018-8-8 08:30
个人感觉,既然说到分子受力,为何还要用DFT?分子力学不行吗?

是这样的,这个过程质子受分子各部位的电子密度的影响很大,所以不能用简单的分子力学来描述。也正如7楼所说,想描述断键过程的话,用DFT免不了的。

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

 楼主| 发表于 2018-8-10 18:37:27 | 显示全部楼层
sobereva 发表于 2018-8-8 09:13
直接在Gaussian中通过ADMP或BOMD关键词做从头算动力学模拟即可

谢谢这个建议。我去Gaussian官网看了一下BOMD,没太明白它这个步长是怎么弄的……它原话是这么写的:
StepSize=n
Sets the dynamic step size to n*0.0001, in the appropriate units. The default step size is 0.25 amu1/2*Bohr except for GradientOnly calculations where it is 0.025 femtoseconds.
请问这个amu开平方再乘以Bohr是描述什么物理量的?如果我选后面那种,就是gradientonly的话,我想把步长设为1飞秒,那我是不是应该写stepsize=10000?

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

 楼主| 发表于 2018-8-10 18:47:33 | 显示全部楼层
liyuanhe211 发表于 2018-8-8 15:56
你想通过“质子转移的轨迹”说明什么问题?

要用“分子模拟”的软件做,分子模拟即使用周期性盒子也可以 ...

是我们头儿丢给我的任务,他就是要看质子运动的路线,让我把每一步的坐标都搞出来。我暂时也不知道他到底想用这个说明什么问题。
你说的从能量极小的几何结构出发受力是0,这个我清楚,我会选择受力不是0的出发点的。
你说的“简单的从头算动力学的问题也不是“电子结构是能量最低的”,” 那它是怎么选择电子结构的呢?我以为软件在每一步都是会在当时的几何结构上优化出能量最低的波函数,如果不是那样的话,软件是怎么做的,它是根据前一步的波函数做初猜?然后尽量让电子结构保持在同一个势能面上是吗?

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

 楼主| 发表于 2018-8-10 18:52:23 | 显示全部楼层
laomalao 发表于 2018-8-8 18:03
IRC分析就可以了

你和楼上都提到IRC了,但是我跟我们头儿商量了一下,他坚持让我做动力学分析,然后把质子运动路线给他看。虽然他也不清楚具体该怎么做。

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

 楼主| 发表于 2018-8-10 18:53:39 | 显示全部楼层
Accelerator 发表于 2018-8-9 06:05
想用MD研究这个过程的话可以从质子转移的过渡态开始吧
虽然不知道这样做能得到什么额外的信息

嗯嗯,我也在考虑从过渡态开始呢

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

 楼主| 发表于 2018-8-10 18:55:00 | 显示全部楼层
Daniel_Arndt 发表于 2018-8-10 05:51
你的描述让我想起来14年的一篇Nature Chemistry用TeraChem这个软件来模拟原始大气产生氨基酸的。它的SI里面 ...

这个信息量好大,好的,我一个一个找找看,我不怕数据多,就怕自己方法用的不好。

61

帖子

0

威望

162

eV
积分
223

Level 3 能力者

 楼主| 发表于 2018-8-10 18:59:16 | 显示全部楼层
谢谢大家热情回复,刚发现每天的评分是有配额的,那我就明天继续弄吧。不过这个评分也不足以表达我的感激。你们提供的各种信息都对我很有帮助,真心谢谢
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-12-10 15:42 , Processed in 0.121548 second(s), 25 queries .

快速回复 返回顶部 返回列表