24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2168  |  回复: 3

木棉和我

木虫 (职业作家)

[求助] 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倍)。请问,哪里出了问题?
回复此楼
已阅   回复此楼   关注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,是否与其沸点较低有关。
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的回帖
相关版块跳转 我要订阅楼主 木棉和我 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 光量子物理方向 博士招生 1人(2026.09) +3 sandyworld 2026-05-15 4/200 2026-05-17 14:38 by sandyworld
[教师之家] 上海大学实验技术岗位非升即走 +6 嘻嘻哈哈乐呵呵 2026-05-15 6/300 2026-05-17 14:35 by Ermito
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 xx7gd5zq4e 2026-05-15 7/350 2026-05-17 07:58 by 11n4dfd8yn
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 7hlccowb3h 2026-05-15 4/200 2026-05-17 07:46 by 11n4dfd8yn
[找工作] 售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,科目齐全,可+急 +6 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:16 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
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:55 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 v9tggjlwd0 2026-05-15 5/250 2026-05-17 00:32 by xiangfeng
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:30 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 l7k6xnh0yc 2026-05-14 4/200 2026-05-16 23:10 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 6/300 2026-05-16 22:55 by ue3ir18jc3
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见