请教下输入文件，结构不收敛

dcckin

我算的是一个蛋白质，算了几何结构一直不收敛，SFC已收敛，希望各位老师给点意见。

! BP SV(P) SV/J opt TightSCF D3 grid5 nofinalgrid VerySlowConv LSHIFT KDIIS scfconv6 xyzfile SmallPrint PDBFILE Pal4
%geom
  optimizeHydrogens true # do not constrain any hydrogen position
end
%basis NewGTO Fe "TZVPP" end
       NewAuxGTO Fe "TZV/J" end
       end
%scf  
          maxiter 500  

    end

* xyz 1 5

  Fe  7.41453763072420      6.36079858202043      6.21309620534936
  Fe  5.25780583407153      6.12931094853904      5.06591368242220
  C   5.57587347803631      4.71035482335136      3.85083367859818
  N   5.85314853155085      3.82477516916566      3.11218529510779
  C   8.24394691987708      7.99132575751955      5.69892767207861
  N   8.80298552604534      9.00494952241834      5.45276604418698
  Fe  9.75128914176056      2.66960320134394     10.01444946985218
  Fe  9.78627021607296      4.80825932222152      8.75455690166511
  Fe  11.99442813563189      3.82694127772671      9.63639960329682
  Fe  10.86522365814514      2.68374517069882      7.50722687113502
  S   11.71796039056521      4.70113038492783      7.62323347380132
  S   11.65781733774154      1.73068428710476      9.80376858145445
  S   8.65968956468069      2.94595239520095      8.23260127357117
  S   10.22667088358964      4.61554945296494     10.94433124589732
  C   4.47999564823760      5.60901716138331      8.46588543191243
  C   4.94280486745931      7.05778312390661      8.40014243616637
  S   5.61774431462538      7.66029631933295      6.79704885503543
  C   5.63385598083710      4.60050247564093      8.41790542618960
  S   6.42985594604542      4.45116208331788      6.78018070963848
  C   8.59262218473907      5.41752132882892      5.28200769012304
  O   9.42321848068625      4.83503024244028      4.71699896600509
  C   4.87044152967654      7.35477042114921      4.04846124658678
  N   19.43203544306628     10.88051104596565      6.40289662302426
  C   19.95073069834854      9.96927810508087      7.50437140886998
  C   18.63890311518140      9.37146561169458      8.12662401178089
  O   17.56605188065926      9.69676938068335      7.66358358938352
  C   20.82508139020403     10.75393464666739      8.50265270782353
  S   22.21220807804552     11.49185771975139      7.55404472804138
  O   18.75884752560874      8.51435272152234      9.14568423266050
  H   21.16799983855277     10.06306251737541      9.31042837796982
  H   20.18960455124462     11.53268068457781      8.99239637095179
  H   22.55778677060632     12.42767095737542      8.49856518058359
  H   20.53669814475475      9.15691724005861      7.01331556684293
  H   18.36126753119828     10.76098871259596      6.47563266320676
  H   19.71629375324035     11.87561537268005      6.55985301400412
  H   19.76162611203678     10.61493352188856      5.45259572786632
  H   3.92038032496450      5.45555643170304      9.41460329366417
  H   3.74479994046524      5.40286388916705      7.66361269186898
  H   4.11191281048323      7.75372203327993      8.65695202549785
  H   5.74507543065950      7.22743248787538      9.16002986092740
  H   6.41355839368573      4.86320155549274      9.17823878598726
  H   5.27947859832760      3.57714070799335      8.67891413285131
  C   9.74267355736337      7.08348371198658      8.66352672337011
  S   8.36191885026351      6.39486487819754      8.19353606476970
  H   9.95983297638126      8.14289093666352      8.87383667561097
  H   19.69176633169075      8.28279602200360      9.37780560318086
  H   9.13374781724131      2.61571057841744     11.41652801650195
  H   12.21870051183093      2.21544697711727      6.77293278149132
  H   13.40985211309591      3.39045304098121      9.17857275119500




*

ggdh

看了下你的结构，
发现结构不合理。有C或者CH去和铁或者硫配位的。现实中蛋白质里面会有碳化铁这种结构？
结构不合理，scf收敛就难，收敛了优化后的结果也没有意义。
另外你用grid5 要知道默认的grid是2 scfconver 却只有6 基组也很小。这样真的好么、

plus

你算的是两个东西结合在一起的稳定构型，如果找不到合适的相互吸引的作用力，他们会越来越远。
建议：
1.优化中心结构。
2.固定中心结构，选择合适的位置加上底物，优化底物的位置，找到稳定的结合方式。
3.全放松优化。

如果实在结合不上去，分开算吧，用H讨论相对能量。

happyrabbit

可以先用mopac2016初步优化一下，看下结构是否合理，然后再用orca计算