Выбор базиса

[Darvin56]

Добрый день!
Уважаемые ученые и умники стоит вот какой вопрос. Необходимо выбрать базис для тяжелых элементов в комплексах с органикой.
Т.е. сами комплексы уже имеются и экспериментально изучаются. Хотим сравнить экспериментальные данные с теоретическими(модельными), а именно посчитать колебательные спектры и т.д.. Но как я понял, нужно сначала молекулу оптимизировать в программе, а уже потом ставить "hessian". Программа "FireFly". Метод DFT( дальше МКССП планирую).
Металлы которые мы используем: Zn, Co, Cd, Pb
Покапавшись в литературе с первыми двумя не возникло вопросов, большинство берут 6-31G/B3LYP.
Для кадмия 6-21G, а вот для свинца проблема.
Скажите пожалуйста, какой базис Вы берете для своих исследований, для каждого металла?

[Shorku]

6-31G- это практически самый минимум, его берут, вероятно, не от хорошей жизни.

[Darvin56]

Siegbahn например очень хорошо остановился на этом базисе. Я наверно ученный 1 уровня) поэтому и спрашиваю тех кто очень хорошо разбирается. Посоветуйте.

[alxyppv]

Идете на сайт https://bse.pnl.gov/bse/portal. Выбираете нужные Вам элементы, нужный формат - и вперед. Лучше если Вы будете использовать один тип базисов для всех элементов (если это возможно, конечно). Например, в Вашем случае я бы остановился на семействе
def2-TZVP ( ну или для начала более "легкий" аналог def2-SVP). Эти базисы есть для всех Ваших элементов (обратите внимание, что для Pb и Cd используюся остовные потенциалы, aka effective core potentials, ECP).

[Гесс]

+1 к alxyppv
6-31G может хватить на описание органики (а может и нехватить). Его наверняка не хватит для металла.
def2-TZVP (особенно с дальнейшей накруткой PP, PD, PPD) может хватить для металла. Но молекулу какого размера ТС вытянет полностью на этом базисе? Атомов 50-60 если у него доступ к хорошему кластеру? Атомов 20 если нет? Это среднепотолочные оценки, но погрешность скажем раза в 2 максимум.
Да, скрещивание базисов иногда приводит к неприятным результатам, но это вполне распространенный подход. LANL2DZ один из наиболее популярных выборов для металла.
Мультиреференсная идея - это хорошо, ДФТ может вообще никак не описывать металлокомплексы (всмысле описывать предельно лажево). Подробнее следует читать посты на форуме и ЖЖ sanya1024 она в мультиреференсных металлокомплексах сьела не одну собаку.

[alxyppv]

Ну в принципе да, при не очень больших ресурсах def2-TZVP может оказаться дороговат, но тогда какой вообще смысл говорить про мультиреференсные методы? там все равно лучше большие базисы использовать...
если def2-TZVP все же слишком дорого, можно взять def2-SVP для органической части, и def2-TZVP для металла. Для качественных расчетов в рамках DFT это уже нормально. Если же железо все-таки слабовато, то можно попровать софт с разного рода ускорялками. Orca с гибридными функционалами очень неплохо работает (опция RIJCOSX ускоряет в разы) + есть свои базисы для почти всех елементов с релятивисткими поправками (ZORA, DKH)

[Darvin56]

def2-SVP(TZVP) это ECP базис? Или только для конкретных элементов используется ECP? Просто мне кажется сложновато ECP для свинца будет делать. Но в целом я Вас услышал.
И еще в общем можете прикинуть на вскидку какой по мощности нужен компьютер для данного базиса на систему из 20 молекул и еще отдельно 80?
Сколько памяти он будет использовать?
Для 6-31 я уже научился примерно вычислять.
И еще один вопрос, это все вы указываете для функционала B3LYP?

[Гесс]

def2-SVP(TZVP) это ECP базис? Или только для конкретных элементов используется ECP? Просто мне кажется сложновато ECP для свинца будет делать. Но в целом я Вас услышал.

ECP начиная с криптона. bse.pnl куда вам дали ссылку выше дает базис для любого из элементов. В файерфлай разве карлсруевские базисы не вшиты?

И еще в общем можете прикинуть на вскидку какой по мощности нужен компьютер для данного базиса на систему из 20 молекул и еще отдельно 80?
Сколько памяти он будет использовать?

Наверное все же атомов, а не молекул. Вы указываете ВСЕ или неводородные? В оценке железа я не мастак, но учтите что стоимость DFT растет как куб от числа функций. От есть ооочень приближенно увеличение системы в 4 раза приведет к удорожанию расчета в 64 раза.

Для 6-31 я уже научился примерно вычислять.

Насколько дороже def2-SVP по сравнению с 6-31G - посмотрите число функций в каждом базисе, среднепотолочно прикинете.

И еще один вопрос, это все вы указываете для функционала B3LYP?

ИМХО гибридные функционалы различаются по стоимости менее чем в разы. "чистые" намного дешевле, "двойные" намного дороже, а так... B3LYP, BHHLYP, PBE0... примерно одни и те же яйца. CAM-B3LYP чуть подороже, но опять же не критично. Миннесота регулярно выпускает функционалы для разных задач, [...]
P.S. а вот он, их свежак http://arxiv.org/ftp/arxiv/papers/1212/1212.0944.pdf (когда нибудь этот линк устареет, это The quest for a universal density functional: The accuracy of density functionals across a broad spectrum of databases in chemistry and physics)

[Darvin56]

У меня в выдаче появилась вот такая ошибка. Расчет идет.

THE POINT GROUP OF THE MOLECULE IS C1
THE ORDER OF THE PRINCIPAL AXIS IS 0
*** WARNING! ATOM 1 SHELL 1 TYPE S HAS NORMALIZATION 1.22471038
*** WARNING! ATOM 1 SHELL 2 TYPE S HAS NORMALIZATION 0.98072600
*** WARNING! ATOM 1 SHELL 3 TYPE S HAS NORMALIZATION 0.96315922
*** WARNING! ATOM 1 SHELL 6 TYPE P HAS NORMALIZATION 0.96635290
*** WARNING! ATOM 1 SHELL 7 TYPE P HAS NORMALIZATION 0.97852202
*** WARNING! ATOM 1 SHELL 9 TYPE D HAS NORMALIZATION 1.19016493

Молекула у меня не симметрична

ECP начиная с криптона. bse.pnl куда вам дали ссылку выше дает базис для любого из элементов. В файерфлай разве карлсруевские базисы не вшиты?

• Безымянный.png
  Это все что в него входит. Но я понял как их добавлять.

Вообщем, для своего расчета я взял def2-TZVP (металл)- def2-SVP (органика)/B3LYP 19 молекул, пока машина тянет. Завтра отпишусь

[Darvin56]

Извините, стыдно спрашивать. А точно def2-SVP вкупе с def2-TZVP будет больше.
Сравнивал их

Код: Выделить всё

!  6-311G  EMSL  Basis Set Exchange Library   27/01/15 3:59
! Elements                             References
! --------                             ----------
! H, Li - Ne: R. Krishnan, J.S. Binkley, R. Seeger and J.A. Pople,
!             J. Chem. Phys. 72, 650 (1980)
! Na - Ar:    A.D. McLean and G.S. Chandler J. Chem. Phys. 72, 5639, (1980).
! K  - Ca:    J-P. Blaudeau, M. P. McGrath, L.A. Curtiss and L. Radom,
!             J. Chem. Phys. 107, 5016 (1997).
! Ga - Kr:    L. A. Curtiss, M. P. McGrath, J-P. Blandeau, N. E. Davis, 
!             R. C. Binning, Jr. L. Radom, J. Chem. Phys. 103, 6104 (1995).
! I      :    M.N. Glukhovstev, A. pross, M.P. McGrath, L. Radom, J. Chem. Phys.
!             103, 1878 (1995)
!   


$DATA
CARBON
S   6
  1   4563.2400000              0.00196665       
  2    682.0240000              0.0152306        
  3    154.9730000              0.0761269        
  4     44.4553000              0.2608010        
  5     13.0290000              0.6164620        
  6      1.8277300              0.2210060        
L   3
  1     20.9642000              0.1146600              0.0402487        
  2      4.8033100              0.9199990              0.2375940        
  3      1.4593300             -0.00303068             0.8158540        
L   1
  1      0.4834560              1.0000000              1.0000000        
L   1
  1      0.1455850              1.0000000              1.0000000        
$END

Код: Выделить всё

