| 查看: 1242 | 回复: 6 | |||
| 当前主题已经存档。 | |||
tfl03木虫 (职业作家)
|
[交流]
【求助】如何使用XMD来进行氧化物的计算
|
||
|
XMD声称能进行陶瓷的计算，陶瓷氧化物居多，都有正负离子，不知道XMD如何进行EVALD求解或其它方法求解库仑能？下面给一个MgO的例子，短程相互作用为Buckingham势参数形式。我已经检验过，如果只计算一个团簇，所得到的结果与GULP没有区别。但对于，晶体，似乎不能处理。 #作为晶体似乎不可能，应为XMD没有EVALD方法来求库仑能 #下面是Buckingham势参数，只有Core，没有shel EUNIT eV POTENTIAL SET PAIR 2 CALC Qmg = 2.0 CALC Qo = -2.0 CALC A1 = 1280.0 CALC ROU1 = 0.3 CALC C1 = 0.0 CALC A2 = 1280.0 CALC ROU2 = 0.3 CALC C2 = 27.88 CALC TIMES =1000 CALC START =1.0 CALC FINAL =8.0 POTENTIAL PAIR 1 1 TIMES START FINAL FORMULA (14.40057*Qmg*Qmg)/r POTENTIAL PAIR 1 2 TIMES START FINAL FORMULA (14.40057*Qmg*Qo)/r+ A1*EXP(-r/ROU1)-C1/(r*r*r*r*r*r) POTENTIAL PAIR 2 2 TIMES START FINAL FORMULA (14.40057*Qo*Qo)/r + A2*EXP(-r/ROU2)-C2/(r*r*r*r*r*r) calc MG=1 calc O=2 calc MassMG=24.0 calc MassO=16.0 calc A0=4.05 Box 15 15 15 fill particle 8 1 0 0 0 1 0 0.5 0.5 1 0.5 0 0.5 1 0.5 0.5 0 2 0.5 0.5 0.5 2 0.5 0 0 2 0 0.5 0 2 0 0 0.5 fill go scale 4.235384 4.235384 4.235384 write particle select all move all 50 50 50 write particle #下面一行将晶格常数保存在txt文件里 erase constant.txt bsave 1 constant.txt erase energy.txt esave 1 energy.txt select type MG mass MassMG select type O mass MassO #下面两行的数字设定温度 clamp 1 itemp 1 dtime 1e-15 #下面三行的数字设定压力，如果考察温度的影响时把后两行关闭，前加#即可 pressure clamp 1.2 repeat 20 cmd 10 write box end #结果无法求得晶格常数 [ Last edited by lei0736 on 2009-11-25 at 12:29 ] |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
2楼2009-05-26 20:53:25
老虎大王
木虫 (著名写手)
- 模拟EPI: 17
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+8,VIP+0):谢谢 解释得非常清楚 5-27 09:24
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+8,VIP+0):谢谢 解释得非常清楚 5-27 09:24
|
一、XMD所称的陶瓷主要指碳化物、氮化物等，因为它可以使用Tersoff势和SW势。至于你所说的氧化物陶瓷，原子间结合有一定的离子键特征，可视为离子晶体来处理。 二、对于你所说的含离子键，或者说原子带有电荷、需要计算库仑力的体系，XMD也能处理，方法就像你所做的那样，把库仑势做为一般的双体势的一部分。但是因为XMD里面没有专门处理库仑力的程序，比如没有用Ewald求和等方法，而是直接算1/r，所以效率是比较低下的。但也确实是可以算的。 三、你的Script我看过了，形式上没有什么问题。我也试着运行了，可以顺利执行。 四、你所说的不能求得晶格常数，应该是指运行过程中元胞体积不断增大。如果你说的是这个意思，那么不是程序问题。程序的任务是输出你要算的东西。现在是可以输出的，说明程序没有问题。双体势对于MD来说总是相对比较粗糙的。对双体势作用下的离子体系，等压MD中元胞尺寸不断长大，这是一个常见的事情。我本人就经常遇见。问题应该在于你的势函数。 五、压力在MD中本来就不是一个容易控制的量。用MD方法来优化晶胞尺寸，不是一个好的选择。要考察特定势函数对晶胞尺寸的预测，还是要用晶格动力学方法，使用Gulp是较好的选择。 Gulp预测和优化晶胞尺寸，基本方法是能量极小化，这与MD并非相同的概念。 [ Last edited by 老虎大王 on 2009-5-27 at 10:25 ] |
3楼2009-05-27 09:19:36
tfl03
木虫 (职业作家)
- 应助: 24 (小学生)
- 金币: 3128.5
- 散金: 1126
- 红花: 42
- 帖子: 4453
- 在线: 1812.2小时
- 虫号: 772313
- 注册: 2009-05-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2009-05-27 19:14:04
5楼2009-07-19 21:04:14
6楼2009-07-19 22:24:58
7楼2009-07-20 09:17:33
15