24小时热门版块排行榜    

查看: 5774  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 294求调剂材料与化工专硕 +8 陌の森林 2026-03-18 8/400 2026-03-19 10:26 by 是鸭鸭呀呀
[考研] 085600材料与化工调剂 324分 +9 llllkkkhh 2026-03-18 9/450 2026-03-19 07:40 by BruceLiu320
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 10/500 2026-03-18 20:41 by Wangjingyue
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 0703化学调剂 +4 pupcoco 2026-03-17 7/350 2026-03-18 12:14 by djl2006
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
信息提示
请填处理意见