!  Def2-SVP  EMSL  Basis Set Exchange Library   27/01/15 4:00
! Elements                             References
! --------                             ----------
! H He Li Be B C N O F Ne Na Mg Al Si P S Cl Ar K Ca Sc Ti V Cr Mn Fe Co Ni Cu Zn Ga Ge As Se Br Kr Rb Sr Y Zr Nb Mo Tc Ru Rh Pd Ag Cd In Sn Sb Te I Xe Cs Ba La Hf Ta W Re Os Ir Pt Au Hg Tl Pb Bi Po At Rn : F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys., Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence
! quality for H to Rn: Design and assessment of accuracy 7, 3297 (2005).
! 



$DATA
CARBON
S   5
  1   1238.4016938              0.54568832082E-02      
  2    186.29004992             0.40638409211E-01      
  3     42.251176346            0.18025593888    
  4     11.676557932            0.46315121755    
  5      3.5930506482           0.44087173314    
S   1
  1      0.40245147363          1.0000000        
S   1
  1      0.13090182668          1.0000000        
P   3
  1      9.4680970621           0.38387871728E-01      
  2      2.0103545142           0.21117025112    
  3      0.54771004707          0.51328172114    
P   1
  1      0.15268613795          1.0000000        
D   1
  1      0.8000000              1.0000000        
$END

Даже моя машина 6-311 долго переваривает, а со вторым как-то больше итераций считает за то же время. Но это пока только предварительные результаты.

[alxyppv]

def2-SVP несколько меньше чем 6-311, но в просто 6-311 нет поляризационных функций (т.е. надо добавлять, чтобы было 6-311G**). Но если Вам надо расширить базис, то от def2-SVP естественый путь к def2-TZVP

[Гесс]

А точно def2-SVP вкупе с def2-TZVP будет больше.
Сравнивал их

Код: Выделить всё

!  6-311G  

Код: Выделить всё

!  Def2-SVP 

Даже моя машина 6-311 долго переваривает, а со вторым как-то больше итераций считает за то же время. Но это пока только предварительные результаты.

Пожалуйста задавайте вопросы конкренее, мне довольно непонятно "больше" чего?, "вторым" - это которым?
6-311G триплзетабазис, Def2-SVP - даблзета, то есть при прочих очень равных 6-311G БОЛЬШЕ чем Def2-SVP.
В то же время надиффуживание поппловского 6-311G до по прежнему триплзета 6-311++G ОЧЕНЬ сильно удорожит базис.
Что учтено в Def2 базисах уже написал alxyppv. Нужно ли вам упираться в увеличение базиса от дабл до трипл или в диффузные и поляризационные функции - я не знаю.

[Darvin56]

Большое спасибо за советы. Мы уже почти окончательно определились. Еще один вопрос. Какой из базисов будет больше по размеру (и точнее), если считать органику 6-31G или def2-SVP, добавлять функции не требуется. 6-311 слишком большой, очень сильно нагружает систему, а 6-31 думаем хватит. Металл мы берем def2-TZVPD

[Гесс]

Большое спасибо за советы. Мы уже почти окончательно определились. Еще один вопрос. Какой из базисов будет больше по размеру (и точнее), если считать органику 6-31G или def2-SVP, добавлять функции не требуется. 6-311 слишком большой, очень сильно нагружает систему, а 6-31 думаем хватит. Металл мы берем def2-TZVPD

Если металлы def2-TZVPD то все остальное базисом того же типа. То есть def2-SVP. Иначе появление проблем от "стыка" вероятнее. Кроме того если вы накручиваете диффузию на металл, возможно неплохо увеличить и остаток до def2-SVPD. Разумеется станет дороже, но мне это кажется корректнее.

[Darvin56]

А вообще есть какие-нибудь источники на скорую руку об использовании def2? Мы сами этот вопрос еще разберем поподробнее. Просто чтобы было с чего начинать, запустить пока расчеты и не сомневаться.
С нас будет огромное спасибо
Еще вопрос. Не вижу разницу в двух базисах
def2-TZVP

Код: Выделить всё

!  Def2-TZVP  EMSL  Basis Set Exchange Library   28/01/15 8:37
! Elements                             References
! --------                             ----------
! H He Li Be B C N O F Ne Na Mg Al Si P S Cl Ar K Ca Sc Ti V Cr Mn Fe Co Ni Cu Zn Ga Ge As Se Br Kr Rb Sr Y Zr Nb Mo Tc Ru Rh Pd Ag Cd In Sn Sb Te I Xe Cs Ba La Hf Ta W Re Os Ir Pt Au Hg Tl Pb Bi Po At Rn: F. Weigend and R. Ahlrichs, Phys. Chem. Chem. Phys., Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence
! quality for H to Rn: Design and assessment of accuracy 7, 3297 (2005).
! 



$DATA
CADMIUM
S   2
  1     20.000000000           -0.17401196927    
  2     16.309051661            0.37137277194    
S   1
  1      6.4512709901           1.0000000        
S   1
  1      1.5350047390           1.0000000        
S   1
  1      0.68983631523          1.0000000        
S   1
  1      0.13475720977          1.0000000        
S   1
  1      0.47416874527E-01            1.0000000        
P   4
  1     14.000681404            0.69333062793E-01      
  2      8.3094019872          -0.25420103681    
  3      2.2020058122           0.49200980368    
  4      1.0779246137           0.49702118131    
P   1
  1      0.50111879846          1.0000000        
P   1
  1      0.16000000000          1.0000000        
P   1
  1      0.51600000000E-01            1.0000000        
D   4
  1     30.380789793            0.32545123835E-02      
  2     11.474551578           -0.14212074843E-01      
  3      3.0507394903           0.24961756436    
  4      1.3622028524           0.44905635176    
D   1
  1      0.56399409888          1.0000000        
D   1
  1      0.20916686535          1.0000000        
F   1
  1      1.5981300              1.0000000       

и

def2-TZVPD

Код: Выделить всё

!  Def2-TZVPD  EMSL  Basis Set Exchange Library   28/01/15 8:37
! Elements                             References
! --------                             ----------
! H He Li Be B C N O F Ne Na Mg Al Si P S Cl Ar K Ca Sc Ti V Cr Mn Fe Co Ni Cu Zn Ga Ge As Se Br Kr Rb Sr Y Zr Nb Mo Tc Ru Rh Pd Ag Cd In Sn Sb Te I Xe Cs Ba La Hf Ta W Re Os Ir Pt Au Hg Tl Pb Bi Po At Rn: Dmitrij Rappoport and F. Furche, Property-optimized Gaussian basis sets for molecular response calculations, J. Chem. Phys. 133, 134105 (2010).
! 



$DATA
CADMIUM
S   2
  1     20.000000000           -0.17401196927    
  2     16.309051661            0.37137277194    
S   1
  1      6.4512709901           1.0000000        
S   1
  1      1.5350047390           1.0000000        
S   1
  1      0.68983631523          1.0000000        
S   1
  1      0.13475720977          1.0000000        
S   1
  1      0.47416874527E-01            1.0000000        
P   4
  1     14.000681404            0.69333062793E-01      
  2      8.3094019872          -0.25420103681    
  3      2.2020058122           0.49200980368    
  4      1.0779246137           0.49702118131    
P   1
  1      0.50111879846          1.0000000        
P   1
  1      0.16000000000          1.0000000        
P   1
  1      0.51600000000E-01            1.0000000        
D   4
  1     30.380789793            0.32545123835E-02      
  2     11.474551578           -0.14212074843E-01      
  3      3.0507394903           0.24961756436    
  4      1.3622028524           0.44905635176    
D   1
  1      0.56399409888          1.0000000        
D   1
  1      0.20916686535          1.0000000        
F   1
  1      1.5981300000           1.0000000     

Такое бывает? Или это означает, что диффузных функций для данного металла в данном базисе нет?

[Гесс]

Thisiscaster, я вынес ваше сообщение в отдельную тему. viewtopic.php?f=71&t=110563
Давайте не будем смешивать мух с немухами.

[Гесс]

Такое бывает? Или это означает, что диффузных функций для данного металла в данном базисе нет?

Навскидку не скажу, думаю надо почитать Property-optimized Gaussian basis sets for molecular response calculations в котором вводились доаугуменченные def2 базисы.

[Darvin56]

базис def2-SVPD(органика)-def2-TZVPD (Co²⁺)M=4 UHF, DFT

