Проблема с рестартом сканирования ППЭ в Gaussian

[David Beckham]

Уважаемые товарищи вот такой вопрос.
При сканировании ППЭ в Гауссиане расчет вылетел. Чтобы продолжить скан с того места, где он закончился, я нашел в Help Гауссиана следующее:
"Restart
Restarts a PES scan calculation. A failed Scan calculation may be restarted from its checkpoint file by simply repeating the route section of the original job, adding the Restart option to the Scan keyword. No other input is required."
Соответственно сделал новый инпут следующего содержания:

%mem=30000000
%chk=E:\gaussian_files\transB_scan3.chk
#p Scan=Restart

К сожалению ничего не получилось, а в *.out
Можно прочитать следующее:
......
Title: transB_scan3 B3LYP /6-31G* on B3LYP Geom
Route: #p opt=modredundant b3lyp/6-31g(d) nosymm
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -507 Len= 170002 IPos= 0 Q= 6119648
.....
С файлом transB_scan3.chk я ничего не делал.
Объясните пожалуйста что я не так делаю и как его все-таки запустить?

[o-oxhem]

Уважаемые товарищи вот такой вопрос.
При сканировании ППЭ в Гауссиане расчет вылетел. Чтобы продолжить скан с того места, где он закончился, я нашел в Help Гауссиана следующее:
"Restart
Restarts a PES scan calculation. A failed Scan calculation may be restarted from its checkpoint file by simply repeating the route section of the original job, adding the Restart option to the Scan keyword. No other input is required."
Соответственно сделал новый инпут следующего содержания:

%mem=30000000
%chk=E:\gaussian_files\transB_scan3.chk
#p Scan=Restart

К сожалению ничего не получилось, а в *.out
Можно прочитать следующее:
......
Title: transB_scan3 B3LYP /6-31G* on B3LYP Geom
Route: #p opt=modredundant b3lyp/6-31g(d) nosymm
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)= -507 Len= 170002 IPos= 0 Q= 6119648
.....
С файлом transB_scan3.chk я ничего не делал.
Объясните пожалуйста что я не так делаю и как его все-таки запустить?

Итак, Ваш роут
#p opt=modredundant b3lyp/6-31g(d) nosymm
restart:
# opt=(RESTART,modredundant) b3lyp/6-31g(d) nosymm

[David Beckham]

Спасибо, уже лучше.
Теперь другая проблема, скан закончился неудачей из-за превышения количества иттераций на одном из шагов (по крайней мере я так понял). Как увеличить количество шагов. Куда добавить MaxCycle=....
Пробовал вот так: #p opt=(restart, modredundant, MaxCycle=999) b3lyp/6-31g(d) no symm
И не получилось.

[o-oxhem]

Спасибо, уже лучше.
Теперь другая проблема, скан закончился неудачей из-за превышения количества иттераций на одном из шагов (по крайней мере я так понял). Как увеличить количество шагов. Куда добавить MaxCycle=....
Пробовал вот так: #p opt=(restart, modredundant, MaxCycle=999) b3lyp/6-31g(d) no symm
И не получилось.

Прошу выкладывать последние 200 строк файла - Ваша секретность не пострадает, а помогать будет проще.

Предел для MaxCycle опредеяется gau$$ian'ом - его Вам не преодолеть.
Во-вторых, если есть keyword restart, то нельзя ничего прибавлять.
Придется начинать ещё раз, учитывая ошибки.

Можно попробовать:
# opt=(calcfc, modredundant, MaxCycle=посмотрите какое максимальное значение указывает в out файле гауссиан) b3lyp/6-31g(d) nosymm

[David Beckham]

