24小时热门版块排行榜    

查看: 4435  |  回复: 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

至尊木虫 (职业作家)

3楼2021-03-27 17:13:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

sanwen2100

至尊木虫 (职业作家)

2楼2021-03-27 17:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +20 医学老男孩 2026-08-13 45/2250 2026-08-14 15:50 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 医学老男孩
[基金申请] filecode +5 cratir 2026-08-14 9/450 2026-08-14 12:25 by 医学老男孩
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +10 Tide man 2026-08-10 11/550 2026-08-14 10:30 by 沁言学术
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
信息提示
请填处理意见