24小时热门版块排行榜    

查看: 2369  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

smartly

金虫 (小有名气)

[交流] SIESTA scf不收敛,怎么办? 已有4人参与

计算的石墨烯,相同设置下,4x4,3x3,都可以,但是计算到5x5,就不收敛了。。。

[ Last edited by smartly on 2011-2-26 at 14:37 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

smartly

金虫 (小有名气)

引用回帖:
Originally posted by wangsong1016 at 2011-02-28 21:50:27:
输出文件贴出来看看啊,也好大家重复一下,看看有问题没有。

SystemName         g55v                                               
SystemLabel        g55v                                               
NumberOfAtoms      49                                                
NumberOfSpecies     1                                                
%block ChemicalSpeciesLabel                                          
1   6   C                                                            
%endblock ChemicalSpeciesLabel                                       
LatticeConstant       1 Ang                                          
%block LatticeParameters                                             
12.30  12.30    10.000   90.0000   90.0000   120.0000                 
%endblock LatticeParameters                                          
                                                                     
%block kgrid_Monkhorst_Pack                                          
    20   0   0  0.5                                                   
    0   20   0  0.5                                                   
    0   0   1  0.5                                                   
%endblock kgrid_Monkhorst_Pack                                       
PAO.BasisSize       DZP                                               
PAO.EnergyShift     50 meV                                            
MeshCutoff          250.0 Ry                                          
MaxSCFIterations    500                                               
DM.MixingWeight      0.01                                             
DM.NumberPulay       9                                                
DM.Tolerance         1.d-4                                            
MD.MaxForceTol      0.02 eV/Ang                                       
MD.MaxStressTol    0.02 GPa                                          
MD.VariableCell     T                                                
SolutionMethod       diagon                                          
SpinPolarized    T                                                   
xc.functional     GGA     # Default value                             
xc.authors        PBE     # Default value                             
                                                                     
MD.TypeOfRun         CG                                               
MD.NumCGsteps         5000                                            
AtomicCoordinatesFormat Ang                                          
%block AtomicCoordinatesAndAtomicSpecies                              
0.013989001    0.007458892    0.088856371   1                        
1.243456925    0.717911479    0.080748034   1                        
2.473989001    0.007458892    0.088856371   1                        
3.703456925    0.717911479    0.080748034   1                        
4.933989001    0.007458892    0.088856371   1                        
6.163456925    0.717911479    0.080748034   1                        
7.393989001    0.007458892    0.088856371   1                        
8.623456925    0.717911479    0.080748034   1                        
9.853989001    0.007458892    0.088856371   1                        
11.083456925    0.717911479    0.080748034   1                        
-1.216010999    2.137881385    0.088856371   1                        
0.013456925    2.848333972    0.080748034   1                        
1.243989001    2.137881385    0.088856371   1                        
2.473456925    2.848333972    0.080748034   1                        
3.703989001    2.137881385    0.088856371   1                        
4.933456925    2.848333972    0.080748034   1                        
6.163989001    2.137881385    0.088856371   1                        
7.393456925    2.848333972    0.080748034   1                        
8.623989001    2.137881385    0.088856371   1                        
9.853456925    2.848333972    0.080748034   1                        
-2.446010999    4.268303878    0.088856371   1                        
-1.216543075    4.978756465    0.080748034   1                        
0.013989001    4.268303878    0.088856371   1                        
1.243456925    4.978756465    0.080748034   1                        
2.473989001    4.268303878    0.088856371   1                        
4.933989001    4.268303878    0.088856371   1                        
6.163456925    4.978756465    0.080748034   1                        
7.393989001    4.268303878    0.088856371   1                        
8.623456925    4.978756465    0.080748034   1                        
-3.676010999    6.398726372    0.088856371   1                        
-2.446543075    7.109178959    0.080748034   1                        
-1.216010999    6.398726372    0.088856371   1                        
0.013456925    7.109178959    0.080748034   1                        
1.243989001    6.398726372    0.088856371   1                        
2.473456925    7.109178959    0.080748034   1                        
3.703989001    6.398726372    0.088856371   1                        
4.933456925    7.109178959    0.080748034   1                        
6.163989001    6.398726372    0.088856371   1                        
7.393456925    7.109178959    0.080748034   1                        
-4.906010999    8.529148865    0.088856371   1                        
-3.676543075    9.239601452    0.080748034   1                        
-2.446010999    8.529148865    0.088856371   1                        
-1.216543075    9.239601452    0.080748034   1                        
0.013989001    8.529148865    0.088856371   1                        
1.243456925    9.239601452    0.080748034   1                        
2.473989001    8.529148865    0.088856371   1                        
3.703456925    9.239601452    0.080748034   1                        
4.933989001    8.529148865    0.088856371   1                        
6.163456925    9.239601452    0.080748034   1                        
%endblock AtomicCoordinatesAndAtomicSpecies
3楼2011-03-03 14:36:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

wangsong1016

铜虫 (初入文坛)


zzy870720z(金币+1): 谢谢参与交流 2011-03-01 08:29:04
引用回帖:
Originally posted by smartly at 2011-02-25 15:59:42:
计算的石墨烯,相同设置下,4x4,3x3,都可以,但是计算到5x5,就不收敛了。。。

[ Last edited by smartly on 2011-2-26 at 14:37 ]

输出文件贴出来看看啊,也好大家重复一下,看看有问题没有。
2楼2011-02-28 21:50:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yashumimi

银虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
不收敛指的是500步还没有算完么?
还是出现错误了?
4楼2011-03-03 15:00:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsong1016

铜虫 (初入文坛)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gzqdyouxia(金币+2): 谢谢提示 2011-03-04 11:04:43
引用回帖:
Originally posted by smartly at 2011-03-03 14:36:17:
SystemName         g55v                                               
SystemLabel        g55v                                               
NumberOfAtoms      49                               ...

50个原子,选择20*20的K点,计算起来很慢啊,你用并行计算吗?感觉应该不用那么多的K点。我没用并行,所以重复你的计算时间太长了,没有算,估计问题出在你的MixingWeight选择的太小了,0.01. 这一项用来保证电子迭代的收敛。我计算一般都用到0.3或者0.4.你可以试试看,我不保证能结果收敛。
5楼2011-03-04 10:55:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见