Код: Выделить всё

 ITER EX   TOTAL ENERGY      E CHANGE  DENSITY CHANGE    DIIS ERROR
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0 -2178.411424014 -2178.411424014**************   0.531860489
   2  1  2500.775790104  4679.187214118  75.3413423612607.522547867
   3  2  2470.972888737   -29.802901367  68.2367728192606.141739198
   4  3  2444.567414749   -26.405473988  28.9861729332605.583262099
   5  4  2443.000153303    -1.567261446  54.3948912902610.745711138
   6  5  2432.882086347   -10.118066955 113.7230437622617.517384484
   7  6  2420.699855180   -12.182231167  22.8703506382618.471155116
   8  7  2418.855597803    -1.844257377  15.7511999362622.162007329
   9  8  2416.706241320    -2.149356483   8.1878913622622.308616512
          * * *    REDUCING  DIIS SUBSPACE    * * *
          * * *    REDUCING  DIIS SUBSPACE    * * *
          * * *    REDUCING  DIIS SUBSPACE    * * *
          * * *    REDUCING  DIIS SUBSPACE    * * *
          * * *    REDUCING  DIIS SUBSPACE    * * *
          * * *    REDUCING  DIIS SUBSPACE    * * *
          * * *    REDUCING  DIIS SUBSPACE    * * *
          * * *    REDUCING  DIIS SUBSPACE    * * *
  10  9  2414.861382970    -1.844858350**************2620.719945188
  11 10-24769.897142692-27184.758525663****************************
  12 11-70579.246416051-45809.3492733592573.5121578069292.582958504
  13 12-80953.431520707-10374.185104656 550.8123285448295.462187263
  14 13-83238.867381448 -2285.4358607401417.6763885768160.818387689
  15 14-89615.419881691 -6376.5525002441517.1936511267840.852064858
  16 15-91964.943979913 -2349.5240982211627.4173909787613.673627429
  17 16-94562.386081067 -2597.4421011542410.1277603296705.419791090
  18 17**************** -7480.4506702053788.0131137965701.550067594
  19 18**************** -5341.485806539**************4381.494662128
  20 19**************** -9534.5813407022397.3674376791757.582810861
  21 20**************** -2688.4368425516527.7773631771556.851531400
  22 21**************** -6806.6657513318022.5262555091033.262453443
  23 22**************** -5272.5809316423250.271850864 636.933027385
  24 23**************** -1319.2067561291681.634363094 254.681803401
  25 24****************  -190.676345407 471.519827797 126.097086194
  26 25****************   -24.389744519 142.550991545  50.098544642
  27 26****************    -4.096689355  36.372615573  12.097820090
  28 27****************    -1.200256277  23.801808197   7.087699469
  29 28****************     0.462879384  14.583723759   3.110192757
  30 29****************    -0.426995922   5.615883611   1.459413624
  31 30****************    -0.191903281   2.033611943   0.879668082
  32 31****************    -0.295552940   1.542255933   0.241674740
  33 32****************    -0.051327003   0.360166782   0.152409889
  34 33****************    -0.001587656   0.218897253   0.050625765
  35 34****************    -0.000859472   0.141524208   0.025485486
  36 35****************     0.000049439   0.103121677   0.018226943
  37 36****************    -0.000053454   0.042755291   0.008630700
  38 37****************    -0.000059113   0.015278196   0.003183913
  39 38****************     0.000002665   0.005007992   0.001201602
  40 39****************    -0.000021128   0.002723502   0.000736499
  41 40****************     0.000006888   0.001612797   0.000460845
  42 41****************     0.000033680   0.001058576   0.000275611
  43 42****************    -0.000043952   0.001420988   0.000162235
  44 43****************    -0.000000507   0.000397790   0.000083051
  45 44****************     0.000002513   0.000181939   0.000042273
  46 45****************     0.000023507   0.000075390   0.000036559
  47 46****************    -0.000019775   0.000065841   0.000036744
  48 47****************     0.000047206   0.000067010   0.000036588
  49 48****************    -0.000028090   0.000013025   0.000036973
  50 49****************     0.000020298   0.000039930   0.000036701
  51 50****************    -0.000007695   0.000021235   0.000036763
  52 51****************    -0.000006274   0.000019042   0.000036690
  53 52****************    -0.000001091   0.000023072   0.000036719
  54 53****************    -0.000007395   0.000018497   0.000036623
  55 54****************     0.000023577   0.000010368   0.000036885
  56 55****************    -0.000056295   0.000013216   0.000037664
  57 56****************     0.000027159   0.000012411   0.000036935
  58 57****************    -0.000002387   0.000014147   0.000037242
  59 58****************     0.000008308   0.000026315   0.000036905
  60 59****************     0.000002628   0.000017480   0.000036891
  61 60****************     0.000030776   0.000038681   0.000036808
  62 61****************    -0.000066013   0.000029073   0.000036965
  63 62****************     0.000039446   0.000021328   0.000036723
  64 63****************     0.000013391   0.000026587   0.000036911
  65 64****************    -0.000004821   0.000015052   0.000036915
  66 65****************    -0.000031603   0.000006045   0.000037046
  67 66****************     0.000036527   0.000010229   0.000036831
  68 67****************    -0.000011560   0.000020830   0.000036864
  69 68****************    -0.000008738   0.000021761   0.000036901
  70 69****************    -0.000040187   0.000002159   0.000037186
  71 70****************     0.000045149   0.000024052   0.000036760
  72 71****************     0.000025964   0.000013543   0.000036807
  73 72****************    -0.000013178   0.000035279   0.000036790
  74 73****************    -0.000031154   0.000013802   0.000036882
  75 74****************     0.000010681   0.000011182   0.000037235
  76 75****************     0.000011775   0.000017273   0.000036975
  77 76****************     0.000015574   0.000012385   0.000036918
  78 77****************    -0.000021428   0.000016178   0.000037188
  79 78****************    -0.000002582   0.000013462   0.000036745
  80 79****************    -0.000006504   0.000009886   0.000037016
  81 80****************    -0.000028643   0.000013499   0.000037055
  82 81****************     0.000044696   0.000013703   0.000036917
  83 82****************    -0.000026122   0.000014807   0.000036847
  84 83****************     0.000048989   0.000017511   0.000036834
  85 84****************    -0.000043197   0.000016578   0.000036716
  86 85****************     0.000009921   0.000010982   0.000037176
  87 86****************     0.000012416   0.000014650   0.000036629
  88 87****************    -0.000020412   0.000013867   0.000037053
  89 88****************    -0.000016240   0.000012657   0.000036774
  90 89****************     0.000019776   0.000009967   0.000036900
  91 90****************     0.000012914   0.000026063   0.000036747
  92 91****************    -0.000004135   0.000018425   0.000036903
  93 92****************    -0.000012991   0.000019006   0.000036616
  94 93****************    -0.000016678   0.000017142   0.000036838
  95 94****************     0.000047327   0.000015620   0.000036773
  96 95****************    -0.000031574   0.000007726   0.000036837
  97 96****************     0.000018589   0.000016432   0.000037018
  98 97****************    -0.000026707   0.000012398   0.000036786
  99 98****************     0.000012122   0.000013939   0.000036668
 100 99****************     0.000027630   0.000013681   0.000036655
 101100****************     0.000010266   0.000014760   0.000036932
 102101****************    -0.000011331   0.000014958   0.000036944
 103102****************    -0.000005782   0.000014075   0.000036841
 104103****************    -0.000029057   0.000009932   0.000037161
 105104****************    -0.000004220   0.000008736   0.000037212
 106105****************     0.000014108   0.000015382   0.000037392
 107106****************    -0.000009230   0.000008182   0.000037054
 108107****************     0.000028056   0.000015545   0.000036779
 109108****************     0.000004715   0.000013017   0.000036996
 110109****************    -0.000029766   0.000007617   0.000036889
 111110****************     0.000030150   0.000010019   0.000036659
 112111****************    -0.000018571   0.000011488   0.000036963
 113112****************     0.000026937   0.000017857   0.000036690
 114113****************     0.000003556   0.000023505   0.000036662
 115114****************    -0.000007280   0.000011217   0.000037290
 116115****************     0.000004178   0.000012583   0.000036897
 117116****************    -0.000005140   0.000013394   0.000036695
 118117****************     0.000005620   0.000013009   0.000036817
 119118****************    -0.000021344   0.000025827   0.000037192
 120119****************    -0.000010052   0.000011366   0.000036798
 121120****************     0.000011578   0.000042293   0.000037024
 122121****************    -0.000016704   0.000014370   0.000037179
 123122****************    -0.000003763   0.000025562   0.000036899

