24小时热门版块排行榜    

查看: 1371  |  回复: 3

凯旋1029

铁虫 (初入文坛)

[求助] 在CP运行中关于Error in routine ortho (1): ortho went bananas报错

我在运行cp.x过程中,遇到了“ortho went bananas”这个报错,弄了很久都没找到问题在哪里,有人能帮我看一下么?下面是我的输入文件:
  &CONTROL                                                                     
  calculation = 'cp'                                   
  dt =   2.0                                          
  iprint = 10                                                   
  isave = 100                                               
  ndr = 52                                               
  ndw = 53                                                
  nstep = 10000                                         
  outdir = './out/'                              
  prefix = 'ABX3wfopt'                                 
  restart_mode = 'restart'                       
  verbosity = 'high'                     
  wf_collect = .false.
  ekin_conv_thr=1.e-5,
  etot_conv_thr=1.e-7,
  forc_conv_thr=1.e-5,
  pseudo_dir='./'                                             
/                                                   
&SYSTEM                                               
  ecutrho =   120.0                                    
  ecutwfc =   30.0                     
  ibrav = 7         
  celldm(1)=33.1948  
  celldm(3)=1.4459                                                      
  nat = 96                                       
  ntyp = 5                                                        
  nr1b = 24                                    
  nr2b = 24
  nr3b = 24                                                              
  /                                          
&ELECTRONS                                               
  electron_damping =   0.1                           
  electron_dynamics = 'damp'                        
  emass =   400                                      
  emass_cutoff =   3.0                              
  electron_temperature = 'not_controlled'
/                                                         
&IONS                                          
  ion_dynamics = 'damp'                       
  ion_damping=0.02                           
  ion_temperature='not_controlled'               
/                                                     
ATOMIC_SPECIES                 
Pb 207.2 Pb.pbe-dn-rrkjus_psl.0.2.2.UPF              
I 126.9 I.pbe-n-rrkjus_psl.0.2.UPF                  
C 12.0 C.pbe-n-rrkjus_psl.0.1.UPF                    
N 14.0 N.pbe-n-rrkjus_psl.0.1.UPF         
H 1.0 H.pbe-rrkjus_psl.0.1.UPF                                                            
ATOMIC_POSITIONS angstrom                          
Pb       4.2403    4.7040    6.2142              
Pb      13.0233    4.7040    6.2142              
Pb       4.2403   13.4870    6.2142              
Pb      13.0233   13.4870    6.2142              
Pb       4.2403    4.7040   18.9132              
Pb      13.0233    4.7040   18.9132              
Pb       4.2403   13.4870   18.9
谢谢啦
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

凯旋1029

铁虫 (初入文坛)

有没有大神看一下
2楼2016-11-03 08:54:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪中谷雨

新虫 (初入文坛)

把dr的值改小些试试,这是不收敛

发自小木虫Android客户端
3楼2017-03-31 10:30:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪中谷雨

新虫 (初入文坛)

dt或者是ortho增加的方向相反,你试试

发自小木虫Android客户端
4楼2017-03-31 10:33:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 凯旋1029 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +5 Kittylucky 2026-08-27 5/250 2026-08-27 12:10 by 姜太公烤鱼
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +5 苏东坡二世 2026-08-23 5/250 2026-08-27 10:39 by 炜雨
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见