24小时热门版块排行榜    

查看: 425  |  回复: 0

ra2ghgzh

木虫 (正式写手)

学痴

[求助] Si的bulk modulus是多少?

我分别使用了murn.f和origin拟合Si的晶格常数出现了不同的结果
使用murn.f拟合时的输入文件是
2
0.25
5.40        5.50        50
11
2
0.25
5.40        5.50        50
11
5.4        -10.79352
5.41        -10.7981
5.42        -10.80327
5.43        -10.80685
5.44        -10.80921
5.45        -10.81107
5.46        -10.81213
5.47        -10.81247
5.48        -10.81241
5.49        -10.8115
5.5        -10.80964
得到的结果是
ALATT0   VOL0 (ULA**3)  B0 (MBAR)     B0PRIME       E0 (EV)      E0 (RY)

alat=   5.47046 b0=    6.0777 b0p=     9.732 E0=  -0.397347
        5.4705     40.92725      6.07770      9.73227    -10.81262     -0.79471
        2.8948
用orign拟合时用的输入是
39.37        -10.79352
39.59        -10.7981
39.81        -10.80327
40.03        -10.80685
40.25        -10.80921
40.47        -10.81107
40.69        -10.81213
40.92        -10.81247
41.14        -10.81241
41.37        -10.8115
41.59        -10.80964
拟合函数是
E=E0+9/16*V0^3*B0*Bprime/V^2-27/16*V0*B0*Bprime*(V0/V)^(4/3)+27/16*V0*B0*Bprime*(V0/V)^(2/3)-9/16*V0*B0*Bprime+63/8*V0*B0*(V0/V)^(4/3)-9/4*V0^3*B0/V^2-9*V0*B0*(V0/V)^(2/3)+27/8*V0*B0
得到的结果是
V0=41.2
B0=0.256
Bprime=6.03
E0=-10.8
其中B0的变化最大,我想知道哪个是正确的,谢谢。
回复此楼

» 猜你喜欢

我永远当不了老板
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ra2ghgzh 的主题更新
信息提示
请填处理意见