24小时热门版块排行榜    

查看: 4437  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 水壶呼呼 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +6 Tide man 2026-08-14 7/350 2026-08-14 18:36 by 正能量1斤
[基金申请] filecode +6 cratir 2026-08-14 10/500 2026-08-14 18:29 by 笑叹辞穷
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见