24小时热门版块排行榜    

查看: 1789  |  回复: 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
  • 2017-02-12 15:22:34, 762.16 K
  • 附件 2 : GB.data
  • 2017-02-12 15:22:35, 43.68 K

» 猜你喜欢

哦哦哦哦哦哦哦哦哦哦哦哦哦哦
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by xxxxzzz at 2017-02-14 09:02:32
远离晶界的时候还可以,不知道为啥,大神能不能帮我调试一下啊,非常感谢了...

我也没算过这个呀,只能提些建议,既然远离晶界算的对,那是不是对晶界原子的处理有些问题?建议你先把某个确定距离上的几个原子形成能算对,再把所有的原子算出来看趋势
4楼2017-02-14 15:09:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
BaiXM那篇文章吧?我没算过这一部分,形成能在远离晶界的时候和完美晶格的结果能对上吗?如果不能应该就是计算方法有点问题
2楼2017-02-13 21:24:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxzzz

金虫 (正式写手)

引用回帖:
2楼: Originally posted by lx_PICO at 2017-02-13 21:24:08
BaiXM那篇文章吧?我没算过这一部分,形成能在远离晶界的时候和完美晶格的结果能对上吗?如果不能应该就是计算方法有点问题

远离晶界的时候还可以,不知道为啥,大神能不能帮我调试一下啊,非常感谢了
哦哦哦哦哦哦哦哦哦哦哦哦哦哦
3楼2017-02-14 09:02:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxzzz

金虫 (正式写手)

引用回帖:
4楼: Originally posted by lx_PICO at 2017-02-14 15:09:57
我也没算过这个呀,只能提些建议,既然远离晶界算的对,那是不是对晶界原子的处理有些问题?建议你先把某个确定距离上的几个原子形成能算对,再把所有的原子算出来看趋势...

能不能加个好友啊?

发自小木虫Android客户端
哦哦哦哦哦哦哦哦哦哦哦哦哦哦
5楼2017-02-16 23:42:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 环境学硕288求调剂 +6 皮皮皮123456 2026-03-22 6/300 2026-03-22 16:52 by i_cooler
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 285求调剂 +6 ytter 2026-03-22 6/300 2026-03-22 12:09 by 星空星月
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 311求调剂 +3 勇敢的小吴 2026-03-20 3/150 2026-03-21 17:40 by ColorlessPI
[考研] 313求调剂 +4 肆叁贰壹22 2026-03-19 4/200 2026-03-21 17:33 by ColorlessPI
[考研] 268求调剂 +9 简单点0 2026-03-17 9/450 2026-03-21 15:37 by lature00
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
信息提示
请填处理意见