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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » Co3O4的LDA+U计算报错求助
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 3242  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 QiMing7 2026-03-25 5/250 2026-03-29 21:10 by 唐沐儿
[考研] 296求调剂 +6 彼岸t 2026-03-29 6/300 2026-03-29 20:22 by 永字号
[考研] 化学工程085602 305分求调剂 +22 RichLi_ 2026-03-25 22/1100 2026-03-29 19:51 by 无际的草原
[考研] 329求调剂 +10 钮恩雪 2026-03-25 10/500 2026-03-29 13:32 by peike
[考研] 070300化学354求调剂 +6 101次希望 2026-03-28 6/300 2026-03-29 12:57 by 无际的草原
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 2026年华南师范大学欢迎化学，化工，生物，生医工等专业优秀学子加入！ +3 llss0711 2026-03-28 6/300 2026-03-29 10:26 by llss0711
[考研] 0856求调剂 +13 zhn03 2026-03-25 14/700 2026-03-29 08:13 by fmesaito
[考研] 330分求调剂 +5 qzenlc 2026-03-29 5/250 2026-03-29 07:37 by 无际的草原
[考研] 322求调剂 +7 宋明欣 2026-03-27 7/350 2026-03-28 21:27 by sanrepian
[考研] 本科新能源科学与工程，一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 085600，材料与化工321分求调剂 +9 大馋小子 2026-03-28 9/450 2026-03-28 14:56 by 神马都不懂
[考研] 322求调剂 +5 旧吢 2026-03-24 5/250 2026-03-28 13:26 by Iveryant
[考研] 一志愿南京航空航天大学材料学硕求调剂 +3 @taotao 2026-03-28 3/150 2026-03-28 10:26 by JourneyLucky
[考研] 085602 307分 求调剂 +7 不知道叫什么！ 2026-03-26 7/350 2026-03-28 09:57 by 神马都不懂
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
[考研] 317求调剂 +7 蛋黄咸肉粽 2026-03-26 7/350 2026-03-27 02:29 by fmesaito
[考研] 281求调剂 +6 Koxui 2026-03-24 7/350 2026-03-26 15:37 by 无际的草原
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
信息提示
关闭
请填处理意见
关闭