1. 第12行定义的变量x为温度,即计算这些温度下体系的晶格常数.
我没有做过非晶物质的建模,可能需要使用一些软件或者编程来产生.
是不是你把脚本的文件名改了? 找不到文件了,就会出现这个提示.
可以通过体积换算,也可以通过修改脚本直接输出晶格常数 可以考虑将温度控制在一个较小的范围内
脚本写的简洁明了,但是我想问下,我运行之后,改变npt的步数,发现温度下对应的体积会发生变化,然后热膨胀系数也会变化,请问那么如何确定正确的npt步数呢?
这个温度差不多是随便给的,是为了将体系从一个比较低的温度开始升温.