24小时热门版块排行榜    

查看: 6312  |  回复: 32

huang998

木虫 (小有名气)

[求助] 如何对nve系综进行弛豫? 已有3人参与

在用lammps模拟裂纹扩展中,如何对nve系综进行弛豫?是用mini_style,和minimize?
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wxrui

新虫 (正式写手)

nve系综直接run不就可以了?min_style和 minimize 能量最小化,不是弛豫

» 本帖已获得的红花(最新10朵)

3楼2017-06-05 21:19:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wxrui

新虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
huang998: 金币+3, 有帮助 2017-07-10 11:02:24
nve系综直接run不就可以了?min_style和 minimize 能量最小化,不是弛豫
4楼2017-06-05 21:20:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wxrui

新虫 (正式写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by huang998 at 2017-06-05 21:54:58
安装包里的例子里没有min _style. minimize这些命令,但是删除原子后怎么驰豫?
...

这样说吧,min_style 和minimize 是用于优化结构的方法,有几类算法可以实现,常用的有共轭梯度法(cg),模拟退火,最速下降法等,目的都是找个尽可能全局的(局域)能量最低点最稳定的结构,这时候的优化出来的结构没有温度。另外,你删除部分原子与弛豫没有必然的联系,也就是说,你删除了部分原子之后你依然可以nve弛豫,或者minimize,如果所有原子都删除,,也就不存在研究体系了,也就没法弛豫或者minimize了

» 本帖已获得的红花(最新10朵)

7楼2017-06-06 12:15:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wxrui

新虫 (正式写手)

引用回帖:
11楼: Originally posted by huang998 at 2017-06-06 16:02:46
我也不太清楚,如果要删除部分原子的话,应该用npt?
...

删除原子是建立模型的方法之一,和后期你运行哪种系综是没有关系的,具体的系综使用需要结合你的研究体系以及研究条件
12楼2017-06-06 16:34:06
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cl1994861

金虫 (小有名气)

引用回帖:
16楼: Originally posted by huang998 at 2017-06-07 11:35:22
如果是先删除原子,然后弛豫,可以像这样吗?控制温度的变化,run一段时间??
# equilibrate

velocity        mobile create 300.0 5812775 temp new3d
fix                1 all nve
fix                2 boundary setforce 0.0 0.0 0.0

f ...

对的,如果你指的是删除原子后,想让系综重新稳定的话,不过你的这个时间有点短

» 本帖已获得的红花(最新10朵)

19楼2017-06-07 16:45:52
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

huang998

木虫 (小有名气)

2楼2017-06-05 20:39:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

huang998

木虫 (小有名气)

送红花一朵
引用回帖:
3楼: Originally posted by wxrui at 2017-06-05 21:19:21
nve系综直接run不就可以了?min_style和 minimize 能量最小化,不是弛豫

安装包里的例子里没有min _style. minimize这些命令,但是删除原子后怎么驰豫?

发自小木虫IOS客户端
5楼2017-06-05 21:54:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cl1994861

金虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
huang998: 金币+2, 有帮助 2017-07-10 11:02:48
引用回帖:
5楼: Originally posted by huang998 at 2017-06-05 21:54:58
安装包里的例子里没有min _style. minimize这些命令,但是删除原子后怎么驰豫?
...

这2个命令是用来能量最小化的,弛豫过程说的就是运行过程中达到所设置系综条件的过程。你说的删除原子后弛豫?原子都没有不存在弛豫一说

» 本帖已获得的红花(最新10朵)

6楼2017-06-06 10:48:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

huang998

木虫 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by wxrui at 2017-06-06 12:15:33
这样说吧,min_style 和minimize 是用于优化结构的方法,有几类算法可以实现,常用的有共轭梯度法(cg),模拟退火,最速下降法等,目的都是找个尽可能全局的(局域)能量最低点最稳定的结构,这时候的优化出来的结 ...

比如我先建一个0 30 0 20 -5 -5的box,然后用delete_atoms删除一些原子后,要先equilibrate吧,如果是nve系综可以用temp/rescale控制温度进行驰豫?如果不是那怎么办?

发自小木虫IOS客户端
8楼2017-06-06 14:30:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

huang998

木虫 (小有名气)

送红花一朵
引用回帖:
6楼: Originally posted by cl1994861 at 2017-06-06 10:48:40
这2个命令是用来能量最小化的,弛豫过程说的就是运行过程中达到所设置系综条件的过程。你说的删除原子后弛豫?原子都没有不存在弛豫一说...

有原子,我一个模型几万个原子中,删除几十个后,要先平衡再进行其他吧?

发自小木虫IOS客户端
9楼2017-06-06 14:33:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wxrui

新虫 (正式写手)

【答案】应助回帖

★ ★ ★
huang998: 金币+3, 有帮助 2017-07-10 11:02:58
引用回帖:
8楼: Originally posted by huang998 at 2017-06-06 14:30:39
比如我先建一个0 30 0 20 -5 -5的box,然后用delete_atoms删除一些原子后,要先equilibrate吧,如果是nve系综可以用temp/rescale控制温度进行驰豫?如果不是那怎么办?
...

删除部分原子,表示你已经建好模型了,然后就优化结构达到你想要的一个状态,run 就可以了。nve系综我没用过,所以不太清楚,这需要结合你的模拟体系,研究目的。温控算法有好几种,nose-hoover算法用的比较多,这也需要根据你的研究体系决定。另外nve系综应该不能控温吧?本身能量就是恒定了,,再标定粒子速度控温,,感觉有些奇怪
10楼2017-06-06 15:28:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huang998 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见