| 查看: 745 | 回复: 1 | ||
jjkzh金虫 (初入文坛)
|
[求助]
关于CPMD的几个困惑，还请高手前辈解惑
|
|
第一个问题是QUENCH BO 的问题，按照手册上的说法，这样是为了让电子波函数接近绝热表面。我的理解 是这一步等同于波函数优化。但是我不明白的是这个是只有第一步优化那？还是模拟过程里的每一步？我理解应该是第一步。但在其后的模拟里，我怎么知道波函数没有偏离BO表面那？看总能量和电子动能就可以了吗？ 第二个问题还是QUENCH BO 的问题，按照手册上的说法，第一步的所有粒子的速度为零，我有个疑问，如果 这样，假定我的温度是1000K,那么第二步开始，系统要给电子，原子加以动能，满足粒子速率在1000K的分布。 这样系统总体能量不变。但动能会增加。换句话说，第一步系统没有动能，但后面的步骤不再受QUENCH BO的影响，开始给予粒子动能。我的理解对吗？ 第三个问题是如果我想做NVT模拟，怎么样才能知道系统达到平衡态？我的做法是第一步做几何优化(geometry optimization)，第二步是quench bo，如果目标温度是1000K,我会设定TEMPRATURE 1100K TEMPCONTROL IONS 1000K. 第三步是RESTART, 但设定是 TEMPERATURE 1100. TEMPCONTROL IONS 1000. 150. 没有QUENCH BO。 第四步是TEMPERATURE 1000 NOSE IONS 1000 3400 NOSE ELECTRONS 0.0003 10000。 我觉得第三步不是很必要。请大伙看看有没有问题，基本错误？另外，我怎么知道系统绝热很好那？第四步只是看电子动能在小范围震荡就够了吗？我怎么知道电子波函数没有偏离绝热面？？ 第四个问题是什么情况下使用NOSE IONS MASSIVE？这个对到达震动态平衡很有用。但我不是很清楚它对计算资源消耗有多厉害，什么情况必须用？？ 多谢各位。 [ Last edited by jjkzh on 2012-11-26 at 19:14 ] |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
2楼2012-12-11 21:14:44
