24小时热门版块排行榜    

查看: 1558  |  回复: 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 的主题更新
信息提示
请填处理意见