File: h2o.efield.out

package info (click to toggle)
espresso 5.1%2Bdfsg-3
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 146,004 kB
  • ctags: 17,245
  • sloc: f90: 253,041; sh: 51,271; ansic: 27,494; tcl: 15,570; xml: 14,508; makefile: 2,958; perl: 2,035; fortran: 1,924; python: 337; cpp: 200; awk: 57
file content (714 lines) | stat: -rw-r--r-- 41,473 bytes parent folder | download | duplicates (6)
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
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
=------------------------------------------------------------------------------=

     CP: variable-cell Car-Parrinello molecular dynamics
     using norm-conserving and ultrasoft Vanderbilt pseudopotentials

     Version: 4.0 - Mon Apr 28 15:32:33 CEST 2008                      
     Authors: Alfredo Pasquarello, Kari Laasonen, Andrea Trave, Roberto Car,
       Paolo Giannozzi, Nicola Marzari, Carlo Cavazzoni, Guido Chiarotti,
       Sandro Scandolo, Paolo Focher, Gerardo Ballabio, and others

=------------------------------------------------------------------------------=

   This run was started on:  16:11:11  28Apr2008               

   Serial Build

   Job Title: MD Simulation


   Atomic Pseudopotentials Parameters
   ----------------------------------

   Reading pseudopotential for specie #  1 from file :
   /home/giannozz/espresso/pseudo/O_US.van
   file type is  1: Vanderbilt US PP

    ============================================================
    |  pseudopotential report for atomic species:  1           |
    |        pseudo potential version   7   3   2              |
    ------------------------------------------------------------
    |  oxygen              PBE                  exchange-corr  |
    |  z =   8.    zv( 1) =   6.    exfact =   5.00000         |
    |  ifpcor =  0           atomic energy = -31.58351 Ry      |
    |  index    orbital      occupation    energy              |
    |    1        200           2.00       -1.76               |
    |    2        210           4.00       -0.66               |
    |  rinner =  0.7000  0.7000  0.7000                        |
    |    new generation scheme:                                |
    |    nbeta =  4     kkbeta =  519     rcloc =    1.0000    |
    |    ibeta    l     epsilon   rcut                         |
    |      1       0      -1.76   1.20                         |
    |      2       0      -0.66   1.20                         |
    |      3       1      -1.76   1.20                         |
    |      4       1      -0.66   1.20                         |
    ============================================================

   Reading pseudopotential for specie #  2 from file :
   /home/giannozz/espresso/pseudo/H_US.van
   file type is  1: Vanderbilt US PP

    ============================================================
    |  pseudopotential report for atomic species:  2           |
    |        pseudo potential version   7   3   2              |
    ------------------------------------------------------------
    |  hydrogen            PBE                  exchange-corr  |
    |  z =   1.    zv( 2) =   1.    exfact =   5.00000         |
    |  ifpcor =  0           atomic energy =  -0.91772 Ry      |
    |  index    orbital      occupation    energy              |
    |    1        100           1.00       -0.48               |
    |  rinner =  0.5000                                        |
    |    new generation scheme:                                |
    |    nbeta =  1     kkbeta =  271     rcloc =    0.6000    |
    |    ibeta    l     epsilon   rcut                         |
    |      1       0      -0.48   0.80                         |
    ============================================================


   Main Simulation Parameters (from input)
   ---------------------------------------
   Restart Mode       =       0   reset_counters 
   Number of MD Steps =     200
   Print out every          100 MD Steps
   Reads from unit    =      50
   Writes to unit     =      51
   MD Simulation time step            =       6.00
   Electronic fictitious mass (emass) =     350.00
   emass cut-off                      =       3.00

   Simulation Cell Parameters (from input)
   external pressure       =            0.00 [GPa]
   wmass (calculated)      =         2770.46 [AU]
   ibrav =    1
   alat  =    20.00000000
   a1    =    20.00000000    0.00000000    0.00000000
   a2    =     0.00000000   20.00000000    0.00000000
   a3    =     0.00000000    0.00000000   20.00000000
 
   b1    =     0.05000000    0.00000000    0.00000000
   b2    =     0.00000000    0.05000000    0.00000000
   b3    =     0.00000000    0.00000000    0.05000000
   omega =  8000.00000000

   Energy Cut-offs
   ---------------
   Ecutwfc =   25.0 Ry,      Ecutrho =  100.0 Ry,      Ecuts =  100.0 Ry
   Gcutwfc =   15.9     ,    Gcutrho =   31.8          Gcuts =   31.8
   NOTA BENE: refg, mmx =   0.050000  2400
   Eigenvalues calculated without the kinetic term contribution
   Orthog. with lagrange multipliers : eps =   0.50E-07,  max = 250
   Electron dynamics with newton equations
   Electron dynamics : the temperature is not controlled

   Electronic states
   -----------------
   Number of Electron =     8, of States =     4
   Occupation numbers :
   2.00 2.00 2.00 2.00


   Exchange and correlations functionals
   -------------------------------------
   Using Local Density Approximation with
     Exchange functional: SLATER                                                      
     Correlation functional: PERDEW AND WANG                                             
   Using Generalized Gradient Corrections with
     Exchange functional: PERDEW BURKE ERNZERHOF                                      
     Correlation functional: PERDEW BURKE ERNZERHOF                                      
     Exchange-correlation      = PBE (1434)


   Ions Simulation Parameters
   --------------------------
   Ions are allowed to move
   Ions dynamics with newton equations
   the temperature is computed for     9 degrees of freedom
   ion dynamics with fricp =  0.0000 and greasp =  1.0000
   Ionic position (from input)
   sorted by specie, and converted to real a.u. coordinates
   Species   1 atoms =    1 mass =     29166.22 (a.u.),        16.00 (amu) rcmax =   0.50 (a.u.)
       10.000000    10.000000    10.000000
   Species   2 atoms =    2 mass =      3645.78 (a.u.),         2.00 (amu) rcmax =   0.50 (a.u.)
       11.732500     9.675700    10.000000
        9.675700    11.732500    10.000000
   Ionic position will be re-read from restart file

   All atoms are allowed to move
   Ionic temperature is not controlled


   Cell Dynamics Parameters (from STDIN)
   -------------------------------------
   Starting cell generated from CELLDM
   Cell parameters will be re-read from restart file
   Constant VOLUME Molecular dynamics
   cell parameters are not allowed to move

   Verbosity: iprsta =  1



   Simulation dimensions initialization
   ------------------------------------

   unit vectors of full simulation cell
   in real space:                         in reciprocal space (units 2pi/alat):
   1    20.0000    0.0000    0.0000              1.0000    0.0000    0.0000
   2     0.0000   20.0000    0.0000              0.0000    1.0000    0.0000
   3     0.0000    0.0000   20.0000              0.0000    0.0000    1.0000
 
   Stick Mesh
   ----------
   nst =  1591,  nstw =   397, nsts =  1591
        PEs    n.st   n.stw   n.sts    n.g    n.gw   n.gs
          1    3181     793    3181  135043   16879  135043
          0    3181     793    3181  135043   16879  135043
 
 
   Real Mesh
   ---------
   Global Dimensions   Local  Dimensions   Processor Grid
   .X.   .Y.   .Z.     .X.   .Y.   .Z.     .X.   .Y.   .Z.
    64    64    64      64    64    64       1     1     1
   Array leading dimensions ( nr1x, nr2x, nr3x )   =     64    64    64
   Local number of cell to store the grid ( nnrx ) =     262144
   Number of x-y planes for each processors: 
   nr3l =    64
 
   Smooth Real Mesh
   ----------------
   Global Dimensions   Local  Dimensions   Processor Grid
   .X.   .Y.   .Z.     .X.   .Y.   .Z.     .X.   .Y.   .Z.
    64    64    64      64    64    64       1     1     1
   Array leading dimensions ( nr1x, nr2x, nr3x )   =     64    64    64
   Local number of cell to store the grid ( nnrx ) =     262144
   Number of x-y planes for each processors: 
   nr3sl =    64
 
   Small Box Real Mesh
   -------------------
   Global Dimensions   Local  Dimensions   Processor Grid
   .X.   .Y.   .Z.     .X.   .Y.   .Z.     .X.   .Y.   .Z.
    10    10    10      10    10    10       1     1     1
   Array leading dimensions ( nr1x, nr2x, nr3x )   =     10    10    10
   Local number of cell to store the grid ( nnrx ) =       1000
 
   unit vectors of box grid cell
   in real space:                         in reciprocal space:
       3.1250    0.0000    0.0000              1.0000    0.0000    0.0000
       0.0000    3.1250    0.0000              0.0000    1.0000    0.0000
       0.0000    0.0000    3.1250              0.0000    0.0000    1.0000
 
   Reciprocal Space Mesh
   ---------------------
                Large Mesh
   PE   Global(ngmt)     Local(ngm) MaxLocal(ngmx)
    1          67522          67522          67522
                Smooth Mesh
   PE   Global(ngst)     Local(ngs) MaxLocal(ngsx)
    1          67522          67522          67522
                Wave function Mesh
   PE   Global(ngwt)     Local(ngw) MaxLocal(ngwx)
    1           8440           8440           8440

                Small box Mesh
   ngb =          243 not distributed to processors


   System geometry initialization
   ------------------------------
 
   Scaled positions from standard input
   O    0.500000E+00  0.500000E+00  0.500000E+00
   H    0.586625E+00  0.483785E+00  0.500000E+00
   H    0.483785E+00  0.586625E+00  0.500000E+00
   ibrav =    1       cell parameters 

      20.00000       0.00000       0.00000
       0.00000      20.00000       0.00000
       0.00000       0.00000      20.00000
 

     !----------------------------------!
     !                                  !
     ! ADIABATIC SWITCHING OF THE FIELD !
     !                                  !
     !----------------------------------!



   Pseudopotentials initialization
   -------------------------------
  nlinit  nh(is), ngb, is, kkbeta, lmaxq =            8         243           1
         519           3
 
                        qqq 
  -0.2814  -0.1752   0.0000   0.0000
  -0.1752  -0.1519   0.0000   0.0000
   0.0000   0.0000   3.4543  -1.8676
   0.0000   0.0000  -1.8676   0.9140
 
  nlinit  nh(is), ngb, is, kkbeta, lmaxq =            1         243           2
         271           1
 
                        qqq 
   9.4958
 


   Common initialization

   Specie:     1
   1  indv=  1   ang. mom=  0
   2  indv=  2   ang. mom=  0
   3  indv=  3   ang. mom=  1
   4  indv=  3   ang. mom=  1
   5  indv=  3   ang. mom=  1
   6  indv=  4   ang. mom=  1
   7  indv=  4   ang. mom=  1
   8  indv=  4   ang. mom=  1
 
                        dion 
   1.1855  -0.0482   0.0000   0.0000
  -0.0482  -0.1194   0.0000   0.0000
   0.0000   0.0000  24.3757 -18.2972
   0.0000   0.0000 -18.2972  13.6557

   Specie:     2
   1  indv=  1   ang. mom=  0
 
                        dion 
   3.3279
 ------------------------DYNAMICS IN THE WANNIER BASIS--------------------------
                              DYNAMICS PARAMETERS 
 FICTITIOUS MASS PARAMETER USED FOR SPREAD OPTIMIZATION            =  1500.0
 TIME STEP USED FOR DAMPED DYNAMICS                                =  5.0000000
 DAMPING COEFFICIENT USED FOR WANNIER FUNCTION SPREAD OPTIMIZATION =  0.3000000
          20 STEPS OF DAMPED MOLECULAR DYNAMICS FOR OPTIMIZATION OF THE SPREAD
 AVERAGE WANNIER FUNCTION SPREAD WRITTEN TO     FORT.24
 INDIVIDUAL WANNIER FUNCTION SPREAD WRITTEN TO  /h2o_mol.spr
 WANNIER CENTERS WRITTEN TO                     /h2o_mol.wfc
 SOME PERTINENT RUN-TIME INFORMATION WRITTEN TO FORT.27
 -------------------------------------------------------------------------------
 
 IBRAV SELECTED:           1
 WANNIER SETUP : check G vectors and weights
