24小时热门版块排行榜    

查看: 1791  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 接收2026硕士调剂(学硕+专硕) +3 allen-yin 2026-03-23 3/150 2026-03-23 11:29 by 奶油草莓.
[考研] 070300化学求调剂 +5 苑豆豆 2026-03-20 5/250 2026-03-23 09:39 by guoliang1816
[考研] 招08考数学 +4 laoshidan 2026-03-20 11/550 2026-03-23 08:47 by laoshidan
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 0854电子信息求调剂 +3 α____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 307求调剂 +11 冷笙123 2026-03-17 11/550 2026-03-22 20:16 by edmund7
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 一志愿华中科技大学071000,求调剂 +4 沿岸有贝壳6 2026-03-21 4/200 2026-03-22 07:21 by ilovexiaobin
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 294求调剂材料与化工专硕 +15 陌の森林 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见