Задача не сходится, считает больше 10 часов.
def2-SVP(органика)-def2-TZVP (Co²⁺)M=4 UHF,DFT без диффузных, сходится. Оптимизация за 7 часов.
def2-SVP(органика)-def2-SVPD (Cd²⁺) M=1, UHF(но не важно каким и так и так пробовал) DFT

Код: Выделить всё

 ITER EX   TOTAL ENERGY      E CHANGE  DENSITY CHANGE    DIIS ERROR
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0 -2019.137800649 -2019.1378006494282.610076830  15.387229747
   2  1 -2169.975790785  -150.837990137 464.686148272   1.176917554
   3  2 -2148.534344213    21.441446572  41.690143246   2.475171782
   4  3 -2230.127603320   -81.593259107 130.461454145   2.452464813
   5  4 -2231.858235987    -1.730632667  50.793094746   1.190732925
   6  5 -2238.392677965    -6.534441978  43.495023002   0.827442614
   7  6 -2243.203894050    -4.811216086  33.300841436   0.410961862
   8  7 -2244.584667558    -1.380773507  23.886671526   0.170853788
   9  8 -2244.882319356    -0.297651799   5.550279533   0.039717380
  10  9 -2244.958853795    -0.076534439   1.637260216   0.017068154
  11 10 -2244.980173527    -0.021319732   1.696542478   0.037020780
  12 11 -2244.989360447    -0.009186920   0.480569871   0.037784797
  13 12 -2244.996739582    -0.007379135   0.872685296   0.037533416
  14 13 -2245.014786119    -0.018046538   1.238778244   0.019386577
  15 14 -2245.028380876    -0.013594757   0.885671384   0.017129611
  16 15 -2245.036465074    -0.008084198   0.291406518   0.008835349
  17 16 -2245.038844010    -0.002378937   0.184299117   0.003994219
  18 17 -2245.039675791    -0.000831781   0.104956852   0.001017126
  19 18 -2245.039808491    -0.000132700   0.076387588   0.001225170
  20 19 -2245.039929766    -0.000121274   0.036080932   0.000929674
  21 20 -2245.040053830    -0.000124064   0.034687949   0.000423125
  22 21 -2245.040127897    -0.000074067   0.024840730   0.000350035
  23 22 -2245.040179817    -0.000051920   0.030835541   0.000322506
  24 23 -2245.040224260    -0.000044443   0.039724270   0.000318987
  25 24 -2245.040258600    -0.000034340   0.038790035   0.000407226
  26 25 -2245.040273404    -0.000014804   0.026134900   0.000377539
  27 26 -2245.040295767    -0.000022363   0.026407556   0.000312546
  28 27 -2245.040314935    -0.000019168   0.031235119   0.000181326
  29 28 -2245.040329942    -0.000015008   0.039725597   0.000140277
  30 29 -2245.040342909    -0.000012967   0.039922686   0.000140596
  31 30 -2245.040353578    -0.000010669   0.030167710   0.000128928
  32 31 -2245.040361170    -0.000007593   0.037085912   0.000080600
  33 32 -2245.040366919    -0.000005749   0.046419606   0.000089297
  34 33 -2245.040373064    -0.000006144   0.025778160   0.000119922
  35 34 -2245.040377027    -0.000003963   0.029112000   0.000132088
  36 35 -2245.040379230    -0.000002203   0.013157538   0.000153647
  37 36 -2245.040381470    -0.000002240   0.004999084   0.000093001
  38 37 -2245.040380939     0.000000531   0.003988600   0.000127543
  39 38 -2245.040380519     0.000000420   0.008799338   0.000153222
  40 39 -2245.040381047    -0.000000528   0.010434415   0.000082924
  41 40 -2245.040381394    -0.000000347   0.002174988   0.000063324
  42 41 -2245.040381987    -0.000000593   0.012113600   0.000082225
  43 42 -2245.040382649    -0.000000662   0.012267016   0.000040777
  44 43 -2245.040383350    -0.000000701   0.006482436   0.000039050
  45 44 -2245.040384018    -0.000000668   0.004934900   0.000034063
  46 45 -2245.040384202    -0.000000184   0.002583154   0.000038198
  47 46 -2245.040384338    -0.000000136   0.003505030   0.000056410
  48 47 -2245.040384284     0.000000054   0.001138100   0.000053650
  49 48 -2245.040384456    -0.000000172   0.001925861   0.000065945
  50 49 -2245.040384557    -0.000000101   0.003615905   0.000060482
  51 50 -2245.040384860    -0.000000304   0.002441995   0.000055651
  52 51 -2245.040384784     0.000000076   0.008433958   0.000053188
  53 52 -2245.040385294    -0.000000509   0.023373332   0.000027923
  54 53 -2245.040385770    -0.000000476   0.004824932   0.000069390
  55 54 -2245.040385941    -0.000000171   0.004254182   0.000027775
  56 55 -2245.040386188    -0.000000247   0.007117619   0.000021237
  57 56 -2245.040386699    -0.000000511   0.004943916   0.000027106
  58 57 -2245.040386655     0.000000044   0.002838485   0.000019896
  59 58 -2245.040387070    -0.000000415   0.001939680   0.000020207
  60 59 -2245.040386956     0.000000114   0.001181434   0.000023836
  61 60 -2245.040386955     0.000000001   0.000316155   0.000028386
  62 61 -2245.040386961    -0.000000006   0.000942754   0.000033843
  63 62 -2245.040386978    -0.000000017   0.000677123   0.000039177
  64 63 -2245.040386862     0.000000116   0.001342903   0.000038965
  65 64 -2245.040387016    -0.000000154   0.000799482   0.000028177
  66 65 -2245.040387076    -0.000000061   0.001120910   0.000031571
  67 66 -2245.040387150    -0.000000074   0.001860744   0.000029254
  68 67 -2245.040386928     0.000000223   0.003025971   0.000025226
  69 68 -2245.040387290    -0.000000363   0.002506289   0.000018774
  70 69 -2245.040387643    -0.000000353   0.002052831   0.000023036
  71 70 -2245.040387261     0.000000382   0.001768351   0.000023021
  72 71 -2245.040387475    -0.000000214   0.001040891   0.000025922
  73 72 -2245.040387393     0.000000082   0.002230619   0.000026549
  74 73 -2245.040387546    -0.000000153   0.001939213   0.000029403
  75 74 -2245.040387652    -0.000000106   0.001968971   0.000033594
  76 75 -2245.040387266     0.000000386   0.001212520   0.000037729
  77 76 -2245.040387319    -0.000000052   0.000826253   0.000041770
  78 77 -2245.040387570    -0.000000251   0.000661134   0.000039795
  79 78 -2245.040387348     0.000000221   0.001089014   0.000033466
  80 79 -2245.040387205     0.000000143   0.004780159   0.000029820
  81 80 -2245.040387497    -0.000000292   0.002088056   0.000018791
  82 81 -2245.040387398     0.000000099   0.002135249   0.000013802
  83 82 -2245.040387579    -0.000000181   0.001209298   0.000014974
  84 83 -2245.040387415     0.000000164   0.001347340   0.000017373
  85 84 -2245.040387537    -0.000000122   0.001287456   0.000020730
  86 85 -2245.040387435     0.000000102   0.001052999   0.000018732
  87 86 -2245.040387589    -0.000000154   0.000344597   0.000020787
  88 87 -2245.040387503     0.000000087   0.000262952   0.000023633
  89 88 -2245.040387469     0.000000033   0.001109103   0.000024258
  90 89 -2245.040387919    -0.000000450   0.000864248   0.000021641
  91 90 -2245.040387798     0.000000121   0.000509246   0.000023078
  92 91 -2245.040387251     0.000000548   0.000440654   0.000024901
  93 92 -2245.040387597    -0.000000346   0.001280211   0.000026388
  94 93 -2245.040387930    -0.000000333   0.000809462   0.000032007
  95 94 -2245.040387989    -0.000000059   0.001272214   0.000029993
  96 95 -2245.040387794     0.000000195   0.001868563   0.000029619
  97 96 -2245.040387911    -0.000000117   0.001495070   0.000026831
  98 97 -2245.040387728     0.000000183   0.002457661   0.000026286
  99 98 -2245.040387745    -0.000000017   0.001128622   0.000024064
 100 99 -2245.040387681     0.000000064   0.001006258   0.000025128
 101100 -2245.040388226    -0.000000545   0.001632072   0.000027341
 102101 -2245.040388044     0.000000182   0.002238575   0.000032964
 103102 -2245.040388266    -0.000000222   0.001106560   0.000036323
 104103 -2245.040388429    -0.000000163   0.000706030   0.000041279
 105104 -2245.040388843    -0.000000414   0.001521438   0.000044852
 106105 -2245.040388146     0.000000697   0.000623322   0.000052501
 107106 -2245.040388458    -0.000000312   0.000369929   0.000052637
 108107 -2245.040388294     0.000000164   0.000422177   0.000050382
 109108 -2245.040387985     0.000000310   0.000157434   0.000051356
 110109 -2245.040388107    -0.000000123   0.000919778   0.000051588
 111110 -2245.040388226    -0.000000118   0.000593167   0.000046477
 112111 -2245.040388262    -0.000000036   0.001271648   0.000046383
 113112 -2245.040388329    -0.000000067   0.004008990   0.000055278
 114113 -2245.040388286     0.000000043   0.001624753   0.000037699
 115114 -2245.040388211     0.000000075   0.001453177   0.000060162
 116115 -2245.040388629    -0.000000418   0.001402726   0.000051388
 117116 -2245.040388767    -0.000000138   0.000650476   0.000052852
 118117 -2245.040388815    -0.000000047   0.000417145   0.000056996
 119118 -2245.040389108    -0.000000293   0.001112006   0.000057490
 120119 -2245.040388667     0.000000440   0.000317860   0.000059595
 121120 -2245.040388746    -0.000000079   0.001295463   0.000062918
 122121 -2245.040388716     0.000000030   0.004577713   0.000059406
 123122 -2245.040389337    -0.000000621   0.005175590   0.000072254
 124123 -2245.040390508    -0.000001171   0.002630615   0.000098224
 125124 -2245.040391373    -0.000000865   0.001384607   0.000098903
 126125 -2245.040391172     0.000000201   0.010883100   0.000102567
 127126 -2245.040389006     0.000002166   0.006012027   0.000114762
 128127 -2245.040388270     0.000000737   0.004458700   0.000092548
 129128 -2245.040388729    -0.000000459   0.002468922   0.000092719
 130129 -2245.040388580     0.000000149   0.000960579   0.000044263
 131130 -2245.040389284    -0.000000704   0.003473906   0.000046234
 132131 -2245.040389633    -0.000000350   0.002255046   0.000073532
 133132 -2245.040390556    -0.000000922   0.000828143   0.000069478
 134133 -2245.040390629    -0.000000073   0.002169568   0.000076225
 135134 -2245.040389957     0.000000672   0.003871336   0.000100072
 136135 -2245.040389823     0.000000134   0.003182780   0.000120271
 137136 -2245.040390491    -0.000000668   0.001831630   0.000127654
 138137 -2245.040390282     0.000000209   0.003023007   0.000130871
 139138 -2245.040390606    -0.000000325   0.004394603   0.000131112
 140139 -2245.040390676    -0.000000069   0.004945456   0.000147577
 141140 -2245.040393216    -0.000002540   0.003482251   0.000105572
 142141 -2245.040395142    -0.000001926   0.009547195   0.000113346
 143142 -2245.040394512     0.000000631   0.006461850   0.000130345
 144143 -2245.040393520     0.000000992   0.008733866   0.000228260
 145144 -2245.040392178     0.000001342   0.008262449   0.000187837
 146145 -2245.040391822     0.000000356   0.007454287   0.000147844
 147146 -2245.040392208    -0.000000386   0.005522403   0.000156471
 148147 -2245.040392603    -0.000000395   0.004658897   0.000149534
 149148 -2245.040393220    -0.000000617   0.002407178   0.000129558
 150149 -2245.040393532    -0.000000312   0.010749854   0.000125529
 151150 -2245.040398607    -0.000005075   0.009089070   0.000140070
 152151 -2245.040402437    -0.000003830   0.004001718   0.000184361
 153152 -2245.040400307     0.000002130   0.011935879   0.000274244
 154153 -2245.040396079     0.000004228   0.014940749   0.000288309
 155154 -2245.040392616     0.000003463   0.013741891   0.000235646
 156155 -2245.040390941     0.000001675   0.005579216   0.000142729
 157156 -2245.040390531     0.000000410   0.007152488   0.000138353
 158157 -2245.040391578    -0.000001047   0.002132727   0.000151511
 159158 -2245.040392307    -0.000000729   0.003351428   0.000112427
 160159 -2245.040393439    -0.000001132   0.004502187   0.000120891
 161160 -2245.040395778    -0.000002339   0.005994412   0.000131280
 162161 -2245.040398860    -0.000003082   0.003637602   0.000159430
 163162 -2245.040401360    -0.000002500   0.003027444   0.000175252
 164163 -2245.040400878     0.000000482   0.004512216   0.000216212
 165164 -2245.040399811     0.000001067   0.007509320   0.000275371
 166165 -2245.040399020     0.000000791   0.004745635   0.000255446
 167166 -2245.040399144    -0.000000124   0.002383034   0.000231923
 168167 -2245.040399564    -0.000000420   0.002409586   0.000237547
 169168 -2245.040400859    -0.000001295   0.000995916   0.000248698
 170169 -2245.040401401    -0.000000542   0.011126114   0.000245798
 171170 -2245.040401845    -0.000000444   0.014851504   0.000207072
 172171 -2245.040410284    -0.000008440   0.014641970   0.000229246
 173172 -2245.040414094    -0.000003810   0.002641468   0.000243069
 174173 -2245.040418629    -0.000004535   0.003987883   0.000259127
 175174 -2245.040419114    -0.000000485   0.006775604   0.000262358
 176175 -2245.040421090    -0.000001976   0.004607869   0.000320057
 177176 -2245.040423378    -0.000002288   0.005571939   0.000302949
 178177 -2245.040424049    -0.000000671   0.003858839   0.000324783
 179178 -2245.040424479    -0.000000429   0.016875361   0.000337018
 180179 -2245.040427996    -0.000003517   0.015122705   0.000304017
 181180 -2245.040432510    -0.000004514   0.012268304   0.000316895
 182181 -2245.040437165    -0.000004655   0.014311912   0.000331792
 183182 -2245.040446972    -0.000009808   0.013241329   0.000355016
 184183 -2245.040447461    -0.000000489   0.014039270   0.000376838
 185184 -2245.040461425    -0.000013964   0.014174879   0.000475897
 186185 -2245.040479209    -0.000017784   0.003712471   0.000447814
 187186 -2245.040487727    -0.000008518   0.008609086   0.000473810
 188187 -2245.040499761    -0.000012034   0.005302941   0.000505826
 189188 -2245.040509211    -0.000009450   0.012242171   0.000525262
 190189 -2245.040506069     0.000003142   0.012246856   0.000516414
 191190 -2245.040507446    -0.000001377   0.017517626   0.000520134
 192191 -2245.040505745     0.000001702   0.037975331   0.000502812
 193192 -2245.040512688    -0.000006943   0.015315945   0.000615084
 194193 -2245.040527065    -0.000014378   0.003020583   0.000536205
 195194 -2245.040535246    -0.000008181   0.011799843   0.000565916
 196195 -2245.040532779     0.000002467   0.073149129   0.000573087
 197196 -2245.040597540    -0.000064761   0.024879747   0.000704211
 198197 -2245.040667212    -0.000069672   0.013678232   0.000799135
 199198 -2245.040682096    -0.000014884   0.027595184   0.000814330
 200199 -2245.040708511    -0.000026415   0.036664023   0.000860010
 201200 -2245.040714284    -0.000005773   0.010680999   0.000860273
 202201 -2245.040718672    -0.000004388   0.022016011   0.000888437
 203202 -2245.040721855    -0.000003183   0.015940244   0.000905732
 204203 -2245.040785303    -0.000063447   0.043954412   0.000924425
 205204 -2245.040813159    -0.000027856   0.013035763   0.000943030
 206205 -2245.040869174    -0.000056015   0.048327201   0.001054255
 207206 -2245.040950390    -0.000081215   0.070156025   0.001159097
 208207 -2245.040985671    -0.000035281   0.042804745   0.001184363
 209208 -2245.041071892    -0.000086221   0.034834635   0.001176509
 210209 -2245.041151873    -0.000079982   0.041237430   0.001217258
 211210 -2245.041208931    -0.000057058   0.079227655   0.001292112
 212211 -2245.041102150     0.000106781   0.035740649   0.001219638
 213212 -2245.041274251    -0.000172101   0.114501811   0.001400416
 214213 -2245.041391329    -0.000117079   0.066506000   0.001508857
 215214 -2245.041717387    -0.000326058   0.023775184   0.001640230
 216215 -2245.042024005    -0.000306618   0.054909934   0.001786854
 217216 -2245.042054597    -0.000030591   0.027030710   0.001756589
 218217 -2245.042000924     0.000053673   0.020204322   0.001777562
 219218 -2245.041887497     0.000113427   0.046126680   0.001732961
 220219 -2245.041429633     0.000457864   0.035572647   0.001723520
 221220 -2245.041350123     0.000079509   0.028292873   0.001631664
 222221 -2245.041285576     0.000064548   0.088367294   0.001438271
 223222 -2245.040927448     0.000358128   0.057751676   0.001156318
 224223 -2245.040858318     0.000069130   0.054275217   0.001362588
 225224 -2245.040858696    -0.000000378   0.104798482   0.001184640
 226225 -2245.040874623    -0.000015927   0.048657663   0.000951771
 227226 -2245.040959034    -0.000084411   0.053170698   0.000899880
 228227 -2245.041121197    -0.000162163   0.045059672   0.001094588
 229228 -2245.041311295    -0.000190099   0.056005101   0.001513732
 230229 -2245.041486954    -0.000175659   0.036762924   0.001558137
 231230 -2245.041328160     0.000158794   0.043575400   0.001752166
 232231 -2245.041255461     0.000072699   0.070316006   0.001473679
 233232 -2245.041118019     0.000137442   0.024931904   0.001728672
 234233 -2245.041061022     0.000056997   0.023685086   0.001720218
 235234 -2245.041018208     0.000042814   0.013766572   0.001721287
 236235 -2245.041070405    -0.000052197   0.020518477   0.001565453
 237236 -2245.041191692    -0.000121287   0.096669984   0.001474588
 238237 -2245.041716632    -0.000524940   0.047361197   0.001363959
 239238 -2245.042103242    -0.000386610   0.061859280   0.001528529
 240239 -2245.042287664    -0.000184422   0.018588928   0.001471890
 241240 -2245.042386929    -0.000099265   0.017433514   0.001544898
 242241 -2245.042488682    -0.000101754   0.015487452   0.001626780
 243242 -2245.042572153    -0.000083471   0.069235578   0.001640664
 244243 -2245.042732935    -0.000160781   0.072706478   0.001859646
 245244 -2245.042983131    -0.000250196   0.043042735   0.001792135
 246245 -2245.043165715    -0.000182584   0.038964676   0.001990522
 247246 -2245.043413487    -0.000247772   0.033992660   0.003297410
 248247 -2245.043365358     0.000048129   0.034123746   0.003383877
 249248 -2245.043287174     0.000078184   0.032696506   0.003388859
 250249 -2245.043168950     0.000118224   0.173506556   0.003167959
 251250 -2245.043843520    -0.000674569   0.047596927   0.003083808
 252251 -2245.044345224    -0.000501704   0.062882916   0.003083787
 253252 -2245.044303255     0.000041969   0.041333167   0.003697350
 254253 -2245.043881416     0.000421839   0.058798470   0.004428430
 255254 -2245.044340903    -0.000459487   0.410025827   0.003693449
 256255 -2245.045573418    -0.001232515   0.155215513   0.003314278
 257256 -2245.046369749    -0.000796331   0.079104782   0.003188930
 258257 -2245.047182323    -0.000812574   0.089964112   0.003391194
 259258 -2245.047056818     0.000125505   0.173249687   0.003618546
 260259 -2245.048861665    -0.001804847   0.217234772   0.005533889
 261260 -2245.050163604    -0.001301939   0.307184796   0.005241067
 262261 -2245.051813941    -0.001650336   0.250024099   0.004013510
 263262 -2245.051466748     0.000347193   0.076263721   0.003797312
 264263 -2245.051514664    -0.000047916   0.126353809   0.003814208
 265264 -2245.051816394    -0.000301730   0.126690440   0.005358170
 266265 -2245.051942314    -0.000125920   0.122297301   0.005779845
 267266 -2245.051465321     0.000476993   0.060854300   0.005632603
 268267 -2245.051434994     0.000030327   0.576392792   0.005791080
 269268 -2245.053633547    -0.002198553   0.260598090   0.004087416
 270269 -2245.055807277    -0.002173730   0.626913528   0.004244813
 271270 -2245.060656375    -0.004849098   0.589601268   0.005933210
 272271 -2245.063576314    -0.002919940   0.406092100   0.005917040
 273272 -2245.065137023    -0.001560709   0.744643056   0.004481066
 274273 -2245.066196221    -0.001059198   0.292016206   0.008908255
 275274 -2245.067686528    -0.001490307   0.154394174   0.005780339
 276275 -2245.067850872    -0.000164343   0.249925264   0.005454407
 277276 -2245.068807031    -0.000956160   1.255604389   0.005104442
 278277 -2245.072256890    -0.003449859   0.627887770   0.003370216
 279278 -2245.074122436    -0.001865546   0.442859827   0.003449145
 280279 -2245.075319471    -0.001197036   0.861669385   0.002924210
 281280 -2245.076650572    -0.001331100   0.226068703   0.002118487
 282281 -2245.077048648    -0.000398076   0.487318368   0.001837157
 283282 -2245.077464975    -0.000416327   0.308136352   0.002699655
 284283 -2245.077810035    -0.000345060   0.130826245   0.000797306
 285284 -2245.077941840    -0.000131805   0.136208827   0.001201883
 286285 -2245.078107167    -0.000165327   0.172359369   0.000791082
 287286 -2245.078218102    -0.000110935   0.039455556   0.000598950
 288287 -2245.078213993     0.000004110   0.057030368   0.000640275
 289288 -2245.078241539    -0.000027547   0.097510898   0.001026310
 290289 -2245.078228700     0.000012839   0.068123904   0.001491373
 291290 -2245.078302005    -0.000073305   0.034560660   0.001124211
 292291 -2245.078318831    -0.000016826   0.267224902   0.001169005
 293292 -2245.078167216     0.000151616   0.185845958   0.001181428
 294293 -2245.078103662     0.000063554   0.207518047   0.001047879
 295294 -2245.078071054     0.000032608   0.224247004   0.000903015
 296295 -2245.078131128    -0.000060074   0.202087494   0.000391715
 297296 -2245.078176777    -0.000045650   0.085014222   0.000343801
 298297 -2245.078232065    -0.000055287   0.037793888   0.000327402
 299298 -2245.078275591    -0.000043526   0.029685478   0.000261088
 300299 -2245.078311406    -0.000035816   0.031625211   0.000267470
 301300 -2245.078364093    -0.000052687   0.043482269   0.000201421
 302301 -2245.078405294    -0.000041201   0.032325591   0.000209330
 303302 -2245.078435514    -0.000030220   0.055543903   0.000189165
 304303 -2245.078476995    -0.000041481   0.030447211   0.000181572
 305304 -2245.078501828    -0.000024833   0.014400959   0.000158371
 306305 -2245.078496662     0.000005166   0.017150864   0.000227343
 307306 -2245.078500848    -0.000004186   0.003709714   0.000225234
 308307 -2245.078498793     0.000002055   0.003150639   0.000234371
 309308 -2245.078497918     0.000000875   0.004912599   0.000241109
 310309 -2245.078498410    -0.000000492   0.006514064   0.000234997
 311310 -2245.078500202    -0.000001792   0.003833909   0.000232819
 312311 -2245.078502232    -0.000002030   0.033976998   0.000235908
 313312 -2245.078517764    -0.000015532   0.020185213   0.000193897
 314313 -2245.078526252    -0.000008488   0.014437168   0.000194316
 315314 -2245.078532540    -0.000006289   0.035714088   0.000210228
 316315 -2245.078518607     0.000013933   0.026989952   0.000222380
 317316 -2245.078509496     0.000009112   0.008491000   0.000274462
 318317 -2245.078506916     0.000002580   0.004188978   0.000195991
 319318 -2245.078508053    -0.000001138   0.007002227   0.000182208
 320319 -2245.078506745     0.000001308   0.007713461   0.000184564
 321320 -2245.078505645     0.000001100   0.008285948   0.000178176
 322321 -2245.078506315    -0.000000670   0.012693617   0.000180098
 323322 -2245.078506034     0.000000281   0.020565267   0.000171599
 324323 -2245.078502969     0.000003065   0.057728327   0.000184353
 325324 -2245.078500741     0.000002228   0.045177503   0.000145204
 326325 -2245.078507248    -0.000006507   0.015317086   0.000113073
 327326 -2245.078514032    -0.000006784   0.019081000   0.000104134
 328327 -2245.078521516    -0.000007484   0.009247509   0.000128205
 329328 -2245.078526998    -0.000005482   0.022341002   0.000103740
 330329 -2245.078532175    -0.000005177   0.012526221   0.000095992
 331330 -2245.078534346    -0.000002171   0.016415915   0.000084618
 332331 -2245.078537466    -0.000003120   0.017359238   0.000095356
 333332 -2245.078541548    -0.000004082   0.009602166   0.000078438
 334333 -2245.078542837    -0.000001289   0.004046590   0.000098610
 335334 -2245.078543671    -0.000000834   0.001863989   0.000099381
 336335 -2245.078543950    -0.000000279   0.006630665   0.000099515
 337336 -2245.078543432     0.000000518   0.002379664   0.000098122
 338337 -2245.078543570    -0.000000138   0.001517638   0.000090990
 339338 -2245.078543478     0.000000092   0.018737539   0.000090100
 340339 -2245.078543701    -0.000000223   0.011107278   0.000086809
 341340 -2245.078543546     0.000000154   0.012274587   0.000090494
 342341 -2245.078545150    -0.000001603   0.009404754   0.000116340
 343342 -2245.078545630    -0.000000480   0.015951761   0.000098639
 344343 -2245.078547335    -0.000001705   0.007852038   0.000131461
 345344 -2245.078547216     0.000000119   0.003050787   0.000142136
 346345 -2245.078547593    -0.000000378   0.005555230   0.000131178
 347346 -2245.078546949     0.000000645   0.008541071   0.000132059
 348347 -2245.078545986     0.000000962   0.019090081   0.000131095
 349348 -2245.078544698     0.000001288   0.029700589   0.000106894
 350349 -2245.078546715    -0.000002017   0.012610818   0.000089418
 351350 -2245.078546053     0.000000662   0.009333408   0.000108679
 352351 -2245.078546583    -0.000000530   0.009100364   0.000133783
 353352 -2245.078548872    -0.000002289   0.027275660   0.000125934
 354353 -2245.078547264     0.000001608   0.006805587   0.000178087
 355354 -2245.078549097    -0.000001833   0.019875368   0.000145665
 356355 -2245.078551460    -0.000002363   0.010835582   0.000142532
 357356 -2245.078552682    -0.000001221   0.040405584   0.000159800
 358357 -2245.078546869     0.000005812   0.049199249   0.000183713
 359358 -2245.078540175     0.000006694   0.032240042   0.000148948
 360359 -2245.078538727     0.000001448   0.012034386   0.000129330
 361360 -2245.078539065    -0.000000338   0.027246028   0.000081616
 362361 -2245.078538093     0.000000972   0.014425229   0.000076592
 363362 -2245.078537604     0.000000489   0.011051583   0.000047398

