24小时热门版块排行榜    

查看: 6135  |  回复: 11

pinwei2014

铜虫 (初入文坛)

[交流] Vasp相比经典分子动力学的优缺点 已有4人参与

刚接触分子动力学,一直有个疑惑,既然Vasp是基于量子力学的分子动力学模拟,那它是不是完胜了经典的分子动力学?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur


小木虫: 金币+0.5, 给个红包,谢谢回帖
经典动力学根本不涉及电子,只是分子骨架的运动。原子间的作用分为两种,键力和非键力。前者根据胡克定律(目前可能不少力场也考虑非简谐项,比如x^3项)定义,把分子拆分成键长振动、键角振动、二面角变化导致的能量变化,不同力场中的不同类型的键、不同原子间的键他们的劲度系数都不同,通过一组经验化的值来描述;后者则是拟合LJ曲线,来描述vdW和静电力。所以绝大多数情况不适合描述成键、断键的过程(反应力场除外),严重牵涉到电子效应的时候肯定解决不了问题。

从头算动力学就要分为CP、BO以及其他的动力学,BOMD和上述经典力学大致操作差不多,但是在计算分子的能量时(此时分子内、分子间的作用都是计算的),每一步产生的构型不是通过力场,而是严格的求解K-S方程,每一步都类似计算单点能那样需要进行scf。计算量可想而知,这样成键、断键就可以描述。但是由于DFT的一些特性,弱作用处理起来不是那么好,好的方法代价又太大,对于那些弱作用主导的大体系(蛋白、胰岛素),经典动力学还是很有必要的。

所以,除了计算每个构形的能量处理可能不一样,对于原子核的处理没什么太大不同,所以那些MSD等性质仍然可以计算,然而一些和电子尺度相关的性质,经典动力学就无能为力了。BOMD建议先做单点能测试看看每个构形计算的时间大概多少,MD准备做多少个step(还需检验平衡),大致乘一下就可以估计计算时间及可行性了。

一点个人观点,望其他老前辈批评指正了。。。
不一定挂在论坛,计算问题问题欢迎留言。
7楼2014-04-03 09:36:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞行鸟

至尊木虫 (文坛精英)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by pinwei2014 at 2014-04-02 14:17:35
vasp是不是只能计算0K的情形,而经典的可以各种算...

当然不是啊,各种压力,温度都可以,电场都行
快乐吧
4楼2014-04-02 15:04:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

飞行鸟

至尊木虫 (文坛精英)


小木虫: 金币+0.5, 给个红包,谢谢回帖
如果你的计算硬件够强大的话,那当然是可以完胜了。 但是如果不是的话,你做能做的体系就不能足够大了,只能使用比较小的,在某些计算上就无法使用完善的模型。
快乐吧
2楼2014-04-02 12:26:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pinwei2014

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by 飞行鸟 at 2014-04-02 12:26:18
如果你的计算硬件够强大的话,那当然是可以完胜了。 但是如果不是的话,你做能做的体系就不能足够大了,只能使用比较小的,在某些计算上就无法使用完善的模型。

vasp是不是只能计算0K的情形,而经典的可以各种算
3楼2014-04-02 14:17:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pinwei2014

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by 飞行鸟 at 2014-04-02 15:04:55
当然不是啊,各种压力,温度都可以,电场都行...

那我想请教一下vasp怎么计算晶体内部原子不同温度下的均方位移?
的确,如果性能足够好的话,能跑更多原子就完胜了。
Vasp相比经典分子动力学的优缺点
1.png

5楼2014-04-02 15:42:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by pinwei2014 at 2014-04-02 15:42:26
那我想请教一下vasp怎么计算晶体内部原子不同温度下的均方位移?
的确,如果性能足够好的话,能跑更多原子就完胜了。

1.png
...

VASP做AIMD的时候会输出各个模拟时间步长的原子坐标,由各个原子的坐标变化情况就可以得到MSD。
6楼2014-04-02 23:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pinwei2014

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by stractor at 2014-04-02 23:03:06
VASP做AIMD的时候会输出各个模拟时间步长的原子坐标,由各个原子的坐标变化情况就可以得到MSD。...

能详细介绍下么,本人小白,不太懂。
8楼2014-04-24 13:37:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pinwei2014

铜虫 (初入文坛)

引用回帖:
6楼: Originally posted by stractor at 2014-04-02 23:03:06
VASP做AIMD的时候会输出各个模拟时间步长的原子坐标,由各个原子的坐标变化情况就可以得到MSD。...

用什么软件来处理计算后的数据呢?
9楼2014-04-24 14:05:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

内容已删除
10楼2014-04-24 15:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pinwei2014 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +4 布布和一二 2026-08-10 4/200 2026-08-10 14:50 by 医学老男孩
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 关于代码变化问题,想知道的进来 +16 且听虎啸 2026-08-07 23/1150 2026-08-10 12:35 by majunge000
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见