24小时热门版块排行榜    

查看: 2044  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 木棉和我 的主题更新
信息提示
请填处理意见