Все так же, энергия то сходится, то расходится.
def2-SVP(органика)-def2-TZVPD (Cd²⁺)он же TZVP по https://bse.pnl.gov/bse/portal M=1 UHF,DFT

Код: Выделить всё

 ITER EX   TOTAL ENERGY      E CHANGE  DENSITY CHANGE    DIIS ERROR
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0 -2464.546196063 -2464.5461960633238.587603870  43.687923881
   2  1 -2519.137347844   -54.5911517816552.860655860   2.157472422
   3  2 -2286.976680841   232.160667003**************   8.120906795
   4  3 -2302.863385151   -15.886704311 592.442945344   9.300456542
   5  4 -2305.713034657    -2.849649506  77.490857074  10.034335904
   6  5 -2306.175420602    -0.462385945  34.607747487  10.321296352
   7  6 -2306.444863822    -0.269443220  83.669476053  10.457012451
   8  7 -2306.553038262    -0.108174440 358.751947289  10.434523905
   9  8 -2339.642091828   -33.089053566  13.903182153  10.691000173
  10  9 -2344.661637221    -5.019545393 205.715235603  10.716295114
  11 10 -2350.284412728    -5.622775507**************  10.621797111

Скачет по бешеному.
def2-SVP(органика)-def2-SVP (Cd²⁺) M=1 UHF,DFT

