| 查看: 1599 | 回复: 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 ] |
回复此楼» 猜你喜欢
291求调剂
已经有11人回复
-
304求调剂
已经有8人回复
-
求调剂推荐
已经有9人回复
-
B区0809 ，数一英一，290 求调剂
已经有5人回复
-
求调剂
已经有12人回复
-
291求调剂
已经有13人回复
-
085404 22408 309分求调剂
已经有11人回复
-
求调剂
已经有7人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有6人回复
-
307中医考研调剂
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算过渡态问题
已经有13人回复
- 巨正则系综蒙特卡洛方法是怎样计算吸附过程的
已经有14人回复
- vasp算法的一些基础问题
已经有10人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- vasp优化结构中出现的问题
已经有10人回复
- VASP优化结构遇到的问题。
已经有9人回复
- 【求助】VASP收敛问题
已经有15人回复
- 【求助】vasp分子动力学模拟求助
已经有15人回复
- 【原创】vasp的分子动力学模拟
已经有44人回复
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 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
