24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1601  |  回复: 5

cgzhang_gg

铜虫 (著名写手)

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

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

» 猜你喜欢

已阅   回复此楼   关注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朵)

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 324求调剂 +6 hanamiko 2026-03-26 8/400 2026-03-27 04:40 by wxiongid
[考研] 287求调剂 +7 land xuxu 2026-03-26 7/350 2026-03-27 00:44 by wxiongid
[考研] 316求调剂 +4 江辞666 2026-03-26 4/200 2026-03-27 00:33 by wxiongid
[考研] 336材料求调剂 +7 陈滢莹 2026-03-26 9/450 2026-03-27 00:20 by wxiongid
[考研] 279求调剂 +6 红衣隐官 2026-03-21 6/300 2026-03-26 18:32 by 不吃魚的貓
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 10/500 2026-03-26 13:11 by 公瑾逍遥
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 一志愿上海交大生物与医药专硕324分,求调剂 +6 jiajunX 2026-03-22 6/300 2026-03-25 23:05 by licg0208
[考研] 07化学303求调剂 +5 睿08 2026-03-25 5/250 2026-03-25 22:46 by 418490947
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 359求调剂 +3 王了个楠 2026-03-25 3/150 2026-03-25 12:50 by Dyhoer
[考研] 285求调剂 +3 AZMK 2026-03-24 3/150 2026-03-25 12:23 by userper
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
信息提示
请填处理意见