24小时热门版块排行榜    

查看: 3441  |  回复: 8

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
回复此楼
已阅   回复此楼   关注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客户端
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:22 by o3EMxHUAQDtS
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
请填处理意见