24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2178  |  回复: 4

xueke227

新虫 (初入文坛)

[求助] SIESTA几何优化不收敛是不是就说明了这个结构是不稳定的?

RT:算得是一个BN sheet  掺一个Ti原子 后 再吸附一个氢气分子的模型
     算了两次  没有考虑自旋没有收敛 老师让加入自旋再算一次 还是没有收敛
     总分子数:55 个  Ti原子位于BN  sheet的 H 位
     上脚本:


SystemName             BNSheetGGA     # Descriptive name of the system

SystemLabel            1H2TiH  # Short name for naming files



NumberOfAtoms       53

NumberOfSpecies     4



%block ChemicalSpeciesLabel

  1  5   B

  2  7   N

  3  22  Ti

  4  1   H   

%endblock ChemicalSpeciesLabel



XC.functional  GGA

XC.authors     PBE



#----------------------------------------------------------------

LatticeConstant    1 Ang   # Lattice constant alat

%block LatticeVectors               

12.55    0.0    0.0

0.0      12.55  0.00

0.0      0.0    15.00    # c=15 A for the vacuum space

%endblock LatticeVectors

#---------------------------------------------------------------

#---------------------------------------------------------------

%block kgrid_Monkhorst_Pack

1 0 0   0.0

0 1 0   0.0

0 0 1   0.0

%endblock kgrid_Monkhorst_Pack

#-----------------------------------------------------------------------

#-----------------------------------------------------------------------

# SCF options

MaxSCFIterations      100      # Maximum number of SCF iter

DM.NumberPulay         5            # One Pulay every 5 iter

DM.MixingWeight       0.25          # New DM amount for next SCF cycle

DM.Tolerance          1.d-4         # Tolerance in maximum difference

#DM.UseSaveDM          true            # to use continuae..ion files

#SolutionMethod        diagon        # OrderN or Diagon



MeshCutoff           150. Ry         # This data is powered by MS

OccupationFunction   FD

ElectronicTemperature  250 meV       # Temp. for Fermi smearing

# MD

MD.TypeOfRun         broyden

MD.NumCGsteps        16          #Please Pay Enough ATTENTION on this!!!

MD.MaxForceTol       0.05 eV/Ang

#Spined

SpinPolarized T

#------------------------------------------------------------------------

#OUTPUT OPTION

LongOutput   T

AtomicCoordinatesFormat    Ang  

%block AtomicCoordinatesAndAtomicSpecies <1H2TiH.xyz
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ifmc1234

木虫 (著名写手)

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-05-02 19:59:46
xueke227(金币+1): 2011-05-03 14:39:43
首先不知道你这个体系是否有周期性,假如有,在这块
%block kgrid_Monkhorst_Pack
1 0 0   0.0
0 1 0   0.0
0 0 1   0.0
%endblock kgrid_Monkhorst_Pack
增加K点。
第二,对于吸附一个氢气分子的模型
个人觉得除非你用vwd泛函(siesta有development版本),否则考虑的方法会有问题,方法没选对,结果就没多大的意义。
2楼2011-05-02 17:11:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangsong1016

铜虫 (初入文坛)

【答案】应助回帖


youzhizhe(金币+1): 谢谢交流。 2011-05-03 14:43:25
有没有尝试过MD.TypeOfRun CG? 我很早之前用broyden做过简单的尝试,好像不是那么容易收敛,CG相比较容易收敛一些。但是两种算法的差别我还是不明白,希望能有高人指教。到底其中的物理思想是什么?
3楼2011-05-03 14:06:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xueke227

新虫 (初入文坛)

引用回帖:
Originally posted by ifmc1234 at 2011-05-02 17:11:47:
首先不知道你这个体系是否有周期性,假如有,在这块
%block kgrid_Monkhorst_Pack
1 0 0   0.0
0 1 0   0.0
0 0 1   0.0
%endblock kgrid_Monkhorst_Pack
增加K点。
第二,对于吸附一个氢气分子的模型
个 ...

1.周期性是有的 1X1X1的k点只是为了先粗略的判定一下收敛情况 这样想对不对?
2.如果说不收敛是不是就可以认为这组模型是不稳定的?
3.VWD泛函  这个我真的不知道
  计算方法我们用的就是GGA+PBE 有问题?
4楼2011-05-03 14:47:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ifmc1234

木虫 (著名写手)

【答案】应助回帖

xueke227(金币+2): 2011-05-04 15:33:11
引用回帖:
Originally posted by xueke227 at 2011-05-03 14:47:02:
1.周期性是有的 1X1X1的k点只是为了先粗略的判定一下收敛情况 这样想对不对?
2.如果说不收敛是不是就可以认为这组模型是不稳定的?
3.VWD泛函  这个我真的不知道
  计算方法我们用的就是GGA+PBE 有问题?

理论上没问题,但是有时k点太少,也影响收敛。
还有3楼说的改为CG,用Pulay方法而不是broyden,虽然很多时候broyden收敛要快些。
假如想先判断收敛情况,也可以用harris_functional这个方法算一下。
考虑到分子之间的范德瓦尔斯作用,gga是不够的。
5楼2011-05-03 17:33:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xueke227 的主题更新
信息提示
请填处理意见