Код: Выделить всё

 ITER EX   TOTAL ENERGY      E CHANGE  DENSITY CHANGE    DIIS ERROR
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0 -2014.609584573 -2014.6095845734262.890655630  15.953713115
   2  1 -2173.006131338  -158.396546765 471.273661265   1.348024473
   3  2 -2152.092363552    20.913767785  91.008104267   2.224598604
   4  3 -2225.586970681   -73.494607129 119.084522036   2.165957823
   5  4 -2230.979221221    -5.392250540  57.166920335   1.162472764
   6  5 -2238.112014649    -7.132793428  39.914651620   0.672438416
   7  6 -2241.546662789    -3.434648140  31.610561116   0.401298840
   8  7 -2242.729630834    -1.182968045  21.491861183   0.176470180
   9  8 -2243.031601966    -0.301971132   4.742745689   0.031339933
  10  9 -2243.094667167    -0.063065201   0.612090412   0.016640410
  11 10 -2243.116866157    -0.022198990   1.057567475   0.017216546
  12 11 -2243.134324799    -0.017458642   0.860295507   0.013612708
  13 12 -2243.146110860    -0.011786061   0.949133395   0.007234971
  14 13 -2243.154418411    -0.008307551   0.141008360   0.003611596
  15 14 -2243.157508684    -0.003090273   0.276285176   0.004467466
  16 15 -2243.159508923    -0.002000239   0.194956189   0.005607623
  17 16 -2243.161043158    -0.001534236   0.286680658   0.004424241
  18 17 -2243.162792510    -0.001749352   0.076985245   0.001866345
  19 18 -2243.163855694    -0.001063184   0.180839843   0.001462432
  20 19 -2243.164592835    -0.000737141   0.059681507   0.001762802
  21 20 -2243.165065215    -0.000472381   0.030308869   0.002117947
  22 21 -2243.165029327     0.000035888   0.045708970   0.003400039
  23 22 -2243.165411323    -0.000381996   0.085468574   0.002623661
  24 23 -2243.166135751    -0.000724428   0.093289521   0.001967237
  25 24 -2243.166697342    -0.000561590   0.038650155   0.001533651
  26 25 -2243.167051173    -0.000353831   0.017233667   0.001432945
  27 26 -2243.167403879    -0.000352706   0.057996873   0.001526066
  28 27 -2243.167474723    -0.000070844   0.026346538   0.001815968
  29 28 -2243.167421422     0.000053301   0.078455885   0.002329821
  30 29 -2243.167048219     0.000373203   0.116200190   0.002340260
  31 30 -2243.167139278    -0.000091059   0.055761949   0.001854584
  32 31 -2243.167102289     0.000036989   0.009481403   0.002638033
  33 32 -2243.167001075     0.000101214   0.023000842   0.002912969
  34 33 -2243.166996247     0.000004827   0.235453821   0.002800243
  35 34 -2243.167868209    -0.000871962   0.027666952   0.001517832
  36 35 -2243.168378869    -0.000510660   0.047175228   0.001495147
  37 36 -2243.169189659    -0.000810790   0.051914710   0.001818588
  38 37 -2243.169839684    -0.000650025   0.023709080   0.002039203
  39 38 -2243.170341366    -0.000501682   0.025317041   0.002360721
  40 39 -2243.170613747    -0.000272381   0.017819604   0.002536132
  41 40 -2243.170519772     0.000093974   0.019547545   0.003237441
  42 41 -2243.169852398     0.000667375   0.017533135   0.004121765
  43 42 -2243.169713033     0.000139365   0.083906412   0.004478482
  44 43 -2243.169976557    -0.000263524   0.040326176   0.003679572
  45 44 -2243.170523487    -0.000546930   0.033276830   0.003932811
  46 45 -2243.171493782    -0.000970295   0.058766993   0.003429390
  47 46 -2243.172206841    -0.000713059   0.044327531   0.003081099
  48 47 -2243.172495467    -0.000288626   0.068073363   0.003787213
  49 48 -2243.172992588    -0.000497120   0.122227981   0.002889032
  50 49 -2243.173590544    -0.000597957   0.170432776   0.002830890
  51 50 -2243.174157864    -0.000567320   0.090490824   0.003695563
  52 51 -2243.174545659    -0.000387795   0.099630736   0.003954507
  53 52 -2243.174917056    -0.000371397   0.068787635   0.003283421
  54 53 -2243.175315579    -0.000398523   0.130979801   0.003429773
  55 54 -2243.175882886    -0.000567307   0.333296474   0.002985224
  56 55 -2243.177090542    -0.001207656   0.608565710   0.002452083
  57 56 -2243.178102758    -0.001012216   0.236241108   0.004565171
  58 57 -2243.178681199    -0.000578442   0.078084388   0.002285077
  59 58 -2243.179004066    -0.000322867   0.035635356   0.002137786
  60 59 -2243.179014603    -0.000010537   0.036639753   0.002312027
  61 60 -2243.178920653     0.000093950   0.035980412   0.002262712
  62 61 -2243.178803674     0.000116979   0.044214289   0.002388076
  63 62 -2243.178984383    -0.000180709   0.048113792   0.001696080
  64 63 -2243.179016414    -0.000032031   0.119494549   0.001146330
  65 64 -2243.179004700     0.000011713   0.091833369   0.001466627
  66 65 -2243.179101682    -0.000096982   0.032838963   0.000875706
  67 66 -2243.179028858     0.000072824   0.068988651   0.001862544
  68 67 -2243.179206107    -0.000177249   0.112227609   0.001183703
  69 68 -2243.179326595    -0.000120487   0.015561179   0.000964324
  70 69 -2243.179395562    -0.000068968   0.016026357   0.000877946
  71 70 -2243.179384923     0.000010640   0.020869695   0.000853063
  72 71 -2243.179376547     0.000008376   0.017807340   0.000803993
  73 72 -2243.179445742    -0.000069195   0.015700152   0.000861105
  74 73 -2243.179492580    -0.000046839   0.009247769   0.000990219
  75 74 -2243.179465558     0.000027023   0.028946809   0.000894256
  76 75 -2243.179581777    -0.000116219   0.019346678   0.001112174
  77 76 -2243.179701501    -0.000119724   0.008976671   0.001306648
  78 77 -2243.179679357     0.000022144   0.012016905   0.001505813
  79 78 -2243.179608012     0.000071345   0.031243261   0.001557722
  80 79 -2243.179608921    -0.000000909   0.024013540   0.001424468
  81 80 -2243.179755050    -0.000146129   0.027541565   0.001531383
  82 81 -2243.180033011    -0.000277961   0.019685556   0.001773689
  83 82 -2243.180223561    -0.000190550   0.027726134   0.002008995
  84 83 -2243.180438841    -0.000215280   0.019525122   0.002188613
  85 84 -2243.180419000     0.000019841   0.024342365   0.002205221
  86 85 -2243.180674758    -0.000255758   0.067583996   0.002219251
  87 86 -2243.180634357     0.000040401   0.074593973   0.002688336
  88 87 -2243.180808925    -0.000174568   0.014408173   0.002552739
  89 88 -2243.180907123    -0.000098198   0.024878030   0.002496634
  90 89 -2243.181154509    -0.000247387   0.019517346   0.002826567
  91 90 -2243.181385155    -0.000230645   0.097414635   0.002859373
  92 91 -2243.182325645    -0.000940490   0.058644081   0.003478443
  93 92 -2243.183481755    -0.001156110   0.068946388   0.003900479
  94 93 -2243.181832270     0.001649486   0.080985014   0.003566333
  95 94 -2243.180425084     0.001407185   0.168056710   0.002947018
  96 95 -2243.179528268     0.000896816   0.080853309   0.001682596
  97 96 -2243.179791093    -0.000262826   0.043944086   0.002169666
  98 97 -2243.180251730    -0.000460637   0.024713651   0.002515465
  99 98 -2243.180619264    -0.000367534   0.015178617   0.002879218
 100 99 -2243.180822163    -0.000202899   0.038451789   0.003260192
 101100 -2243.180747835     0.000074328   0.089137930   0.003921616
 102101 -2243.181090131    -0.000342295   0.049132631   0.002475889
 103102 -2243.181807644    -0.000717513   0.020352229   0.002813174
 104103 -2243.181389776     0.000417868   0.043883069   0.002629285
 105104 -2243.180148538     0.001241238   0.034649347   0.002107100
 106105 -2243.179679458     0.000469080   0.032557489   0.001214125
 107106 -2243.179584862     0.000094596   0.031227037   0.001119191
 108107 -2243.179823917    -0.000239054   0.021294925   0.001239908
 109108 -2243.179944060    -0.000120143   0.015856305   0.001529041
 110109 -2243.180046854    -0.000102794   0.024523255   0.001581325
 111110 -2243.180207078    -0.000160224   0.023754988   0.001463192
 112111 -2243.180424986    -0.000217908   0.013457797   0.001616761
 113112 -2243.180397860     0.000027126   0.037338662   0.001707942
 114113 -2243.179952895     0.000444965   0.058271997   0.001369171
 115114 -2243.179916272     0.000036624   0.027376844   0.001289674
 116115 -2243.180012109    -0.000095837   0.021821349   0.001217391
 117116 -2243.179988001     0.000024108   0.056368505   0.001328283
 118117 -2243.179987919     0.000000082   0.030545083   0.002563389
 119118 -2243.179804482     0.000183437   0.027174468   0.001413452
 120119 -2243.180314279    -0.000509797   0.048946311   0.002099810
 121120 -2243.180599845    -0.000285565   0.014095691   0.001681567
 122121 -2243.180726678    -0.000126834   0.024638561   0.001818613
 123122 -2243.180746128    -0.000019450   0.041319918   0.002092503

