计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9053|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] MD模拟过程疑问

[复制链接 Copy URL]

113

帖子

0

威望

903

eV
积分
1016

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
我在看文献的时候,看到一些MD模拟,开始时是在NPT系综跑2ns,接着继续在NVT系综跑。我想请教为什么要跑NPT系综的2ns?不能一直跑NVT么?我是MD小白,没做过MD的

294

帖子

0

威望

3528

eV
积分
3822

Level 5 (御坂)

2#
发表于 Post on 2016-5-10 12:36:14 | 只看该作者 Only view this author
简单点说就是  在NPT系综下进行体系的平衡,在NVT系综下进行数据收集和分析。  有些情况下,是先进行能量最小化(EM),然后直接进行NVT模拟,或者是EM,然后在NPT系综下进行数据的收集。

一般的套路:建模-EM-NPT平衡-NVT数据收集分析

有些情况比较特殊就会出现:建模-EM-NPT模拟 或者 建模-EM-NVT模拟

具体的操作方式和你所研究的体系相关,一般文献在simulation methods部分会写出来,多看点文献就会了解一点!

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sobereva + 2

查看全部评分 View all ratings

336

帖子

2

威望

4531

eV
积分
4907

Level 6 (一方通行)

3#
发表于 Post on 2016-5-10 14:10:33 | 只看该作者 Only view this author
ruanyang 发表于 2016-5-10 12:36
简单点说就是  在NPT系综下进行体系的平衡,在NVT系综下进行数据收集和分析。  有些情况下,是先进行能量最 ...

什么体系在NVT下收集数据?

294

帖子

0

威望

3528

eV
积分
3822

Level 5 (御坂)

4#
发表于 Post on 2016-5-10 14:25:22 | 只看该作者 Only view this author
很多啊,用gromacs的都是这样

294

帖子

0

威望

3528

eV
积分
3822

Level 5 (御坂)

5#
发表于 Post on 2016-5-10 14:25:35 | 只看该作者 Only view this author
也不能说都是,只是大部分

1096

帖子

6

威望

6271

eV
积分
7487

Level 6 (一方通行)

6#
发表于 Post on 2016-5-10 22:15:08 | 只看该作者 Only view this author
用npt系综先跑是因为要让体系的压力合理,自建的体系压力很难恰好是1atm

用nvt续跑大多数情况下没什么优势,一般来说直接用npt系综跑就行了,特殊情况大多是和边界有关。比如无限长的碳纳米管,在周期性盒子之间是有化学键相连的,这个时候通常我们不希望盒子的大小(尤其是有化学键的方向)变化,一般就是先在np(z)t系综下面把压力跑合理,再在nvt系综下进行平衡模拟

评分 Rate

参与人数
Participants 2
eV +4 收起 理由
Reason
stecue + 3 赞!
sobereva + 1

查看全部评分 View all ratings

113

帖子

0

威望

903

eV
积分
1016

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2016-5-11 08:30:03 | 只看该作者 Only view this author
ruanyang 发表于 2016-5-10 12:36
简单点说就是  在NPT系综下进行体系的平衡,在NVT系综下进行数据收集和分析。  有些情况下,是先进行能量最 ...

好的,谢谢你的回答

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 18:16 , Processed in 1.509211 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list