inw = 1:   1   0   0   1.000000
inw = 2:   0   1   0   1.000000
inw = 3:   0   0   1   1.000000
 Translations to be done           3
 ibrav selected:           1
           1 :           1           0           0
           2 :           0           1           0
           3 :           0           0           1
 Translation           1 for        8440 G vectors
 Translation           2 for        8440 G vectors
 Translation           3 for        8440 G vectors
 
 GRID MAPPING DONE
 DYNAMICS IN THE PRESENCE OF AN EXTERNAL ELECTRIC FIELD
 
 POLARIZATION CONTRIBUTION OUTPUT TO FORT.28 IN THE FOLLOWING FORMAT
 
 EFX, EFY, EFZ, ELECTRIC ENTHALPY(ELECTRONIC), ELECTRIC ENTHALPY(IONIC)
 
 E0(x) =  0.0000000
 E0(y) =  0.0000000
 E0(z) =  0.0000000
 E1(x) =  0.0000000
 E1(y) =  0.0000000
 E1(z) =  0.0050000
 FIELD WILL BE TURNED ON ADIBATICALLY OVER   100 STEPS
 
 
   Short Legend and Physical Units in the Output
   ---------------------------------------------
   NFI    [int]          - step index
   EKINC  [HARTREE A.U.] - kinetic energy of the fictitious electronic dynamics
   TEMPH  [K]            - Temperature of the fictitious cell dynamics
   TEMP   [K]            - Ionic temperature
   ETOT   [HARTREE A.U.] - Scf total energy (Kohn-Sham hamiltonian)
   ENTHAL [HARTREE A.U.] - Enthalpy ( ETOT + P * V )
   ECONS  [HARTREE A.U.] - Enthalpy + kinetic energy of ions and cell
   ECONT  [HARTREE A.U.] - Constant of motion for the CP lagrangian
 

   reading restart file: /h2o_mol_50.save
   restart file read in    0.019 sec.


   formf: eself=  30.31961
   formf:     vps(g=0)=  -0.0002325     rhops(g=0)=  -0.0007500
   formf: sum_g vps(g)=  -0.4490836 sum_g rhops(g)=  -4.2852036
   formf:     vps(g=0)=  -0.0000709     rhops(g=0)=  -0.0001250
   formf: sum_g vps(g)=  -0.7398336 sum_g rhops(g)=  -0.7142006
   Delta V(G=0):   0.000785Ry,    0.021372eV
 
  nfi    ekinc  temph  tempp        etot      enthal       econs       econt    vnhh   xnhh0    vnhp   xnhp0
    1  0.00000    0.0    0.7   -17.20048   -17.20048   -17.20047   -17.20047  0.0000  0.0000  0.0000  0.0000
    2  0.00001    0.0    6.0   -17.20016   -17.20016   -17.20008   -17.20007  0.0000  0.0000  0.0000  0.0000
    3  0.00006    0.0   16.0   -17.20033   -17.20033   -17.20010   -17.20004  0.0000  0.0000  0.0000  0.0000
    4  0.00007    0.0   30.2   -17.20055   -17.20055   -17.20012   -17.20006  0.0000  0.0000  0.0000  0.0000
    5  0.00006    0.0   47.6   -17.20082   -17.20082   -17.20014   -17.20008  0.0000  0.0000  0.0000  0.0000
    6  0.00006    0.0   67.7   -17.20113   -17.20116   -17.20019   -17.20014  0.0000  0.0000  0.0000  0.0000
    7  0.00011    0.0   89.9   -17.20148   -17.20148   -17.20020   -17.20009  0.0000  0.0000  0.0000  0.0000
    8  0.00019    0.0  113.9   -17.20188   -17.20188   -17.20025   -17.20007  0.0000  0.0000  0.0000  0.0000
    9  0.00026    0.0  139.6   -17.20229   -17.20229   -17.20030   -17.20005  0.0000  0.0000  0.0000  0.0000
   10  0.00027    0.0  166.7   -17.20271   -17.20235   -17.19998   -17.19971  0.0000  0.0000  0.0000  0.0000
   11  0.00021    0.0  194.7   -17.20313   -17.20258   -17.19980   -17.19959  0.0000  0.0000  0.0000  0.0000
   12  0.00020    0.0  222.8   -17.20355   -17.20355   -17.20038   -17.20017  0.0000  0.0000  0.0000  0.0000
   13  0.00023    0.0  250.1   -17.20397   -17.20397   -17.20040   -17.20018  0.0000  0.0000  0.0000  0.0000
   14  0.00025    0.0  275.3   -17.20435   -17.20435   -17.20043   -17.20018  0.0000  0.0000  0.0000  0.0000
   15  0.00025    0.0  297.7   -17.20469   -17.20469   -17.20045   -17.20019  0.0000  0.0000  0.0000  0.0000
   16  0.00026    0.0  316.4   -17.20497   -17.20198   -17.19747   -17.19721  0.0000  0.0000  0.0000  0.0000
   17  0.00031    0.0  331.2   -17.20522   -17.20132   -17.19660   -17.19629  0.0000  0.0000  0.0000  0.0000
   18  0.00039    0.0  342.1   -17.20543   -17.20043   -17.19555   -17.19516  0.0000  0.0000  0.0000  0.0000
   19  0.00047    0.0  349.5   -17.20558   -17.20558   -17.20060   -17.20013  0.0000  0.0000  0.0000  0.0000
   20  0.00050    0.0  353.8   -17.20566   -17.21349   -17.20845   -17.20795  0.0000  0.0000  0.0000  0.0000
   21  0.00037    0.0  355.3   -17.20568   -17.21529   -17.21023   -17.20985  0.0000  0.0000  0.0000  0.0000
   22  0.00029    0.0  354.1   -17.20565   -17.20565   -17.20061   -17.20032  0.0000  0.0000  0.0000  0.0000
   23  0.00028    0.0  350.3   -17.20559   -17.20559   -17.20060   -17.20032  0.0000  0.0000  0.0000  0.0000
   24  0.00025    0.0  343.6   -17.20547   -17.20547   -17.20057   -17.20032  0.0000  0.0000  0.0000  0.0000
   25  0.00022    0.0  333.8   -17.20530   -17.20530   -17.20054   -17.20032  0.0000  0.0000  0.0000  0.0000
   26  0.00020    0.0  320.7   -17.20508   -17.22794   -17.22336   -17.22317  0.0000  0.0000  0.0000  0.0000
   27  0.00019    0.0  304.5   -17.20484   -17.23132   -17.22698   -17.22679  0.0000  0.0000  0.0000  0.0000
   28  0.00026    0.0  285.6   -17.20458   -17.20458   -17.20051   -17.20025  0.0000  0.0000  0.0000  0.0000
   29  0.00026    0.0  264.8   -17.20429   -17.20430   -17.20052   -17.20026  0.0000  0.0000  0.0000  0.0000
   30  0.00022    0.0  243.4   -17.20262   -17.20262   -17.19915   -17.19893  0.0000  0.0000  0.0000  0.0000
   31  0.00018    0.0  222.1   -17.20230   -17.24692   -17.24375   -17.24357  0.0000  0.0000  0.0000  0.0000
   32  0.00021    0.0  201.4   -17.20200   -17.20200   -17.19913   -17.19892  0.0000  0.0000  0.0000  0.0000
   33  0.00023    0.0  181.6   -17.20170   -17.14579   -17.14320   -17.14297  0.0000  0.0000  0.0000  0.0000
   34  0.00015    0.0  162.9   -17.20141   -17.13930   -17.13698   -17.13682  0.0000  0.0000  0.0000  0.0000
   35  0.00011    0.0  145.1   -17.20112   -17.20113   -17.19906   -17.19895  0.0000  0.0000  0.0000  0.0000
   36  0.00008    0.0  128.4   -17.20084   -17.12523   -17.12340   -17.12332  0.0000  0.0000  0.0000  0.0000
   37  0.00009    0.0  112.6   -17.20059   -17.20061   -17.19900   -17.19891  0.0000  0.0000  0.0000  0.0000
   38  0.00008    0.0   98.0   -17.20037   -17.08727   -17.08588   -17.08580  0.0000  0.0000  0.0000  0.0000
   39  0.00010    0.0   84.8   -17.20018   -17.20020   -17.19900   -17.19890  0.0000  0.0000  0.0000  0.0000
   40  0.00006    0.0   73.5   -17.20001   -17.09339   -17.09235   -17.09229  0.0000  0.0000  0.0000  0.0000
   41  0.00007    0.0   64.3   -17.19987   -17.22867   -17.22776   -17.22769  0.0000  0.0000  0.0000  0.0000
   42  0.00005    0.0   57.6   -17.19976   -17.07591   -17.07509   -17.07503  0.0000  0.0000  0.0000  0.0000
   43  0.00007    0.0   53.5   -17.19971   -17.06688   -17.06612   -17.06605  0.0000  0.0000  0.0000  0.0000
   44  0.00007    0.0   52.0   -17.19969   -17.05767   -17.05693   -17.05686  0.0000  0.0000  0.0000  0.0000
   45  0.00006    0.0   53.1   -17.19970   -17.19974   -17.19898   -17.19892  0.0000  0.0000  0.0000  0.0000
   46  0.00005    0.0   56.5   -17.19972   -17.03880   -17.03800   -17.03795  0.0000  0.0000  0.0000  0.0000
   47  0.00005    0.0   62.2   -17.19980   -17.02924   -17.02836   -17.02830  0.0000  0.0000  0.0000  0.0000
   48  0.00006    0.0   70.1   -17.19992   -17.01962   -17.01862   -17.01856  0.0000  0.0000  0.0000  0.0000
   49  0.00006    0.0   80.2   -17.20007   -17.00997   -17.00883   -17.00876  0.0000  0.0000  0.0000  0.0000
   50  0.00006    0.0   92.5   -17.20024   -16.95031   -16.94899   -16.94893  0.0000  0.0000  0.0000  0.0000
   51  0.00006    0.0  106.8   -17.20044   -17.25298   -17.25146   -17.25139  0.0000  0.0000  0.0000  0.0000
   52  0.00008    0.0  123.0   -17.20067   -16.98109   -16.97934   -16.97926  0.0000  0.0000  0.0000  0.0000
   53  0.00010    0.0  140.8   -17.20094   -16.97164   -16.96963   -16.96954  0.0000  0.0000  0.0000  0.0000
   54  0.00011    0.0  159.9   -17.20123   -16.96229   -16.96001   -16.95989  0.0000  0.0000  0.0000  0.0000
   55  0.00012    0.0  180.0   -17.20152   -16.89092   -16.88835   -16.88823  0.0000  0.0000  0.0000  0.0000
   56  0.00013    0.0  200.7   -17.20183   -17.26643   -17.26357   -17.26344  0.0000  0.0000  0.0000  0.0000
   57  0.00016    0.0  221.9   -17.20214   -17.26903   -17.26587   -17.26571  0.0000  0.0000  0.0000  0.0000
   58  0.00021    0.0  243.2   -17.20246   -16.85742   -16.85395   -16.85374  0.0000  0.0000  0.0000  0.0000
   59  0.00020    0.0  264.5   -17.20278   -16.91804   -16.91428   -16.91407  0.0000  0.0000  0.0000  0.0000
   60  0.00020    0.0  285.6   -17.20308   -16.90986   -16.90579   -16.90559  0.0000  0.0000  0.0000  0.0000
   61  0.00022    0.0  306.0   -17.20338   -16.82655   -16.82219   -16.82196  0.0000  0.0000  0.0000  0.0000
   62  0.00026    0.0  326.0   -17.20522   -16.81850   -16.81386   -16.81360  0.0000  0.0000  0.0000  0.0000
   63  0.00027    0.0  344.8   -17.20552   -17.28496   -17.28005   -17.27978  0.0000  0.0000  0.0000  0.0000
   64  0.00023    0.0  361.5   -17.20575   -16.80043   -16.79527   -16.79505  0.0000  0.0000  0.0000  0.0000
   65  0.00022    0.0  375.4   -17.20594   -16.79199   -16.78664   -16.78642  0.0000  0.0000  0.0000  0.0000
   66  0.00022    0.0  386.2   -17.20609   -16.78393   -16.77842   -16.77820  0.0000  0.0000  0.0000  0.0000
   67  0.00023    0.0  393.6   -17.20619   -16.77628   -16.77067   -16.77044  0.0000  0.0000  0.0000  0.0000
   68  0.00025    0.0  397.5   -17.20626   -16.76906   -16.76340   -16.76315  0.0000  0.0000  0.0000  0.0000
   69  0.00026    0.0  397.8   -17.20628   -16.76225   -16.75658   -16.75633  0.0000  0.0000  0.0000  0.0000
   70  0.00025    0.0  394.6   -17.20623   -17.29657   -17.29095   -17.29070  0.0000  0.0000  0.0000  0.0000
   71  0.00023    0.0  387.7   -17.20610   -16.74982   -16.74429   -16.74406  0.0000  0.0000  0.0000  0.0000
   72  0.00023    0.0  377.1   -17.20594   -16.74425   -16.73887   -16.73865  0.0000  0.0000  0.0000  0.0000
   73  0.00023    0.0  362.8   -17.20573   -16.73906   -16.73389   -16.73367  0.0000  0.0000  0.0000  0.0000
   74  0.00022    0.0  345.0   -17.20546   -16.73427   -16.72935   -16.72914  0.0000  0.0000  0.0000  0.0000
   75  0.00020    0.0  324.0   -17.20516   -17.30050   -17.29589   -17.29569  0.0000  0.0000  0.0000  0.0000
   76  0.00021    0.0  300.2   -17.20478   -16.72577   -16.72150   -16.72129  0.0000  0.0000  0.0000  0.0000
   77  0.00018    0.0  274.0   -17.20438   -16.72212   -16.71821   -16.71803  0.0000  0.0000  0.0000  0.0000
   78  0.00018    0.0  246.3   -17.20398   -16.71880   -16.71529   -16.71511  0.0000  0.0000  0.0000  0.0000
   79  0.00017    0.0  217.5   -17.20356   -16.71582   -16.71272   -16.71255  0.0000  0.0000  0.0000  0.0000
   80  0.00015    0.0  188.4   -17.20313   -17.30142   -17.29873   -17.29858  0.0000  0.0000  0.0000  0.0000
   81  0.00011    0.0  159.6   -17.20268   -16.71078   -16.70850   -16.70839  0.0000  0.0000  0.0000  0.0000
   82  0.00010    0.0  131.8   -17.20226   -16.70875   -16.70688   -16.70677  0.0000  0.0000  0.0000  0.0000
   83  0.00010    0.0  105.5   -17.20187   -16.70700   -16.70550   -16.70540  0.0000  0.0000  0.0000  0.0000
   84  0.00008    0.0   81.5   -17.20153   -17.30103   -17.29987   -17.29979  0.0000  0.0000  0.0000  0.0000
   85  0.00006    0.0   60.4   -17.20120   -16.70424   -16.70338   -16.70332  0.0000  0.0000  0.0000  0.0000
   86  0.00004    0.0   42.7   -17.20093   -16.70325   -16.70264   -16.70260  0.0000  0.0000  0.0000  0.0000
   87  0.00003    0.0   28.9   -17.20072   -16.70245   -16.70204   -16.70201  0.0000  0.0000  0.0000  0.0000
   88  0.00004    0.0   19.5   -17.20059   -16.70187   -16.70160   -16.70156  0.0000  0.0000  0.0000  0.0000
   89  0.00004    0.0   14.4   -17.20053   -16.70147   -16.70127   -16.70122  0.0000  0.0000  0.0000  0.0000
   90  0.00004    0.0   13.9   -17.20053   -16.70122   -16.70103   -16.70099  0.0000  0.0000  0.0000  0.0000
   91  0.00004    0.0   17.8   -17.20059   -16.70111   -16.70086   -16.70082  0.0000  0.0000  0.0000  0.0000
   92  0.00006    0.0   25.9   -17.20072   -16.70113   -16.70076   -16.70070  0.0000  0.0000  0.0000  0.0000
   93  0.00009    0.0   37.9   -17.20092   -16.70126   -16.70072   -16.70063  0.0000  0.0000  0.0000  0.0000
   94  0.00011    0.0   53.6   -17.20118   -16.70148   -16.70071   -16.70060  0.0000  0.0000  0.0000  0.0000
   95  0.00012    0.0   72.5   -17.20147   -16.70175   -16.70072   -16.70060  0.0000  0.0000  0.0000  0.0000
   96  0.00013    0.0   94.1   -17.20180   -16.70207   -16.70073   -16.70060  0.0000  0.0000  0.0000  0.0000
   97  0.00015    0.0  117.8   -17.20217   -16.70244   -16.70076   -16.70061  0.0000  0.0000  0.0000  0.0000
   98  0.00017    0.0  143.0   -17.20258   -17.30284   -17.30081   -17.30063  0.0000  0.0000  0.0000  0.0000
   99  0.00019    0.0  168.9   -17.20297   -16.70321   -16.70080   -16.70062  0.0000  0.0000  0.0000  0.0000

 * Physical Quantities at step:   100

   from rhoofr: total integrated electronic density
   in g-space =    8.000000   in r-space =   8.000000
 


                total energy =      -17.20337 Hartree a.u.
              kinetic energy =        7.89488 Hartree a.u.
        electrostatic energy =      -23.63568 Hartree a.u.
                         esr =        0.00199 Hartree a.u.
                       eself =       30.31961 Hartree a.u.
      pseudopotential energy =       -1.77286 Hartree a.u.
  n-l pseudopotential energy =        4.52955 Hartree a.u.
 exchange-correlation energy =       -4.21926 Hartree a.u.
           average potential =        0.00000 Hartree a.u.



   Eigenvalues (eV), kp =   1 , spin =  1

  -25.50  -12.80   -9.50   -7.19
 
   Allocated memory (kb) =     94136
 
   CELL_PARAMETERS
   20.00000000    0.00000000    0.00000000
    0.00000000   20.00000000    0.00000000
    0.00000000    0.00000000   20.00000000

   System Density [g/cm^3] :     0.0280


   Center of mass square displacement (a.u.):   0.000000

   ATOMIC_POSITIONS
   O    0.997384E+01  0.997384E+01  0.100000E+02
   H    0.117766E+02  0.984081E+01  0.100003E+02
   H    0.984085E+01  0.117767E+02  0.999973E+01
 
   ATOMIC_VELOCITIES
   O   -0.914210E-04 -0.913682E-04  0.233414E-06
   H    0.827931E-03 -0.969419E-04  0.842079E-06
   H   -0.965636E-04  0.827887E-03 -0.270939E-05
 


   Partial temperatures (for each ionic specie) 
   Species  Temp (K)   Mean Square Displacement (a.u.)
        1      51.29     0.0014
        2     266.63     0.0292
 
  nfi    ekinc  temph  tempp        etot      enthal       econs       econt    vnhh   xnhh0    vnhp   xnhp0
  100  0.00019    0.0  194.8   -17.20337   -17.30364   -17.30086   -17.30067  0.0000  0.0000  0.0000  0.0000

   writing restart file: /h2o_mol_51.save
   restart file written in    0.071 sec.

  101  0.00028    0.0  220.0   -17.20373   -16.70398   -16.70084   -16.70057  0.0000  0.0000  0.0000  0.0000
  102  0.00029    0.0  243.8   -17.20410   -16.70436   -16.70089   -16.70060  0.0000  0.0000  0.0000  0.0000
  103  0.00028    0.0  265.6   -17.20444   -16.70471   -16.70092   -16.70064  0.0000  0.0000  0.0000  0.0000
  104  0.00027    0.0  285.0   -17.20476   -16.70502   -16.70096   -16.70069  0.0000  0.0000  0.0000  0.0000
  105  0.00024    0.0  301.6   -17.20503   -17.30529   -17.30099   -17.30075  0.0000  0.0000  0.0000  0.0000
  106  0.00027    0.0  315.3   -17.20522   -16.70546   -16.70096   -16.70069  0.0000  0.0000  0.0000  0.0000
  107  0.00024    0.0  325.9   -17.20539   -16.70564   -16.70100   -16.70076  0.0000  0.0000  0.0000  0.0000
  108  0.00023    0.0  333.1   -17.20551   -16.70576   -16.70102   -16.70079  0.0000  0.0000  0.0000  0.0000
  109  0.00024    0.0  336.9   -17.20557   -16.70583   -16.70103   -16.70079  0.0000  0.0000  0.0000  0.0000
  110  0.00026    0.0  337.3   -17.20557   -16.70583   -16.70103   -16.70077  0.0000  0.0000  0.0000  0.0000
  111  0.00028    0.0  334.1   -17.20552   -16.70578   -16.70102   -16.70074  0.0000  0.0000  0.0000  0.0000
  112  0.00029    0.0  327.6   -17.20542   -16.70568   -16.70102   -16.70072  0.0000  0.0000  0.0000  0.0000
  113  0.00031    0.0  318.0   -17.20529   -16.70555   -16.70102   -16.70072  0.0000  0.0000  0.0000  0.0000
  114  0.00031    0.0  305.6   -17.20512   -16.70538   -16.70103   -16.70072  0.0000  0.0000  0.0000  0.0000
  115  0.00028    0.0  290.9   -17.20490   -16.70516   -16.70102   -16.70074  0.0000  0.0000  0.0000  0.0000
  116  0.00024    0.0  274.3   -17.20464   -16.70491   -16.70100   -16.70075  0.0000  0.0000  0.0000  0.0000
  117  0.00021    0.0  256.9   -17.20280   -16.70307   -16.69941   -16.69920  0.0000  0.0000  0.0000  0.0000
  118  0.00019    0.0  238.9   -17.20253   -16.70280   -16.69939   -16.69921  0.0000  0.0000  0.0000  0.0000
  119  0.00017    0.0  220.4   -17.20225   -16.70252   -16.69938   -16.69921  0.0000  0.0000  0.0000  0.0000
  120  0.00016    0.0  201.4   -17.20196   -16.70223   -16.69936   -16.69920  0.0000  0.0000  0.0000  0.0000
  121  0.00014    0.0  182.4   -17.20166   -16.70193   -16.69933   -16.69919  0.0000  0.0000  0.0000  0.0000
  122  0.00013    0.0  163.5   -17.20136   -16.70163   -16.69930   -16.69918  0.0000  0.0000  0.0000  0.0000
  123  0.00012    0.0  144.9   -17.20109   -16.70136   -16.69929   -16.69917  0.0000  0.0000  0.0000  0.0000
  124  0.00012    0.0  127.1   -17.20083   -16.70110   -16.69929   -16.69917  0.0000  0.0000  0.0000  0.0000
  125  0.00010    0.0  110.2   -17.20057   -16.70084   -16.69927   -16.69917  0.0000  0.0000  0.0000  0.0000
  126  0.00008    0.0   94.6   -17.20033   -16.70060   -16.69925   -16.69917  0.0000  0.0000  0.0000  0.0000
  127  0.00007    0.0   80.5   -17.20011   -16.70039   -16.69924   -16.69917  0.0000  0.0000  0.0000  0.0000
  128  0.00007    0.0   68.2   -17.19994   -16.70021   -16.69924   -16.69917  0.0000  0.0000  0.0000  0.0000
  129  0.00006    0.0   57.8   -17.19979   -16.70006   -16.69924   -16.69918  0.0000  0.0000  0.0000  0.0000
  130  0.00005    0.0   49.7   -17.19967   -16.69994   -16.69923   -16.69918  0.0000  0.0000  0.0000  0.0000
  131  0.00004    0.0   43.8   -17.19957   -16.69984   -16.69922   -16.69918  0.0000  0.0000  0.0000  0.0000
  132  0.00004    0.0   40.2   -17.19951   -16.69979   -16.69921   -16.69918  0.0000  0.0000  0.0000  0.0000
  133  0.00005    0.0   38.8   -17.19950   -16.69978   -16.69922   -16.69918  0.0000  0.0000  0.0000  0.0000
  134  0.00006    0.0   39.8   -17.19953   -16.69980   -16.69924   -16.69918  0.0000  0.0000  0.0000  0.0000
  135  0.00006    0.0   43.0   -17.19958   -16.69986   -16.69924   -16.69919  0.0000  0.0000  0.0000  0.0000
  136  0.00005    0.0   48.6   -17.19966   -16.69994   -16.69924   -16.69919  0.0000  0.0000  0.0000  0.0000
  137  0.00006    0.0   56.5   -17.19978   -16.70006   -16.69925   -16.69919  0.0000  0.0000  0.0000  0.0000
  138  0.00007    0.0   66.7   -17.19994   -16.70021   -16.69926   -16.69919  0.0000  0.0000  0.0000  0.0000
  139  0.00008    0.0   79.1   -17.20014   -16.70041   -16.69928   -16.69920  0.0000  0.0000  0.0000  0.0000
  140  0.00009    0.0   93.6   -17.20035   -16.70062   -16.69929   -16.69920  0.0000  0.0000  0.0000  0.0000
  141  0.00009    0.0  109.8   -17.20058   -16.70086   -16.69929   -16.69920  0.0000  0.0000  0.0000  0.0000
  142  0.00011    0.0  127.6   -17.20085   -16.70112   -16.69930   -16.69919  0.0000  0.0000  0.0000  0.0000
  143  0.00014    0.0  146.7   -17.20115   -16.70142   -16.69933   -16.69919  0.0000  0.0000  0.0000  0.0000
  144  0.00016    0.0  166.7   -17.20146   -16.70173   -16.69936   -16.69920  0.0000  0.0000  0.0000  0.0000
  145  0.00017    0.0  187.5   -17.20178   -16.70205   -16.69938   -16.69921  0.0000  0.0000  0.0000  0.0000
  146  0.00017    0.0  208.8   -17.20210   -16.70237   -16.69939   -16.69922  0.0000  0.0000  0.0000  0.0000
  147  0.00019    0.0  230.3   -17.20243   -16.70269   -16.69941   -16.69922  0.0000  0.0000  0.0000  0.0000
  148  0.00021    0.0  251.8   -17.20276   -16.70302   -16.69944   -16.69923  0.0000  0.0000  0.0000  0.0000
  149  0.00022    0.0  272.9   -17.20309   -16.70335   -16.69946   -16.69924  0.0000  0.0000  0.0000  0.0000
  150  0.00023    0.0  293.3   -17.20340   -16.70366   -16.69948   -16.69926  0.0000  0.0000  0.0000  0.0000
  151  0.00022    0.0  313.0   -17.20524   -16.70551   -16.70105   -16.70082  0.0000  0.0000  0.0000  0.0000
  152  0.00023    0.0  331.5   -17.20551   -16.70578   -16.70106   -16.70082  0.0000  0.0000  0.0000  0.0000
  153  0.00025    0.0  348.1   -17.20576   -16.70603   -16.70107   -16.70082  0.0000  0.0000  0.0000  0.0000
  154  0.00027    0.0  362.3   -17.20598   -16.70624   -16.70108   -16.70081  0.0000  0.0000  0.0000  0.0000
  155  0.00028    0.0  373.6   -17.20614   -16.70641   -16.70109   -16.70081  0.0000  0.0000  0.0000  0.0000
  156  0.00028    0.0  381.8   -17.20625   -16.70652   -16.70108   -16.70079  0.0000  0.0000  0.0000  0.0000
  157  0.00030    0.0  386.6   -17.20631   -16.70658   -16.70107   -16.70077  0.0000  0.0000  0.0000  0.0000
  158  0.00032    0.0  387.8   -17.20633   -16.70660   -16.70108   -16.70076  0.0000  0.0000  0.0000  0.0000
  159  0.00033    0.0  385.5   -17.20630   -16.70657   -16.70107   -16.70075  0.0000  0.0000  0.0000  0.0000
  160  0.00032    0.0  379.6   -17.20620   -16.70647   -16.70106   -16.70074  0.0000  0.0000  0.0000  0.0000
  161  0.00030    0.0  370.3   -17.20604   -16.70631   -16.70103   -16.70073  0.0000  0.0000  0.0000  0.0000
  162  0.00029    0.0  357.7   -17.20583   -16.70610   -16.70100   -16.70071  0.0000  0.0000  0.0000  0.0000
  163  0.00027    0.0  342.0   -17.20559   -17.30586   -17.30098   -17.30071  0.0000  0.0000  0.0000  0.0000
  164  0.00029    0.0  323.4   -17.20529   -17.30553   -17.30092   -17.30063  0.0000  0.0000  0.0000  0.0000
  165  0.00033    0.0  302.2   -17.20494   -17.30518   -17.30088   -17.30054  0.0000  0.0000  0.0000  0.0000
  166  0.00037    0.0  278.8   -17.20453   -16.70478   -16.70081   -16.70043  0.0000  0.0000  0.0000  0.0000
  167  0.00026    0.0  253.7   -17.20412   -16.70439   -16.70078   -16.70052  0.0000  0.0000  0.0000  0.0000
  168  0.00018    0.0  227.5   -17.20371   -16.70398   -16.70074   -16.70055  0.0000  0.0000  0.0000  0.0000
  169  0.00015    0.0  200.6   -17.20330   -16.70356   -16.70070   -16.70055  0.0000  0.0000  0.0000  0.0000
  170  0.00015    0.0  173.9   -17.20288   -16.70314   -16.70066   -16.70051  0.0000  0.0000  0.0000  0.0000
  171  0.00018    0.0  147.9   -17.20249   -17.30274   -17.30063   -17.30045  0.0000  0.0000  0.0000  0.0000
  172  0.00012    0.0  123.5   -17.20209   -16.70232   -16.70056   -16.70043  0.0000  0.0000  0.0000  0.0000
  173  0.00015    0.0  101.4   -17.20175   -16.70200   -16.70055   -16.70040  0.0000  0.0000  0.0000  0.0000
  174  0.00015    0.0   82.2   -17.20145   -16.90170   -16.90053   -16.90038  0.0000  0.0000  0.0000  0.0000
  175  0.00013    0.0   66.6   -17.20120   -16.90145   -16.90050   -16.90037  0.0000  0.0000  0.0000  0.0000
  176  0.00012    0.0   54.8   -17.20099   -16.90124   -16.90046   -16.90034  0.0000  0.0000  0.0000  0.0000
  177  0.00011    0.0   47.2   -17.20085   -16.90110   -16.90043   -16.90032  0.0000  0.0000  0.0000  0.0000
  178  0.00008    0.0   43.9   -17.20078   -16.90104   -16.90041   -16.90033  0.0000  0.0000  0.0000  0.0000
  179  0.00007    0.0   44.7   -17.20079   -16.90104   -16.90041   -16.90034  0.0000  0.0000  0.0000  0.0000
  180  0.00005    0.0   49.6   -17.20085   -16.90111   -16.90040   -16.90035  0.0000  0.0000  0.0000  0.0000
  181  0.00005    0.0   58.4   -17.20099   -17.50124   -17.50041   -17.50036  0.0000  0.0000  0.0000  0.0000
  182  0.00011    0.0   71.0   -17.20116   -17.50140   -17.50039   -17.50028  0.0000  0.0000  0.0000  0.0000
  183  0.00026    0.0   86.8   -17.20139   -16.70163   -16.70039   -16.70013  0.0000  0.0000  0.0000  0.0000
  184  0.00027    0.0  105.6   -17.20159   -16.90185   -16.90035   -16.90008  0.0000  0.0000  0.0000  0.0000
  185  0.00023    0.0  126.6   -17.20190   -16.70216   -16.70036   -16.70013  0.0000  0.0000  0.0000  0.0000
  186  0.00018    0.0  149.4   -17.20222   -16.70249   -16.70036   -16.70019  0.0000  0.0000  0.0000  0.0000
  187  0.00015    0.0  173.1   -17.20257   -16.70284   -16.70038   -16.70023  0.0000  0.0000  0.0000  0.0000
  188  0.00015    0.0  197.0   -17.20294   -16.70321   -16.70040   -16.70025  0.0000  0.0000  0.0000  0.0000
  189  0.00017    0.0  220.5   -17.20331   -16.70358   -16.70043   -16.70026  0.0000  0.0000  0.0000  0.0000
  190  0.00020    0.0  243.0   -17.20367   -17.30393   -17.30047   -17.30028  0.0000  0.0000  0.0000  0.0000
  191  0.00019    0.0  263.9   -17.20399   -17.30422   -17.30046   -17.30027  0.0000  0.0000  0.0000  0.0000
  192  0.00022    0.0  282.8   -17.20426   -16.70450   -16.70047   -16.70025  0.0000  0.0000  0.0000  0.0000
  193  0.00022    0.0  299.3   -17.20451   -16.70478   -16.70051   -16.70029  0.0000  0.0000  0.0000  0.0000
  194  0.00024    0.0  313.0   -17.20473   -16.70500   -16.70054   -16.70030  0.0000  0.0000  0.0000  0.0000
  195  0.00026    0.0  323.7   -17.20490   -16.70517   -16.70055   -16.70030  0.0000  0.0000  0.0000  0.0000
  196  0.00026    0.0  331.1   -17.20501   -16.70528   -16.70056   -16.70030  0.0000  0.0000  0.0000  0.0000
  197  0.00026    0.0  335.0   -17.20507   -16.70533   -16.70056   -16.70030  0.0000  0.0000  0.0000  0.0000
  198  0.00026    0.0  335.4   -17.20508   -16.70535   -16.70057   -16.70030  0.0000  0.0000  0.0000  0.0000
  199  0.00026    0.0  332.3   -17.20506   -17.30532   -17.30059   -17.30032  0.0000  0.0000  0.0000  0.0000

 * Physical Quantities at step:   200

   from rhoofr: total integrated electronic density
   in g-space =    8.000000   in r-space =   8.000000
 


                total energy =      -17.20497 Hartree a.u.
              kinetic energy =        7.80633 Hartree a.u.
        electrostatic energy =      -23.59375 Hartree a.u.
                         esr =        0.00118 Hartree a.u.
                       eself =       30.31961 Hartree a.u.
      pseudopotential energy =       -1.74434 Hartree a.u.
  n-l pseudopotential energy =        4.51370 Hartree a.u.
 exchange-correlation energy =       -4.18692 Hartree a.u.
           average potential =        0.00000 Hartree a.u.



   Eigenvalues (eV), kp =   1 , spin =  1

  -24.90  -13.13   -8.82   -7.07
 
   Allocated memory (kb) =     94136
 
   CELL_PARAMETERS
   20.00000000    0.00000000    0.00000000
    0.00000000   20.00000000    0.00000000
    0.00000000    0.00000000   20.00000000

   System Density [g/cm^3] :     0.0280


   Center of mass square displacement (a.u.):   0.000000

   ATOMIC_POSITIONS
   O    0.999211E+01  0.999216E+01  0.100000E+02
   H    0.118289E+02  0.964196E+01  0.100024E+02
   H    0.964243E+01  0.118290E+02  0.999756E+01
 
   ATOMIC_VELOCITIES
   O   -0.823051E-04 -0.822689E-04  0.643465E-06
   H    0.103717E-02 -0.379607E-03  0.377913E-05
   H   -0.378725E-03  0.103776E-02 -0.892684E-05
 


   Partial temperatures (for each ionic specie) 
   Species  Temp (K)   Mean Square Displacement (a.u.)
        1      41.58     0.0001
        2     468.23     0.0104
 
  nfi    ekinc  temph  tempp        etot      enthal       econs       econt    vnhh   xnhh0    vnhp   xnhp0
  200  0.00025    0.0  326.0   -17.20497   -17.30521   -17.30057   -17.30032  0.0000  0.0000  0.0000  0.0000

   writing restart file: /h2o_mol_51.save
   restart file written in    0.089 sec.



   Averaged Physical Quantities
                      accomulated      this run
   ekinc         :        0.00018       0.00018 (AU)
   ekin          :        7.82053       7.82053 (AU)
   epot          :      -29.53578     -29.53578 (AU)
   total energy  :      -17.20285     -17.20285 (AU)
   temperature   :      201.05257     201.05257 (K )
   enthalpy      :      -16.90456     -16.90456 (AU)
   econs         :      -16.90170     -16.90170 (AU)
   pressure      :        0.00000       0.00000 (Gpa)
   volume        :     8000.00000    8000.00000 (AU)


     initialize   :     1.06s CPU
     total_time   :   106.08s CPU (     200 calls,   0.530 s avg)
     formf        :     0.05s CPU
     rhoofr       :    19.19s CPU (     200 calls,   0.096 s avg)
     vofrho       :    64.63s CPU (     200 calls,   0.323 s avg)
     dforce       :    10.87s CPU (     400 calls,   0.027 s avg)
     calphi       :     0.19s CPU (     200 calls,   0.001 s avg)
     ortho        :     0.62s CPU (     200 calls,   0.003 s avg)
     ortho_iter   :     0.01s CPU (     200 calls,   0.000 s avg)
     rsg          :     0.00s CPU (     200 calls,   0.000 s avg)
     rhoset       :     0.05s CPU (     200 calls,   0.000 s avg)
     updatc       :     0.03s CPU (     200 calls,   0.000 s avg)
     newd         :     0.24s CPU (     200 calls,   0.001 s avg)
     calbec       :     0.00s CPU (     201 calls,   0.000 s avg)
     prefor       :     0.30s CPU (     401 calls,   0.001 s avg)
     strucf       :     0.54s CPU (     201 calls,   0.003 s avg)
     nlfl         :     0.00s CPU (     200 calls,   0.000 s avg)
     nlfq         :     0.67s CPU (     200 calls,   0.003 s avg)
     rhov         :     4.55s CPU (     200 calls,   0.023 s avg)
     nlsm1        :     0.47s CPU (     601 calls,   0.001 s avg)
     nlsm2        :     0.67s CPU (     200 calls,   0.003 s avg)
     fft          :    29.04s CPU (    1800 calls,   0.016 s avg)
     ffts         :     6.48s CPU (     400 calls,   0.016 s avg)
     fftw         :    12.11s CPU (    1200 calls,   0.010 s avg)
     fftb         :     0.41s CPU (   31914 calls,   0.000 s avg)
 
 
 
     CP           :  1m48.85s CPU time,     1m58.66s wall time

 
   This run was terminated on:  16:13:10  28Apr2008            

=------------------------------------------------------------------------------=
   JOB DONE.
=------------------------------------------------------------------------------=