Тут получше, но все равно не нравится мне это.

Получается что эти базисы не подходят? Я попробовал посчитать кадмий 3-21+G* сходится отлично. За достоверность последнего я не берусь.

Если металлы def2-TZVPD то все остальное базисом того же типа. То есть def2-SVP. Иначе появление проблем от "стыка" вероятнее. Кроме того если вы накручиваете диффузию на металл, возможно неплохо увеличить и остаток до def2-SVPD. Разумеется станет дороже, но мне это кажется корректнее.

Вы говорите про суперпозиционную ошибку?

[Darvin56]

Делаю вывод, за все эти дни так сходимости и не добился. Все считает нормально, но никак не может сойтись то требуемого уровня. Ставил все вот инпут

Код: Выделить всё

 $CONTRL SCFTYP=UHF RUNTYP=optimize MULT=1 icut=6 icharg=0 $END
 $CONTRL DFTTYP=B3LYP $END
 $CONTRL MAXIT=3000 NZVAR=1 $END
 $SCF    COUPLE=.TRUE.
         F(1)=1.0, 0.5
         ALPHA(1)=  0.0,  1.0,  1.0
         BETA(1)=   10.0,  0.0,  0.0
 $END
 $SYSTEM TIMLIM=30000 MWORDS=260 $END
 $STATPT OPTTOL=1.0E-3 NSTEP=2000 METHOD=GDIIS hssend=.t. $END
 $ZMAT dlc=.t. auto=.t. $END
