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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps 求 MgO 弹性常数 C44, 非正交盒子不能采用 Ewald 求和，怎么解？？
61/1返回列表
查看: 1662  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cgzhang_gg

铜虫 (著名写手)

[求助] lammps 求 MgO 弹性常数 C44, 非正交盒子不能采用 Ewald 求和，怎么解？？ 已有1人参与

我现在想用 lammps 求 MgO 的弹性常数，所采用经验势要用到库伦长程作用，可以用 Ewald 求和或者 PPPM 求解。采用正交模拟盒子，通过改变模拟盒子某一（多）个方向的长度求能量差，可以得到C11和C12的值，和实验相符。现在想求 C44, 需要改变盒子的形状（不是正交盒子了），问题就来了，lammps 不允许在非正交盒子内 Ewald 求和 或者 PPPM 求库伦长程作用，有谁碰到过这个问题，怎么解决的呢？非常感谢帮助！
回复此楼

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

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
cgzhang_gg: 金币+60, ★★★很有帮助, 谢谢回复，我尝试一下，遇到什么问题再向您请教。 2016-06-17 08:15:49
可以尝试用Wolf的直接截断法。我的印象是lammps支持这种方法。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
2楼2016-06-17 02:03:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ayl11112

新虫 (初入文坛)
不过对于建立不同的斜方晶系还是用正交的超晶胞来算试试吧，不然看看gromacs，说不定gromacs支持斜方的

发自小木虫Android客户端
赞一下 回复此楼
3楼2016-06-17 19:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cgzhang_gg

铜虫 (著名写手)
引用回帖:
3楼: Originally posted by ayl11112 at 2016-06-17 19:38:30
不过对于建立不同的斜方晶系还是用正交的超晶胞来算试试吧，不然看看gromacs，说不定gromacs支持斜方的

我试过正交的超晶胞，只是把晶胞内的原子按照斜方位置周期排列，但是采用CG方法优化后就变成完美的MgO结构（fcc结构），我只好没有优化，直接求能量，但这样计算出的弹性常数与实验有差别。
赞一下 回复此楼
4楼2016-06-18 09:01:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ayl11112

新虫 (初入文坛)
我后天帮你问问我导师吧，他什么都做应该知道吧。我的想法是min中固定晶胞中某个位置的原子，这样可以强迫保持斜方的晶系.

发自小木虫Android客户端
赞一下 回复此楼

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

cgzhang_gg
5楼2016-06-18 20:30:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cgzhang_gg

铜虫 (著名写手)
送红花一朵
引用回帖:
5楼: Originally posted by ayl11112 at 2016-06-18 20:30:17
我后天帮你问问我导师吧，他什么都做应该知道吧。我的想法是min中固定晶胞中某个位置的原子，这样可以强迫保持斜方的晶系.

我感觉只固定一个原子达不到效果，太谢谢你了，那你帮我问问吧
赞一下 回复此楼
6楼2016-06-19 11:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cgzhang_gg 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然上会要求 +5 无名者登山 2026-05-18 9/450 2026-05-18 17:50 by BlakeReary
[考博] 博士申请 +4 星…… 2026-05-18 5/250 2026-05-18 17:34 by 炎甲00
[教师之家] 上海大学实验技术岗位非升即走 +8 嘻嘻哈哈乐呵呵 2026-05-15 8/400 2026-05-18 11:53 by ning312
[硕博家园] 考博自荐 +5 科研狗111 2026-05-13 6/300 2026-05-18 11:22 by 糊糊涂涂好
[基金申请] 重磅！青年科学基金项目（C类）资助增幅预计超过50% +7 水和泥不是水泥 2026-05-13 10/500 2026-05-18 07:50 by 水和泥不是水泥
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 ky2p12rrjj 2026-05-15 4/200 2026-05-17 19:47 by Equinoxhua
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 7/350 2026-05-17 19:42 by Equinoxhua
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-17 19:39 by Equinoxhua
[公派出国] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 5/250 2026-05-17 18:39 by Equinoxhua
[考博] 找博士生导师 +4 小代想上岸 2026-05-15 5/250 2026-05-17 17:44 by 脱颖而出
[考博] 2026博士还有哪些学校有名额 +6 小王求读研 2026-05-15 7/350 2026-05-17 16:54 by 知音湖畔
[基金申请] 这年头没有找到涵评专家，还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[找工作] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[博后之家] 售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 k37jurhrau 2026-05-16 4/200 2026-05-17 01:35 by ue3ir18jc3
[有机交流] 求有机合成大神指点三硫酸乙烯酯（CAS：2793408-99-6）的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[教师之家] 教学课件你会给同学吗 +8 硕士研究生吗 2026-05-13 8/400 2026-05-14 22:23 by 常规沥青
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
关闭
请填处理意见
关闭