24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1797  |  回复: 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的回帖

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的回帖
查看全部 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +6 HanBeiNingZC 2026-03-24 6/300 2026-03-28 07:55 by baoball
[考研] 275求调剂 +10 jjjjjjjjjjl 2026-03-27 10/500 2026-03-27 23:47 by barnett0632
[考研] 085602 307分 求调剂 +6 不知道叫什么! 2026-03-26 6/300 2026-03-27 22:00 by 催化大白
[考研] 070300求调剂306分 +3 26要上岸 2026-03-27 3/150 2026-03-27 17:57 by arrow8852
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +9 vv迷 2026-03-22 9/450 2026-03-27 15:59 by 不吃魚的貓
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +10 哈哈157349 2026-03-21 10/500 2026-03-27 13:06 by zzll406
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 0703化学338求调剂! +6 Zuhui0306 2026-03-26 7/350 2026-03-27 10:35 by shangxh
[考研] 求调剂,一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +4 @taotao 2026-03-26 5/250 2026-03-27 08:10 by hypershenger
[考研] 321求调剂 +6 Ymlll 2026-03-24 6/300 2026-03-26 20:50 by 不吃魚的貓
[考研] 化学工程085602 305分求调剂 +17 RichLi_ 2026-03-25 17/850 2026-03-26 19:44 by plmuchong
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 一志愿吉林大学材料与化工303分求调剂 +4 为学666 2026-03-24 4/200 2026-03-25 11:27 by BruceLiu320
[考研] 318求调剂 +3 plum李子 2026-03-23 3/150 2026-03-25 09:42 by 雾散后相遇lc
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 考研调剂 +3 呼呼?~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
信息提示
请填处理意见