| 查看: 9764 | 回复: 13 | ||||||||||||||||
| 本帖产生 2 个 模拟EPI ，点击这里进行查看 | ||||||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||||||||||||
ChemiAndy木虫 (正式写手)
|
[交流]
第一原理动力学的系综选择，和各种系综的输入文件Example 已有9人参与
|
|||||||||||||||
|
做Car-Paranello MD，无论是用VASP, DMOL还是CPMD/cp2k， 都面临着如何选择系综的问题。那该怎么选择呢？我的看法是，应该尽可能使用NVE系综，而不是NVT/NPT。为什么呢？ 因为NVE系综是真正体系能量守恒的。 像使用Nose-Hoover chain 的NVT系综固然也能保证能量守恒,但是它使用一个或多个虚拟的储热罐来缓存由温度波动引起的能量波动，从而模仿一个NVE系综。这种扩展拉格朗日体系方法对体系性质的热力学平均和采样的影响不大，但是，对具体的动力学性质的细节会产生一定影响，因为虚拟的储热罐的动力学行为随着其本身耦合频率，幅度等参数的调整，对体系的演进有着一定影响，具体就是对分子的扩散，转动，化学反应的时间特征影响。另外，任何的拉格朗日虚拟扩展量都会增加实际的运动自由度，所以会略微慢点。 NVE则没有这些缺点，它是完美的隔离体系，完全遵循能量守恒定律的条件，也没有虚拟的动力学变量。 但NVE有个缺点妨害了它的使用。什么缺点呢？就是的它的能量守恒于起始能量。即守恒能量等于起始构型的势能+起始设定的动能。很可能就是你前面预平衡阶段最后一帧的速度和构型。NVE开始后，能量在动能、势能之间波动式分配。如果一开始势能高，能量就向动能转化，平衡后温度就高；反之，动能向势能转化，平衡后温度就低。显然，这更接近于真实体系的动能-势能的转化模式，但是这就使得其最终的温度其平衡之前很难预料，你本来想模拟300K，结果最终温度却稳定在了250K, 或者400K！你怎么办？只好使用NVT, NPT去控制温度。 从上面分析看，使用NVT/NPT是不得以而强为之。还是NVE好。 NVE的温度问题也不是不能解决。一般用Velocity scaling积分算法跑一个短时间的NVT模拟，比如100fs，取最后一帧Restart一个NVE，一般都能把平衡后的温度控制在一定的范围内。 附， CPMD各种系综的设定。 1. NVE呢，不指定任何控温关键词的，又没有设定任何控压关键词的，即为NVE。For example: &CPMD MOLECULAR DYNAMICS CP RESTART WAVEFUNCTION COORDINATES VELOCITIES CELL LATEST EMASS 500 TIMESTEP // 0.096fs/步 4 MAXSTEP 50000 &END 如果不想使用前面模拟最后一帧的速度，可以加关键词： TEMPERATURE 300 它的意思是在NVE的开始前，重新为每个原子/动力学变量设定动能，使得总的温度为300K。NVE开始后，并不控制体系守在这个温度。 2. NVT 有多种多种控制温度的关键词。 (1) Velocity Scaling控温可以帮助体系快速预平衡到设定温度/势能附近，但是不满足能量守恒，一般不用于采样阶段的模拟。关键词： TEMPCONTROL &CPMD MOLECULAR DYNAMICS CP RESTART WAVEFUNCTION COORDINATES VELOCITIES CELL LATEST TEMPCONTROL IONS, ELECTRONS 300 0.006 EMASS 500 TIMESTEP 4 MAXSTEP 50000 &END 其中，IONS 300设定原子的温度为300K; ELECTRON 0.006(au), 是设定电子，即波函数作为运动量(或叫扩展拉格朗日变量)的温度，需预先跑个100步的NVE，取最后几步的EKINC即电子动能，作为这里的控温参数。可比之稍大一点。 (2) NOSE-HOOVER Chain控温是遵循能量守恒的NVT系综，或者说它通过增加变量的方法让体系模拟一个隔离的NVE体系。第一原理模拟中广泛使用。 &CPMD MOLECULAR DYNAMICS CP RESTART WAVEFUNCTION COORDINATES VELOCITIES CELL ACCUMULATORS NOSEE NOSEP LATEST // 重接上次模拟继续 TIMESTEP 4 MAXSTEP 10000 TRAJECTORY XYZ SAMPLE //每隔10步采个样存到TRAJECTORY中 10 DIPOLE DYNAMICS WANNIER SAMPLE //每隔10步存储Dipole moment到DIPOLE文件中 10 STORE //每隔500步写一次RESTART.1文件 500 EMASS 500 NOSE IONS // NOSE控制原子温度和频率(单位波数，即cm-1) 300 1600 NOSE ELECTRONS // NOSE控制电子动能和频率(单位波数，即cm-1) 0.020 25000 &END (3) NPT的设定：NPT的目的是最终使得盒子在每个方向上的受力与外部的压力相等。压力控制要在&CPMD部分使用PARRINELLO-RAHMAN NPT关键词，并在&SYSTEM中设定PRESSURE的大小。此外，还有一些其它的关键词控制压力的各向异性，请查阅手册。 &CPMD MOLECULAR DYNAMICS CP RESTART WAVEFUNCTION COORDINATES VELOCITIES CELL ACCUMULATORS NOSEE NOSEP LATEST // 重接上次模拟继续 TIMESTEP 4 MAXSTEP 10000 TRAJECTORY XYZ SAMPLE //每隔10步采个样存到TRAJECTORY中 10 DIPOLE DYNAMICS WANNIER SAMPLE //每隔10步存储Dipole moment到DIPOLE文件中 10 STORE //每隔500步写一次RESTART.1文件 500 EMASS 500 NOSE IONS // NOSE控制原子温度和频率(单位波数，即cm-1) 300 1600 NOSE ELECTRONS // NOSE控制电子动能和频率(单位波数，即cm-1) 0.020 25000 PARRINELLO-RAHMAN NPT //设定NPT系综和控压方法 STRESS TENSOR //每隔100步输出一次压力值 100 &END &SYSTEM ANGSTROM PRESSURE // 设定外部压力为100Kbar 100 CELL VECTORS // NPT模拟需要使用CELL VECTOR来设置盒子 13.080 0.000 0.000 0.000 11.050 0.000 -9.805 0.000 14.374 CUTOFF 100.0 &END [ Last edited by ChemiAndy on 2012-2-15 at 05:23 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 | 仿真建模与计算 | 分子动力学 | 晶体模型计算分析 |
| 第一性原理计算-tg | 基础知识 | 关于Lammps | 我学习计算的一些帖子 |
| 量化,第一性,MD笔记 | 动力学或热力学计算 | 量子化学计算与使用经验 | 分子模拟 |
| vasp 学习贴 |
» 本帖已获得的红花（最新10朵）
qphll» 猜你喜欢
化学工程085602 305分求调剂
已经有10人回复
-
289求调剂
已经有15人回复
-
291 求调剂
已经有7人回复
-
274求调剂
已经有14人回复
-
北京林业大学硕导招生广告
已经有3人回复
-
309求调剂
已经有5人回复
-
292求调剂
已经有8人回复
-
求调剂
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有3人回复
-
调剂310
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- MD系综选择的问题？
已经有5人回复
- nvt下模拟之后获得压力的问题
已经有5人回复
- 分子动力学能量守恒的相关问题
已经有18人回复
- 【英美经典教材】《园艺学的原理和实践》(英文彩色高清)【已搜索，无重复】
已经有198人回复
- 自学求助输入文件，谢谢
已经有8人回复
- 关于MS动力学精度设置问题
已经有4人回复
- LAMMPS分子动力学模拟学习用书推荐
已经有6人回复
- 复旦大学quantum espresso培训及其它
已经有33人回复
- NPT系综下判断是否平衡的标志
已经有12人回复
- K-point的选择问题
已经有8人回复
- MS 系综平均 分子动力学
已经有5人回复
- NPT系综动力学
已经有19人回复
- 化学反应能量真的守恒吗
已经有42人回复
- 分子动力学模拟中的系综问题
已经有5人回复
- 模拟轧制，过程简单，但是结果一直不对，问了很多人也没找到原因，求助
已经有9人回复
- 求助：谁有计算双原子分子能级跃迁Franck-condon因子的程序，最好是fortran的，谢谢！
已经有8人回复
- 怎么把Discover动力学完成后的轨迹文件提取出来？
已经有10人回复
- 做材料力学性能研究是用vasp好还是ab initio 好啊？！
已经有20人回复
- NPT系综的平衡问题
已经有10人回复
- 【讨论】系综问题
已经有32人回复
- 【求助】请教分子动力学模拟
已经有8人回复
yoghurt117
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 8018.3
- 散金: 56
- 红花: 8
- 帖子: 576
- 在线: 134.9小时
- 虫号: 761174
- 注册: 2009-05-01
- 性别: GG
- 专业: 理论和计算化学
7楼2012-03-05 10:17:27
safina.din
金虫 (初入文坛)
- 应助: 7 (幼儿园)
- 金币: 1021.9
- 红花: 1
- 帖子: 31
- 在线: 534.2小时
- 虫号: 740379
- 注册: 2009-04-04
- 专业: 理论和计算化学
2楼2012-02-16 02:54:06
meixinqing
金虫 (著名写手)
梅心晴｜一只蘑菇
- 应助: 46 (小学生)
- 贵宾: 0.122
- 金币: 5691.4
- 散金: 423
- 红花: 35
- 帖子: 1923
- 在线: 369.4小时
- 虫号: 280674
- 注册: 2006-09-17
- 性别: MM
- 专业: 实验地球化学和计算地球化
- 管辖: 签证指南
4楼2012-02-16 13:36:04
ChemiAndy
木虫 (正式写手)
- 模拟EPI: 40
- 应助: 52 (初中生)
- 金币: 2573.3
- 散金: 18
- 红花: 104
- 帖子: 590
- 在线: 553.8小时
- 虫号: 836907
- 注册: 2009-08-30
- 性别: GG
- 专业: 理论和计算化学
御剑江湖(金币+5): 谢谢 2012-02-17 21:42:18
|
手册建议是大于3倍,以免它们之间形成耦合. 这里选择NOSE ELECTRONS的耦合频率25000cm-1，大概是每1个fs就让电子储热罐与体系交换一次能量，即耦合一次。这个频率的确很高，目的是希望电子动能的波动小一些。不过该值的选取对模拟结果的影响不大。使用5000-30000都是允许的。5000基本上对应大约每6fs耦合一次。 选取小的值，计算会略微快一点点，但是容易出问题。电子动能可能在耦合前就波动到灾难性的结果。因此，如果碰到NaN了，可把NOSE ELECTRONS调大一些试试。 [ Last edited by ChemiAndy on 2012-2-16 at 07:10 ] |
5楼2012-02-17 03:07:12