!$BASIS  GBASIS=N31 NGAUSS=6 NPFUNC=2 NDFUNC=2 $END
!$BASIS GBASIS=TZV $END
!$GUESS GUESS=moread $END 
 $GUESS GUESS=HUCKEL $END
!$SCF  SORS=.TRUE. $END DIRSCF=.True.  CICOEF=.True. MVOQ=.TRUE. $END
 $SCF  DAMP=.TRUE. DMPCUT=0.6 $END

Не знаю что делать. Прошу знающих людей опять, либо в подборе базиса, либо оптимизации настроек.

исходный инпут

Код: Выделить всё

 $CONTRL SCFTYP=UHF RUNTYP=optimize MULT=1 icut=6 icharg=0 $END
 $CONTRL DFTTYP=B3LYP $END
 $CONTRL MAXIT=3000 NZVAR=0 $END
 $SYSTEM TIMLIM=30000 MWORDS=260 $END
 $STATPT OPTTOL=1.0E-4 NSTEP=3000  $END
!$BASIS  GBASIS=N31 NGAUSS=6 NPFUNC=2 NDFUNC=2 $END
!$BASIS GBASIS=TZV $END
!$GUESS GUESS=moread $END 
 $GUESS GUESS=HUCKEL $END
 $SCF  DAMP=.TRUE. DMPCUT=0.6 $END

[Гесс]

Возможно я был неправ отделив viewtopic.php?f=71&t=110563.
Там есть интересный набор рекомендаций от alien308. Ну и как обычно ответ почти столь же универсальный как "42" - "а попробуйте мультиреференс".
Заселенность орбиталей хоть каким-то боком известна? (скажем из синглпоинта даже если и не в минимуме, но хоть как то).