24小时热门版块排行榜

返回列表

本帖产生 1 个模拟EPI ，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ChemiAndy

木虫 (正式写手)

模拟EPI: 40
应助: 52 (初中生)
金币: 2573.3
散金: 18
红花: 104
帖子: 590
在线: 553.8小时
虫号: 836907
注册: 2009-08-30
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhang_jaj: 金币+15, ★★★★★最佳答案, 谢谢指教 2012-07-09 07:27:16
chaizhm: 金币+5, 谢谢~ 2012-07-09 09:19:18
zh1987hs: 模拟EPI+1, 谢谢！ 2012-07-09 09:36:15
zhang_jaj: 金币+5, ★★★★★最佳答案, 谢谢 2012-07-09 19:40:23

你这个H2， D2， T2的气体力场应该很好做，可以通过气相量子化学计算拟合全部参数。

因为只有一个键，键参数可以通过量子力学力常数计算直接得到。LJ势可以通过做H2-H2, H2-D2, H2-T2, D2-T2的分子间势能面拟合得到。单质气体还无需考虑静电参数。所以十分容易开发。

一般的开发力场步骤的讨论可以看这里：http://muchong.com/bbs/viewthread.php?tid=3744217

储氢曾经很热，应该很容易找到大量的现有H2力场参数。你看看它们怎么做的参数，照样做个D2和T2的即可。你这个模拟的难点不在势参数，而在于量子效应，需要特殊的蒙特卡洛方法，对每一步模拟的坐标不确定性进行计算。