24小时热门版块排行榜

返回列表

jjkzh

金虫 (初入文坛)

应助: 1 (幼儿园)
金币: 1050.4
散金: 120
帖子: 43
在线: 129.4小时
虫号: 1628245
注册: 2012-02-20
性别: GG
专业: 理论和计算化学

[求助] 关于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 ]

回复此楼