关于Δρ计算的一些注意点

一个虫友学作Δρ，选择了HCl来练手，结果作出来了，发现与文献上看到的大不一样，问我为什么。
因为不知道他是如何具体操作的，便拟写了几点让他对照自查，不想写下来还不少，虽然是简单一个HCl例子，但还是涉及到方方面面，麻雀虽小，肝胆具全！
现贴在这里供学作Δρ虫友参考。

据我的经验基本原则是：
1、HCl、H、Cl的坐标都取HCl 中的，并用-1和nosymm将它们固定在某一个坐标平面。
2、对HCl、H、Cl统一自定义网格点数、网格的起始点和步长，不能让程序自选。
3、HCl分子取向的主轴是z还是x？Cl呢？指的是Cl的那个单电子轨道的取向要与HCl一致。Cl的那个单电子轨道的取向可以在multiwfn中看到，如果需要可以在这里可以改变取向。
深入的考虑是：
4、关于Cl的计算，如果用ROHF，单电子将放在孤对电子之上，是HOMO，能级比H还高。如果用HF计算，单电子将放在孤对电子之下，能级比H低得多。Cl用HF计算的总能量也比ROHF的低得多。故Cl应用HF计算的数据。
5、对于单原子Cl的计算，还不能加nosymm，因为加了，不作对称性操作，Cl的单电子取向不在某一个主轴方向，不能与HCl对应。这时是先算Cl不加nosymm，计算自动会将Cl放在坐标原点，单电子取向在主轴z方向。再安排HCl在主轴z方向，Cl放在坐标原点，并在输入坐标数据中在Cl和H之后加-1固定，并加nosymm以免改变了坐标系的主轴取向和挪动了电荷重心。
6、HCl的Δρ，既可作全电子的，又可以只作相关轨道的，如果只看成键，后者更明晰一些。因为HCl的成键，就只是Cl的那个单电子和H之间的，Cl上的其他电子没有参加与H的成键，MO的安排充分体现了这点。这时可取HCl的成键MO与Cl的那个单电子MO、H的即成。如果是作全电子的，Δρ还会包含Cl因吸引了电荷过来，使其上的孤对电子能级被抬升，也会反映到Δρ的图形上。

[ Last edited by yjcmwgk on 2011-2-10 at 10:45 ] 返回小木虫查看更多