应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » CPMD/CP2K » cp2k与vasp在分子动力学模拟方面的比较
121/1返回列表
查看: 10893  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

bluejh

金虫 (小有名气)

[交流] cp2k与vasp在分子动力学模拟方面的比较

最近在用vasp做abmd的模拟工作，结果很令人郁闷，15x15x15大小的盒子加入80个水左右，计算速度特别慢，基本没法做。每天也就是迭代200步。基本参数设置如下
IBRION = 0
  NSW = 3000
  POTIM = 1
  SMASS = 2
  TEBEG = 298.0
  TEEND = 298.0

# ISIF = 0
  ISYM = 0
  NELMIN=8
  ALGO =normal
  IALGO=38  
  ISPIN = 2
  LWAVE = .FALSE.
  LCHARG = .FALSE.
  ISMEAR = -1
  SIGMA = 0.05
  LREAL= Auto
  NELECT=613

  APACO=10
  NPACO=200
  NBLOCK = 1
  KBLOCK = 50

想请教用过cp2k的大牛们，cp2k处理类似体系是否能快很多？ 实在不行就换软件啦。。。
回复此楼

» 猜你喜欢

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

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2015-04-14 11:37:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bluejh

金虫 (小有名气)
引用回帖:
2楼: Originally posted by lsloneil at 2015-04-14 11:55:15
AIMD就是很慢，cp2k是专门做这个的，应该会快一点，但像你这样的体系在16个核的机器上cp2k大概一天也只能跑1-2ps， 一般常规的采样也至少要20-30 ps，所以一般要跑一个月，如果没有足够的计算资源还是放弃吧。

To Isloneil，
请问20-30ps是针对CP2K的么？vasp是否可以用10-20ps？从原理上来说vasp是基于BOMD，弛豫过程中似乎结果更接近于基态！
赞一下(1人) 回复此楼
6楼2015-04-15 18:41:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lsloneil

专家顾问 (正式写手)
★ ★
bluejh(金币+1): 谢谢参与
abinitio: 金币+1, 鼓励交流 2015-04-15 07:33:20
AIMD就是很慢，cp2k是专门做这个的，应该会快一点，但像你这样的体系在16个核的机器上cp2k大概一天也只能跑1-2ps， 一般常规的采样也至少要20-30 ps，所以一般要跑一个月，如果没有足够的计算资源还是放弃吧。
赞一下 回复此楼
2楼2015-04-14 11:55:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)
★ ★
bluejh(金币+1): 谢谢参与
abinitio: 金币+1, 鼓励交流 2015-04-15 07:33:28
加快速度建议如下：
ALGO =normal 改为Fast or very Fast
NELMIN = 6
ISPIN = 1 关掉ISPIN，没必要
SIGMA = 0.05 再该大一点 0.1 or 0.2？（这个需要test）
EDIFF = 1E-4再小点 改为 1E-3
另外， 你的NSW=3000 偏小了。
赞一下 回复此楼
3楼2015-04-14 13:54:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluejh

金虫 (小有名气)
引用回帖:
3楼: Originally posted by emilyoyang at 2015-04-14 13:54:53
加快速度建议如下：
ALGO =normal 改为Fast or very Fast
NELMIN = 6
ISPIN = 1 关掉ISPIN，没必要
SIGMA = 0.05 再该大一点 0.1 or 0.2？（这个需要test）
EDIFF = 1E-4再小点 改为 1E-3
另外， 你的NSW=300 ...

谢谢emilyoang，我的体系设置ALGO=48开始时老报错，所以才使用normal先预优化，后面会改用48. EDIFF = 1E-4改为 1E-3会不会有点小么？ 看文献上似乎都是E-4居多
赞一下 回复此楼
5楼2015-04-15 09:14:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
为何不用CPMD,做这种体系的话快得多吧
赞一下 回复此楼
10楼2015-12-03 19:00:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingweilu

新虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
做类似楼主这样的体系，gromacs是强项，4核8线程的电脑一周就能完成几十个ns

发自小木虫Android客户端
赞一下 回复此楼
11楼2017-01-20 13:29:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingweilu

新虫 (著名写手)
cp2k也是基于BOMD的

发自小木虫Android客户端
赞一下 回复此楼
12楼2017-01-20 13:31:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
magicmonk4楼
2015-04-15 09:06   回复  
bluejh(金币+1): 谢谢参与
天涯20127楼
2015-04-18 08:20   回复  
bluejh(金币+1): 谢谢参与
shenrenren8楼
2015-04-18 08:26   回复  
bluejh(金币+1): 谢谢参与
miaojiabing9楼
2015-04-18 08:27   回复  
相关版块跳转 我要订阅楼主 bluejh 的主题更新
121/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