source: ThirdParty/mpqc_open/src/bin/mpqc/validate/ref/mbpt_mp2r12_ne2_tz.out

Candidate_v1.7.1 stable v1.7.1
Last change on this file was 860145, checked in by Frederik Heber <heber@…>, 9 years ago

Merge commit '0b990dfaa8c6007a996d030163a25f7f5fc8a7e7' as 'ThirdParty/mpqc_open'

  • Property mode set to 100644
File size: 19.0 KB
Line 
1
2 MPQC: Massively Parallel Quantum Chemistry
3 Version 2.2.0-alpha
4
5 Machine: i686-pc-linux-gnu
6 User: cljanss@aros.ca.sandia.gov
7 Start Time: Thu Aug 7 17:30:18 2003
8
9 Using ProcMessageGrp for message passing (number of nodes = 1).
10 Using PthreadThreadGrp for threading (number of threads = 1).
11 Using ProcMemoryGrp for distributed shared memory.
12 Total number of processors = 1
13
14 Using IntegralV3 by default for molecular integrals evaluation
15
16 Reading file /home/cljanss/src/SC/lib/atominfo.kv.
17 Molecule: setting point group to d2h
18 WARNING: two unbound groups of atoms
19 consider using extra_bonds input
20
21 adding bond between 1 and 2
22
23 IntCoorGen: generated 1 coordinates.
24 Forming optimization coordinates:
25 SymmMolecularCoor::form_variable_coordinates()
26 expected 0 coordinates
27 found 1 variable coordinates
28 found 0 constant coordinates
29 Reading file /home/cljanss/src/SC/lib/basis/cc-pvtz.kv.
30 Reading file /home/cljanss/src/SC/lib/basis/dz_LdunningR.kv.
31
32 CLSCF::init: total charge = 0
33
34 Starting from core Hamiltonian guess
35
36 Using symmetric orthogonalization.
37 n(basis): 6 0 2 2 0 6 2 2
38 Maximum orthogonalization residual = 2.48511
39 Minimum orthogonalization residual = 0.051287
40 docc = [ 3 0 1 1 0 3 1 1 ]
41 nbasis = 20
42
43 CLSCF::init: total charge = 0
44
45 Projecting guess wavefunction into the present basis set
46
47 SCF::compute: energy accuracy = 1.0000000e-06
48
49 integral intermediate storage = 959288 bytes
50 integral cache = 31037352 bytes
51 nuclear repulsion energy = 13.2294312250
52
53 13009 integrals
54 iter 1 energy = -252.7634629761 delta = 2.97136e-01
55 13009 integrals
56 iter 2 energy = -257.0370876104 delta = 1.09021e-01
57 13009 integrals
58 iter 3 energy = -257.0432373325 delta = 1.12627e-02
59 13009 integrals
60 iter 4 energy = -257.0445247364 delta = 4.83949e-03
61 13009 integrals
62 iter 5 energy = -257.0447095061 delta = 1.26965e-03
63 13009 integrals
64 iter 6 energy = -257.0447095249 delta = 1.58571e-05
65 13009 integrals
66 iter 7 energy = -257.0447095260 delta = 3.65216e-06
67
68 HOMO is 3 B1u = -0.843668
69 LUMO is 4 Ag = 1.256348
70
71 total scf energy = -257.0447095260
72
73 Projecting the guess density.
74
75 The number of electrons in the guess density = 20
76 Using symmetric orthogonalization.
77 n(basis): 13 3 7 7 3 13 7 7
78 Maximum orthogonalization residual = 2.21515
79 Minimum orthogonalization residual = 0.0426604
80 The number of electrons in the projected density = 19.9994
81
82 docc = [ 3 0 1 1 0 3 1 1 ]
83 nbasis = 60
84 Reading file /home/cljanss/src/SC/lib/basis/aug-cc-pvtz.kv.
85
86 Molecular formula Ne2
87
88 MPQC options:
89 matrixkit = <ReplSCMatrixKit>
90 filename = mbpt_mp2r12_ne2_tz
91 restart_file = mbpt_mp2r12_ne2_tz.ckpt
92 restart = no
93 checkpoint = no
94 savestate = no
95 do_energy = yes
96 do_gradient = no
97 optimize = no
98 write_pdb = no
99 print_mole = yes
100 print_timings = yes
101
102
103 SCF::compute: energy accuracy = 1.0000000e-08
104
105 integral intermediate storage = 2028280 bytes
106 integral cache = 29942440 bytes
107 nuclear repulsion energy = 13.2294312250
108
109 956892 integrals
110 iter 1 energy = -257.0616414261 delta = 8.49687e-02
111 956892 integrals
112 iter 2 energy = -257.0637140245 delta = 1.84855e-03
113 956892 integrals
114 iter 3 energy = -257.0637239355 delta = 1.24417e-04
115 956892 integrals
116 iter 4 energy = -257.0637246470 delta = 2.23547e-05
117 956892 integrals
118 iter 5 energy = -257.0637247553 delta = 5.95508e-06
119 956892 integrals
120 iter 6 energy = -257.0637247556 delta = 5.35576e-07
121
122 HOMO is 3 B1u = -0.845157
123 LUMO is 4 Ag = 1.092828
124
125 total scf energy = -257.0637247556
126 NOTE: There are degenerate orbitals within an irrep. This will make
127 some diagnostics, such as the largest amplitude, nonunique.
128 NOTE: There are degenerate orbitals within an irrep. This will make
129 some diagnostics, such as the largest amplitude, nonunique.
130
131 Entered SBS A intermediates evaluator
132 nproc = 1
133 Memory available per node: 32000000 Bytes
134 Static memory used per node: 3015568 Bytes
135 Total memory used per node: 10399888 Bytes
136 Memory required for one pass: 10399888 Bytes
137 Minimum memory required: 3938608 Bytes
138 Batch size: 8
139 npass rest nbasis nshell nfuncmax
140 1 0 60 18 7
141 nocc nvir nfzc nfzv
142 10 50 2 0
143 Memory used for integral storage: 2906768 Bytes
144 Size of global distributed array: 5529600 Bytes
145 Will hold transformed integrals in memory
146 Beginning pass 1
147 Begin loop over shells (grt, 1.+2. q.t.)
148 working on shell pair ( 0 0), 0.6% complete
149 working on shell pair ( 5 2), 10.5% complete
150 working on shell pair ( 7 6), 20.5% complete
151 working on shell pair ( 9 6), 30.4% complete
152 working on shell pair ( 11 2), 40.4% complete
153 working on shell pair ( 12 7), 50.3% complete
154 working on shell pair ( 13 11), 60.2% complete
155 working on shell pair ( 14 14), 70.2% complete
156 working on shell pair ( 16 0), 80.1% complete
157 working on shell pair ( 17 0), 90.1% complete
158 working on shell pair ( 17 17), 100.0% complete
159 End of loop over shells
160 Begin third q.t.
161 End of third q.t.
162 Begin fourth q.t.
163 End of fourth q.t.
164
165 Entered ABS A intermediates evaluator
166 nproc = 1
167 Memory available per node: 32000000 Bytes
168 Static memory used per node: 3094800 Bytes
169 Total memory used per node: 5864960 Bytes
170 Memory required for one pass: 5864960 Bytes
171 Minimum memory required: 3448000 Bytes
172 Batch size: 8
173 npass rest nbasis nshell nfuncmax nbasis(ABS) nshell(ABS) nfuncmax(ABS)
174 1 0 60 18 7 92 26 7
175 nocc nvir nfzc nfzv
176 10 50 2 0
177 Using canonical orthogonalization.
178 n(basis): 19 5 11 11 5 19 11 11
179 Maximum orthogonalization residual = 2.67573
180 Minimum orthogonalization residual = 0.0271252
181 Memory used for integral storage: 2918288 Bytes
182 Size of global distributed array: 1884160 Bytes
183 Will hold transformed integrals in memory
184 Beginning pass 1
185 Begin loop over shells (grt, 1.+2. q.t.)
186 working on shell pair ( 0 0), 0.2% complete
187 working on shell pair ( 1 20), 10.0% complete
188 working on shell pair ( 3 14), 19.9% complete
189 working on shell pair ( 5 8), 29.7% complete
190 working on shell pair ( 7 2), 39.5% complete
191 working on shell pair ( 8 22), 49.4% complete
192 working on shell pair ( 10 16), 59.2% complete
193 working on shell pair ( 12 10), 69.0% complete
194 working on shell pair ( 14 4), 78.8% complete
195 working on shell pair ( 15 24), 88.7% complete
196 working on shell pair ( 17 18), 98.5% complete
197 End of loop over shells
198 Begin fourth q.t.
199 End of fourth q.t.
200
201 Basis Set completeness diagnostics:
202 -Tr(V)/Tr(B) for alpha-alpha pairs: 0.048017
203 -Tr(V)/Tr(B) for alpha-beta pairs: 0.109657
204
205 Alpha-alpha MBPT2-R12/A pair energies:
206 i j mp2(ij) r12(ij) mp2-r12(ij)
207 ----- ----- ------------ ------------ ------------
208 2 1 -0.000000025 -0.000000046 -0.000000071
209 3 1 -0.002521679 -0.000317732 -0.002839412
210 3 2 -0.002521878 -0.000317755 -0.002839632
211 4 1 -0.002523304 -0.000319910 -0.002843214
212 4 2 -0.002523437 -0.000319796 -0.002843233
213 4 3 -0.009022240 -0.000312787 -0.009335026
214 5 1 -0.002523304 -0.000319910 -0.002843214
215 5 2 -0.002523437 -0.000319796 -0.002843233
216 5 3 -0.009022240 -0.000312787 -0.009335026
217 5 4 -0.009027521 -0.000312618 -0.009340139
218 6 1 -0.002523536 -0.000319182 -0.002842719
219 6 2 -0.002523661 -0.000319043 -0.002842703
220 6 3 -0.009023040 -0.000311690 -0.009334730
221 6 4 -0.000000483 -0.000001581 -0.000002063
222 6 5 -0.009028282 -0.000311733 -0.009340015
223 7 1 -0.002523536 -0.000319182 -0.002842719
224 7 2 -0.002523661 -0.000319043 -0.002842703
225 7 3 -0.009023040 -0.000311690 -0.009334730
226 7 4 -0.009028282 -0.000311733 -0.009340015
227 7 5 -0.000000483 -0.000001581 -0.000002063
228 7 6 -0.009029013 -0.000310867 -0.009339880
229 8 1 -0.002525223 -0.000316973 -0.002842196
230 8 2 -0.002525318 -0.000316667 -0.002841985
231 8 3 -0.000001269 -0.000006401 -0.000007670
232 8 4 -0.009033817 -0.000310655 -0.009344472
233 8 5 -0.009033817 -0.000310655 -0.009344472
234 8 6 -0.009034547 -0.000309773 -0.009344321
235 8 7 -0.009034547 -0.000309773 -0.009344321
236
237 Alpha-beta MBPT2-R12/A pair energies:
238 i j mp2(ij) r12(ij) mp2-r12(ij)
239 ----- ----- ------------ ------------ ------------
240 1 1 -0.004782085 -0.000832435 -0.005614520
241 1 2 -0.004782175 -0.000832351 -0.005614527
242 1 3 -0.004198950 -0.000969198 -0.005168148
243 1 4 -0.004201346 -0.000968661 -0.005170007
244 1 5 -0.004201346 -0.000968661 -0.005170007
245 1 6 -0.004201682 -0.000968079 -0.005169761
246 1 7 -0.004201682 -0.000968079 -0.005169761
247 1 8 -0.004204120 -0.000969099 -0.005173218
248 2 1 -0.004782175 -0.000832351 -0.005614527
249 2 2 -0.004782272 -0.000832365 -0.005614637
250 2 3 -0.004199078 -0.000969067 -0.005168145
251 2 4 -0.004201477 -0.000968598 -0.005170075
252 2 5 -0.004201477 -0.000968598 -0.005170075
253 2 6 -0.004201813 -0.000968017 -0.005169831
254 2 7 -0.004201813 -0.000968017 -0.005169831
255 2 8 -0.004204251 -0.000968971 -0.005173222
256 3 1 -0.004198950 -0.000969198 -0.005168148
257 3 2 -0.004199078 -0.000969067 -0.005168145
258 3 3 -0.011069281 -0.001266825 -0.012336106
259 3 4 -0.005708562 -0.000635445 -0.006344008
260 3 5 -0.005708562 -0.000635445 -0.006344008
261 3 6 -0.005709013 -0.000634715 -0.006343728
262 3 7 -0.005709013 -0.000634715 -0.006343728
263 3 8 -0.011080973 -0.001264046 -0.012345018
264 4 1 -0.004201346 -0.000968661 -0.005170007
265 4 2 -0.004201477 -0.000968598 -0.005170075
266 4 3 -0.005708562 -0.000635445 -0.006344008
267 4 4 -0.011080036 -0.001264410 -0.012344446
268 4 5 -0.005711858 -0.000634587 -0.006346445
269 4 6 -0.011080811 -0.001263589 -0.012344400
270 4 7 -0.005712325 -0.000633957 -0.006346282
271 4 8 -0.005715771 -0.000633930 -0.006349701
272 5 1 -0.004201346 -0.000968661 -0.005170007
273 5 2 -0.004201477 -0.000968598 -0.005170075
274 5 3 -0.005708562 -0.000635445 -0.006344008
275 5 4 -0.005711858 -0.000634587 -0.006346445
276 5 5 -0.011080036 -0.001264410 -0.012344446
277 5 6 -0.005712325 -0.000633957 -0.006346282
278 5 7 -0.011080811 -0.001263589 -0.012344400
279 5 8 -0.005715771 -0.000633930 -0.006349701
280 6 1 -0.004201682 -0.000968079 -0.005169761
281 6 2 -0.004201813 -0.000968017 -0.005169831
282 6 3 -0.005709013 -0.000634715 -0.006343728
283 6 4 -0.011080811 -0.001263589 -0.012344400
284 6 5 -0.005712325 -0.000633957 -0.006346282
285 6 6 -0.011081589 -0.001262782 -0.012344371
286 6 7 -0.005712792 -0.000633345 -0.006346137
287 6 8 -0.005716230 -0.000633203 -0.006349434
288 7 1 -0.004201682 -0.000968079 -0.005169761
289 7 2 -0.004201813 -0.000968017 -0.005169831
290 7 3 -0.005709013 -0.000634715 -0.006343728
291 7 4 -0.005712325 -0.000633957 -0.006346282
292 7 5 -0.011080811 -0.001263589 -0.012344400
293 7 6 -0.005712792 -0.000633345 -0.006346137
294 7 7 -0.011081589 -0.001262782 -0.012344371
295 7 8 -0.005716230 -0.000633203 -0.006349434
296 8 1 -0.004204120 -0.000969099 -0.005173218
297 8 2 -0.004204251 -0.000968971 -0.005173222
298 8 3 -0.011080973 -0.001264046 -0.012345018
299 8 4 -0.005715771 -0.000633930 -0.006349701
300 8 5 -0.005715771 -0.000633930 -0.006349701
301 8 6 -0.005716230 -0.000633203 -0.006349434
302 8 7 -0.005716230 -0.000633203 -0.006349434
303 8 8 -0.011092865 -0.001261267 -0.012354131
304
305 RHF energy [au]: -257.063724755551
306 MP2 correlation energy [au]: -0.528658888169
307 (MBPT2)-R12/ A correlation energy [au]: -0.064532741399
308 MBPT2-R12/ A correlation energy [au]: -0.593191629568
309 MBPT2-R12/ A energy [au]: -257.656916385119
310
311Value of the MolecularEnergy: -257.6569163851
312
313 MBPT2_R12:
314 Standard Approximation: A
315 Spin-adapted algorithm: false
316 How to Store Transformed Integrals: mem_posix
317
318 Transformed Integrals file: mbpt_mp2r12_ne2_tz./r12ints.dat
319
320 Auxiliary Basis:
321 GaussianBasisSet:
322 nbasis = 92
323 nshell = 26
324 nprim = 44
325 name = "aug-cc-pVTZ"
326
327 MBPT2:
328 Function Parameters:
329 value_accuracy = 3.945593e-07 (1.000000e-06) (computed)
330 gradient_accuracy = 0.000000e+00 (1.000000e-06)
331 hessian_accuracy = 0.000000e+00 (1.000000e-04)
332
333 Molecular Coordinates:
334 IntMolecularCoor Parameters:
335 update_bmat = no
336 scale_bonds = 1
337 scale_bends = 1
338 scale_tors = 1
339 scale_outs = 1
340 symmetry_tolerance = 1.000000e-05
341 simple_tolerance = 1.000000e-03
342 coordinate_tolerance = 1.000000e-07
343 have_fixed_values = 0
344 max_update_steps = 100
345 max_update_disp = 0.500000
346 have_fixed_values = 0
347
348 Molecular formula: Ne2
349 molecule<Molecule>: (
350 symmetry = d2h
351 unit = "angstrom"
352 { n atoms geometry }={
353 1 Ne [ 0.0000000000 0.0000000000 2.0000000000]
354 2 Ne [ 0.0000000000 0.0000000000 -2.0000000000]
355 }
356 )
357 Atomic Masses:
358 19.99244 19.99244
359
360 Bonds:
361 STRE s1 4.00000 1 2 Ne-Ne
362
363 SymmMolecularCoor Parameters:
364 change_coordinates = no
365 transform_hessian = yes
366 max_kappa2 = 10.000000
367
368 GaussianBasisSet:
369 nbasis = 60
370 nshell = 18
371 nprim = 36
372 name = "cc-pVTZ"
373 Reference Wavefunction:
374 Function Parameters:
375 value_accuracy = 3.945593e-09 (1.000000e-08) (computed)
376 gradient_accuracy = 0.000000e+00 (1.000000e-06)
377 hessian_accuracy = 0.000000e+00 (1.000000e-04)
378
379 Molecule:
380 Molecular formula: Ne2
381 molecule<Molecule>: (
382 symmetry = d2h
383 unit = "angstrom"
384 { n atoms geometry }={
385 1 Ne [ 0.0000000000 0.0000000000 2.0000000000]
386 2 Ne [ 0.0000000000 0.0000000000 -2.0000000000]
387 }
388 )
389 Atomic Masses:
390 19.99244 19.99244
391
392 GaussianBasisSet:
393 nbasis = 60
394 nshell = 18
395 nprim = 36
396 name = "cc-pVTZ"
397 SCF Parameters:
398 maxiter = 40
399 density_reset_frequency = 10
400 savestate_iter = 0
401 savestate_frequency = 1
402 level_shift = 0.000000
403
404 CLSCF Parameters:
405 charge = 0
406 ndocc = 10
407 docc = [ 3 0 1 1 0 3 1 1 ]
408
409
410 The following keywords in "mbpt_mp2r12_ne2_tz.in" were ignored:
411 mpqc:mole:reference:guess_wavefunction:multiplicity
412 mpqc:mole:reference:multiplicity
413
414 CPU Wall
415mpqc: 36.06 73.63
416 calc: 35.53 72.63
417 mp2-r12/a energy: 35.53 72.63
418 mp2-r12/a pair energies: 0.02 0.02
419 r12a-abs-mem: 22.00 45.14
420 mp2-r12/a passes: 21.86 45.00
421 4. q.t.: 0.14 0.29
422 MO ints store: 0.00 0.00
423 grt+1.qt+2.qt: 21.72 44.71
424 mp2-r12a intermeds: 0.06 0.06
425 MO ints contraction: 0.05 0.05
426 MO ints retrieve: 0.00 0.00
427 r12a-sbs-mem: 9.06 18.36
428 mp2-r12/a passes: 8.75 17.75
429 3. q.t.: 0.10 0.26
430 4. q.t.: 0.11 0.27
431 MO ints store: 0.00 0.00
432 compute emp2: 0.01 0.01
433 grt+1.qt+2.qt: 8.50 17.19
434 mp2-r12a intermeds: 0.30 0.60
435 MO ints contraction: 0.29 0.59
436 MO ints retrieve: 0.00 0.00
437 vector: 4.44 9.10
438 density: 0.00 0.00
439 evals: 0.01 0.00
440 extrap: 0.01 0.01
441 fock: 4.42 9.08
442 accum: 0.00 0.00
443 ao_gmat: 4.27 8.77
444 start thread: 4.26 8.77
445 stop thread: 0.00 0.00
446 init pmax: 0.00 0.00
447 local data: 0.00 0.00
448 setup: 0.07 0.07
449 sum: 0.00 0.00
450 symm: 0.08 0.23
451 input: 0.53 1.00
452 vector: 0.37 0.84
453 density: 0.00 0.00
454 evals: 0.01 0.00
455 extrap: 0.00 0.01
456 fock: 0.36 0.82
457 accum: 0.00 0.00
458 ao_gmat: 0.33 0.79
459 start thread: 0.33 0.79
460 stop thread: 0.00 0.00
461 init pmax: 0.00 0.00
462 local data: 0.00 0.00
463 setup: 0.01 0.01
464 sum: 0.00 0.00
465 symm: 0.02 0.02
466
467 End Time: Thu Aug 7 17:31:32 2003
468
Note: See TracBrowser for help on using the repository browser.