Файл не прикрепляется, поэтому скопировал последние строки
...................
Item Value Threshold Converged?
Maximum Force 0.000200 0.000450 YES
RMS Force 0.000047 0.000300 YES
Maximum Displacement 0.000663 0.001800 YES
RMS Displacement 0.000184 0.001200 YES
Predicted change in Energy=-1.281002D-07
Optimization completed.
-- Stationary point found.
----------------------------
! Optimized Parameters !
! (Angstroms and Degrees) !
-------------------------- --------------------------
! Name Definition Value Derivative Info. !
--------------------------------------------------------------------------------
! R1 R(1,2) 1.3416 -DE/DX = 0.0 !
! R2 R(1,5) 1.3496 -DE/DX = -0.0001 !
! R3 R(1,14) 1.0123 -DE/DX = 0.0 !
! R4 R(2,3) 1.305 -DE/DX = 0.0001 !
! R5 R(3,4) 1.3523 -DE/DX = -0.0001 !
! R6 R(4,5) 1.3245 -DE/DX = 0.0001 !
! R7 R(5,6) 1.4071 -DE/DX = 0.0 !
! R8 R(6,7) 1.2498 -DE/DX = 0.0 !
! R9 R(7,12) 1.3095 -DE/DX = 0.0 !
! R10 R(8,9) 1.3273 -DE/DX = 0.0 !
! R11 R(8,12) 1.3491 -DE/DX = -0.0001 !
! R12 R(9,10) 1.3331 -DE/DX = 0.0 !
! R13 R(10,11) 1.3134 -DE/DX = -0.0001 !
! R14 R(11,12) 1.3772 -DE/DX = -0.0002 !
! R15 R(11,13) 1.0148 -DE/DX = 0.0 !
! A1 A(2,1,5) 108.8981 -DE/DX = -0.0001 !
! A2 A(2,1,14) 121.5885 -DE/DX = 0.0001 !
! A3 A(5,1,14) 129.5134 -DE/DX = 0.0 !
! A4 A(1,2,3) 105.7875 -DE/DX = 0.0001 !
! A5 A(2,3,4) 111.4713 -DE/DX = -0.0001 !
! A6 A(3,4,5) 105.7164 -DE/DX = 0.0 !
! A7 A(1,5,4) 108.1267 -DE/DX = 0.0001 !
! A8 A(1,5,6) 119.833 -DE/DX = -0.0001 !
! A9 A(4,5,6) 132.0403 -DE/DX = 0.0 !
! A10 A(5,6,7) 111.9782 -DE/DX = 0.0 !
! A11 A(6,7,12) 160.0 -DE/DX = -0.0587 !
! A12 A(9,8,12) 107.0843 -DE/DX = 0.0 !
! A13 A(8,9,10) 111.188 -DE/DX = 0.0 !
! A14 A(9,10,11) 105.9109 -DE/DX = 0.0 !
! A15 A(10,11,12) 110.1191 -DE/DX = 0.0001 !
! A16 A(10,11,13) 122.4985 -DE/DX = 0.0 !
! A17 A(12,11,13) 127.3824 -DE/DX = 0.0 !
! A18 A(7,12,8) 133.7112 -DE/DX = 0.0 !
! A19 A(7,12,11) 120.5912 -DE/DX = 0.0 !
! A20 A(8,12,11) 105.6977 -DE/DX = 0.0 !
! D1 D(5,1,2,3) -0.0054 -DE/DX = 0.0 !
! D2 D(14,1,2,3) -180.0029 -DE/DX = 0.0 !
! D3 D(2,1,5,4) 0.0024 -DE/DX = 0.0 !
! D4 D(2,1,5,6) 180.0025 -DE/DX = 0.0 !
! D5 D(14,1,5,4) -180.0003 -DE/DX = 0.0 !
! D6 D(14,1,5,6) -0.0002 -DE/DX = 0.0 !
! D7 D(1,2,3,4) 0.0065 -DE/DX = 0.0 !
! D8 D(2,3,4,5) -0.0051 -DE/DX = 0.0 !
! D9 D(3,4,5,1) 0.0015 -DE/DX = 0.0 !
! D10 D(3,4,5,6) 180.0013 -DE/DX = 0.0 !
! D11 D(1,5,6,7) 180.0007 -DE/DX = 0.0 !
! D12 D(4,5,6,7) 0.0009 -DE/DX = 0.0 !
! D13 D(6,7,12,8) -0.0003 -DE/DX = 0.0 !
! D14 D(6,7,12,11) 180.0001 -DE/DX = 0.0 !
! D15 D(12,8,9,10) 0.0016 -DE/DX = 0.0 !
! D16 D(9,8,12,7) -180.0006 -DE/DX = 0.0 !
! D17 D(9,8,12,11) -0.001 -DE/DX = 0.0 !
! D18 D(8,9,10,11) -0.0015 -DE/DX = 0.0 !
! D19 D(9,10,11,12) 0.0008 -DE/DX = 0.0 !
! D20 D(9,10,11,13) 180.0006 -DE/DX = 0.0 !
! D21 D(10,11,12,7) -180.0003 -DE/DX = 0.0 !
! D22 D(10,11,12,8) 0.0001 -DE/DX = 0.0 !
! D23 D(13,11,12,7) 0.0 -DE/DX = 0.0 !
! D24 D(13,11,12,8) 180.0004 -DE/DX = 0.0 !
--------------------------------------------------------------------------------
Iteration 1 RMS(Cart)= 0.13970365 RMS(Int)= 0.02167806
Iteration 2 RMS(Cart)= 0.00017393 RMS(Int)= 0.02164926
Iteration 3 RMS(Cart)= 0.00017326 RMS(Int)= 0.02162058
Iteration 4 RMS(Cart)= 0.00017260 RMS(Int)= 0.02159200
Iteration 5 RMS(Cart)= 0.00017194 RMS(Int)= 0.02156354
Iteration 6 RMS(Cart)= 0.00017128 RMS(Int)= 0.02153519
Iteration 7 RMS(Cart)= 0.00017063 RMS(Int)= 0.02150695
Iteration 8 RMS(Cart)= 0.00016999 RMS(Int)= 0.02147882
Iteration 9 RMS(Cart)= 0.00016934 RMS(Int)= 0.02145079
Iteration 10 RMS(Cart)= 0.00016870 RMS(Int)= 0.02142288
Iteration 11 RMS(Cart)= 0.00016807 RMS(Int)= 0.02139507
Iteration 12 RMS(Cart)= 0.00016744 RMS(Int)= 0.02136736
Iteration 13 RMS(Cart)= 0.00016681 RMS(Int)= 0.02133977
Iteration 14 RMS(Cart)= 0.00016619 RMS(Int)= 0.02131227
Iteration 15 RMS(Cart)= 0.00016557 RMS(Int)= 0.02128488
Iteration 16 RMS(Cart)= 0.00016495 RMS(Int)= 0.02125760
Iteration 17 RMS(Cart)= 0.00016434 RMS(Int)= 0.02123042
Iteration 18 RMS(Cart)= 0.00016373 RMS(Int)= 0.02120333
Iteration 19 RMS(Cart)= 0.00016313 RMS(Int)= 0.02117636
Iteration 20 RMS(Cart)= 0.00016253 RMS(Int)= 0.02114948
Iteration 21 RMS(Cart)= 0.00016193 RMS(Int)= 0.02112270
Iteration 22 RMS(Cart)= 0.00016134 RMS(Int)= 0.02109602
Iteration 23 RMS(Cart)= 0.00016075 RMS(Int)= 0.02106944
Iteration 24 RMS(Cart)= 0.00016016 RMS(Int)= 0.02104296
Iteration 25 RMS(Cart)= 0.00015958 RMS(Int)= 0.02101658
Iteration 26 RMS(Cart)= 0.00015900 RMS(Int)= 0.02099029
Iteration 27 RMS(Cart)= 0.00015843 RMS(Int)= 0.02096411
Iteration 28 RMS(Cart)= 0.00015785 RMS(Int)= 0.02093801
Iteration 29 RMS(Cart)= 0.00015729 RMS(Int)= 0.02091202
Iteration 30 RMS(Cart)= 0.00015672 RMS(Int)= 0.02088611
Iteration 31 RMS(Cart)= 0.00015616 RMS(Int)= 0.02086031
Iteration 32 RMS(Cart)= 0.00015560 RMS(Int)= 0.02083459
Iteration 33 RMS(Cart)= 0.00015505 RMS(Int)= 0.02080897
Iteration 34 RMS(Cart)= 0.00015449 RMS(Int)= 0.02078345
Iteration 35 RMS(Cart)= 0.00015394 RMS(Int)= 0.02075801
Iteration 36 RMS(Cart)= 0.00015340 RMS(Int)= 0.02073267
Iteration 37 RMS(Cart)= 0.00015286 RMS(Int)= 0.02070741
Iteration 38 RMS(Cart)= 0.00015232 RMS(Int)= 0.02068225
Iteration 39 RMS(Cart)= 0.00015178 RMS(Int)= 0.02065718
Iteration 40 RMS(Cart)= 0.00015125 RMS(Int)= 0.02063219
Iteration 41 RMS(Cart)= 0.00015072 RMS(Int)= 0.02060730
Iteration 42 RMS(Cart)= 0.00015019 RMS(Int)= 0.02058249
Iteration 43 RMS(Cart)= 0.00014967 RMS(Int)= 0.02055778
Iteration 44 RMS(Cart)= 0.00014915 RMS(Int)= 0.02053315
Iteration 45 RMS(Cart)= 0.00014863 RMS(Int)= 0.02050860
Iteration 46 RMS(Cart)= 0.00014812 RMS(Int)= 0.02048414
Iteration 47 RMS(Cart)= 0.00014761 RMS(Int)= 0.02045977
Iteration 48 RMS(Cart)= 0.00014710 RMS(Int)= 0.02043549
Iteration 49 RMS(Cart)= 0.00014659 RMS(Int)= 0.02041129
Iteration 50 RMS(Cart)= 0.00014609 RMS(Int)= 0.02038717
Iteration 51 RMS(Cart)= 0.00014559 RMS(Int)= 0.02036313
Iteration 52 RMS(Cart)= 0.00014509 RMS(Int)= 0.02033918
Iteration 53 RMS(Cart)= 0.00014460 RMS(Int)= 0.02031532
Iteration 54 RMS(Cart)= 0.00014410 RMS(Int)= 0.02029153
Iteration 55 RMS(Cart)= 0.00014362 RMS(Int)= 0.02026783
Iteration 56 RMS(Cart)= 0.00014313 RMS(Int)= 0.02024421
Iteration 57 RMS(Cart)= 0.00014265 RMS(Int)= 0.02022067
Iteration 58 RMS(Cart)= 0.00014217 RMS(Int)= 0.02019721
Iteration 59 RMS(Cart)= 0.00014169 RMS(Int)= 0.02017383
Iteration 60 RMS(Cart)= 0.00014121 RMS(Int)= 0.02015053
Iteration 61 RMS(Cart)= 0.00014074 RMS(Int)= 0.02012731
Iteration 62 RMS(Cart)= 0.00014027 RMS(Int)= 0.02010417
Iteration 63 RMS(Cart)= 0.00013980 RMS(Int)= 0.02008111
Iteration 64 RMS(Cart)= 0.00013934 RMS(Int)= 0.02005812
Iteration 65 RMS(Cart)= 0.00013888 RMS(Int)= 0.02003521
Iteration 66 RMS(Cart)= 0.00013842 RMS(Int)= 0.02001238
Iteration 67 RMS(Cart)= 0.00013796 RMS(Int)= 0.01998963
Iteration 68 RMS(Cart)= 0.00013750 RMS(Int)= 0.01996695
Iteration 69 RMS(Cart)= 0.00013705 RMS(Int)= 0.01994435
Iteration 70 RMS(Cart)= 0.00013660 RMS(Int)= 0.01992182
Iteration 71 RMS(Cart)= 0.00013615 RMS(Int)= 0.01989937
Iteration 72 RMS(Cart)= 0.00013571 RMS(Int)= 0.01987699
Iteration 73 RMS(Cart)= 0.00013527 RMS(Int)= 0.01985469
Iteration 74 RMS(Cart)= 0.00013483 RMS(Int)= 0.01983246
Iteration 75 RMS(Cart)= 0.00013439 RMS(Int)= 0.01981030
Iteration 76 RMS(Cart)= 0.00013395 RMS(Int)= 0.01978822
Iteration 77 RMS(Cart)= 0.00013352 RMS(Int)= 0.01976621
Iteration 78 RMS(Cart)= 0.00013309 RMS(Int)= 0.01974427
Iteration 79 RMS(Cart)= 0.00013266 RMS(Int)= 0.01972240
Iteration 80 RMS(Cart)= 0.00013224 RMS(Int)= 0.01970060
Iteration 81 RMS(Cart)= 0.00013181 RMS(Int)= 0.01967888
Iteration 82 RMS(Cart)= 0.00013139 RMS(Int)= 0.01965722
Iteration 83 RMS(Cart)= 0.00013097 RMS(Int)= 0.01963564
Iteration 84 RMS(Cart)= 0.00013055 RMS(Int)= 0.01961413
Iteration 85 RMS(Cart)= 0.00013014 RMS(Int)= 0.01959268
Iteration 86 RMS(Cart)= 0.00012973 RMS(Int)= 0.01957130
Iteration 87 RMS(Cart)= 0.00012932 RMS(Int)= 0.01955000
Iteration 88 RMS(Cart)= 0.00012891 RMS(Int)= 0.01952876
Iteration 89 RMS(Cart)= 0.00012850 RMS(Int)= 0.01950759
Iteration 90 RMS(Cart)= 0.00012810 RMS(Int)= 0.01948648
Iteration 91 RMS(Cart)= 0.00012770 RMS(Int)= 0.01946545
Iteration 92 RMS(Cart)= 0.00012730 RMS(Int)= 0.01944448
Iteration 93 RMS(Cart)= 0.00012690 RMS(Int)= 0.01942358
Iteration 94 RMS(Cart)= 0.00012650 RMS(Int)= 0.01940274
Iteration 95 RMS(Cart)= 0.00012611 RMS(Int)= 0.01938197
Iteration 96 RMS(Cart)= 0.00012572 RMS(Int)= 0.01936126
Iteration 97 RMS(Cart)= 0.00012533 RMS(Int)= 0.01934062
Iteration 98 RMS(Cart)= 0.00012494 RMS(Int)= 0.01932005
Iteration 99 RMS(Cart)= 0.00012455 RMS(Int)= 0.01929954
Iteration100 RMS(Cart)= 0.00012417 RMS(Int)= 0.01927909
New curvilinear step not converged.
RedCar failed.
Error termination via Lnk1e in e:\G03W\l103.exe at Wed Dec 09 16:52:06 2009.
Job cpu time: 0 days 0 hours 55 minutes 44.0 seconds.
File lengths (MBytes): RWF= 30 Int= 0 D2E= 0 Chk= 6 Scr= 1

[o-oxhem]

Геометрия предыдущего шага сошлась. НО gau$$ian не может сдалать следующий шаг. В добавок к calcfc можно написать gdiis.
Ещё раз посмотрите, что и куда вращается/движется. К сообщению об ошибке прошу выкладывать route (всё для наглядности, а Вам для более конкретных ответов).

scf(xqc,maxcycle=512,tight)

[David Beckham]

Прикрепляю *.out файл моего расчета.

[o-oxhem]

Прикрепляю *.out файл моего расчета.

аттачмент работает для архивов и .txt. Для .out не работает. Либо архивируете, либо меняете расширение.

[David Beckham]

В архиве инпут и аут.

[o-oxhem]

В архиве инпут и аут.

Вы указали релаксированный скан по диэд. углу D 5 6 7 12 ,
но D * 6 7 * R
При сканировании по валентному углу - достижение величины 180 и более вызывает сообщение об ошибке при сканировании.

[David Beckham]

С удалением координат (R) я уже разобрался, с углом тоже, более 180 не получится - это я знал, не знал что 180 также не прокатит.
Спасибо за помощь!

[VTur]

При сканировании по валентному углу - достижение величины 180 и более вызывает сообщение об ошибке при сканировании.

Да вроде катит.