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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » SIESTA scf不收敛，怎么办？
51/1返回列表
查看: 2370  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

smartly

金虫 (小有名气)

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

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

[ Last edited by smartly on 2011-2-26 at 14:37 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-02-25 15:59:42
已阅   回复此楼   关注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.你可以试试看，我不保证能结果收敛。
赞一下(2人) 回复此楼
5楼2011-03-04 10:55:33
已阅   回复此楼   关注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 ]

输出文件贴出来看看啊，也好大家重复一下，看看有问题没有。
赞一下(1人) 回复此楼
2楼2011-02-28 21:50:27
已阅   回复此楼   关注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的回帖

yashumimi

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
不收敛指的是500步还没有算完么？
还是出现错误了？
赞一下(1人) 回复此楼
4楼2011-03-03 15:00:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