| 查看: 1748 | 回复: 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 |
回复此楼» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急！在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