| 查看: 9635 | 回复: 23 | |||||||||||||||||
| 本帖产生 2 个 模拟EPI ，点击这里进行查看 | |||||||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||||||||
[交流]
关于势函数的求助
|
|||||||||||||||||
|
相信做分子模拟亲们都会用到势函数？ 但是请问大家的势函数参数如何取得的呢？ 我也相信有部分童鞋是查得文献，同时肯定也有部分童鞋是自己拟合的。 请问是通过什么方法或软件拟合呢？希望有详细的步骤和例子。 感谢大家！感谢虫子们！ |
回复此楼» 收录本帖的淘帖专辑推荐
 资源收集 | 分子动力学 | 分子模拟 | 模拟知识汇总 |
| MS 吸附模拟 | 模型 - 模拟 - 统计 - 热力学 | 我学习计算的一些帖子 | 动力学模拟 |
| 模拟软件安装精华淘 | 动力学强帖 | 我的分享 | MS学习 |
| 力场 | 分子动力学 |
» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于吸附MD方面，跪求指导。
已经有7人回复
- 关于周期信号的收敛性的问题
已经有4人回复
- 空位迁移能的计算
已经有10人回复
- 请教吸附动力学拟合
已经有11人回复
- 势函数截断距离如何取？
已经有11人回复
- 请问势函数在lammps官网中怎么下载？
已经有6人回复
- 分子动力学势函数
已经有6人回复
- 离子迁移数的测定
已经有5人回复
- 拟和Cu-Al势函数
已经有6人回复
- 【求助】量子化学方法求amber力场参数
已经有10人回复
- 【求助】求 嵌入势函数 的物理解释
已经有4人回复
- 【求助】ms中如何假如势函数，比如EAM，比如自己修改过的势函数。
已经有12人回复
- 【求助】塑性与弹性势函数
已经有5人回复
- 【求助】关于ISIF的选择
已经有8人回复
- 【求助】MS建模导入lammps问题求救
已经有14人回复
- 【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置
已经有4人回复
» 抢金币啦！回帖就可以得到:
-
基金祈福
+1/964
-
面上第三年申请，改了一遍又一遍，希望能中啊！散金祈福！！！
+5/850
-
何不食肉糜
+1/477
-
每年这时候我都再仔细读一遍《范进中举》
+1/472
-
华西生物制药研究院招聘科研助理
+2/400
-
基金申请小调查！预祝申请必定中中中！！祈福！！！
+5/260
-
坐标北京诚征男友
+1/112
-
美国University of Texas at Arlington土木结构工程方向2027年春季全奖博士招生
+1/88
-
中国药科大学博士，博后招聘（高分子，金属有机方向），欢迎大家交流推荐
+1/86
-
Polymers 期刊特刊 组稿
+1/76
-
那年春暖花开，今朝如梦初醒 —— 北漂十五载，相亲十年，归乡孑然一身
+1/69
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
+1/32
-
上海大学王亮课题组诚招博士后（碳材料、光催化、电催化方向）
+2/28
-
***大叔找对象（沪深等地）
+1/20
-
上海交通大学化学化工学院倪伟焱课题组招收2027年申请考核制博士生（电催化+高分子）
+1/12
-
pcold 原核表达载体
+1/8
-
2027环境工程博士申请
+1/8
-
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/6
-
出全新蔡司 Plan?APOCHROMAT 20×/0.8 平场复消色差物镜
+1/2
-
找女友
+1/1
16楼2011-10-29 13:34:57
3楼2011-10-26 15:17:52
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+10, 模拟EPI+1): 谢谢 2011-10-27 17:16:14
小木虫(金币+0.5):给个红包，谢谢回帖
御剑江湖(金币+10, 模拟EPI+1): 谢谢 2011-10-27 17:16:14
|
MS中的Gulp模块适合你。 看看手册吧。我也没用过。 晶体结构(a,b,c)4%的误差是稍微有点大。拟合力场晶体结构方面的偏差一般能达到2%以内，密度方面小于5%, 晶格能方面小于10%。因为必须考虑实验值也是有误差的。晶格能，蒸发焓这些实验值的误差也在10%左右，刻意弄成跟实验值完全一致没有必要，只要跟类似的物质比，结合能定性大小顺序符合即可； 另外，你还要问自己，即你要研究的问题是什么？这个问题需要多精确的晶格常数值？因为有时候你调了参数，看上去与实验值更接近了，可是一跑动力学，一升高温度，得到的结果反而更差。比如你现在要预测体积膨胀系数，那我建议你先用现在的参数跑一些模拟，得到一些结果，看看跟实验值的误差有多大。如果在5%以内是能够接受的。 初学者往往会产生错觉，以为分子模拟能真实再现所有原子过程并给出精确的实验结果。但其实经典力场本身是个十分粗糙的模型，它忽略了电子的运动，把原子看成是带电荷的球，分子间的作用通过一个简化的数学模型（Buckingham+coulomb）来表达，没有多体作用(可达~5%)，没有极化(可达1%)，请问你觉得这在多大程度上符合真实的原子结构和相互作用？一般来说， 我们使用分子力学模型并非是为了完全精确预测实验结果，而是来定性地研究某一性质的变化规律，或者现象背后的原子过程。即，关注相对值，而非绝对值。 我并不是不建议你调参数，相反我认为做模拟的都要尝试验证参数，可能的话通过拟合更多的实验数据来改进参数。但是在调参数的时候，要认识到力场模型的局限性，以及它的误差范围和精度范围。想像你拿的是一把最小刻度为cm的尺子，而两根棍子的长度差别是几个mm。你能得到精确的结果吗? [ Last edited by ChemiAndy on 2011-10-26 at 11:32 ] |
4楼2011-10-26 23:52:22
|
你好！你的回答让我茅塞顿开，非常感谢！ 我研究的是材料的缺陷，是化学缺陷，不是物理缺陷，举个例子吧：1、研究晶体中某个原子的取代，某个原子的空缺，某个原子离开自己的位置跑到间隙去了等等一系列的问题，然后计算这个过程产生的能量，再举个例子吧：2、当我向这个材料体系中去掺杂的时候，我要掺杂的离子（非常微量、基本不影响原材料体系的晶体结构）可能会占据原体系材料的那个位置，通过模拟计算得出能量值，然后就可以判断杂质离子可能会进入哪个位置，然后来指导我的实验。 在这种情况下，我在选取势参数（原材料体系的）的时候，验证参数的合理性，应该考虑它的哪些因素？比如晶格常数、弹性模量、密度、结合能等等，你觉得哪个最重要，其实我也知道研究不同的内容，要求的不一样。也就是说选取的势参数不可能同时保证晶格常数、弹性模量、密度、结合能都和实验值想符合，在这种情况下，你觉得我在判断势参数合理性时，最应该关注哪一个参考量呢？目前我只考虑了晶格常数这一个参考量。 |
5楼2011-10-27 09:53:02
