| 查看: 1244 | 回复: 8 | |||
[交流]
【求助】dmol3的结果怎么拟合力场？ 已有6人参与
|
|
DMOL3得到的结果怎么能拟合得到想要的力场参数？比如说先优化CH2CH2CH2NH2和CO2之间的相互作用能，然后改怎么做，拟合得到这些不同原子的力场参数？请高手指教。。。 [ Last edited by qli0218 on 2010-4-15 at 12:54 ] |
回复此楼» 猜你喜欢
不自信的我
已经有10人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有10人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
-
要不要辞职读博？
已经有3人回复
-
论文终于录用啦！满足毕业条件了
已经有26人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
自荐读博
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教有关DMOL3的设置问题
已经有10人回复
- 【求助】Dmol3 晶体计算
已经有15人回复
- 【求助】dmol3中外加电场的单位a.u.怎么换算成国际单位
已经有11人回复
- 【求助】Gulp中的几何优化，MD等和Dmol3的中的有什么区别？
已经有11人回复
- 【求助】谁知道Dmol3 加电场实现的细节？？
已经有8人回复
- 【求助】如何根据势能曲线拟合力场参数
已经有19人回复
- 【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】
已经有15人回复
- 【求助】DMOl3中电场怎么加?
已经有4人回复
2楼2010-04-15 09:15:40
3楼2010-04-15 11:45:22
4楼2010-04-16 08:41:46
wangtingyun
银虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 186.8
- 散金: 12
- 红花: 2
- 帖子: 147
- 在线: 65.1小时
- 虫号: 494160
- 注册: 2008-01-12
- 专业: 金属有机化学
5楼2010-04-21 14:35:31
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
zyj8119(金币+1):谢谢参与
zyj8119(金币+5): 2010-04-21 19:02
gavinliu7390(金币+2):谢谢交流！ 2010-04-23 00:58
zyj8119(金币+1):谢谢参与
zyj8119(金币+5): 2010-04-21 19:02
gavinliu7390(金币+2):谢谢交流！ 2010-04-23 00:58
|
这个要求太强大了吧，还真不知道MS软件自身有没有拟合力场这个功能。 如果是自己手动拟合的话，至少要对力场中势函数的形式非常熟悉吧。一般函数都与原子间距有关，形如E=f(r)之类的，而其中函数f的形式是已知的，其中的一些参数才是真正需要拟合的。那么根据不同r值下dmol的计算E的结果应该是可以拟合出函数f中的参数的。但应该是很复杂的工作，而且如何提高势函数的迁移性等等，都是很难解决的问题。 如果更夸张的想自己设计f的形式，那将是一个很大的项目了。 另外楼主似乎要处理分子间的相互作用，提醒一下，我记得用DFT的描述分子间作用力的时候是需要特别注意的，相关文章应该有不少。 |
6楼2010-04-21 17:48:15
ym23
铁杆木虫 (著名写手)
- 应助: 4 (幼儿园)
- 贵宾: 0.186
- 金币: 9769.9
- 散金: 690
- 红花: 10
- 帖子: 1462
- 在线: 258.8小时
- 虫号: 837382
- 注册: 2009-08-31
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
zyj8119(金币+1):谢谢参与
zyj8119(金币+5): 2010-04-21 22:11
zyj8119(金币+4): 2010-04-23 00:07
gavinliu7390(金币+2):谢谢交流！ 2010-04-23 00:58
zyj8119(金币+1):谢谢参与
zyj8119(金币+5): 2010-04-21 22:11
zyj8119(金币+4): 2010-04-23 00:07
gavinliu7390(金币+2):谢谢交流！ 2010-04-23 00:58
|
正如6楼所说，我觉得体材料可以这样做。 1. 用MS的castep模块，计算不同晶格常数下材料的总能。即晶格常数每变化0.1或0.2A，就算一个能量，这个只计算能量不用几何优化。至少算30组数据。 2. 将这几十组数据输入到orgin里，横坐标用晶格常数，纵坐标用能量值，用描点法画条曲线，这条曲线应该是像“对号”一样的曲线。 3. 找一个函数，即6楼说的形式已知的函数。我们通常用morse或rose势，来拟合你描出的曲线，orgin拟合方法网上有视频，请楼主自己找找吧。 4. 拟合的结果就是你想要的势参数，即力场参数。 希望我们做体材料力场参数的方法对你有帮助。 |
7楼2010-04-21 21:43:57
8楼2010-04-23 00:07:13
qzhaosdu
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1401.6
- 散金: 2052
- 红花: 7
- 帖子: 1969
- 在线: 263.5小时
- 虫号: 493351
- 注册: 2008-01-10
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 量子化学
9楼2010-04-23 02:12:34