24小时热门版块排行榜    

查看: 3221  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wangwen552

铜虫 (小有名气)

[求助] 如何利用LAMMPS产生层错 已有2人参与

请各位大侠帮忙看一下下面的内容如何用LAMMPS实现:
我想用LAMMPS计算层错能,看了几篇文献都是说把一个晶体分成两部分(比如上下两部分),一部分固定,另一部分相对做一个柏氏矢量的刚性移动,通过LAMMPS计算势能以此计算层错能,试了几次都不成功,哪位大侠知道到底怎样在LAMMPS中实现一部分固定,另一部分刚性移动啊?跪求帮忙!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leeve

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by leoluode at 2016-11-25 16:19:43
同学,请问你现在是采用什么命令做刚性位移的啊?我现在采用fix setforce 刚性固定原子,并且加diaplace_atoms实现位移,但是感觉计算的不正确啊,不知道你现在哟想到什么好的方法没有呢?

displace_atoms 是在计算(simulations)前对原子进行位移(参手册),与计算层错时的刚性位移应该不是一个概念
9楼2018-03-24 16:43:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

badercao

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wangwen552: 金币+10, ★★★很有帮助 2016-09-13 16:27:21
本帖内容被屏蔽

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

2楼2016-09-09 19:29:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangwen552

铜虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by badercao at 2016-09-09 19:29:37
fix rigid
fix setforce
fix addforce

大神你好,能否详细说明一下,这几个命令我都试过,但是总是不对。我是建了一个单晶,沿着Y方向分为上下两个region(block),也是两个group,我是想固定下半部分,将上半部分沿着X方向有一个柏氏矢量的刚性位移,那个fix rigid命令是只针对上半部分吗?还有后面两条命令分别是针对哪一部分呢?还有如果用了fix rigid命令,那系综是不是默认就是nve了,可我需要LAMMPS计算的是势能,如果用了nve那最后系统的能量都相等了这不就矛盾了吗?纠结死了。。。
3楼2016-09-11 10:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wangwen552: 金币+10, ★有帮助 2016-09-13 16:27:05
http://wenku.baidu.com/view/9053 ... a4.html?from=search
这个上海交大引入柏氏矢量的写的还可以,但是稍微有的繁杂,你参考参考吧
4楼2016-09-11 13:12:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-15 00:35 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-15 00:24 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:51 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 11:04 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 16:55 by zzy933268
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:49 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-12 10:38 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:27 by HajghKOMynmz
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:16 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:14 by HajghKOMynmz
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cdqGD6ftHmbl 2026-11-11 3/150 2026-11-12 07:53 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 zSbvVX9JFl9J 2026-11-09 8/400 2026-11-12 04:57 by 9fOTqj074tIg
信息提示
请填处理意见