| 查看: 763 | 回复: 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 ] |
回复此楼» 猜你喜欢
大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。
已经有11人回复
-
两块石头
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
2楼2012-12-11 21:14:44
100