应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 如何利用LAMMPS产生层错
101/1返回列表
查看: 3227  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wangwen552

铜虫 (小有名气)

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

请各位大侠帮忙看一下下面的内容如何用LAMMPS实现：
我想用LAMMPS计算层错能，看了几篇文献都是说把一个晶体分成两部分（比如上下两部分），一部分固定，另一部分相对做一个柏氏矢量的刚性移动，通过LAMMPS计算势能以此计算层错能，试了几次都不成功，哪位大侠知道到底怎样在LAMMPS中实现一部分固定，另一部分刚性移动啊？跪求帮忙！
回复此楼

» 猜你喜欢
1楼 2016-09-09 16:10:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

badercao

禁虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wangwen552: 金币+10, ★★★很有帮助 2016-09-13 16:27:21
本帖内容被屏蔽

» 本帖已获得的红花（最新10朵）

wangwen552
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的回帖

wangwen552

铜虫 (小有名气)
引用回帖:
4楼: Originally posted by lx_PICO at 2016-09-11 13:12:42
http://wenku.baidu.com/view/90536b3bee06eff9aef807a4.html?from=search
这个上海交大引入柏氏矢量的写的还可以，但是稍微有的繁杂，你参考参考吧

多谢！我去看一下
回复此楼
5楼2016-09-11 13:17:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

badercao

禁虫 (正式写手)
本帖内容被屏蔽
6楼2016-09-11 16:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangwen552

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by badercao at 2016-09-11 16:29:12
大神谈不上，只是针对你讲的，我觉得用上面的命令可能会实现。（1）fix rigid 命令针对group。（2）后面两条命令是加力的另外一种方法。（3）使用fix rigid也是可以计算势能的，具体check一下manual。 如果还有疑问 ...

恩，谢谢！
回复此楼
7楼2016-09-12 09:39:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leoluode

银虫 (小有名气)
同学，请问你现在是采用什么命令做刚性位移的啊？我现在采用fix setforce 刚性固定原子，并且加diaplace_atoms实现位移，但是感觉计算的不正确啊，不知道你现在哟想到什么好的方法没有呢？
赞一下(1人) 回复此楼
8楼2016-11-25 16:19:43
已阅   回复此楼   关注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的回帖

Stoudemire29

金虫 (小有名气)
8楼同学在吗？fix setforce 和displace_atoms的用法能介绍一下吗？新手一枚
赞一下 回复此楼
共同进步
10楼2018-10-29 14:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangwen552 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:43 by xSrQNSSjcgoX
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 03:32 by xSrQNSSjcgoX
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-14 13:35 by KXypaGgBBlw7
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-13 23:37 by 6PrO1XyZ5URF
[考博] 考博求助 10+3 社会实践对别 2026-11-08 3/150 2026-11-13 14:04 by 沐小虫252
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1Nb0IVlmV6v0 2026-11-11 4/200 2026-11-13 11:18 by Equinoxhua
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 04:21 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:27 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:16 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一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:46 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 09:21 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:03 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 05:19 by 9fOTqj074tIg
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
信息提示
关闭
请填处理意见
关闭