| 查看: 1408 | 回复: 4 | |||
stractor金虫 (著名写手)
|
[交流]
【求助】VASP正则系综模拟问题
|
|
我从来没做过第一性原理分子动力学模拟。最近想在这方面尝试一下。我模拟的是NVT系综。现在对SMASS这个关键字的设置不是很清楚。 手册上说了，对于SMASS>=0，就是选择正则系综，而且模拟采用Nose质量控制温度波动。SMASS的值就对应Nose质量。 手册上这样描述的：The Nose massshould be set so that the induced temperature fluctuation show approximately the same frequencies as the typical phono frequencies …… The frequency of the approximate temperature fluctuations induced by the Nose thermostat is written to the OUTCAR file. 我的问题是，这个质量应该设置多少？ 是不是需要尝试以些SMASS的值，再在OUTCAR中找温度波动的频率？ OUTCAR中这个波动频率在什么位置？ [ Last edited by stractor on 2010-9-2 at 14:33 ] |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有13人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
vasp计算过渡态问题
已经有13人回复- 巨正则系综蒙特卡洛方法是怎样计算吸附过程的
已经有14人回复
- vasp算法的一些基础问题
已经有10人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- vasp优化结构中出现的问题
已经有10人回复
- VASP优化结构遇到的问题。
已经有9人回复
- 【求助】VASP收敛问题
已经有15人回复
- 【求助】vasp分子动力学模拟求助
已经有15人回复
- 【原创】vasp的分子动力学模拟
已经有44人回复
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
2楼2010-09-02 14:29:53
freer007
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 1 (幼儿园)
- 金币: 2421.7
- 红花: 1
- 帖子: 676
- 在线: 146小时
- 虫号: 707811
- 注册: 2009-02-24
- 性别: GG
- 专业: 半导体材料
3楼2010-09-02 15:27:49
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2010-09-03 08:11:42
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
★ ★ ★
aylayl08(金币+3):谢谢分享 2010-09-08 08:38:29
aylayl08(金币+3):谢谢分享 2010-09-08 08:38:29
|
我发邮件问了我所看的一篇做BOMD模拟的老外，很快他就回答了我的问题， 下面是他的回答： The nose mass is related to the average vibrational frequency in the system (at least the way that VASP implements it). The reason for this is you are trying to make the 'frequency' of the thermostat match that of the system you are simulating. If the frequency is too large, you will get unphysical excitations, if it is too low your system will take a very long time to equilibrate. Unfortunately, in the case of something like Li4BN3H10, there is a wide range of vibrational frequencies (from the heavy nitrogen to the very light hydrogen). So, what I ended up doing is taking a system average vibrational frequency based on some earlier phonon calculations (I forget what this frequency was, though) then doing a series of shortruns with VASP to see what the nose frequency was based on its output (it prints out the nose frequency). Since the mass is related to the frequency by f ~ sqrt(1/m), I was able to zero in on a 'close enough' value (it ended up being 5.0 for my calculations). You can probably also go into the source code to see how they set up the default SMASS to make this process even easier. Hope that helps, |
5楼2010-09-08 08:21:48