LJ参数ε和σ的量子化学计算

wzkchem5

吞吞12138 发表于 2023-10-17 14:50
老师请问这两个参数对于某个分子来说不是常数吗，仿真中只是因为材料库中没有这种气体，需要输入这两个参 ...

我的意思是说某个分子的σ和ε不是唯一的，适合用来描述总相互作用能的σ和ε，和适合描述范德华相互作用的σ和ε，一般而言并不一样。需要知道你们材料库里的σ和ε是描述什么作用能（总作用能？Keesom+Debye+London+Pauli排斥？只有London+Pauli排斥？）用的

吞吞12138

wzkchem5 发表于 2023-10-17 22:44
我的意思是说某个分子的σ和ε不是唯一的，适合用来描述总相互作用能的σ和ε，和适合描述范德华相互作用 ...

是用来描述范德华相互作用的σ和ε，不考虑Keesom+Debye，请问老师这个数据需要学习软件自行计算吗，这个数据在分子动力学仿真中是不是属于基础呀。

wzkchem5

吞吞12138 发表于 2023-10-18 05:54
是用来描述范德华相互作用的σ和ε，不考虑Keesom+Debye，请问老师这个数据需要学习软件自行计算吗，这个 ...

那可以这样做：随机产生一些由两个C4F7N分子组成的体系，其中两个分子的距离涵盖从略小于范德华半径之和到远大于范德华半径之和的各种值，分子的初始构象、取向是任意的。然后在固定两个分子质心距离的情况下，做结构优化至收敛。再用任意一种能给出London项和Pauli排斥项的能量分解方法（例如LED）得到当前质心距离下的London项和Pauli排斥项的大小。最后把London项和Pauli排斥项之和对质心距离作图，用Lennard-Jones函数拟合即可。
但前提是你确定你确实不要Keesom和Debye项，这两项对于你的体系或许小，但是也是可以计算的。而且大部分实验只能测到总作用能，对气体的很多宏观性质起作用的也是总作用能，而不仅仅是London项和Pauli排斥项。