24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1822  |  回复: 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

金虫 (正式写手)

引用回帖:
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的回帖
查看全部 5 个回答

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考研] 一志愿沪9,326求生物学调剂 +11 刘墨墨 2026-04-13 11/550 2026-04-17 15:35 by wanganpin
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 7/350 2026-04-16 14:31 by zhouxiaoyu
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 346分,工科0854求调剂,专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 调剂结束 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
信息提示
请填处理意见