| 查看: 829 | 回复: 3 | |||
[求助]
Abinit结构优化时使用自动优化给出的不是最优解？ 已有1人参与
|
|
先谢过 初学abinit 真心问题多。。。 按照论坛里面台湾人的教程 把经验值acell c/a带入optcell=0 ionmov=3 得到一组xred 把xred带入做optcell=2 ionmov=3 计算 给输入文件： ------------------- dilatmx 1.05 tolmxf 1.0d-05 ionmov 3 optcell 2 ntime 20 ecutsm 0.5 getwfk -1 prtden 0 prteig 0 chksymbreak 0 kptopt 1 ngkpt 6 6 6 toldfe 1.0d-6 nstep 30 diemac 12 ecut 100 enunit 2 ntypat 3 natom 7 znucl 49 80 34 typat 1 1 2 3 3 3 3 acell 10.87 10.87 10.87 rprim 0.5 -0.5 1.02 -0.5 -0.5 1.02 -0.5 -0.5 -1.02 xred 0.0000000000E+00 -5.5511151231E-17 5.5511151231E-17 5.0000000000E-01 2.5000000000E-01 -7.5000000000E-01 0.0000000000E+00 5.0000000000E-01 -5.0000000000E-01 5.1408395117E-01 7.8573180120E-02 -1.2857754903E-01 4.8591604883E-01 6.4266150021E-01 -5.9265713129E-01 5.0004368913E-02 8.7142245097E-01 -3.5733849979E-01 -5.0004368913E-02 4.0734286871E-01 -9.2142681988E-01 ------------------------------------------------------ 给出out结果： acell 1.3257791925E+01 1.3257791925E+01 1.3257791925E+01 rprim 4.0680634832E-01 -4.0680634832E-01 8.1793471007E-01 -4.0680634832E-01 -4.0680634832E-01 8.1793471007E-01 -4.0680634832E-01 -4.0680634832E-01 -8.1793471007E-01 虽然这组rprim有点奇怪。。。但是还是验证下 都保持不变 给acell+ 3*0.001 acell1 1.3257791925E+01 1.3257791925E+01 1.3257791925E+01 Bohr acell2 1.3267791925E+01 1.3267791925E+01 1.3267791925E+01 Bohr acell3 1.3277791925E+01 1.3277791925E+01 1.3277791925E+01 Bohr acell4 1.3287791925E+01 1.3287791925E+01 1.3287791925E+01 Bohr acell5 1.3297791925E+01 1.3297791925E+01 1.3297791925E+01 Bohr etotal1 -1.3142002925E+02 etotal2 -1.3142007513E+02 etotal3 -1.3142010670E+02 etotal4 -1.3142012419E+02 etotal5 -1.3142012794E+02  不知道怎么个情况了 一般来说 程序寻值不是应该比手工靠谱+精确吗？。。。另外 这个13.25Bohr*0.407*2的结果倒是和经验值很接近呢。。。 |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有6人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有12人回复
-
filecode
已经有16人回复
-
重要来源：本周末出结果
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有11人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有6人回复
-
不应该看fileCode
已经有9人回复
-
Filecode 又变了，巨变
已经有4人回复
-
分享一下我之前已中青C的计划书的filecode
已经有5人回复
-
FileCode能看出啥？
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这样的优化问题有哪些方法可以做？
已经有11人回复
- CASTEP几何结构优化的问题 虽然收敛了 感觉不是很对 请大家帮忙找一下原因
已经有11人回复
- 进行结构优化时就不动了 完全是菜鸟 不懂阿 求助！！！！
已经有12人回复
- 收敛测试和结构优化
已经有10人回复
- 用castep进行结构优化的遇到的一些问题
已经有17人回复
- 对irc得到的结构是否可以（or 需要）作进一步结构优化？
已经有7人回复
- vasp结构优化中EDIFF以及KPOINTS设置，以及氢吸附结合能问题！
已经有17人回复
- 结构优化的结果竟然重复不出来，每次优化后的结构都不一样，求高手解答……
已经有9人回复
- 请教一下那款免费软件的结构优化功能比较好用？
已经有7人回复
- VASP结构优化时出现的POTIM问题
已经有4人回复
- 晶体结构优化和画能带图的方法是否保持一致？？？
已经有12人回复
- vasp 结构优化问题
已经有14人回复
- vasp进行固定某些原子的某些坐标的结构优化
已经有5人回复
- 结构优化中问题
已经有6人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- 关于结构优化中的对称性问题
已经有9人回复
- 【求助】结构优化不收敛
已经有9人回复
- 【求助】用Gaussian进行结构优化
已经有19人回复
- 【求助】晶体结构优化，计算出问题，请大家帮忙！
已经有7人回复
- 【求助】Castep中结构优化的问题
已经有11人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
2楼2014-03-30 16:40:49
3楼2014-04-19 19:25:59
4楼2014-04-19 19:31:12
