| 查看: 1737 | 回复: 4 | |||
[交流]
【求助】用murn.f来拟合晶格常数 已有3人参与
|
|
各位虫友，你们好！我用网上的murn.f来拟合晶格常数，得到的out.m文件如下，请问那个是晶格常数？另外侯老师给出的g/alat=是什么意思？ # unit of length = 0.529177 angstroem # energies input in ELECTRONVOLTS # unit cell volume = a**3 * 0.250000 # a(min)= 3.5000000 # a(max)= 4.3000000 # points: 50 # Zero point lattice energy by Moruzzi-Janak-Williams: # Grueneisen constant = 3.500 # T(Debye) = -5.73 # Volume/at for T(Debye) = 3.60 # Tested only for fcc lattice (Al) # number of data points: n = 3 #---------- # i alatt volume E # 1 3.800000 13.718000 -7.7110664890 -7.7110663000 # 2 3.900000 14.829750 -7.8155879439 -7.8155878000 # 3 4.000000 16.000000 -7.7645597103 -7.7645596000 #---------- # alatt0 vol0 (ula**3) b0 (mbar) b0prime e0 # 3.9000 14.82975 1.00000 4.00000 0.00000 # min: 1.9500 1.85372 0.01000 0.01000 # max: 5.8500 50.05041 10.00000 10.00000 #---------- # fit by murnaghan equation equation of state #---------- #iter, alatt0, vol0 (ula**3), b0 (mbar), b0prime, e0, sq sum, ierr # 20 3.9006 14.83620 17.31752 14.43339 -7.8155902 0.00000 0 #---------- # a (a.u.), E #---------- 3.50000 1.262708 9.078298 3.51633 -0.372188 7.443402 3.53265 -1.718709 6.096881 3.54898 -2.828148 4.987442 3.56531 -3.742495 4.073095 3.58163 -4.496152 3.319438 3.59796 -5.117333 2.698257 3.61429 -5.629198 2.186392 3.63061 -6.050774 1.764816 3.64694 -6.397705 1.417885 3.66327 -6.682869 1.132721 3.67959 -6.916873 0.898717 3.69592 -7.108463 0.707127 3.71224 -7.264853 0.550737 3.72857 -7.392004 0.423586 3.74490 -7.494842 0.320749 3.76122 -7.577440 0.238151 3.77755 -7.643172 0.172418 3.79388 -7.694834 0.120756 3.81020 -7.734744 0.080846 3.82653 -7.764825 0.050766 3.84286 -7.786674 0.028916 3.85918 -7.801621 0.013969 3.87551 -7.810769 0.004821 3.89184 -7.815038 0.000552 3.90816 -7.815195 0.000395 3.92449 -7.811878 0.003712 3.94082 -7.805620 0.009970 3.95714 -7.796866 0.018725 3.97347 -7.785987 0.029603 3.98980 -7.773296 0.042295 4.00612 -7.759050 0.056540 4.02245 -7.743469 0.072121 4.03878 -7.726735 0.088855 4.05510 -7.709000 0.106591 4.07143 -7.690392 0.125198 4.08776 -7.671019 0.144571 4.10408 -7.650972 0.164618 4.12041 -7.630326 0.185264 4.13673 -7.609145 0.206446 4.15306 -7.587482 0.228108 4.16939 -7.565382 0.250208 4.18571 -7.542884 0.272706 4.20204 -7.520019 0.295571 4.21837 -7.496814 0.318776 4.23469 -7.473291 0.342299 4.25102 -7.449470 0.366120 4.26735 -7.425366 0.390224 4.28367 -7.400993 0.414598 4.30000 -7.376361 0.439229 |
回复此楼» 猜你喜欢
我面上完蛋了
已经有11人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急！在vasp中使用murn.f程序拟合晶格常数的问题
已经有18人回复
- 用LDA赝势优化出来的晶格常数一定比实验值小吗？谢谢！
已经有10人回复
- 采用murn.f计算体积弹性模量不正确？
已经有17人回复
- 晶格常数优化
已经有7人回复
- 【求助】为什么用castep几何优化晶格常数下不来？
已经有22人回复
- 【求助成功】晶格优化中，murn.f的使用方法
已经有7人回复
- 【求助】拟合晶格常数
已经有14人回复
- 【求助】求教关于优化
已经有6人回复
condensed
木虫 (著名写手)
- 应助: 21 (小学生)
- 金币: 3223.5
- 散金: 44
- 红花: 3
- 帖子: 1002
- 在线: 208.4小时
- 虫号: 292567
- 注册: 2006-11-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-06-23 09:58:15
3楼2010-06-23 11:15:33
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 1ST强帖: 5
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7441.4
- 散金: 9654
- 红花: 119
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
hedaors(金币+1):Thanks, 其实用 isif=3 还是更保险点 2010-06-23 13:20:18
小木虫(金币+0.5):给个红包，谢谢回帖交流
hedaors(金币+1):Thanks, 其实用 isif=3 还是更保险点 2010-06-23 13:20:18
|
晶格常数已经优化出来了,就没有必要用ISIF=3俩,2就好. 其它随便,收敛就好. 呵呵:认真想了想,结构优化我还是没掌握到位. 对于原胞,我直接优化晶格参数,然后就用保持形状和体积来优化原子位置了。这个对于简单的晶体结构在能量和结合结构上也能满足一定精度要求.如果是复杂的结构,可能有必要做ISIF=3的细致优化. 这种全面优化在能量和几何参数上能在精度上提高多少(比如受力从 0.02->0.01),我没有概念.希望有人可以提供参考. [ Last edited by cenwanglai on 2010-6-23 at 13:48 ] |
4楼2010-06-23 13:07:42
5楼2010-06-23 15:04:23
20