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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】边界能量处理
31/1返回列表
查看: 778  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ly67255776

新虫 (初入文坛)

[交流] 【求助】边界能量处理

求教，我在一个box里面设置上下两个region不同的晶格取向，左边的晶格是右边绕Z方向旋转45度所得，我本意是研究界面的能量，但是结果边界处的能量异常高，而界面处的能量却无法突出，求教如何解决这个问题！！！！in文件如下
units         metal
boundary    p p p
atom_style     atomic

lattice         fcc 3.61
region          box  block 0 6 0 6 0 36
create_box      1 box
create_atoms    1 box
region          up   block INF INF INF INF 18 INF
delete_atoms    region up
lattice         none
lattice         fcc 3.61 orient x 1 1 0 orient y -1 1 0
create_atoms    1 region up

pair_style     eam/alloy
pair_coeff     * * jin_copper_lammps.setfl Cu

neighbor     0.5 bin
neigh_modify     every 5 delay 0 check yes
timestep    0.01
thermo         10
variable      N equal atoms
variable     E equal pe

compute     pe all pe/atom
compute     ke all ke/atom
compute     coord all coord/atom 3.0

dump            1 all cfg 1 a.*.cfg id type xs ys zs c_pe c_ke c_coord
dump_modify     1 element Cu
run             0
undump          1

min_style     cg
minimize     1.0e-30 1.0e-15 1000 1000
print        "^---------fcc_boundary------------- N=$N  E=$E ------------------------"
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-01-09 11:36:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

he66yuan88

金虫 (正式写手)
注册之后4天才能发帖，先顶一个，再问问题。
赞一下 回复此楼
2楼2011-01-09 16:24:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beastx

银虫 (小有名气)
ly67255776(金币+50): 2011-01-14 09:56:19
这确实是一个值得研究的问题，我回去研究研究。界面原子配位能算全了没有
赞一下 回复此楼
3楼2011-01-13 10:06:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ly67255776 的主题更新
31/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:42 by XiYP04It3Evb
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:20 by XiYP04It3Evb
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:28 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:18 by XiYP04It3Evb
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 13:15 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 7xHjyB3bAiwg 2026-10-24 6/300 2026-10-25 13:11 by XiYP04It3Evb
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 9/450 2026-10-25 13:02 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-25 12:49 by XiYP04It3Evb
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:49 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-25 12:45 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:36 by XiYP04It3Evb
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:32 by XiYP04It3Evb
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 GVGRDgNfv7P4 2026-10-23 4/200 2026-10-25 12:10 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 01:22 by XcOYG53d3x4p
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-24 23:35 by XcOYG53d3x4p
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-24 22:56 by XcOYG53d3x4p
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 18:23 by WqtG0Aem8Ufl
[考博] 27秋申博自荐 +4 且听风吟时 2026-10-21 4/200 2026-10-23 13:26 by 大大大肥虫
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
关闭
请填处理意见
关闭