| 查看: 1688 | 回复: 5 | ||
[求助]
小程序Murn.f拟合晶格常数结果出现NaN如何解决？
|
|
初学vasp，在用程序murn.f进行拟合晶格常数时，有的能拟合出来，有的拟合不出来，拟合不出来的情况是输出文件out.m中出现NaN，不知道是什么原因？下面是我的输入和输出文件： 输入文件in.m 2 0.25 5.20 5.70 40 11 5.20 -11.805731 5.25 -11.874634 5.30 -11.920540 5.35 -11.945297 5.40 -11.950657 5.45 -11.938253 5.50 -11.909634 5.55 -11.866236 5.60 -11.809424 5.65 -11.740470 5.70 -11.660566 输出文件out.m least-squares fit of etot(v) by Murnaghan eq. UNIT OF LENGTH: ULA = 0.529177 ANGSTROEM UNIT OF ENERGY: IUNIT = 2: (ENERGIES INPUT IN ELECTRONVOLTS) IN THE STRUCTURE CONSIDERED, THE UNIT CELL VOLUME = A**3 * 0.250000 min and max alat(bohr) to compute Etot for at 40 points = 5.2000000 and 5.7000000 NUMBER OF DATA POINTS: N = 11 ------------------------------------------------------------------------------- I ALATT VOLUME ENERGY (EV) ENERGY (RY) 1 5.200000 35.152000 -11.8057310000 -11.8057310000 2 5.250000 36.175781 -11.8746340000 -11.8746340000 3 5.300000 37.219250 -11.9205400000 -11.9205400000 4 5.350000 38.282594 -11.9452970000 -11.9452970000 5 5.400000 39.366000 -11.9506570000 -11.9506570000 6 5.450000 40.469656 -11.9382530000 -11.9382530000 7 5.500000 41.593750 -11.9096340000 -11.9096340000 8 5.550000 42.738469 -11.8662360000 -11.8662360000 9 5.600000 43.904000 -11.8094240000 -11.8094240000 10 5.650000 45.090531 -11.7404700000 -11.7404700000 11 5.700000 46.298250 -11.6605660000 -11.6605660000 ------------------------------------------------------------------------------- ALATT0 VOL0 (ULA**3) B0 (MBAR) B0PRIME E0 (EV) E0 (RY) 5.4000 39.36600 1.00000 4.00000 -11.95066 0.00000 (=INITIAL GUESS) MIN: 2.7000 4.92075 0.01000 0.01000 MAX: 8.1000 132.86025 10.00000 10.00000 ------------------------------------------------------------------------------- FIT BY MURNAGHAN EQUATION EQUATION OF STATE ------------------------------------------------------------------------------- ITER ALATT0 VOL0 (ULA**3) B0 (MBAR) B0PRIME E0 (EV) SUM OF SQUARES IERR 20 NaN NaN NaN NaN NaN NaN 10 40 NaN NaN NaN NaN NaN NaN 0 40 NaN NaN NaN NaN NaN NaN 0 ------------------------------------------------------------------------------- ALATT0 VOL0 (ULA**3) B0 (MBAR) B0PRIME E0 (EV) E0 (RY) alat= NaN b0= NaN b0p= NaN E0= NaN NaN NaN NaN NaN NaN NaN NaN ------------------------------------------------------------------------------- alat(A) energy(Ryd) alat(bohr) energy(H) ------------------------------------------------------------------------------- 2.75172 NaN 5.20000 NaN 2.75851 NaN 5.21282 NaN 2.76529 NaN 5.22564 NaN 2.77207 NaN 5.23846 NaN 2.77886 NaN 5.25128 NaN 2.78564 NaN 5.26410 NaN 2.79243 NaN 5.27692 NaN 2.79921 NaN 5.28974 NaN 2.80600 NaN 5.30256 NaN 2.81278 NaN 5.31538 NaN 2.81956 NaN 5.32821 NaN 2.82635 NaN 5.34103 NaN 2.83313 NaN 5.35385 NaN 2.83992 NaN 5.36667 NaN 2.84670 NaN 5.37949 NaN 2.85349 NaN 5.39231 NaN 2.86027 NaN 5.40513 NaN 2.86705 NaN 5.41795 NaN 2.87384 NaN 5.43077 NaN 2.88062 NaN 5.44359 NaN 2.88741 NaN 5.45641 NaN 2.89419 NaN 5.46923 NaN 2.90098 NaN 5.48205 NaN 2.90776 NaN 5.49487 NaN 2.91454 NaN 5.50769 NaN 2.92133 NaN 5.52051 NaN 2.92811 NaN 5.53333 NaN 2.93490 NaN 5.54615 NaN 2.94168 NaN 5.55897 NaN 2.94847 NaN 5.57179 NaN 2.95525 NaN 5.58462 NaN 2.96204 NaN 5.59744 NaN 2.96882 NaN 5.61026 NaN 2.97560 NaN 5.62308 NaN 2.98239 NaN 5.63590 NaN 2.98917 NaN 5.64872 NaN 2.99596 NaN 5.66154 NaN 3.00274 NaN 5.67436 NaN 3.00953 NaN 5.68718 NaN 3.01631 NaN 5.70000 NaN |
回复此楼» 猜你喜欢
336求调剂
已经有5人回复
-
277材料科学与工程080500求调剂
已经有6人回复
-
333求调剂
已经有4人回复
-
286求调剂
已经有8人回复
-
298求调剂
已经有4人回复
-
306求0703调剂一志愿华中师范
已经有5人回复
-
299求调剂
已经有7人回复
-
26考研一志愿中国石油大学(华东)305分求调剂
已经有6人回复
-
初试 317
已经有5人回复
-
296求调剂
已经有3人回复
tsglss
铁杆木虫 (著名写手)
- 应助: 41 (小学生)
- 金币: 6014.2
- 红花: 19
- 帖子: 1875
- 在线: 940小时
- 虫号: 321641
- 注册: 2007-03-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2017-04-10 15:05:16
|
小程序murn.f是采用Murn状态方程来拟合晶格常数的。使用方法参考侯柱锋博士写的vasp手册。 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : murn.tar.gz
- 附件 2 : VASP的个人经验手册-侯柱峰2004版.pdf
2017-04-10 15:43:30, 11.79 K
2017-04-10 15:49:19, 245.92 K
3楼2017-04-10 15:49:49
小虫虫LG
木虫 (正式写手)
- 应助: 114 (高中生)
- 金币: 1775.9
- 散金: 1871
- 红花: 25
- 帖子: 975
- 在线: 678.5小时
- 虫号: 2128430
- 注册: 2012-11-15
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2017-04-10 15:51:56
5楼2017-04-10 16:12:55
tsglss
铁杆木虫 (著名写手)
- 应助: 41 (小学生)
- 金币: 6014.2
- 红花: 19
- 帖子: 1875
- 在线: 940小时
- 虫号: 321641
- 注册: 2007-03-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2017-04-10 19:48:18
