应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » Co3O4的LDA+U计算报错求助
91/1返回列表
查看: 3472  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

c_chi

木虫 (正式写手)

[求助] Co3O4的LDA+U计算报错求助

我计算Co3O4报错如下，不知是什么问题，尝试换赝势也不行，in文件贴在了后面。
谢谢！！！
atomic mag. moment =     3.000000
N of occupied +U levels =  168.000000
--- exit write_ns ---
Atomic wfc used for LDA+U Projector are NOT orthogonalized
     Starting wfc are  440 randomized atomic wfcs

     total cpu time spent up to now is       20.5 secs

     per-process dynamical memory:    89.3 Mb

     Self-consistent Calculation

     iteration #  1     ecut=    30.00 Ry     beta=0.40
     Davidson diagonalization with overlap

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     Error in routine cdiaghg (1516):
     S matrix not positive definite
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...


#######################################
我的in文件
&CONTROL
title = "MnO, multiplicity=6, U=2.5"
calculation = "scf",
restart_mode = "from_scratch",
prefix = "MnO.c0.s6.u2.5",
outdir = ".//a0.08",
pseudo_dir = "./",
wf_collect=.TRUE.
/
&SYSTEM
ibrav=1,celldm(1)=8.022446,ecutwfc=30, ecutrho=300,
starting_magnetization(1)=0.25,starting_magnetization(2)=0.25,
nspin=2,occupations="fixed",
nat=56,ntyp=2,tot_charge=0.000000,tot_magnetization=5.000000
lda_plus_U=.true.,Hubbard_U(1)=2.5,Hubbard_alpha(1)=0.08
/
&ELECTRONS
startingpot="file"
mixing_beta=0.4,conv_thr=1D-9,diago_thr_init=1D-11,electron_maxstep=150
/
ATOMIC_SPECIES
Co 58.933 Co.pbe-sp-van.UPF
O 15.999 O.PBE.US.RRKJ3
ATOMIC_POSITIONS (angstrom)
Co 1.003 1.003 1.003
Co  7.020   7.020   7.020
Co  5.014   1.003   5.014
Co  3.008   7.020   3.008
Co  1.003   5.014   5.014
Co  7.020   3.008   3.008
Co  5.014   5.014   1.003
Co  3.008   3.008   7.020
Co  4.011   4.011   4.011
Co  2.006   6.017   0.000
Co  6.017   2.006   0.000
Co  6.017   0.000   2.006
Co  2.006   0.000   6.017
Co  0.000   2.006   6.017
Co  0.000   6.017   2.006
Co  4.011   0.000   0.000
Co  2.006   2.006   4.011
Co  6.017   6.017   4.011
Co  6.017   4.011   6.017
Co  2.006   4.011   2.006
Co  0.000   4.011   0.000
Co  4.011   2.006   2.006
Co  4.011   6.017   6.017
Co  0.000   0.000   4.011
O  2.082   2.082   2.082
O  5.941   5.941   5.941
O  3.935   7.946   6.093
O  4.087   0.076   1.929
O  7.946   6.093   3.935
O  0.076   1.929   4.087
O  6.093   3.935   7.946
O  1.929   4.087   0.076
O  0.076   4.087   1.929
O  7.946   3.935   6.093
O  4.087   1.929   0.076
O  3.935   6.093   7.946
O  1.929   0.076   4.087
O  6.093   7.946   3.935
O  2.082   6.093   6.093
O  5.941   1.929   1.929
O  3.935   3.935   2.082
O  4.087   4.087   5.941
O  7.946   2.082   7.946
O  0.076   5.941   0.076
O  0.076   0.076   5.941
O  7.946   7.946   2.082
O  4.087   5.941   4.087
O  3.935   2.082   3.935
O  6.093   2.082   6.093
O  1.929   5.941   1.929
O  2.082   3.935   3.935
O  5.941   4.087   4.087
O  5.941   0.076   0.076
O  2.082   7.946   7.946
O  6.093   6.093   2.082
O  1.929   1.929   5.941
K_POINTS (automatic)
3 3 3 0 0 0
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-03-16 23:04:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pangrui1985

铜虫 (小有名气)
你的截断能够吗
回复此楼
2楼2015-03-21 23:26:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
请问楼主你这个问题解决了吗？我也碰到了这个错误。
赞一下 回复此楼
未来就是现在
3楼2015-06-01 08:45:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

c_chi

木虫 (正式写手)
引用回帖:
3楼: Originally posted by future_wl at 2015-06-01 08:45:02
请问楼主你这个问题解决了吗？我也碰到了这个错误。

当时没解决，我就用了文献中的数值
赞一下 回复此楼
4楼2015-06-01 19:43:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
4楼: Originally posted by c_chi at 2015-06-01 19:43:56
当时没解决，我就用了文献中的数值...

我的刚解决，是晶格常数单位设置错了，celldm(1)的单位应该是a.u,我当时用成bohr了，改过来之后就好了！
赞一下 回复此楼
未来就是现在
5楼2015-06-01 22:02:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zheng8386

新虫 (初入文坛)
能请教一下你的Co3O4如何建模的，
赞一下 回复此楼
6楼2015-09-24 09:50:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

迷迷物理

新虫 (初入文坛)
请教我算了一个162个原子的LDA+U体系，可是显示atomic wfc used for LDA+U projector are not orthogonalized是什么原因，怎么修改？

发自小木虫Android客户端
赞一下(1人) 回复此楼
7楼2018-04-01 09:06:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lastzealot

新虫 (著名写手)

【乱世的奸雄】
请问tot_magnetization的取值是如何来选取的？

发自小木虫IOS客户端
赞一下 回复此楼
8楼2019-09-05 19:37:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by 迷迷物理 at 2018-04-01 09:06:38
请教我算了一个162个原子的LDA+U体系，可是显示atomic wfc used for LDA+U projector are not orthogonalized是什么原因，怎么修改？

这个不用管，这个指的是做DFT+U需要将能带投影到一组局域轨道，这组局域轨道可以正交也可以不正交。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
9楼2019-09-07 05:31:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 c_chi 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