| 查看: 1932 | 回复: 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
» 猜你喜欢
关于代码变化问题，想知道的进来
已经有12人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有3人回复
-
基金中了
已经有13人回复
-
关于filecode
已经有6人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有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
100