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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps运行时nvt下对应压力10391145bar，npt下对应体积扩大12倍，哪里出了问题？
41/1返回列表
查看: 2347  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

木棉和我

木虫 (职业作家)

[求助] lammps运行时nvt下对应压力10391145bar，npt下对应体积扩大12倍，哪里出了问题？ 已有1人参与

诸位。采用lammps计算MgO-TiO2-B2O3体系的玻璃化转变或结晶转变，关注体积和能量变化。势函数采用Materials Explorer5.0软件势函数库中的Born-Mayer-Huggins potential parameters，文献来源为Hirao K, Kawamura K. Material design using personal computer[J]. Shokabo, Tokyo, 1994, 52。in文件的主要语句为：unit metal；timestep 0.001；velocity all create 3000 4928459 dist gaussian；fix 0 all nve；fix 1 all temp/rescale 1 3000 3000 5.0 1； run 8000；unfix 0；unfix 1；fix 2 all npt temp 3000 300 0.1 iso 10378167 10378167 1；run 27000；unfix 2。密度按照实验测定设置为4.00g/cm3。nvt下运行，显示压力为10391145bar，体积为39136.177埃3(对应于之前设置的密度)。npt下运行，压力设置为1bar，体积则为482515.33  埃3（与之前比扩大了12倍）。请问，哪里出了问题？
回复此楼

» 猜你喜欢
1楼 2018-06-13 13:04:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
实验密度是否为3000开的密度？

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

木棉和我

木虫 (职业作家)
引用回帖:
2楼: Originally posted by brucefan at 2018-06-13 13:51:23
实验密度是否为3000开的密度？

密度为常温下的实测密度值。建立的模型为非晶态模型，即原子随机嵌在立方盒子内。B2O3沸点为2133K，是否与其沸点较低有关。
赞一下(1人) 回复此楼
3楼2018-06-13 14:37:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
木棉和我: 金币+20, ★★★★★最佳答案 2018-06-19 10:15:14
引用回帖:
3楼: Originally posted by 木棉和我 at 2018-06-13 14:37:23
密度为常温下的实测密度值。建立的模型为非晶态模型，即原子随机嵌在立方盒子内。B2O3沸点为2133K，是否与其沸点较低有关。...

你用室温密度做3000开的NVT系综模拟，压强大很正常。做模拟时，如果想要模拟常压条件，就用NPT。只要运用合理，控压算法会给出正确的密度或者盒子大小。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
4楼2018-06-16 19:37:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 木棉和我 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-15 15:14 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:50 by VMDqgnlKimtZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 14:34 by VMDqgnlKimtZ
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:58 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:26 by C79jjtjAKjEn
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-15 06:14 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 23:13 by vZfe6xYu34yj
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