24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2314  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博自荐 +6 Linxia林夏 2026-04-13 6/300 2026-04-20 19:26 by 想申博!
[考研] 320求调剂 +6 深郊akm 2026-04-17 6/300 2026-04-20 18:57 by fs26jie
[教师之家] 又一批高校组建人工智能学院 师资行吗 不是骗人吗 +4 yexuqing 2026-04-19 4/200 2026-04-20 14:47 by brantleo
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 5/250 2026-04-20 10:47 by YuY66
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 通信工程求调剂!!! +7 zlb770521 2026-04-14 7/350 2026-04-19 20:56 by Equinoxhua
[考研] 297,工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 085404 22408 309分求调剂 +10 lzmk 2026-04-14 11/550 2026-04-19 16:42 by 中豫男
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[考研] 一志愿中科大材料与化工,353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9,生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
信息提示
请填处理意见