24小时热门版块排行榜    

查看: 4483  |  回复: 2

水壶呼呼

新虫 (小有名气)

[交流] lammps学习记录----restart write_restart read_restart 已有1人参与

Sample Text
在网上看到一个总结很好的,挂下链接~侵删
https://www.360doc.com/content/12/1120/08/8275160_249016551.shtml
restart命令是lammps提供的重启动命令。比如驰豫各个计算都要用,可以先单独驰豫一下,生成restart文件,供其他计算直接调用,从而避免每个计算都进行驰豫。
中间写一个restart文件,在文件中记录某个状态的绝大部分信息,restart文件是二进制文件。
restart2data可以将restart的二进制文件转化为文本文件,可以由read_data读入。
restart和write_restart不同点:前者周期性将系统状态写入不同的restart文件,后者写出该命令使用前的系统状态。使用位置也不同:

restart files文件必须是N的整数倍才可以输出,与dump不同,dump会输出刚开始的还有如果最后是能量最小化的也会输出,所以可以先用write_restart文件在run之前写一个e二进制文件,对于最后一步是能量最小化,并且收敛的话(N>0),restart会输出最后一步
  read_restart是读前两者的文件,开始新的计算
需要注意的是,并不是所有的模拟信息都保存在restart.equil文件中。
比如fix命令就需要在新的模拟文件中重新设置,官方手册已经列举的很全面,有需要可以查查手册。

关于read_restart:
什么信息包括:重建模拟盒子,原子,和他们的分布,以及一些相关的全局设定
根据当前模拟中的处理器数量和处理器命令的设置,将模拟盒子划分为矩形砖的常规3d网格,每个处理器一个。以后可以通过balance或fix balance命令更改分区。
units

newton bond (see discussion of newton command below)

atom style and atom_modify settings id, map, sort

comm style and comm_modify settings mode, cutoff, vel

timestep

simulation box size and shape and boundary settings

atom group definitions

per-type atom settings such as mass

per-atom attributes including their group assignments and molecular topology attributes (bonds, angles, etc)

force field styles (pair, bond, angle, etc)

force field coefficients (pair, bond, angle, etc) in some cases (see below)

pair_modify settings, except the compute option

special_bonds settings


什么没包括:某些fix命令,包括fix shake  fix langevin。某些 pair styles,因为计算力的时候依赖于速度,只用half-step的值,而重新读入时,刚开始计算时用的时full-step的速度,这些fix命令包括 granular pair styles, pair dpd, and pair lubricate  
(个人建议:将fix命令重写)
newton pair (see discussion of newton command below)

fix commands (see below)

compute commands (see below)

variable commands

region commands

neighbor list criteria including neigh_modify settings

kspace_style and kspace_modify settings

info for thermodynamic, dump, or restart output

many-body和tabulated势不能用,需要重新写pair_coeff  bond_coeff 等
回复此楼

» 猜你喜欢

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

sanwen2100

至尊木虫 (职业作家)

2楼2021-03-27 17:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sanwen2100

至尊木虫 (职业作家)

3楼2021-03-27 17:13:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 水壶呼呼 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见