6楼: Originally posted by ChemiAndy at 2012-07-13 23:38:58
力场不能这样简单改。因为通常Bond项的能量是这样的：
E=k(b-b0)^2
这个参数k就是振动力常数，由振动频率决定；很明显，H2和D2由于质量不同，振动力常数不同。

VASP中可以直接改，是因为它是要通过电子结构计算 ...

2楼: Originally posted by ChemiAndy at 2012-07-09 04:18:25
你这个H2， D2， T2的气体力场应该很好做，可以通过气相量子化学计算拟合全部参数。

因为只有一个键，键参数可以通过量子力学力常数计算直接得到。LJ势可以通过做H2-H2, H2-D2, H2-T2, D2-T2的分子间势能面拟合得 ...

2楼: Originally posted by ChemiAndy at 2012-07-09 04:18:25
你这个H2， D2， T2的气体力场应该很好做，可以通过气相量子化学计算拟合全部参数。

因为只有一个键，键参数可以通过量子力学力常数计算直接得到。LJ势可以通过做H2-H2, H2-D2, H2-T2, D2-T2的分子间势能面拟合得 ...