| 查看: 1972 | 回复: 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
» 猜你喜欢
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
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
xxxxzzz
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1682
- 散金: 200
- 帖子: 303
- 在线: 83小时
- 虫号: 2018899
- 注册: 2012-09-21
- 性别: GG
- 专业: 核技术及其应用
5楼2017-02-16 23:42:16
5
