24小时热门版块排行榜    

查看: 4432  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 水壶呼呼 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +14 医学老男孩 2026-08-13 34/1700 2026-08-14 09:21 by xiaruohan
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 30/1500 2026-08-14 07:55 by 医学老男孩
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +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 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
信息提示
请填处理意见