24小时热门版块排行榜    

查看: 6307  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 关于代码变化问题,想知道的进来 +9 且听虎啸 2026-08-07 12/600 2026-08-07 21:53 by 理想三旬2025
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了,论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +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
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见