24小时热门版块排行榜    

查看: 346  |  回复: 0
【悬赏金币】回答本帖问题,作者水壶呼呼将赠送您 5 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

水壶呼呼

新虫 (小有名气)

[求助] 关于restart的问题

最近在使用restart时出现了一个比较问题,一直没想明白。有没有哪位大佬路过能帮忙看一下
有一组对比,分为A组B组。
A组计算过程:在一个in文件里包含弛豫和后续的计算。
B组计算过程:先弛豫,存一个restart二进制文件,再在一个in文件里将这个文件作为初始模型直接进行后续的计算。
具体代码如下:
B组:
存restart文件
shell mkdir initial
#--------------------------DEFINITION-------------------
units                 metal
dimension         3       
boundary        s s p
atom_style        atomic

read_data Cu.data

pair_style eam/alloy
pair_coeff * * Cu01.eam.alloy Cu
#------------------------equilibrium---------------------
timestep 0.002
velocity all create 300.0 12345

min_style cg
minimize 1e-18 1e-18 100000 100000
fix 1 all nve
fix 2 all temp/rescale 10 300.0 300.0 10.0 1.0

dump 1 all custom 8000 ./initial/relax.*.dump id type x y z

run 8000
undump 1
write_restart ./initial/restart.equil


计算过程
shell mkdir config
read_restart restart.equili
pair_style eam/alloy
pair_coeff * * Cu01.eam.alloy Cu
reset_timestep 0
timestep 0.002
dump 2 all custom 1 ./config/shear.*.dump id type x y z  #c_eng c_1

  • dump_modify 2 delay 8875
    fix 3 all deform 1 z scale 0.92 remap x units box
    run 9001

    -----------------------------------------------------------------------------------------------
    A组
    shell mkdir initial
    shell mkdir config

    #--------------------------DEFINITION-------------------
    units                 metal
    dimension         3       
    boundary        s s p
    atom_style        atomic

    read_data Cu.data

    pair_style eam/alloy
    pair_coeff * * Cu01.eam.alloy Cu

    #compute eng all pe/atom
    #compute 1 all stress/atom NULL
    #------------------------equilibrium---------------------
    timestep 0.002
    velocity all create 300.0 12345

    min_style cg
    minimize 1e-18 1e-18 100000 100000
    fix 1 all nve
    fix 2 all temp/rescale 10 300.0 300.0 10.0 1.0

    dump 1 all custom 1000 ./initial/relax.*.dump id type x y z

    thermo 2000
    thermo_style custom step lx ly lz press pe temp

    run 8000
    undump 1

    #-------------------------SHEAR---------------------------------
    reset_timestep 0
    timestep 0.002
    dump 2 all custom 1 ./config/shear.*.dump id type x y z  #c_eng c_1

  • dump_modify 2 delay 8875

    #thermo         100000
    #thermo_style        custom step temp ke pe press

    fix 3 all deform 1 z scale 0.92 remap x units box

    run 9001
  • 回复此楼
    已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
    相关版块跳转 我要订阅楼主 水壶呼呼 的主题更新
    不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
    最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
    [基金申请] 离职后国自然项目不能变更单位,在新单位还能申请新的国自然项目吗 5+3 宜兰 2024-05-31 6/300 2024-05-31 12:05 by felicity6056
    [教师之家] 研究生在毕业答辩时挂了,遗憾 +9 spacr2022 2024-05-27 12/600 2024-05-31 10:17 by zzq770204
    [硕博家园] 求问endnote上Wiley引用格式在哪下载 +3 Elon. 2024-05-30 3/150 2024-05-31 10:13 by crazy peng
    [考博] 广东以理材料系碳点与功能材料课题组 — 2博士名额 / 科研助理 +3 小城夜很美 2024-05-27 9/450 2024-05-31 09:28 by 小城夜很美
    [文学芳草园] 对对子啊 +5 天若孤独 2024-05-29 7/350 2024-05-31 09:00 by wjykycg
    [硕博家园] 好奇博士每天学习的有效时间 +10 hahamyid 2024-05-25 11/550 2024-05-30 18:46 by lyfbangong
    [电化学] 有锂金属软包电池代工公司代做软包嘛 +5 yyywj995 2024-05-26 6/300 2024-05-30 08:33 by 那年的旧城
    [硕博家园] 论大家对6070后普通教授导师的看法 +5 SNaiL1995 2024-05-28 9/450 2024-05-29 21:54 by SNaiL1995
    [博后之家] 2024公派博后申请 +4 326lhpqk 2024-05-27 5/250 2024-05-29 20:03 by @古月胡
    [考博] 24年博士招生 +7 abinit432 2024-05-27 9/450 2024-05-29 17:16 by 中国银河
    [论文投稿] 有没有老师需要发表论文 +4 金老师论文助理- 2024-05-29 4/200 2024-05-29 16:51 by liuyupu132
    [考研] 研0 +3 徐小七七 2024-05-25 5/250 2024-05-29 13:49 by 1158057902
    [论文投稿] 核心初审被拒,理由是“选题的意义不明确,文章写得不像是科技论文”,怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
    [论文投稿] EI学报,一审返修后,为啥不再送审,直接终审中? +4 qweasd12345 2024-05-27 6/300 2024-05-29 00:02 by dut_ameng
    [基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
    [材料综合] 求助,刚入门MOF合成,实验总是失败 20+3 张月半 2024-05-25 8/400 2024-05-28 09:20 by 我想更懂你29
    [硕博家园] 我是很理想化一人 +6 hahamyid 2024-05-26 6/300 2024-05-27 18:13 by 大飞鱼鱼鱼
    [基金申请] 感觉自然基金限制通过比例就是有点扯,学学B口,化学学部,不限制比例。 +10 wsjing 2024-05-26 14/700 2024-05-27 11:57 by kanmiaolucky
    [硕博家园] 求助 +3 单增李斯特21633 2024-05-25 3/150 2024-05-27 10:33 by hahamyid
    [硕博家园] 研0 +6 控制调剂yl 2024-05-25 11/550 2024-05-26 23:37 by sakuraai
    信息提示
    请填处理意见