| 查看: 1816 | 回复: 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 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
