| 查看: 870 | 回复: 5 | ||||
[求助]
The eigenvector number1 has zero norm
|
|
在計算graphene時 不知道為何會在這一部分算不過去 由於我是由單層BN的input改寫的 應該只需要變更lattice constant吧 但卻出現了虛數無法計算的部分 可以請問是哪裡出問題了呢 para.err為 0: 0: fxphas_seq.F90:137:ERROR 0: fxphas_seq: 0: The eigenvector number1 has zero norm. tnlo_2.in檔如下 #Definition of lattice parameters #******************************** acell 4.6686384981 4.6686384981 30 rprim 1 0 0 -0.5 0.8660254038 0.0 0 0 1 #Definition of atomic positions and types #************************************* natom 2 ntypat 1 znucl 6 typat 1 1 xred 0.33333 0.66667 0.5 0.66667 0.33333 0.5 ixc 7 #Parameters of the SCF cycles #**************************** iscf 7 nstep 200 #Plane wave basis and k-point grid #********************************* ecut 50 # ecutsm 0.5 dilatmx 1.05 ngkpt 12 12 2 nshiftk 1 shiftk 0.0 0.0 0.5 ndtset 5 jdtset 1 2 3 4 5 #DATASET1 : scf calculation: GS WF in the BZ #******************************************** prtden1 1 kptopt1 1 toldfe1 1.0d-12 #DATASET2 : non scf calculation: GS WF in the whole BZ #***************************************************** getden2 1 kptopt2 2 iscf2 -2 getwfk2 1 tolwfr2 1.0d-22 nband2 4 #DATASET3 : ddk on the MP grid of k-points #***************************************** getwfk3 2 rfdir3 1 1 1 rfelfd3 2 tolwfr3 1.0d-22 nband3 4 kptopt3 2 #DATASET4 : ddE on the MP grid of k-points #***************************************** prtden4 1 getwfk4 2 getddk4 3 rfdir4 1 1 1 rfphon4 1 rfatpol4 1 2 tolvrs4 1.0d-12 rfelfd4 3 nband4 4 kptopt4 2 prepanl4 1 #DATASET5 : 3DTE calculation #*************************** getden5 1 get1den5 4 getwfk5 2 get1wf5 4 nband5 4 kptopt5 2 optdriver5 5 rf1elfd5 1 rf1phon5 1 rf1atpol5 1 2 rf1dir5 1 1 1 rf2elfd5 1 rf2dir5 1 1 1 rf3elfd5 1 rf3dir5 1 1 1 nbdbuf 0 |
回复此楼» 猜你喜欢
化学专业申博
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
评审有感
已经有27人回复
-
论文撤稿了
已经有7人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
河北省自然科学基金
已经有9人回复
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
2楼2013-09-18 11:31:30
|
DATASET 4 計算開始會NAN 我附上了output檔 感謝你的回應 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : tnlo_2.outA
2013-09-18 14:11:02, 332.38 K
3楼2013-09-18 14:12:20
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
★ ★
franch: 金币+2, 谢谢回帖交流 2013-09-23 20:44:08
franch: 金币+2, 谢谢回帖交流 2013-09-23 20:44:08
 尝试多次，附件下载不下来～LZ的DataSet 1有没有很好的收敛？我测试了一下楼主的input，发现DataSet 1的时候就非常难以收敛。 仔细观察了一下LZ的input，发现至少有两处不妥。 1，对于xred的有效数字问题。为了更有效准确地找到对称性，Abinit一般要求至少10位有效数字，而LZ明显只给了6位。比较简单的解决方法显然是写为分数： xred 1/3 2/3 1/2 2/3 1/3 1/2； 2，对于Graphene类的东西，由于是模拟2D的物质，所以k点在z方向上为1，也即：ngkpt 12 12 1。 经过此改动后，DataSet 1和DataSet 2的收敛是没有问题了。 但是DataSet 3的收敛还是有问题的(准确说就是无法收敛)，现在在尝试解决RF的收敛问题。 如果DataSet 3无法有效收敛，在达到额定步数后，程序依然会认为计算完成，并把结果用于下一步的计算，自然在DataSet 4会出现NaN～ |
4楼2013-09-22 23:50:39
5楼2013-09-23 10:17:40
6楼2013-09-23 10:45:03
10