QE开smearing后SCF报错

myzpku

最近需要计算一个加了smearing才能收敛的体系，但发现加了以后QE算任何体系都会出错。
一个简单的水分子的输入文件如下：

1. &CONTROL                                                                                         
   
2.    calculation     = 'scf'                                                                       
   
3.    restart_mode    = 'from_scratch'                                                              
   
4.    outdir          = '.'                                                                        
   
5.    pseudo_dir      = '/home/dveracruz/SSSP'                                                      
   
6.    prefix          = 'H2O'                                                                       
   
7.    etot_conv_thr   = 5.0d-5                                                                     
   
8.    forc_conv_thr   = 5.0d-4                                                                     
   
9.    nstep           = 1                                                                           
   
10.    tstress         = .false.                                                                     
    
11.    tprnfor         = .false.                                                                     
    
12. /                                                                                               
    
13. &SYSTEM                                                                                         
    
14.    ibrav           = 1                                                                           
    
15.    celldm(1)       = 5                                                                           
    
16.    nat             = 3                                                                           
    
17.    ntyp            = 2                                                                           
    
18.    ecutwfc         = 30                                                                          
    
19.    ecutrho         = 240                                                                        
    
20.    occupations     = 'smearing'                                                                  
    
21.    degauss         = 0.01                                                                        
    
22.    smearing        = 'gaussian'                                                                  
    
23.                                                                                                 
    
24. /                                                                                                
    
25. &ELECTRONS                                                                                      
    
26.    electron_maxstep = 200                                                                        
    
27.    conv_thr        = 1.0d-5                                                                     
    
28.    mixing_mode     = 'plain'                                                                     
    
29.    mixing_beta     = 0.7                                                                        
    
30.    scf_must_converge = .false.                                                                  
    
31.    startingpot     = 'atomic'                                                                    
    
32. /                                                                                               
    
33. ATOMIC_SPECIES                                                                                   
    
34.     O   15.999    O.pbe-n-kjpaw_psl.0.1.UPF                                                      
    
35.     H    1.00784  H.pbe-rrkjus_psl.1.0.0.UPF                                                     
    
36. ATOMIC_POSITIONS angstrom                                                                        
    
37. O       0.0000000000     0.0000000000     0.1165561912                                          
    
38. H      -0.7635241831     0.0000000000    -0.4687701356                                          
    
39. H       0.7635241831     0.0000000000    -0.4687701356                                          
    
40. K_POINTS automatic                                                                              
    
41.   1 1 1   1 1 1  

复制代码

输出文件末出现错误如下：

1. ===================================================================================
   
2. =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
   
3. =   RANK 0 PID 1312620 RUNNING AT COS-CH-E058916
   
4. =   KILLED BY SIGNAL: 9 (Killed)
   
5. ===================================================================================
   
6. 
   
7. ===================================================================================
   
8. =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
   
9. =   RANK 1 PID 1312621 RUNNING AT COS-CH-E058916
   
10. =   KILLED BY SIGNAL: 9 (Killed)
    
11. ===================================================================================
    
12. 
    
13. ===================================================================================
    
14. =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
    
15. =   RANK 2 PID 1312622 RUNNING AT COS-CH-E058916
    
16. =   KILLED BY SIGNAL: 11 (Segmentation fault)
    
17. ===================================================================================
    
18. 
    
19. ===================================================================================
    
20. =   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
    
21. =   RANK 3 PID 1312623 RUNNING AT COS-CH-E058916
    
22. =   KILLED BY SIGNAL: 9 (Killed)
    
23. ===================================================================================

复制代码

如果不用smearing, 计算没有出现问题，不知道是不是输入文件中smearing的设置有误。