求助：如何解释吸附能为正值

wzkchem5

阿瞳 发表于 2021-6-24 14:38
老师，我已经懂这个原理了，但是具体操作还是很模糊，我用Multiwfn的0功能输出了所有轨道的占据数，上面 ...

要把前线轨道的图作出来，从轨道图判断吸附前的哪个轨道对应吸附后的哪个轨道。从文字上是看不出来的

wzkchem5

阿瞳 发表于 2021-6-24 15:28
老师，您说的是这种图吗，但我发现一个问题，甲烷分子和金属Pd之间没有电子云，这是怎么回事

对。
甲烷和Pd之间没有轨道分布是正常现象，因为本来就应该是比较弱的吸附。

阿瞳

wzkchem5 发表于 2021-6-25 00:45
对。
甲烷和Pd之间没有轨道分布是正常现象，因为本来就应该是比较弱的吸附。

老师，从这种前线轨道图还是对比不出来啊，具体应该操作啊，我用Multiwfn去看每个轨道对应的图，但是发现吸附前后差别好大，我看论坛里有人可以看出电子排布式这是怎么看的

wzkchem5

阿瞳 发表于 2021-6-25 02:19
老师，从这种前线轨道图还是对比不出来啊，具体应该操作啊，我用Multiwfn去看每个轨道对应的图，但是发现 ...

这样的话，说明吸附前后的波函数差别太大了，轨道不仅占据数发生了很大变化，形状也发生了很大变化。这时用片段组合波函数的方法可能会方便一点

阿瞳

wzkchem5 发表于 2021-6-25 14:02
这样的话，说明吸附前后的波函数差别太大了，轨道不仅占据数发生了很大变化，形状也发生了很大变化。这时 ...

老师，我把吸附前后的fchk文件上传下，您能帮我看看是不是有其他问题吗
链接：https://pan.baidu.com/s/1syQZKBoIWEG6UatS91Ke4Q

wzkchem5

阿瞳 发表于 2021-6-25 07:40
老师，我把吸附前后的fchk文件上传下，您能帮我看看是不是有其他问题吗
链接：https://pan.baidu.com/s/ ...

不存在问题不问题，两个波函数的轨道不能一一对应是常事，不代表你计算有问题。可能也怪我，我之前的措辞可能显得好像两个波函数的轨道理应一一对应。
这时通过简单设定占据数是不行的，或者至少成功率低，可以改用片段组合波函数方法，我已经说过了。

阿瞳

wzkchem5 发表于 2021-6-25 18:44
不存在问题不问题，两个波函数的轨道不能一一对应是常事，不代表你计算有问题。可能也怪我，我之前的措辞 ...

老师，您好，按您说的采用片段组合波函数的方法去计算得到的结果还是正值，于是我索性把这个过渡金属Pd原子去掉，直接计算甲烷分子与碳纳米管之间的相互作用能，发现用甲烷分子和碳纳米管的复合物能量减去甲烷分子和碳纳米管的电子能量就是正的，用之前含Pd体系算出来的吸附能减去现在算甲烷与碳纳米管的相互作用能即为正常的吸附能，但还是想问下复合物的能量为什么会大于两个片段能量之和

wzkchem5

阿瞳 发表于 2021-6-29 02:33
老师，您好，按您说的采用片段组合波函数的方法去计算得到的结果还是正值，于是我索性把这个过渡金属Pd原 ...

你的吸附了甲烷的结构是怎么得到的？结构初猜是怎么摆的？
如果初始结构时，甲烷离碳纳米管很远，那么理应优化得到相互作用能为负的结构

wzkchem5

阿瞳 发表于 2021-6-29 08:33
初猜就是如图这样摆的，直接把一个甲烷分子放进含过渡金属的管腔内，然后进行优化

有没有考虑过可能这个碳纳米管太细了，甲烷塞进去能量是升高的？
你试试把甲烷放在管口，但是离管口至少3个埃，再优化。你应该能得到吸附能为负的结构，但是甲烷可能仍然在管口，不在碳管里面。
当碳管很细的时候，虽然最稳定的吸附结构的吸附能仍然是负的，但是在最稳定结构下，甲烷未必在碳管里面。

阿瞳

wzkchem5 发表于 2021-6-29 20:50
有没有考虑过可能这个碳纳米管太细了，甲烷塞进去能量是升高的？
你试试把甲烷放在管口，但是离管口至少 ...

老师，您的意思是碳纳米管的C原子对甲烷存在一定的排斥作用吗，我现在已经在算管口增大的了，估计明天就有数据了，谢谢老师

wzkchem5

阿瞳 发表于 2021-6-29 14:03
老师，您的意思是碳纳米管的C原子对甲烷存在一定的排斥作用吗，我现在已经在算管口增大的了，估计明天就 ...

对。简单的判断方法就是，看分子结构的时候把球棍模型改成space filling模型，也就是用半径等于范德华半径的球代表原子，如果甲烷的那些原子和碳管的原子之间有明显重叠，就说明斥力很大