| 查看: 1796 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
xxxxzzz金虫 (正式写手)
|
[求助]
晶界点缺陷形成能计算求助 已有1人参与
|
|
|
大家好，本人最近试着参照文献（Efficient Annealing of Radiation Damage Near Grain Boundaries via Interstitial EmissionScience 327, 1631 (2010)）在计算Cu晶界附近点缺陷（间隙子）的形成能，但是计算的结果总是乱七八糟的，与文献的结果相差很大，有没有大神能够指导下，给些详细的模拟in文件参考下，谢谢了！联系方式：1572111245@qq.com ########################### units metal # 单位类型 dimension 3 newton on on boundary p p p atom_style atomic # 原子类型 neighbor 2.0 bin neigh_modify delay 5 check yes ##################construction############# read_data GB.data group cu type 1 #############################potential energy################## pair_style eam/alloy pair_coeff * * Cu01.eam.alloy Cu ########################equilibriation phase###################### reset_timestep 0 compute eng all pe/atom compute eatoms all reduce sum c_eng thermo 10 thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms minimize 1.0e-15 1.0e-15 5000 5000 run 0 variable N equal count(all) variable No equal $N variable E equal "c_eatoms" variable Eo equal $E ################################# create_atoms 1 single ${xi} ${yi} ${zi} units box reset_timestep 0 thermo 10 thermo_style custom step pe lx ly lz press pxx pyy pzz c_eatoms min_style cg minimize 1e-15 1e-15 5000 5000 run 0 undump 1 variable Ef equal "c_eatoms" variable Ecoh equal -3.53999996838068 variable Ei equal (${Ef}-${Eo}-${Ecoh}) print "Initial energy of atoms = ${Eo}" print "Final energy of atoms = ${Ef}" print "Interstitial formation energy" print "${Ei}" |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Cu01.eam.alloy
- 附件 2 : GB.data
2017-02-12 15:22:34, 762.16 K
2017-02-12 15:22:35, 43.68 K
» 猜你喜欢
xxxxzzz
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1682
- 散金: 200
- 帖子: 303
- 在线: 83小时
- 虫号: 2018899
- 注册: 2012-09-21
- 性别: GG
- 专业: 核技术及其应用
5楼2017-02-16 23:42:16
2楼2017-02-13 21:24:08
xxxxzzz
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1682
- 散金: 200
- 帖子: 303
- 在线: 83小时
- 虫号: 2018899
- 注册: 2012-09-21
- 性别: GG
- 专业: 核技术及其应用
3楼2017-02-14 09:02:32
4楼2017-02-14 15:09:57
