今年発売された Gaussian 16 ! 新機能も気になるところですが、既存の Gaussian 09 とどれほどの差があるのかという点も気になります。
今回は、構造最適化、遷移状態構造の最適化、振動計算について2つを比較してみました。
今回行った比較から log ファイルの中身や座標、エネルギー、振動など全てにおいて差が見られましたので、以下簡単にまとめました。
今回の記事では、たった 12 個の計算しか行っていませんが、今後暇を見てベンチマークを作成していきたいと考えています。
最近大変忙しく、ソースコードを読んでいる暇がありませんので、、、なぜ差があるのかについては調べていません。。。コンパイルの過程にも問題がある可能性かもしれません、、、。
【2017.04.22 追記】
g16 と g09 ではデフォルトの grid が違うため、正しく評価できてないとの指摘がありました。そのせいで、計算値にズレが生じているとのことです。
そこで、grid 値を揃えて再度検証してみました。
また、より詳細なベンチマークが HPC ソリューションさんから出ているようです。
http://www.hpc.co.jp/benchmark20170307.html
テスト環境
CPU: Intel CPU Core i7-7700 3.6GHz 8Mキャッシュ 4コア/8スレッド
m/b: ASUS PRIME H270-PRO
メモリ(64 GB): Crucial DDR4 2400 MHz 16GB 4枚
SSD: M.2 SSD_250 GB Samsung SSD 960 EVO M.2 Type2280 NVMe1.2
input ファイル
通常のベンチマーク測定で扱う化合物は、面白みがありませんので、実践的な系として実際に論文になっているものを input として用いることにしました。
“(2 + 1)-Cycloaddition Reactions Give Further Evidence of the Nitrenium-like Character of 1-Aza-2-azoniaallene Salts”
J. Org. Chem., 2017, 82 (7), pp 4001–4005. DOI: 10.1021/acs.joc.7b00407
Nezar Al-Bataineh, Kendall N. Houk, Matthias Brewer and Xin Hong
今回は、ケンダール・ハウク Kendall N. Houk らが報告した (2 + 1)-Cycloaddition Reactions を取り上げました。SI の 14, TS15, 16, TS17 の座標を利用しました。
本記事では実際に座標を載せることはしませんので、ご自身で論文と SI をダウンロードしてください。
ベンチマーク測定は 1.構造最適化、2.遷移状態構造の最適化、3.振動計算化合物の3種類について行いました。IRC 計算は今回は行いませんでした。論文では、M06-2X/6-311+G(d,p)//B3LYP/6-31G(d)で計算されていますが、今回は B3LYP/6-31+G(d,p) で計算しました。
まず、14, 16 の構造を用いて構造最適化を行いました。キーワードはシンプルにするため、細かい設定は一切行いませんでした。
1 2 3 4 5 6 7 |
%mem=10000MB %nprocshared=8 %nosave # b3lyp/6-31+g(d,p) opt geom=connectivity |
次に TS15, TS17 の構造を用いて遷移状態構造最適化を行いました。キーワードの記載は省略します。
その後、構造最適化された遷移状態の座標を用いて振動計算を行いました。いつもは、chk ファイルから情報を読み込んで計算時間を短縮させていますが、今回はベンチマーク測定なので、そのようなことはしませんでした。そのため、今回行った全ての計算でchk ファイルを保存しませんでした。
1 2 3 4 5 6 7 |
%mem=10000MB %nprocshared=8 %nosave #p b3lyp/6-31+g(d,p) freq geom=connectivity |
計算時間について(grid を揃えた比較)
g16 では、デフォルトで ultrafine の精度ですが、g09 では デフォルトで fine です。そのため、両方に int=grid=ultrafine というキーワードを追加して比較しました。
計算時間の比較結果を以下にまとめました。縦軸が時間(min.)を表しています。
計算速度は gaussian16 の方が約 1.2 倍早かったです。また、grid を揃えたところ、計算値が一致しました。
管理人の持っている g16 は AVX2 版ですが、 g09 の方は覚えていません。。。
計算時間について(grid を揃えていないデータ)
g09 と g16 の計算時間の比較を以下にまとめました。縦軸が時間(min.)を表しています。grid などキーワード指定しない状態で計算させた場合、両者にあまり大きな差はありませんでした。言い換えれば、同じ計算時間の場合 gaussian 16 の方が精度が高いということでしょうか。
g16 は log ファイルの最後に cpu 時間と実際にかかった時間の両方が表示されます。今回は 8 スレッドで計算しました。実際に計算にかかった時間は cpu 時間 の 7.8 〜 8.0 倍でした。並列化効率もまずまずです。
計算誤差(grid を揃えていないデータ)
grid を揃えない場合、二つの log ファイルを比べると計算値の小数点第4位以下に差があります。a.u. で 0.0001 差があると、kcal/mol では 0.06 kcal/mol 程度差があることになります。
ここで実際に TS15 の構造を用いた遷移状態構造最適化の log ファイルについて比べてみましょう。左側が g16 右側が g09 で計算したものです。まずは、最初の Fock 行列の計算です。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 |
Cycle 1 Pass 1 IDiag 1: Cycle 1 Pass 1 IDiag 1: FoFJK: IHMeth= 1 ICntrl= 0 DoSepK=F KAlg= 0 I1Cent= | FoFJK: IHMeth= 1 ICntrl= 0 DoSepK=F KAlg= 0 I1Cent= IRaf= 880000000 NMat= 1 IRICut= 1 DoRegI=T DoRaf | IRaf= 850000000 NMat= 1 IRICut= 1 DoRegI=T DoRafI=F FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1 FoFCou: FMM=F IPFlag= 0 FMFlag= 100000 FMFlg1 NFxFlg= 0 DoJE=F BraDBF=F KetDBF=F FulRan= NFxFlg= 0 DoJE=F BraDBF=F KetDBF=F FulRan= wScrn= 0.000000 ICntrl= 0 IOpCl= 0 I1Cent= | wScrn= 0.000000 ICntrl= 0 IOpCl= 0 I1Cent= NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 NMat0= 1 NMatS0= 1 NMatT0= 0 NMatD0= 1 Symmetry not used in FoFCou. Symmetry not used in FoFCou. E= -613.283415401633 | E= -613.283384131735 DIIS: error= 1.80D-02 at cycle 1 NSaved= 1. DIIS: error= 1.80D-02 at cycle 1 NSaved= 1. NSaved= 1 IEnMin= 1 EnMin= -613.283415401633 IErMin= 1 E | NSaved= 1 IEnMin= 1 EnMin= -613.283384131735 IErMin= 1 E ErrMax= 1.80D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.18D-01 B ErrMax= 1.80D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.18D-01 B IDIUse=3 WtCom= 8.20D-01 WtEn= 1.80D-01 IDIUse=3 WtCom= 8.20D-01 WtEn= 1.80D-01 Coeff-Com: 0.100D+01 Coeff-Com: 0.100D+01 Coeff-En: 0.100D+01 Coeff-En: 0.100D+01 Coeff: 0.100D+01 Coeff: 0.100D+01 Gap= 0.127 Goal= None Shift= 0.000 Gap= 0.127 Goal= None Shift= 0.000 GapD= 0.127 DampG=1.000 DampE=0.500 DampFc=0.5000 IDamp=- GapD= 0.127 DampG=1.000 DampE=0.500 DampFc=0.5000 IDamp=- Damping current iteration by 5.00D-01 Damping current iteration by 5.00D-01 RMSDP=7.16D-02 MaxDP=7.84D+00 OVMax= 5.01D-01 | RMSDP=7.16D-02 MaxDP=7.83D+00 OVMax= 5.01D-01 Cycle 2 Pass 1 IDiag 1: Cycle 2 Pass 1 IDiag 1: RMSU= 2.47D-02 CP: 1.60D+00 RMSU= 2.47D-02 CP: 1.60D+00 E= -613.312707141472 Delta-E= -0.029291739839 Rise | E= -613.312722347193 Delta-E= -0.029338215457 Rise DIIS: error= 1.17D-02 at cycle 2 NSaved= 2. DIIS: error= 1.17D-02 at cycle 2 NSaved= 2. NSaved= 2 IEnMin= 2 EnMin= -613.312707141472 IErMin= 2 E | NSaved= 2 IEnMin= 2 EnMin= -613.312722347193 IErMin= 2 E ErrMax= 1.17D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 9.33D-02 B ErrMax= 1.17D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 9.33D-02 B IDIUse=3 WtCom= 8.83D-01 WtEn= 1.17D-01 IDIUse=3 WtCom= 8.83D-01 WtEn= 1.17D-01 Coeff-Com: 0.281D+00 0.719D+00 Coeff-Com: 0.281D+00 0.719D+00 Coeff-En: 0.461D+00 0.539D+00 Coeff-En: 0.461D+00 0.539D+00 Coeff: 0.302D+00 0.698D+00 Coeff: 0.302D+00 0.698D+00 Gap= 0.081 Goal= None Shift= 0.000 Gap= 0.081 Goal= None Shift= 0.000 RMSDP=3.32D-02 MaxDP=3.48D+00 DE=-2.93D-02 OVMax= 3.77D-01 | RMSDP=3.33D-02 MaxDP=3.48D+00 DE=-2.93D-02 OVMax= 3.77D-01 Cycle 3 Pass 1 IDiag 1: Cycle 3 Pass 1 IDiag 1: RMSU= 1.20D-02 CP: 1.04D+00 2.17D-01 | RMSU= 1.21D-02 CP: 1.04D+00 2.15D-01 E= -613.488153943089 Delta-E= -0.175446801617 Rise | E= -613.488173199106 Delta-E= -0.175450851913 Rise DIIS: error= 1.25D-02 at cycle 3 NSaved= 3. DIIS: error= 1.25D-02 at cycle 3 NSaved= 3. NSaved= 3 IEnMin= 3 EnMin= -613.488153943089 IErMin= 2 E | NSaved= 3 IEnMin= 3 EnMin= -613.488173199106 IErMin= 2 E ErrMax= 1.25D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.16D-01 B ErrMax= 1.25D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.16D-01 B IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00 IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00 Coeff-En: 0.000D+00 0.218D+00 0.782D+00 Coeff-En: 0.000D+00 0.218D+00 0.782D+00 Coeff: 0.000D+00 0.218D+00 0.782D+00 Coeff: 0.000D+00 0.218D+00 0.782D+00 Gap= 0.112 Goal= None Shift= 0.000 Gap= 0.112 Goal= None Shift= 0.000 RMSDP=1.54D-02 MaxDP=9.49D-01 DE=-1.75D-01 OVMax= 4.55D-01 | RMSDP=1.55D-02 MaxDP=9.53D-01 DE=-1.75D-01 OVMax= 4.54D-01 Cycle 4 Pass 1 IDiag 1: Cycle 4 Pass 1 IDiag 1: RMSU= 8.61D-03 CP: 1.13D+00 3.24D-01 1.42D-02 | RMSU= 8.61D-03 CP: 1.13D+00 3.25D-01 1.33D-02 E= -613.204012690037 Delta-E= 0.284141253052 Rise | E= -613.203964781579 Delta-E= 0.284208417527 Rise DIIS: error= 2.44D-02 at cycle 4 NSaved= 4. DIIS: error= 2.44D-02 at cycle 4 NSaved= 4. NSaved= 4 IEnMin= 3 EnMin= -613.488153943089 IErMin= 2 E | NSaved= 4 IEnMin= 3 EnMin= -613.488173199106 IErMin= 2 E ErrMax= 2.44D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.64D-01 B ErrMax= 2.44D-02 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.64D-01 B IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00 IDIUse=2 WtCom= 0.00D+00 WtEn= 1.00D+00 Coeff-En: 0.000D+00 0.000D+00 0.690D+00 0.310D+00 Coeff-En: 0.000D+00 0.000D+00 0.690D+00 0.310D+00 Coeff: 0.000D+00 0.000D+00 0.690D+00 0.310D+00 Coeff: 0.000D+00 0.000D+00 0.690D+00 0.310D+00 Gap= 0.082 Goal= None Shift= 0.000 Gap= 0.082 Goal= None Shift= 0.000 RMSDP=8.22D-03 MaxDP=6.95D-01 DE= 2.84D-01 OVMax= 2.80D-01 | RMSDP=8.24D-03 MaxDP=6.97D-01 DE= 2.84D-01 OVMax= 2.80D-01 Cycle 5 Pass 1 IDiag 1: Cycle 5 Pass 1 IDiag 1: RMSU= 2.36D-03 CP: 1.04D+00 2.34D-01 4.33D-01 5.57D- | RMSU= 2.37D-03 CP: 1.04D+00 2.34D-01 4.31D-01 5.57D- E= -613.604895118188 Delta-E= -0.400882428152 Rise | E= -613.604926592184 Delta-E= -0.400961810605 Rise DIIS: error= 5.17D-03 at cycle 5 NSaved= 5. DIIS: error= 5.17D-03 at cycle 5 NSaved= 5. NSaved= 5 IEnMin= 5 EnMin= -613.604895118188 IErMin= 5 E | NSaved= 5 IEnMin= 5 EnMin= -613.604926592184 IErMin= 5 E ErrMax= 5.17D-03 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.47D-02 B ErrMax= 5.17D-03 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.47D-02 B IDIUse=3 WtCom= 9.48D-01 WtEn= 5.17D-02 IDIUse=3 WtCom= 9.48D-01 WtEn= 5.17D-02 Coeff-Com: -0.301D-02 0.228D-01 0.248D+00-0.184D-01 0.751D+0 | Coeff-Com: -0.299D-02 0.228D-01 0.248D+00-0.184D-01 0.751D+0 Coeff-En: 0.000D+00 0.000D+00 0.158D+00 0.000D+00 0.842D+0 Coeff-En: 0.000D+00 0.000D+00 0.158D+00 0.000D+00 0.842D+0 Coeff: -0.285D-02 0.216D-01 0.243D+00-0.175D-01 0.756D+0 | Coeff: -0.284D-02 0.216D-01 0.243D+00-0.174D-01 0.756D+0 Gap= 0.104 Goal= None Shift= 0.000 Gap= 0.104 Goal= None Shift= 0.000 RMSDP=3.95D-03 MaxDP=4.05D-01 DE=-4.01D-01 OVMax= 6.90D-02 | RMSDP=3.96D-03 MaxDP=4.07D-01 DE=-4.01D-01 OVMax= 6.90D-02 Cycle 6 Pass 1 IDiag 1: Cycle 6 Pass 1 IDiag 1: RMSU= 1.58D-03 CP: 9.83D-01 1.44D-01 5.47D-01 4.78D- | RMSU= 1.59D-03 CP: 9.83D-01 1.43D-01 5.46D-01 4.78D- E= -613.619086877294 Delta-E= -0.014191759106 Rise | E= -613.619123980574 Delta-E= -0.014197388390 Rise DIIS: error= 1.49D-03 at cycle 6 NSaved= 6. DIIS: error= 1.49D-03 at cycle 6 NSaved= 6. NSaved= 6 IEnMin= 6 EnMin= -613.619086877294 IErMin= 6 E | NSaved= 6 IEnMin= 6 EnMin= -613.619123980574 IErMin= 6 E ErrMax= 1.49D-03 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.80D-03 B ErrMax= 1.49D-03 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.80D-03 B IDIUse=3 WtCom= 9.85D-01 WtEn= 1.49D-02 IDIUse=3 WtCom= 9.85D-01 WtEn= 1.49D-02 Coeff-Com: -0.893D-02 0.201D-01 0.105D+00-0.527D-01 0.511D+0 Coeff-Com: -0.893D-02 0.201D-01 0.105D+00-0.527D-01 0.511D+0 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.124D-0 | Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.125D-0 Coeff: -0.880D-02 0.198D-01 0.103D+00-0.519D-01 0.504D+0 | Coeff: -0.879D-02 0.198D-01 0.103D+00-0.519D-01 0.504D+0 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=1.84D-03 MaxDP=2.01D-01 DE=-1.42D-02 OVMax= 3.21D-02 | RMSDP=1.85D-03 MaxDP=2.01D-01 DE=-1.42D-02 OVMax= 3.22D-02 Cycle 7 Pass 1 IDiag 1: Cycle 7 Pass 1 IDiag 1: RMSU= 2.93D-04 CP: 1.01D+00 1.90D-01 5.64D-01 4.53D- | RMSU= 2.94D-04 CP: 1.01D+00 1.89D-01 5.63D-01 4.53D- CP: 5.88D-01 | CP: 5.89D-01 E= -613.621301122314 Delta-E= -0.002214245020 Rise | E= -613.621340180721 Delta-E= -0.002216200147 Rise DIIS: error= 3.21D-04 at cycle 7 NSaved= 7. | DIIS: error= 3.20D-04 at cycle 7 NSaved= 7. NSaved= 7 IEnMin= 7 EnMin= -613.621301122314 IErMin= 7 E | NSaved= 7 IEnMin= 7 EnMin= -613.621340180721 IErMin= 7 E ErrMax= 3.21D-04 0.00D+00 EMaxC= 1.00D-01 BMatC= 7.52D-05 B | ErrMax= 3.20D-04 0.00D+00 EMaxC= 1.00D-01 BMatC= 7.52D-05 B IDIUse=3 WtCom= 9.97D-01 WtEn= 3.21D-03 | IDIUse=3 WtCom= 9.97D-01 WtEn= 3.20D-03 Coeff-Com: -0.438D-02 0.984D-02 0.462D-01-0.277D-01 0.234D+0 | Coeff-Com: -0.438D-02 0.983D-02 0.462D-01-0.277D-01 0.234D+0 Coeff-Com: 0.517D+00 Coeff-Com: 0.517D+00 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+0 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+0 Coeff-En: 0.100D+01 Coeff-En: 0.100D+01 Coeff: -0.437D-02 0.981D-02 0.460D-01-0.276D-01 0.234D+0 | Coeff: -0.437D-02 0.980D-02 0.460D-01-0.276D-01 0.234D+0 Coeff: 0.518D+00 Coeff: 0.518D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=2.32D-04 MaxDP=1.87D-02 DE=-2.21D-03 OVMax= 4.87D-03 | RMSDP=2.32D-04 MaxDP=1.87D-02 DE=-2.22D-03 OVMax= 4.87D-03 Cycle 8 Pass 1 IDiag 1: Cycle 8 Pass 1 IDiag 1: RMSU= 1.18D-04 CP: 1.01D+00 1.88D-01 5.55D-01 4.56D- | RMSU= 1.19D-04 CP: 1.01D+00 1.87D-01 5.54D-01 4.56D- CP: 5.85D-01 5.99D-01 | CP: 5.86D-01 5.99D-01 E= -613.621363784131 Delta-E= -0.000062661817 Rise | E= -613.621402882576 Delta-E= -0.000062701855 Rise DIIS: error= 1.42D-04 at cycle 8 NSaved= 8. DIIS: error= 1.42D-04 at cycle 8 NSaved= 8. NSaved= 8 IEnMin= 8 EnMin= -613.621363784131 IErMin= 8 E | NSaved= 8 IEnMin= 8 EnMin= -613.621402882576 IErMin= 8 E ErrMax= 1.42D-04 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.00D-05 B ErrMax= 1.42D-04 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.00D-05 B IDIUse=3 WtCom= 9.99D-01 WtEn= 1.42D-03 IDIUse=3 WtCom= 9.99D-01 WtEn= 1.42D-03 Coeff-Com: -0.217D-03 0.139D-02 0.510D-02-0.525D-02 0.226D-0 | Coeff-Com: -0.217D-03 0.139D-02 0.511D-02-0.525D-02 0.227D-0 Coeff-Com: 0.338D+00 0.612D+00 Coeff-Com: 0.338D+00 0.612D+00 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+0 Coeff-En: 0.000D+00 0.000D+00 0.000D+00 0.000D+00 0.000D+0 Coeff-En: 0.217D+00 0.783D+00 Coeff-En: 0.217D+00 0.783D+00 Coeff: -0.217D-03 0.139D-02 0.510D-02-0.524D-02 0.226D-0 Coeff: -0.217D-03 0.139D-02 0.510D-02-0.524D-02 0.226D-0 Coeff: 0.338D+00 0.612D+00 Coeff: 0.338D+00 0.612D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=1.13D-04 MaxDP=8.66D-03 DE=-6.27D-05 OVMax= 2.06D-03 | RMSDP=1.13D-04 MaxDP=8.68D-03 DE=-6.27D-05 OVMax= 2.06D-03 Cycle 9 Pass 1 IDiag 1: Cycle 9 Pass 1 IDiag 1: RMSU= 4.02D-05 CP: 1.01D+00 1.90D-01 5.59D-01 4.53D- | RMSU= 4.02D-05 CP: 1.01D+00 1.89D-01 5.58D-01 4.53D- CP: 5.75D-01 7.36D-01 7.28D-01 | CP: 5.76D-01 7.37D-01 7.29D-01 E= -613.621384289805 Delta-E= -0.000020505674 Rise | E= -613.621423396384 Delta-E= -0.000020513808 Rise DIIS: error= 4.24D-05 at cycle 9 NSaved= 9. | DIIS: error= 4.25D-05 at cycle 9 NSaved= 9. NSaved= 9 IEnMin= 9 EnMin= -613.621384289805 IErMin= 9 E | NSaved= 9 IEnMin= 9 EnMin= -613.621423396384 IErMin= 9 E ErrMax= 4.24D-05 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.21D-06 B | ErrMax= 4.25D-05 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.21D-06 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: 0.315D-03-0.465D-03-0.669D-03-0.526D-03-0.709D-0 | Coeff-Com: 0.315D-03-0.465D-03-0.667D-03-0.527D-03-0.708D-0 Coeff-Com: 0.182D+00 0.391D+00 0.439D+00 Coeff-Com: 0.182D+00 0.391D+00 0.439D+00 Coeff: 0.315D-03-0.465D-03-0.669D-03-0.526D-03-0.709D-0 | Coeff: 0.315D-03-0.465D-03-0.667D-03-0.527D-03-0.708D-0 Coeff: 0.182D+00 0.391D+00 0.439D+00 Coeff: 0.182D+00 0.391D+00 0.439D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=2.16D-05 MaxDP=8.92D-04 DE=-2.05D-05 OVMax= 5.12D-04 | RMSDP=2.16D-05 MaxDP=8.93D-04 DE=-2.05D-05 OVMax= 5.12D-04 Cycle 10 Pass 1 IDiag 1: Cycle 10 Pass 1 IDiag 1: RMSU= 1.48D-05 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 1.48D-05 CP: 1.01D+00 1.89D-01 5.57D-01 4.53D- CP: 5.74D-01 7.48D-01 7.74D-01 7.28D- | CP: 5.74D-01 7.48D-01 7.74D-01 7.27D- E= -613.621386188079 Delta-E= -0.000001898274 Rise | E= -613.621425297075 Delta-E= -0.000001900691 Rise DIIS: error= 8.42D-06 at cycle 10 NSaved= 10. DIIS: error= 8.42D-06 at cycle 10 NSaved= 10. NSaved=10 IEnMin=10 EnMin= -613.621386188079 IErMin=10 E | NSaved=10 IEnMin=10 EnMin= -613.621425297075 IErMin=10 E ErrMax= 8.42D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.22D-07 B ErrMax= 8.42D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.22D-07 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: 0.177D-03-0.430D-03-0.105D-02 0.397D-03-0.611D-0 | Coeff-Com: 0.177D-03-0.430D-03-0.105D-02 0.397D-03-0.610D-0 Coeff-Com: 0.471D-01 0.106D+00 0.224D+00 0.633D+00 Coeff-Com: 0.471D-01 0.106D+00 0.224D+00 0.633D+00 Coeff: 0.177D-03-0.430D-03-0.105D-02 0.397D-03-0.611D-0 | Coeff: 0.177D-03-0.430D-03-0.105D-02 0.397D-03-0.610D-0 Coeff: 0.471D-01 0.106D+00 0.224D+00 0.633D+00 Coeff: 0.471D-01 0.106D+00 0.224D+00 0.633D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=1.36D-05 MaxDP=1.31D-03 DE=-1.90D-06 OVMax= 2.23D-04 RMSDP=1.36D-05 MaxDP=1.31D-03 DE=-1.90D-06 OVMax= 2.23D-04 Cycle 11 Pass 1 IDiag 1: Cycle 11 Pass 1 IDiag 1: RMSU= 4.24D-06 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 4.25D-06 CP: 1.01D+00 1.88D-01 5.58D-01 4.53D- CP: 5.77D-01 7.51D-01 7.71D-01 7.06D- | CP: 5.78D-01 7.51D-01 7.72D-01 7.06D- E= -613.621386322536 Delta-E= -0.000000134457 Rise | E= -613.621425431487 Delta-E= -0.000000134411 Rise DIIS: error= 3.37D-06 at cycle 11 NSaved= 11. DIIS: error= 3.37D-06 at cycle 11 NSaved= 11. NSaved=11 IEnMin=11 EnMin= -613.621386322536 IErMin=11 E | NSaved=11 IEnMin=11 EnMin= -613.621425431487 IErMin=11 E ErrMax= 3.37D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.18D-08 B ErrMax= 3.37D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.18D-08 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: 0.513D-04-0.122D-03-0.601D-03 0.244D-03-0.264D-0 Coeff-Com: 0.513D-04-0.122D-03-0.601D-03 0.244D-03-0.264D-0 Coeff-Com: 0.861D-02 0.170D-01 0.692D-01 0.277D+00 0.633D+0 | Coeff-Com: 0.862D-02 0.171D-01 0.692D-01 0.277D+00 0.633D+0 Coeff: 0.513D-04-0.122D-03-0.601D-03 0.244D-03-0.264D-0 Coeff: 0.513D-04-0.122D-03-0.601D-03 0.244D-03-0.264D-0 Coeff: 0.861D-02 0.170D-01 0.692D-01 0.277D+00 0.633D+0 | Coeff: 0.862D-02 0.171D-01 0.692D-01 0.277D+00 0.633D+0 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=3.16D-06 MaxDP=2.06D-04 DE=-1.34D-07 OVMax= 4.74D-05 | RMSDP=3.17D-06 MaxDP=2.07D-04 DE=-1.34D-07 OVMax= 4.74D-05 Cycle 12 Pass 1 IDiag 1: Cycle 12 Pass 1 IDiag 1: RMSU= 1.88D-06 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 1.89D-06 CP: 1.01D+00 1.88D-01 5.57D-01 4.53D- CP: 5.77D-01 7.50D-01 7.75D-01 7.22D- | CP: 5.77D-01 7.50D-01 7.76D-01 7.21D- CP: 1.05D+00 CP: 1.05D+00 E= -613.621386333573 Delta-E= -0.000000011037 Rise | E= -613.621425442544 Delta-E= -0.000000011057 Rise DIIS: error= 2.48D-06 at cycle 12 NSaved= 12. DIIS: error= 2.48D-06 at cycle 12 NSaved= 12. NSaved=12 IEnMin=12 EnMin= -613.621386333573 IErMin=12 E | NSaved=12 IEnMin=12 EnMin= -613.621425442544 IErMin=12 E ErrMax= 2.48D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.75D-09 B ErrMax= 2.48D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.75D-09 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: -0.162D-04 0.396D-04-0.138D-03 0.385D-04-0.150D-0 | Coeff-Com: -0.162D-04 0.396D-04-0.139D-03 0.386D-04-0.151D-0 Coeff-Com: -0.290D-02-0.101D-01-0.363D-02 0.308D-01 0.329D+0 | Coeff-Com: -0.290D-02-0.101D-01-0.362D-02 0.309D-01 0.329D+0 Coeff: -0.162D-04 0.396D-04-0.138D-03 0.385D-04-0.150D-0 | Coeff: -0.162D-04 0.396D-04-0.139D-03 0.386D-04-0.151D-0 Coeff: -0.290D-02-0.101D-01-0.363D-02 0.308D-01 0.329D+0 | Coeff: -0.290D-02-0.101D-01-0.362D-02 0.309D-01 0.329D+0 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=1.67D-06 MaxDP=1.85D-04 DE=-1.10D-08 OVMax= 2.96D-05 | RMSDP=1.68D-06 MaxDP=1.85D-04 DE=-1.11D-08 OVMax= 2.96D-05 Cycle 13 Pass 1 IDiag 1: Cycle 13 Pass 1 IDiag 1: RMSU= 8.93D-07 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 8.93D-07 CP: 1.01D+00 1.88D-01 5.57D-01 4.53D- CP: 5.77D-01 7.52D-01 7.73D-01 7.30D- | CP: 5.77D-01 7.52D-01 7.73D-01 7.29D- CP: 1.02D+00 7.52D-01 CP: 1.02D+00 7.52D-01 E= -613.621386336648 Delta-E= -0.000000003075 Rise | E= -613.621425445608 Delta-E= -0.000000003064 Rise DIIS: error= 1.25D-06 at cycle 13 NSaved= 13. DIIS: error= 1.25D-06 at cycle 13 NSaved= 13. NSaved=13 IEnMin=13 EnMin= -613.621386336648 IErMin=13 E | NSaved=13 IEnMin=13 EnMin= -613.621425445608 IErMin=13 E ErrMax= 1.25D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.48D-10 B ErrMax= 1.25D-06 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.48D-10 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: -0.166D-04 0.382D-04-0.153D-04-0.575D-05 0.207D-0 | Coeff-Com: -0.166D-04 0.381D-04-0.154D-04-0.573D-05 0.207D-0 Coeff-Com: -0.228D-02-0.731D-02-0.924D-02-0.971D-02 0.107D+0 Coeff-Com: -0.228D-02-0.731D-02-0.924D-02-0.971D-02 0.107D+0 Coeff-Com: 0.550D+00 Coeff-Com: 0.550D+00 Coeff: -0.166D-04 0.382D-04-0.153D-04-0.575D-05 0.207D-0 | Coeff: -0.166D-04 0.381D-04-0.154D-04-0.573D-05 0.207D-0 Coeff: -0.228D-02-0.731D-02-0.924D-02-0.971D-02 0.107D+0 Coeff: -0.228D-02-0.731D-02-0.924D-02-0.971D-02 0.107D+0 Coeff: 0.550D+00 Coeff: 0.550D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=5.28D-07 MaxDP=3.89D-05 DE=-3.08D-09 OVMax= 1.00D-05 | RMSDP=5.29D-07 MaxDP=3.91D-05 DE=-3.06D-09 OVMax= 1.00D-05 Cycle 14 Pass 1 IDiag 1: Cycle 14 Pass 1 IDiag 1: RMSU= 2.40D-07 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 2.41D-07 CP: 1.01D+00 1.88D-01 5.57D-01 4.53D- CP: 5.77D-01 7.52D-01 7.74D-01 7.28D- | CP: 5.77D-01 7.52D-01 7.74D-01 7.27D- CP: 1.06D+00 8.74D-01 7.06D-01 | CP: 1.06D+00 8.75D-01 7.07D-01 E= -613.621386337023 Delta-E= -0.000000000375 Rise | E= -613.621425445998 Delta-E= -0.000000000390 Rise DIIS: error= 3.62D-07 at cycle 14 NSaved= 14. DIIS: error= 3.62D-07 at cycle 14 NSaved= 14. NSaved=14 IEnMin=14 EnMin= -613.621386337023 IErMin=14 E | NSaved=14 IEnMin=14 EnMin= -613.621425445998 IErMin=14 E ErrMax= 3.62D-07 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.66D-11 B ErrMax= 3.62D-07 0.00D+00 EMaxC= 1.00D-01 BMatC= 3.66D-11 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: -0.681D-05 0.157D-04 0.744D-05-0.735D-05 0.114D-0 | Coeff-Com: -0.681D-05 0.157D-04 0.743D-05-0.734D-05 0.114D-0 Coeff-Com: -0.743D-03-0.250D-02-0.418D-02-0.731D-02 0.934D-0 | Coeff-Com: -0.743D-03-0.249D-02-0.418D-02-0.731D-02 0.935D-0 Coeff-Com: 0.293D+00 0.605D+00 Coeff-Com: 0.293D+00 0.605D+00 Coeff: -0.681D-05 0.157D-04 0.744D-05-0.735D-05 0.114D-0 | Coeff: -0.681D-05 0.157D-04 0.743D-05-0.734D-05 0.114D-0 Coeff: -0.743D-03-0.250D-02-0.418D-02-0.731D-02 0.934D-0 | Coeff: -0.743D-03-0.249D-02-0.418D-02-0.731D-02 0.935D-0 Coeff: 0.293D+00 0.605D+00 Coeff: 0.293D+00 0.605D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=3.43D-07 MaxDP=4.54D-05 DE=-3.75D-10 OVMax= 2.64D-06 | RMSDP=3.42D-07 MaxDP=4.54D-05 DE=-3.90D-10 OVMax= 2.64D-06 Cycle 15 Pass 1 IDiag 1: Cycle 15 Pass 1 IDiag 1: RMSU= 9.32D-08 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 9.32D-08 CP: 1.01D+00 1.88D-01 5.57D-01 4.53D- CP: 5.77D-01 7.52D-01 7.74D-01 7.27D- CP: 5.77D-01 7.52D-01 7.74D-01 7.27D- CP: 1.07D+00 8.47D-01 7.49D-01 9.51D- | CP: 1.07D+00 8.48D-01 7.50D-01 9.51D- E= -613.621386337062 Delta-E= -0.000000000038 Rise | E= -613.621425446035 Delta-E= -0.000000000037 Rise DIIS: error= 5.04D-08 at cycle 15 NSaved= 15. DIIS: error= 5.04D-08 at cycle 15 NSaved= 15. NSaved=15 IEnMin=15 EnMin= -613.621386337062 IErMin=15 E | NSaved=15 IEnMin=15 EnMin= -613.621425446035 IErMin=15 E ErrMax= 5.04D-08 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.53D-12 B ErrMax= 5.04D-08 0.00D+00 EMaxC= 1.00D-01 BMatC= 2.53D-12 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: -0.100D-05 0.258D-05 0.516D-05-0.255D-05 0.254D-0 | Coeff-Com: -0.101D-05 0.259D-05 0.516D-05-0.255D-05 0.254D-0 Coeff-Com: -0.926D-04-0.392D-03-0.790D-03-0.142D-02-0.771D-0 | Coeff-Com: -0.928D-04-0.393D-03-0.792D-03-0.142D-02-0.770D-0 Coeff-Com: 0.553D-01 0.215D+00 0.735D+00 | Coeff-Com: 0.554D-01 0.215D+00 0.735D+00 Coeff: -0.100D-05 0.258D-05 0.516D-05-0.255D-05 0.254D-0 | Coeff: -0.101D-05 0.259D-05 0.516D-05-0.255D-05 0.254D-0 Coeff: -0.926D-04-0.392D-03-0.790D-03-0.142D-02-0.771D-0 | Coeff: -0.928D-04-0.393D-03-0.792D-03-0.142D-02-0.770D-0 Coeff: 0.553D-01 0.215D+00 0.735D+00 | Coeff: 0.554D-01 0.215D+00 0.735D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=5.85D-08 MaxDP=3.27D-06 DE=-3.84D-11 OVMax= 1.01D-06 | RMSDP=5.87D-08 MaxDP=3.30D-06 DE=-3.71D-11 OVMax= 1.01D-06 Cycle 16 Pass 1 IDiag 1: Cycle 16 Pass 1 IDiag 1: RMSU= 2.64D-08 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 2.65D-08 CP: 1.01D+00 1.88D-01 5.57D-01 4.53D- CP: 5.77D-01 7.52D-01 7.74D-01 7.27D- CP: 5.77D-01 7.52D-01 7.74D-01 7.27D- CP: 1.07D+00 8.57D-01 7.75D-01 9.72D- | CP: 1.07D+00 8.58D-01 7.76D-01 9.71D- E= -613.621386337060 Delta-E= 0.000000000001 Rise | E= -613.621425446027 Delta-E= 0.000000000008 Rise DIIS: error= 2.03D-08 at cycle 16 NSaved= 16. DIIS: error= 2.03D-08 at cycle 16 NSaved= 16. NSaved=16 IEnMin=15 EnMin= -613.621386337062 IErMin=16 E | NSaved=16 IEnMin=15 EnMin= -613.621425446035 IErMin=16 E ErrMax= 2.03D-08 0.00D+00 EMaxC= 1.00D-01 BMatC= 5.19D-13 B ErrMax= 2.03D-08 0.00D+00 EMaxC= 1.00D-01 BMatC= 5.19D-13 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: 0.349D-06-0.925D-06 0.220D-05-0.307D-06 0.380D-0 | Coeff-Com: 0.348D-06-0.922D-06 0.220D-05-0.309D-06 0.399D-0 Coeff-Com: 0.452D-04 0.657D-04 0.944D-04 0.384D-03-0.393D-0 | Coeff-Com: 0.451D-04 0.654D-04 0.940D-04 0.383D-03-0.393D-0 Coeff-Com: -0.576D-02 0.260D-01 0.326D+00 0.664D+00 | Coeff-Com: -0.572D-02 0.261D-01 0.327D+00 0.664D+00 Coeff: 0.349D-06-0.925D-06 0.220D-05-0.307D-06 0.380D-0 | Coeff: 0.348D-06-0.922D-06 0.220D-05-0.309D-06 0.399D-0 Coeff: 0.452D-04 0.657D-04 0.944D-04 0.384D-03-0.393D-0 | Coeff: 0.451D-04 0.654D-04 0.940D-04 0.383D-03-0.393D-0 Coeff: -0.576D-02 0.260D-01 0.326D+00 0.664D+00 | Coeff: -0.572D-02 0.261D-01 0.327D+00 0.664D+00 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=3.67D-08 MaxDP=3.84D-06 DE= 1.14D-12 OVMax= 3.26D-07 | RMSDP=3.67D-08 MaxDP=3.84D-06 DE= 7.50D-12 OVMax= 3.26D-07 Cycle 17 Pass 1 IDiag 1: Cycle 17 Pass 1 IDiag 1: RMSU= 9.59D-09 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 9.61D-09 CP: 1.01D+00 1.88D-01 5.57D-01 4.53D- CP: 5.77D-01 7.52D-01 7.73D-01 7.27D- | CP: 5.77D-01 7.52D-01 7.74D-01 7.27D- CP: 1.07D+00 8.53D-01 7.76D-01 9.72D- | CP: 1.07D+00 8.54D-01 7.76D-01 9.71D- CP: 1.05D+00 CP: 1.05D+00 E= -613.621386337075 Delta-E= -0.000000000015 Rise | E= -613.621425446053 Delta-E= -0.000000000025 Rise DIIS: error= 8.84D-09 at cycle 17 NSaved= 17. | DIIS: error= 8.85D-09 at cycle 17 NSaved= 17. NSaved=17 IEnMin=17 EnMin= -613.621386337075 IErMin=17 E | NSaved=17 IEnMin=17 EnMin= -613.621425446053 IErMin=17 E ErrMax= 8.84D-09 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.02D-13 B | ErrMax= 8.85D-09 0.00D+00 EMaxC= 1.00D-01 BMatC= 1.02D-13 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: 0.419D-06-0.877D-06 0.355D-06 0.223D-06-0.436D-0 Coeff-Com: 0.419D-06-0.877D-06 0.355D-06 0.223D-06-0.436D-0 Coeff-Com: 0.290D-04 0.761D-04 0.141D-03 0.404D-03-0.831D-0 | Coeff-Com: 0.291D-04 0.761D-04 0.141D-03 0.404D-03-0.832D-0 Coeff-Com: -0.122D-01-0.229D-01 0.401D-01 0.296D+00 0.704D+0 | Coeff-Com: -0.122D-01-0.230D-01 0.401D-01 0.296D+00 0.705D+0 Coeff: 0.419D-06-0.877D-06 0.355D-06 0.223D-06-0.436D-0 Coeff: 0.419D-06-0.877D-06 0.355D-06 0.223D-06-0.436D-0 Coeff: 0.290D-04 0.761D-04 0.141D-03 0.404D-03-0.831D-0 | Coeff: 0.291D-04 0.761D-04 0.141D-03 0.404D-03-0.832D-0 Coeff: -0.122D-01-0.229D-01 0.401D-01 0.296D+00 0.704D+0 | Coeff: -0.122D-01-0.230D-01 0.401D-01 0.296D+00 0.705D+0 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=1.03D-08 MaxDP=8.91D-07 DE=-1.48D-11 OVMax= 1.57D-07 | RMSDP=1.03D-08 MaxDP=8.90D-07 DE=-2.55D-11 OVMax= 1.58D-07 Cycle 18 Pass 1 IDiag 1: Cycle 18 Pass 1 IDiag 1: RMSU= 2.76D-09 CP: 1.01D+00 1.89D-01 5.59D-01 4.53D- | RMSU= 2.77D-09 CP: 1.01D+00 1.88D-01 5.57D-01 4.53D- CP: 5.77D-01 7.52D-01 7.73D-01 7.27D- | CP: 5.77D-01 7.52D-01 7.74D-01 7.27D- CP: 1.07D+00 8.53D-01 7.77D-01 9.80D- | CP: 1.07D+00 8.54D-01 7.78D-01 9.79D- CP: 1.10D+00 1.03D+00 CP: 1.10D+00 1.03D+00 E= -613.621386337061 Delta-E= 0.000000000015 Rise | E= -613.621425446031 Delta-E= 0.000000000021 Rise DIIS: error= 2.57D-09 at cycle 18 NSaved= 18. DIIS: error= 2.57D-09 at cycle 18 NSaved= 18. NSaved=18 IEnMin=17 EnMin= -613.621386337075 IErMin=18 E | NSaved=18 IEnMin=17 EnMin= -613.621425446053 IErMin=18 E ErrMax= 2.57D-09 0.00D+00 EMaxC= 1.00D-01 BMatC= 6.48D-15 B | ErrMax= 2.57D-09 0.00D+00 EMaxC= 1.00D-01 BMatC= 6.50D-15 B IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00 Coeff-Com: 0.139D-06-0.268D-06 0.487D-07 0.907D-07-0.141D-0 | Coeff-Com: 0.139D-06-0.268D-06 0.489D-07 0.907D-07-0.141D-0 Coeff-Com: 0.764D-05 0.191D-04 0.393D-04 0.144D-03 0.111D-0 | Coeff-Com: 0.765D-05 0.191D-04 0.393D-04 0.144D-03 0.110D-0 Coeff-Com: -0.339D-02-0.920D-02-0.655D-02 0.585D-01 0.296D+0 | Coeff-Com: -0.340D-02-0.921D-02-0.653D-02 0.585D-01 0.297D+0 Coeff: 0.139D-06-0.268D-06 0.487D-07 0.907D-07-0.141D-0 | Coeff: 0.139D-06-0.268D-06 0.489D-07 0.907D-07-0.141D-0 Coeff: 0.764D-05 0.191D-04 0.393D-04 0.144D-03 0.111D-0 | Coeff: 0.765D-05 0.191D-04 0.393D-04 0.144D-03 0.110D-0 Coeff: -0.339D-02-0.920D-02-0.655D-02 0.585D-01 0.296D+0 | Coeff: -0.340D-02-0.921D-02-0.653D-02 0.585D-01 0.297D+0 Gap= 0.100 Goal= None Shift= 0.000 Gap= 0.100 Goal= None Shift= 0.000 RMSDP=2.51D-09 MaxDP=1.90D-07 DE= 1.46D-11 OVMax= 2.65D-08 | RMSDP=2.48D-09 MaxDP=1.81D-07 DE= 2.14D-11 OVMax= 2.65D-08 SCF Done: E(RB3LYP) = -613.621386337 A.U. after 18 c | SCF Done: E(RB3LYP) = -613.621425446 A.U. after 18 c NFock= 18 Conv=0.25D-08 -V/T= 2.0098 NFock= 18 Conv=0.25D-08 -V/T= 2.0098 KE= 6.076802317978D+02 PE=-3.304611718755D+03 EE= 1.13327067 | KE= 6.076801358457D+02 PE=-3.304611298029D+03 EE= 1.13327031 Leave Link 502 at Wed Apr 19 20:05:11 2017, MaxMem= 131072 | Leave Link 502 at Wed Apr 19 23:52:06 2017, MaxMem= 131072 (Enter /usr/g16/l801.exe) | (Enter /usr/g09/l801.exe) DoSCS=F DFT=T ScalE2(SS,OS)= 1.000000 1.000000 DoSCS=F DFT=T ScalE2(SS,OS)= 1.000000 1.000000 Range of M.O.s used for correlation: 1 360 Range of M.O.s used for correlation: 1 360 NBasis= 360 NAE= 53 NBE= 53 NFC= 0 NFV= 0 NBasis= 360 NAE= 53 NBE= 53 NFC= 0 NFV= 0 NROrb= 360 NOA= 53 NOB= 53 NVA= 307 NVB= 307 NROrb= 360 NOA= 53 NOB= 53 NVA= 307 NVB= 307 **** Warning!!: The largest alpha MO coefficient is 0.14193 | **** Warning!!: The largest alpha MO coefficient is 0.14191 |
また、構造最適化に要したサイクル数は、g16 と g09 で ±1 cycle くらいの差なので、ほぼ同じと見なして良さそうです。どちらが早いかはその時により異なります。
これだけ見ると、g16 になってもあまり高速化されてなさそうな印象を受けますが、最初の vector の計算は g16 の方が圧倒的に収束が早かったです。(今回はcalcfc を指定しました。)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 |
There are 93 degrees of freedom in the 1st order There are 93 degrees of freedom in the 1st order 87 vectors produced by pass 0 Test12= 1.75D-14 1.08D-09 | 89 vectors produced by pass 0 Test12= 1.75D-14 1.08D-09 AX will form 87 AO Fock derivatives at one time. | AX will form 89 AO Fock derivatives at one time. 87 vectors produced by pass 1 Test12= 1.75D-14 1.08D-09 | 87 vectors produced by pass 1 Test12= 1.75D-14 1.08D-09 87 vectors produced by pass 2 Test12= 1.75D-14 1.08D-09 | 87 vectors produced by pass 2 Test12= 1.75D-14 1.08D-09 87 vectors produced by pass 3 Test12= 1.75D-14 1.08D-09 | 87 vectors produced by pass 3 Test12= 1.75D-14 1.08D-09 87 vectors produced by pass 4 Test12= 1.75D-14 1.08D-09 | 87 vectors produced by pass 4 Test12= 1.75D-14 1.08D-09 81 vectors produced by pass 5 Test12= 1.75D-14 1.08D-09 | 82 vectors produced by pass 5 Test12= 1.75D-14 1.08D-09 17 vectors produced by pass 6 Test12= 1.75D-14 1.08D-09 | 55 vectors produced by pass 6 Test12= 1.75D-14 1.08D-09 5 vectors produced by pass 7 Test12= 1.75D-14 1.08D-09 | 21 vectors produced by pass 7 Test12= 1.75D-14 1.08D-09 3 vectors produced by pass 8 Test12= 1.75D-14 1.08D-09 | 16 vectors produced by pass 8 Test12= 1.75D-14 1.08D-09 3 vectors produced by pass 9 Test12= 1.75D-14 1.08D-09 | 16 vectors produced by pass 9 Test12= 1.75D-14 1.08D-09 2 vectors produced by pass 10 Test12= 1.75D-14 1.08D-09 | 15 vectors produced by pass 10 Test12= 1.75D-14 1.08D-09 InvSVY: IOpt=1 It= 1 EMax= 7.77D-16 | 15 vectors produced by pass 11 Test12= 1.75D-14 1.08D-09 Solved reduced A of dimension 546 with 93 vectors. | 15 vectors produced by pass 12 Test12= 1.75D-14 1.08D-09 > 12 vectors produced by pass 13 Test12= 1.75D-14 1.08D-09 > 12 vectors produced by pass 14 Test12= 1.75D-14 1.08D-09 > 10 vectors produced by pass 15 Test12= 1.75D-14 1.08D-09 > 9 vectors produced by pass 16 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 17 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 18 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 19 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 20 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 21 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 22 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 23 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 24 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 25 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 26 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 27 Test12= 1.75D-14 1.08D-09 > 7 vectors produced by pass 28 Test12= 1.75D-14 1.08D-09 > 6 vectors produced by pass 29 Test12= 1.75D-14 1.08D-09 > 6 vectors produced by pass 30 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 31 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 32 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 33 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 34 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 35 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 36 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 37 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 38 Test12= 1.75D-14 1.08D-09 > 4 vectors produced by pass 39 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 40 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 41 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 42 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 43 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 44 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 45 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 46 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 47 Test12= 1.75D-14 1.08D-09 > 3 vectors produced by pass 48 Test12= 1.75D-14 1.08D-09 > 2 vectors produced by pass 49 Test12= 1.75D-14 1.08D-09 > InvSVY: IOpt=1 It= 1 EMax= 1.05D-15 > Solved reduced A of dimension 876 with 93 vectors. End of Minotr F.D. properties file 721 does not exist. End of Minotr F.D. properties file 721 does not exist. End of Minotr F.D. properties file 722 does not exist. End of Minotr F.D. properties file 722 does not exist. End of Minotr F.D. properties file 788 does not exist. End of Minotr F.D. properties file 788 does not exist. Leave Link 1002 at Wed Apr 19 20:19:54 2017, MaxMem= 131072 | Leave Link 1002 at Thu Apr 20 00:10:57 2017, MaxMem= 131072 (Enter /usr/g16/l601.exe) | (Enter /usr/g09/l601.exe) Copying SCF densities to generalized density rwf, IOpCl= 0 I Copying SCF densities to generalized density rwf, IOpCl= 0 I |
g09 で 49 サイクルもかかっているのに、g16 では 10 サイクルで収束しました。
最後に、Optimized Parameter を比較すると g16 と g09 で広範囲にわたり差があります。小数点以下の誤差ですが、少し気になります。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 |
---------------------------- ---------------------------- ! Optimized Parameters ! ! Optimized Parameters ! ! (Angstroms and Degrees) ! ! (Angstroms and Degrees) ! -------------------------- ------ -------------------------- ------ ! Name Definition Value Derivative In ! Name Definition Value Derivative In ------------------------------------------------------------ ------------------------------------------------------------ ! R1 R(1,2) 1.2781 -DE/DX = 0 | ! R1 R(1,2) 1.278 -DE/DX = 0 ! R2 R(1,15) 1.5465 -DE/DX = 0 | ! R2 R(1,15) 1.5469 -DE/DX = 0 ! R3 R(1,16) 1.4965 -DE/DX = 0 ! R3 R(1,16) 1.4965 -DE/DX = 0 ! R4 R(2,3) 1.3138 -DE/DX = 0 ! R4 R(2,3) 1.3138 -DE/DX = 0 ! R5 R(3,4) 1.4073 -DE/DX = 0 ! R5 R(3,4) 1.4073 -DE/DX = 0 ! R6 R(3,21) 1.2551 -DE/DX = 0 | ! R6 R(3,21) 1.2548 -DE/DX = 0 ! R7 R(4,5) 1.4095 -DE/DX = 0 | ! R7 R(4,5) 1.4094 -DE/DX = 0 ! R8 R(4,6) 1.4099 -DE/DX = 0 ! R8 R(4,6) 1.4099 -DE/DX = 0 ! R9 R(5,7) 1.3901 -DE/DX = 0 | ! R9 R(5,7) 1.3902 -DE/DX = 0 ! R10 R(5,8) 1.0844 -DE/DX = 0 ! R10 R(5,8) 1.0844 -DE/DX = 0 ! R11 R(6,9) 1.3894 -DE/DX = 0 | ! R11 R(6,9) 1.3893 -DE/DX = 0 ! R12 R(6,10) 1.0857 -DE/DX = 0 ! R12 R(6,10) 1.0857 -DE/DX = 0 ! R13 R(7,11) 1.4023 -DE/DX = 0 | ! R13 R(7,11) 1.4022 -DE/DX = 0 ! R14 R(7,12) 1.085 -DE/DX = 0 ! R14 R(7,12) 1.085 -DE/DX = 0 ! R15 R(9,11) 1.4018 -DE/DX = 0 ! R15 R(9,11) 1.4018 -DE/DX = 0 ! R16 R(9,13) 1.0848 -DE/DX = 0 ! R16 R(9,13) 1.0848 -DE/DX = 0 ! R17 R(11,14) 1.0855 -DE/DX = 0 ! R17 R(11,14) 1.0855 -DE/DX = 0 ! R18 R(15,20) 1.509 -DE/DX = 0 | ! R18 R(15,20) 1.5092 -DE/DX = 0 ! R19 R(15,26) 1.0976 -DE/DX = 0 ! R19 R(15,26) 1.0976 -DE/DX = 0 ! R20 R(15,28) 1.5626 -DE/DX = 0 | ! R20 R(15,28) 1.5628 -DE/DX = 0 ! R21 R(16,17) 1.0956 -DE/DX = 0 | ! R21 R(16,17) 1.0955 -DE/DX = 0 ! R22 R(16,18) 1.0966 -DE/DX = 0 | ! R22 R(16,18) 1.0967 -DE/DX = 0 ! R23 R(16,19) 1.0905 -DE/DX = 0 ! R23 R(16,19) 1.0905 -DE/DX = 0 ! R24 R(20,21) 1.3369 -DE/DX = 0 | ! R24 R(20,21) 1.3365 -DE/DX = 0 ! R25 R(20,22) 1.4548 -DE/DX = 0 ! R25 R(20,22) 1.4548 -DE/DX = 0 ! R26 R(20,27) 1.0932 -DE/DX = 0 ! R26 R(20,27) 1.0932 -DE/DX = 0 ! R27 R(22,23) 1.3554 -DE/DX = 0 ! R27 R(22,23) 1.3554 -DE/DX = 0 ! R28 R(22,24) 1.0838 -DE/DX = 0 ! R28 R(22,24) 1.0838 -DE/DX = 0 ! R29 R(23,25) 1.0853 -DE/DX = 0 ! R29 R(23,25) 1.0853 -DE/DX = 0 ! R30 R(23,28) 1.5013 -DE/DX = 0 | ! R30 R(23,28) 1.5012 -DE/DX = 0 ! R31 R(28,29) 1.0957 -DE/DX = 0 | ! R31 R(28,29) 1.0956 -DE/DX = 0 ! R32 R(28,30) 1.0988 -DE/DX = 0 | ! R32 R(28,30) 1.0989 -DE/DX = 0 ! A1 A(2,1,15) 120.3159 -DE/DX = 0 | ! A1 A(2,1,15) 120.3319 -DE/DX = 0 ! A2 A(2,1,16) 120.7007 -DE/DX = 0 | ! A2 A(2,1,16) 120.7074 -DE/DX = 0 ! A3 A(15,1,16) 118.983 -DE/DX = 0 | ! A3 A(15,1,16) 118.9603 -DE/DX = 0 ! A4 A(1,2,3) 126.076 -DE/DX = 0 | ! A4 A(1,2,3) 126.081 -DE/DX = 0 ! A5 A(2,3,4) 118.2699 -DE/DX = 0 | ! A5 A(2,3,4) 118.2819 -DE/DX = 0 ! A6 A(2,3,21) 106.4974 -DE/DX = 0 | ! A6 A(2,3,21) 106.4607 -DE/DX = 0 ! A7 A(4,3,21) 110.9882 -DE/DX = 0 | ! A7 A(4,3,21) 110.8722 -DE/DX = 0 ! A8 A(3,4,5) 123.1933 -DE/DX = 0 | ! A8 A(3,4,5) 123.2064 -DE/DX = 0 ! A9 A(3,4,6) 115.9886 -DE/DX = 0 | ! A9 A(3,4,6) 115.9761 -DE/DX = 0 ! A10 A(5,4,6) 120.8119 -DE/DX = 0 | ! A10 A(5,4,6) 120.8135 -DE/DX = 0 ! A11 A(4,5,7) 119.0668 -DE/DX = 0 | ! A11 A(4,5,7) 119.0755 -DE/DX = 0 ! A12 A(4,5,8) 120.0833 -DE/DX = 0 | ! A12 A(4,5,8) 120.0714 -DE/DX = 0 ! A13 A(7,5,8) 120.8085 -DE/DX = 0 | ! A13 A(7,5,8) 120.8131 -DE/DX = 0 ! A14 A(4,6,9) 119.446 -DE/DX = 0 | ! A14 A(4,6,9) 119.4374 -DE/DX = 0 ! A15 A(4,6,10) 119.4206 -DE/DX = 0 | ! A15 A(4,6,10) 119.425 -DE/DX = 0 ! A16 A(9,6,10) 121.1275 -DE/DX = 0 | ! A16 A(9,6,10) 121.1321 -DE/DX = 0 ! A17 A(5,7,11) 120.2207 -DE/DX = 0 | ! A17 A(5,7,11) 120.2121 -DE/DX = 0 ! A18 A(5,7,12) 119.8001 -DE/DX = 0 | ! A18 A(5,7,12) 119.8013 -DE/DX = 0 ! A19 A(11,7,12) 119.9783 -DE/DX = 0 | ! A19 A(11,7,12) 119.9857 -DE/DX = 0 ! A20 A(6,9,11) 119.8766 -DE/DX = 0 | ! A20 A(6,9,11) 119.8871 -DE/DX = 0 ! A21 A(6,9,13) 119.9575 -DE/DX = 0 | ! A21 A(6,9,13) 119.9522 -DE/DX = 0 ! A22 A(11,9,13) 120.1576 -DE/DX = 0 | ! A22 A(11,9,13) 120.1527 -DE/DX = 0 ! A23 A(7,11,9) 120.5685 -DE/DX = 0 | ! A23 A(7,11,9) 120.5651 -DE/DX = 0 ! A24 A(7,11,14) 119.6904 -DE/DX = 0 | ! A24 A(7,11,14) 119.6993 -DE/DX = 0 ! A25 A(9,11,14) 119.7377 -DE/DX = 0 | ! A25 A(9,11,14) 119.7322 -DE/DX = 0 ! A26 A(1,15,20) 113.7738 -DE/DX = 0 | ! A26 A(1,15,20) 113.7364 -DE/DX = 0 ! A27 A(1,15,26) 104.2248 -DE/DX = 0 | ! A27 A(1,15,26) 104.2642 -DE/DX = 0 ! A28 A(1,15,28) 113.1423 -DE/DX = 0 | ! A28 A(1,15,28) 113.0294 -DE/DX = 0 ! A29 A(20,15,26) 109.7928 -DE/DX = 0 | ! A29 A(20,15,26) 109.8724 -DE/DX = 0 ! A30 A(20,15,28) 104.8651 -DE/DX = 0 | ! A30 A(20,15,28) 104.8486 -DE/DX = 0 ! A31 A(26,15,28) 111.1628 -DE/DX = 0 | ! A31 A(26,15,28) 111.2185 -DE/DX = 0 ! A32 A(1,16,17) 110.7372 -DE/DX = 0 | ! A32 A(1,16,17) 110.7309 -DE/DX = 0 ! A33 A(1,16,18) 109.9466 -DE/DX = 0 | ! A33 A(1,16,18) 109.9596 -DE/DX = 0 ! A34 A(1,16,19) 110.814 -DE/DX = 0 | ! A34 A(1,16,19) 110.8096 -DE/DX = 0 ! A35 A(17,16,18) 107.1434 -DE/DX = 0 | ! A35 A(17,16,18) 107.1384 -DE/DX = 0 ! A36 A(17,16,19) 109.5142 -DE/DX = 0 | ! A36 A(17,16,19) 109.5255 -DE/DX = 0 ! A37 A(18,16,19) 108.5901 -DE/DX = 0 | ! A37 A(18,16,19) 108.5815 -DE/DX = 0 ! A38 A(15,20,21) 102.4969 -DE/DX = 0 | ! A38 A(15,20,21) 102.526 -DE/DX = 0 ! A39 A(15,20,22) 108.5083 -DE/DX = 0 | ! A39 A(15,20,22) 108.4941 -DE/DX = 0 ! A40 A(15,20,27) 119.5588 -DE/DX = 0 | ! A40 A(15,20,27) 119.5797 -DE/DX = 0 ! A41 A(21,20,22) 109.2657 -DE/DX = 0 | ! A41 A(21,20,22) 109.2077 -DE/DX = 0 ! A42 A(21,20,27) 92.8288 -DE/DX = 0 | ! A42 A(21,20,27) 92.84 -DE/DX = 0 ! A43 A(22,20,27) 120.6577 -DE/DX = 0 | ! A43 A(22,20,27) 120.6641 -DE/DX = 0 ! A44 A(3,21,20) 148.4114 -DE/DX = 0 | ! A44 A(3,21,20) 148.3471 -DE/DX = 0 ! A45 A(20,22,23) 110.3174 -DE/DX = 0 | ! A45 A(20,22,23) 110.3201 -DE/DX = 0 ! A46 A(20,22,24) 124.0298 -DE/DX = 0 | ! A46 A(20,22,24) 124.0322 -DE/DX = 0 ! A47 A(23,22,24) 125.6426 -DE/DX = 0 | ! A47 A(23,22,24) 125.6366 -DE/DX = 0 ! A48 A(22,23,25) 124.9266 -DE/DX = 0 | ! A48 A(22,23,25) 124.9171 -DE/DX = 0 ! A49 A(22,23,28) 112.4542 -DE/DX = 0 | ! A49 A(22,23,28) 112.4516 -DE/DX = 0 ! A50 A(25,23,28) 122.4935 -DE/DX = 0 | ! A50 A(25,23,28) 122.4959 -DE/DX = 0 ! A51 A(15,28,23) 103.7045 -DE/DX = 0 | ! A51 A(15,28,23) 103.6839 -DE/DX = 0 ! A52 A(15,28,29) 113.7512 -DE/DX = 0 | ! A52 A(15,28,29) 113.8345 -DE/DX = 0 ! A53 A(15,28,30) 110.9382 -DE/DX = 0 | ! A53 A(15,28,30) 110.8388 -DE/DX = 0 ! A54 A(23,28,29) 112.3816 -DE/DX = 0 | ! A54 A(23,28,29) 112.5099 -DE/DX = 0 ! A55 A(23,28,30) 110.1396 -DE/DX = 0 | ! A55 A(23,28,30) 110.0476 -DE/DX = 0 ! A56 A(29,28,30) 106.0115 -DE/DX = 0 | ! A56 A(29,28,30) 106.0038 -DE/DX = 0 ! D1 D(15,1,2,3) -12.6809 -DE/DX = 0 | ! D1 D(15,1,2,3) -12.6907 -DE/DX = 0 ! D2 D(16,1,2,3) 167.5453 -DE/DX = 0 | ! D2 D(16,1,2,3) 167.5421 -DE/DX = 0 ! D3 D(2,1,15,20) 15.6449 -DE/DX = 0 | ! D3 D(2,1,15,20) 15.2889 -DE/DX = 0 ! D4 D(2,1,15,26) 135.2085 -DE/DX = 0 | ! D4 D(2,1,15,26) 134.955 -DE/DX = 0 ! D5 D(2,1,15,28) -103.9023 -DE/DX = 0 | ! D5 D(2,1,15,28) -104.1247 -DE/DX = 0 ! D6 D(16,1,15,20) -164.5775 -DE/DX = 0 | ! D6 D(16,1,15,20) -164.9398 -DE/DX = 0 ! D7 D(16,1,15,26) -45.0138 -DE/DX = 0 | ! D7 D(16,1,15,26) -45.2737 -DE/DX = 0 ! D8 D(16,1,15,28) 75.8753 -DE/DX = 0 | ! D8 D(16,1,15,28) 75.6466 -DE/DX = 0 ! D9 D(2,1,16,17) 129.4266 -DE/DX = 0 | ! D9 D(2,1,16,17) 129.646 -DE/DX = 0 ! D10 D(2,1,16,18) -112.3636 -DE/DX = 0 | ! D10 D(2,1,16,18) -112.1461 -DE/DX = 0 ! D11 D(2,1,16,19) 7.688 -DE/DX = 0 | ! D11 D(2,1,16,19) 7.9003 -DE/DX = 0 ! D12 D(15,1,16,17) -50.3502 -DE/DX = 0 | ! D12 D(15,1,16,17) -50.1244 -DE/DX = 0 ! D13 D(15,1,16,18) 67.8596 -DE/DX = 0 | ! D13 D(15,1,16,18) 68.0835 -DE/DX = 0 ! D14 D(15,1,16,19) -172.0888 -DE/DX = 0 | ! D14 D(15,1,16,19) -171.8701 -DE/DX = 0 ! D15 D(1,2,3,4) 125.9498 -DE/DX = 0 | ! D15 D(1,2,3,4) 125.6643 -DE/DX = 0 ! D16 D(1,2,3,21) 0.259 -DE/DX = 0 | ! D16 D(1,2,3,21) 0.1543 -DE/DX = 0 ! D17 D(2,3,4,5) -9.8323 -DE/DX = 0 | ! D17 D(2,3,4,5) -10.0869 -DE/DX = 0 ! D18 D(2,3,4,6) 169.2711 -DE/DX = 0 | ! D18 D(2,3,4,6) 169.1973 -DE/DX = 0 ! D19 D(21,3,4,5) 113.6476 -DE/DX = 0 | ! D19 D(21,3,4,5) 113.2475 -DE/DX = 0 ! D20 D(21,3,4,6) -67.249 -DE/DX = 0 | ! D20 D(21,3,4,6) -67.4683 -DE/DX = 0 ! D21 D(2,3,21,20) 14.9888 -DE/DX = 0 | ! D21 D(2,3,21,20) 16.0814 -DE/DX = 0 ! D22 D(4,3,21,20) -115.0031 -DE/DX = 0 | ! D22 D(4,3,21,20) -113.8165 -DE/DX = 0 ! D23 D(3,4,5,7) 178.7041 -DE/DX = 0 | ! D23 D(3,4,5,7) 178.9154 -DE/DX = 0 ! D24 D(3,4,5,8) 1.0255 -DE/DX = 0 | ! D24 D(3,4,5,8) 1.1995 -DE/DX = 0 ! D25 D(6,4,5,7) -0.3576 -DE/DX = 0 | ! D25 D(6,4,5,7) -0.3353 -DE/DX = 0 ! D26 D(6,4,5,8) -178.0361 -DE/DX = 0 | ! D26 D(6,4,5,8) -178.0512 -DE/DX = 0 ! D27 D(3,4,6,9) -179.6622 -DE/DX = 0 | ! D27 D(3,4,6,9) -179.8493 -DE/DX = 0 ! D28 D(3,4,6,10) -0.5322 -DE/DX = 0 | ! D28 D(3,4,6,10) -0.6914 -DE/DX = 0 ! D29 D(5,4,6,9) -0.5357 -DE/DX = 0 | ! D29 D(5,4,6,9) -0.5467 -DE/DX = 0 ! D30 D(5,4,6,10) 178.5942 -DE/DX = 0 | ! D30 D(5,4,6,10) 178.6112 -DE/DX = 0 ! D31 D(4,5,7,11) 0.9491 -DE/DX = 0 | ! D31 D(4,5,7,11) 0.936 -DE/DX = 0 ! D32 D(4,5,7,12) -178.7087 -DE/DX = 0 | ! D32 D(4,5,7,12) -178.7244 -DE/DX = 0 ! D33 D(8,5,7,11) 178.6103 -DE/DX = 0 | ! D33 D(8,5,7,11) 178.6345 -DE/DX = 0 ! D34 D(8,5,7,12) -1.0475 -DE/DX = 0 | ! D34 D(8,5,7,12) -1.026 -DE/DX = 0 ! D35 D(4,6,9,11) 0.8307 -DE/DX = 0 | ! D35 D(4,6,9,11) 0.8212 -DE/DX = 0 ! D36 D(4,6,9,13) 179.7862 -DE/DX = 0 | ! D36 D(4,6,9,13) 179.7923 -DE/DX = 0 ! D37 D(10,6,9,11) -178.284 -DE/DX = 0 | ! D37 D(10,6,9,11) -178.3219 -DE/DX = 0 ! D38 D(10,6,9,13) 0.6715 -DE/DX = 0 | ! D38 D(10,6,9,13) 0.6492 -DE/DX = 0 ! D39 D(5,7,11,9) -0.6617 -DE/DX = 0 | ! D39 D(5,7,11,9) -0.6689 -DE/DX = 0 ! D40 D(5,7,11,14) -179.9898 -DE/DX = 0 | ! D40 D(5,7,11,14) 179.9985 -DE/DX = 0 ! D41 D(12,7,11,9) 178.9954 -DE/DX = 0 | ! D41 D(12,7,11,9) 178.9909 -DE/DX = 0 ! D42 D(12,7,11,14) -0.3326 -DE/DX = 0 | ! D42 D(12,7,11,14) -0.3417 -DE/DX = 0 ! D43 D(6,9,11,7) -0.245 -DE/DX = 0 | ! D43 D(6,9,11,7) -0.2264 -DE/DX = 0 ! D44 D(6,9,11,14) 179.0827 -DE/DX = 0 | ! D44 D(6,9,11,14) 179.106 -DE/DX = 0 ! D45 D(13,9,11,7) -179.1984 -DE/DX = 0 | ! D45 D(13,9,11,7) -179.1954 -DE/DX = 0 ! D46 D(13,9,11,14) 0.1293 -DE/DX = 0 | ! D46 D(13,9,11,14) 0.137 -DE/DX = 0 ! D47 D(1,15,20,21) -6.3146 -DE/DX = 0 | ! D47 D(1,15,20,21) -5.6165 -DE/DX = 0 ! D48 D(1,15,20,22) -121.8193 -DE/DX = 0 | ! D48 D(1,15,20,22) -121.0638 -DE/DX = 0 ! D49 D(1,15,20,27) 94.2772 -DE/DX = 0 | ! D49 D(1,15,20,27) 95.0135 -DE/DX = 0 ! D50 D(26,15,20,21) -122.6676 -DE/DX = 0 | ! D50 D(26,15,20,21) -122.0492 -DE/DX = 0 ! D51 D(26,15,20,22) 121.8277 -DE/DX = 0 | ! D51 D(26,15,20,22) 122.5036 -DE/DX = 0 ! D52 D(26,15,20,27) -22.0758 -DE/DX = 0 | ! D52 D(26,15,20,27) -21.4191 -DE/DX = 0 ! D53 D(28,15,20,21) 117.8278 -DE/DX = 0 | ! D53 D(28,15,20,21) 118.3494 -DE/DX = 0 ! D54 D(28,15,20,22) 2.3231 -DE/DX = 0 | ! D54 D(28,15,20,22) 2.9021 -DE/DX = 0 ! D55 D(28,15,20,27) -141.5804 -DE/DX = 0 | ! D55 D(28,15,20,27) -141.0206 -DE/DX = 0 ! D56 D(1,15,28,23) 120.9703 -DE/DX = 0 | ! D56 D(1,15,28,23) 120.1999 -DE/DX = 0 ! D57 D(1,15,28,29) -1.4079 -DE/DX = 0 | ! D57 D(1,15,28,29) -2.3698 -DE/DX = 0 ! D58 D(1,15,28,30) -120.8161 -DE/DX = 0 | ! D58 D(1,15,28,30) -121.7536 -DE/DX = 0 ! D59 D(20,15,28,23) -3.5732 -DE/DX = 0 | ! D59 D(20,15,28,23) -4.216 -DE/DX = 0 ! D60 D(20,15,28,29) -125.9514 -DE/DX = 0 | ! D60 D(20,15,28,29) -126.7857 -DE/DX = 0 ! D61 D(20,15,28,30) 114.6404 -DE/DX = 0 | ! D61 D(20,15,28,30) 113.8305 -DE/DX = 0 ! D62 D(26,15,28,23) -122.1565 -DE/DX = 0 | ! D62 D(26,15,28,23) -122.9128 -DE/DX = 0 ! D63 D(26,15,28,29) 115.4654 -DE/DX = 0 | ! D63 D(26,15,28,29) 114.5175 -DE/DX = 0 ! D64 D(26,15,28,30) -3.9428 -DE/DX = 0 | ! D64 D(26,15,28,30) -4.8664 -DE/DX = 0 ! D65 D(15,20,21,3) -10.6211 -DE/DX = 0 | ! D65 D(15,20,21,3) -12.0221 -DE/DX = 0 ! D66 D(22,20,21,3) 104.3336 -DE/DX = 0 | ! D66 D(22,20,21,3) 102.9071 -DE/DX = 0 ! D67 D(27,20,21,3) -131.7427 -DE/DX = 0 | ! D67 D(27,20,21,3) -133.1743 -DE/DX = 0 ! D68 D(15,20,22,23) 0.0242 -DE/DX = 0 | ! D68 D(15,20,22,23) -0.2543 -DE/DX = 0 ! D69 D(15,20,22,24) -178.8692 -DE/DX = 0 | ! D69 D(15,20,22,24) -179.1023 -DE/DX = 0 ! D70 D(21,20,22,23) -110.9981 -DE/DX = 0 | ! D70 D(21,20,22,23) -111.2749 -DE/DX = 0 ! D71 D(21,20,22,24) 70.1085 -DE/DX = 0 | ! D71 D(21,20,22,24) 69.8771 -DE/DX = 0 ! D72 D(27,20,22,23) 143.4591 -DE/DX = 0 | ! D72 D(27,20,22,23) 143.2061 -DE/DX = 0 ! D73 D(27,20,22,24) -35.4343 -DE/DX = 0 | ! D73 D(27,20,22,24) -35.642 -DE/DX = 0 ! D74 D(20,22,23,25) -178.5654 -DE/DX = 0 | ! D74 D(20,22,23,25) -178.5846 -DE/DX = 0 ! D75 D(20,22,23,28) -2.5689 -DE/DX = 0 | ! D75 D(20,22,23,28) -2.7398 -DE/DX = 0 ! D76 D(24,22,23,25) 0.306 -DE/DX = 0 | ! D76 D(24,22,23,25) 0.2408 -DE/DX = 0 ! D77 D(24,22,23,28) 176.3026 -DE/DX = 0 | ! D77 D(24,22,23,28) 176.0856 -DE/DX = 0 ! D78 D(22,23,28,15) 3.8884 -DE/DX = 0 | ! D78 D(22,23,28,15) 4.4106 -DE/DX = 0 ! D79 D(22,23,28,29) 127.1705 -DE/DX = 0 | ! D79 D(22,23,28,29) 127.8532 -DE/DX = 0 ! D80 D(22,23,28,30) -114.8766 -DE/DX = 0 | ! D80 D(22,23,28,30) -114.1831 -DE/DX = 0 ! D81 D(25,23,28,15) 179.997 -DE/DX = 0 | ! D81 D(25,23,28,15) -179.6288 -DE/DX = 0 ! D82 D(25,23,28,29) -56.7209 -DE/DX = 0 | ! D82 D(25,23,28,29) -56.1863 -DE/DX = 0 ! D83 D(25,23,28,30) 61.232 -DE/DX = 0 | ! D83 D(25,23,28,30) 61.7774 -DE/DX = 0 ------------------------------------------------------------ ------------ |