请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 525|回复: 18

[综合交流] 求助:关于CPMD与AIMD的区别

[复制链接]

11

帖子

0

威望

77

eV
积分
88

Level 2 能力者

发表于 2019-6-22 17:36:44 | 显示全部楼层 |阅读模式
各位老师好,本人在看到一篇文献时,里面在计算气相反应自由能前做了一个PMF采样,同时提到了做AIMD和CPMD模拟,我想问一下各位老师这两个模拟是什么关系?

171

帖子

0

威望

2640

eV
积分
2811

Level 5 (御坂)

跳跳猪

发表于 2019-6-22 21:00:11 | 显示全部楼层
MD计算中的能量和梯度直接通过ab-initio的方式计算的话都可以称为ab-initio molecular dynamics,CPMD也属于AIMD范畴。

跟CPMD对应的是BOMD,一下截取自wikipedia, https://en.wikipedia.org/wiki/Car%E2%80%93Parrinello_molecular_dynamics
The CPMD method is related to the more common Born–Oppenheimer molecular dynamics (BOMD) method in that the quantum mechanical effect of the electrons is included in the calculation of energy and forces for the classical motion of the nuclei. However, whereas BOMD treats the electronic structure problem within the time-independent Schrödinger equation, CPMD explicitly includes the electrons as active degrees of freedom, via (fictitious) dynamical variables.


更详细的可以看一下"Ab Initio Molecular Dynamics: Theory and Implementation", in "Modern Methods and Algorithms of Quantum Chemistry" (p. 301-449), Editor: J. Grotendorst, (NIC, FZ Jülich 2000).  https://juser.fz-juelich.de/record/44687/files/NIC-Band-1.pdf



流年似水,浮生如梦。

2万

帖子

25

威望

2万

eV
积分
46654

管理员

公社社长

发表于 2019-6-23 00:10:24 | 显示全部楼层
AIMD分为BOMD(最主流的实现)、CPMD(相对小众)、ADMP(仅Gaussian支持)
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

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

11

帖子

0

威望

77

eV
积分
88

Level 2 能力者

 楼主| 发表于 2019-6-24 09:05:45 | 显示全部楼层
sobereva 发表于 2019-6-23 00:10
AIMD分为BOMD(最主流的实现)、CPMD(相对小众)、ADMP(仅Gaussian支持)

多谢老师解答

11

帖子

0

威望

77

eV
积分
88

Level 2 能力者

 楼主| 发表于 2019-6-24 09:06:18 | 显示全部楼层
Mikasa 发表于 2019-6-22 21:00
MD计算中的能量和梯度直接通过ab-initio的方式计算的话都可以称为ab-initio molecular dynamics,CPMD也属 ...

多谢老师的帮助

24

帖子

0

威望

111

eV
积分
135

Level 2 能力者

发表于 2019-9-30 15:02:31 | 显示全部楼层
sobereva 发表于 2019-6-23 00:10
AIMD分为BOMD(最主流的实现)、CPMD(相对小众)、ADMP(仅Gaussian支持)

sobereva老师,请问用AIMD跑出来的反应可以认为是近乎无垒的反应吗?

2万

帖子

25

威望

2万

eV
积分
46654

管理员

公社社长

发表于 2019-10-1 07:43:54 | 显示全部楼层
ZXW1993 发表于 2019-9-30 15:02
sobereva老师,请问用AIMD跑出来的反应可以认为是近乎无垒的反应吗?

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

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

24

帖子

0

威望

111

eV
积分
135

Level 2 能力者

发表于 5 天前 | 显示全部楼层
sobereva 发表于 2019-10-1 07:43
否。
有没有势垒的反应都可以做AIMD

那AIMD可以跑出任意能垒的反应吗?还是说接近无垒的反应才能在较短的时间内跑出来?

565

帖子

0

威望

1018

eV
积分
1583

Level 5 (御坂)

发表于 5 天前 来自手机 | 显示全部楼层
ZXW1993 发表于 2019-10-9 09:50
那AIMD可以跑出任意能垒的反应吗?还是说接近无垒的反应才能在较短的时间内跑出来?

取决于你输入的能量

24

帖子

0

威望

111

eV
积分
135

Level 2 能力者

发表于 5 天前 | 显示全部楼层
granvia 发表于 2019-10-9 12:10
取决于你输入的能量

做AIMD的时候如何输入能量呢?是初始构型的能量吗?如果是初始构型的能量的话,需要通过什么方法降低初始构型的能量呢?

99

帖子

1

威望

313

eV
积分
432

Level 3 能力者

发表于 5 天前 | 显示全部楼层
能量就是溫度

565

帖子

0

威望

1018

eV
积分
1583

Level 5 (御坂)

发表于 5 天前 来自手机 | 显示全部楼层
ZXW1993 发表于 2019-10-9 14:06
做AIMD的时候如何输入能量呢?是初始构型的能量吗?如果是初始构型的能量的话,需要通过什么方法降低初始 ...

能量包括振动能与势能。一般通过分配振动能来设置起始能量,势能可设为零(即起始几何构型为分子平衡构型)。振动能的分配是通过随机分配每个原子起始速度来实现(限制条件是体系整体无平动和转动速率),有的MD程序具有自动分配功能(通过指定温度来实现)。当然也可以自写脚本实现。

利用上述方法单独对各反应物分子进行热力学平衡(需要各自设置起始温度),然后各取最后一帧放在一起模拟反应,此时可以额外加上各反应物间的相对碰撞速率。

24

帖子

0

威望

111

eV
积分
135

Level 2 能力者

发表于 5 天前 | 显示全部楼层
granvia 发表于 2019-10-9 17:23
能量包括振动能与势能。一般通过分配振动能来设置起始能量,势能可设为零(即起始几何构型为分子平衡构型 ...

哦哦,我大概明白您的意思了。我在300K温度下对反应物进行模拟,如果有产物生成,这样可以说明这个反应接近无垒吗?

24

帖子

0

威望

111

eV
积分
135

Level 2 能力者

发表于 5 天前 | 显示全部楼层

哦哦,如果我在300K温度下对反应物进行模拟,有产物生成,能够说明该反应接近无垒吗?

565

帖子

0

威望

1018

eV
积分
1583

Level 5 (御坂)

发表于 5 天前 来自手机 | 显示全部楼层
ZXW1993 发表于 2019-10-9 19:16
哦哦,我大概明白您的意思了。我在300K温度下对反应物进行模拟,如果有产物生成,这样可以说明这个反应接 ...

Nope
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2019-10-14 22:23 , Processed in 0.156699 second(s), 24 queries .

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