24小时热门版块排行榜    

查看: 2090  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 木棉和我 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研化学学硕调剂,一志愿985 +3 张vvvv 2026-03-15 3/150 2026-03-16 16:36 by houyaoxu
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 材料与化工(0856)304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 274求调剂0856材料化工 +12 z2839474511 2026-03-11 13/650 2026-03-13 10:39 by peike
[考研] 270求调剂 085600材料与化工专硕 +3 YXCT 2026-03-11 3/150 2026-03-13 10:13 by houyaoxu
[考研] 0857环境调剂 +5 熠熠_11 2026-03-10 5/250 2026-03-11 10:59 by wang_dand
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
信息提示
请填处理意见