argument 1 = s2-example5.nw argument 2 = -np argument 3 = 2 ============================== echo of input deck ============================== title "Born-Oppenheimer AIMD simulation of s2-dimer" echo scratch_dir ./scratch permanent_dir ./perm start s2-bomd geometry S 0.0 0.0 0.0 S 0.0 0.0 1.95 end nwpw simulation_cell SC 20.0 end cutoff 15.0 mult 3 xc lda lmbfgs end task pspw energy nwpw bo_steps 1 500 bo_time_step 10.0 end task pspw born-oppenheimer ================================================================================ Northwest Computational Chemistry Package (NWChem) 5.1.1 -------------------------------------------------------- Environmental Molecular Sciences Laboratory Pacific Northwest National Laboratory Richland, WA 99352 COPYRIGHT (C) 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010 Pacific Northwest National Laboratory, Battelle Memorial Institute. >>> All Rights Reserved <<< DISCLAIMER ---------- This material was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the United States Department of Energy, nor Battelle, nor any of their employees, MAKES ANY WARRANTY, EXPRESS OR IMPLIED, OR ASSUMES ANY LEGAL LIABILITY OR RESPONSIBILITY FOR THE ACCURACY, COMPLETENESS, OR USEFULNESS OF ANY INFORMATION, APPARATUS, PRODUCT, SOFTWARE, OR PROCESS DISCLOSED, OR REPRESENTS THAT ITS USE WOULD NOT INFRINGE PRIVATELY OWNED RIGHTS. LIMITED USE ----------- This software (including any documentation) is being made available to you for your internal use only, solely for use in performance of work directly for the U.S. Federal Government or work under contracts with the U.S. Department of Energy or other U.S. Federal Government agencies. This software is a version which has not yet been evaluated and cleared for commercialization. Adherence to this notice may be necessary for the author, Battelle Memorial Institute, to successfully assert copyright in and commercialize this software. This software is not intended for duplication or distribution to third parties without the permission of the Manager of Software Products at Pacific Northwest National Laboratory, Richland, Washington, 99352. ACKNOWLEDGMENT -------------- This software and its documentation were produced with Government support under Contract Number DE-AC05-76RL01830 awarded by the United States Department of Energy. The Government retains a paid-up non-exclusive, irrevocable worldwide license to reproduce, prepare derivative works, perform publicly and display publicly by or for the Government, including the right to distribute to other Government contractors. Job information --------------- hostname = seattle-1 program = nwchem date = Tue Aug 10 15:41:34 2010 compiled = Wed_Jul_21_12:26:22_2010 source = /home/bylaska/nwchem-releases/nwchem nwchem branch = Development input = s2-example5.nw prefix = s2-bomd. data base = ./perm/s2-bomd.db status = startup nproc = 2 time left = -1s Memory information ------------------ heap = 13107201 doubles = 100.0 Mbytes stack = 13107201 doubles = 100.0 Mbytes global = 26214400 doubles = 200.0 Mbytes (distinct from heap & stack) total = 52428802 doubles = 400.0 Mbytes verify = yes hardfail = no Directory information --------------------- 0 permanent = ./perm 0 scratch = ./scratch NWChem Input Module ------------------- Born-Oppenheimer AIMD simulation of s2-dimer -------------------------------------------- Scaling coordinates for geometry "geometry" by 1.889725989 (inverse scale = 0.529177249) ORDER OF PRIMARY AXIS IS BEING SET TO 4 D4H symmetry detected ------ auto-z ------ Geometry "geometry" -> "" ------------------------- Output coordinates in angstroms (scale by 1.889725989 to convert to a.u.) No. Tag Charge X Y Z ---- ---------------- ---------- -------------- -------------- -------------- 1 S 16.0000 0.00000000 0.00000000 -0.97500000 2 S 16.0000 0.00000000 0.00000000 0.97500000 Atomic Mass ----------- S 31.972070 Effective nuclear repulsion energy (a.u.) 69.4714747405 Nuclear Dipole moment (a.u.) ---------------------------- X Y Z ---------------- ---------------- ---------------- 0.0000000000 0.0000000000 0.0000000000 Symmetry information -------------------- Group name D4h Group number 28 Group order 16 No. of unique centers 1 Symmetry unique atoms 1 Z-matrix (autoz) -------- Units are Angstrom for bonds and degrees for angles Type Name I J K L M Value ----------- -------- ----- ----- ----- ----- ----- ---------- 1 Stretch 1 2 1.95000 XYZ format geometry ------------------- 2 geometry S 0.00000000 0.00000000 -0.97500000 S 0.00000000 0.00000000 0.97500000 ============================================================================== internuclear distances ------------------------------------------------------------------------------ center one | center two | atomic units | angstroms ------------------------------------------------------------------------------ 2 S | 1 S | 3.68497 | 1.95000 ------------------------------------------------------------------------------ number of included internuclear distances: 1 ============================================================================== **************************************************** * * * NWPW PSPW Calculation * * * * [ (Grassman/Stiefel manifold implementation) ] * * * * [ NorthWest Chemistry implementation ] * * * * version #5.10 06/12/02 * * * * This code was developed by Eric J. Bylaska, * * and was based upon algorithms and code * * developed by the group of Prof. John H. Weare * * * **************************************************** >>> JOB STARTED AT Tue Aug 10 15:41:34 2010 <<< ================ input data ======================== library name resolved from: compiled reference NWCHEM_NWPW_LIBRARY set to: < /home/bylaska/nwchem-releases/nwchem/src/nwpw/libraryps/> Generating 1d pseudopotential for S Generated formatted_filename: ./perm/S.vpp random planewave guess, initial psi: s2-bomd.movecs - spin, nalpha, nbeta: 2 7 5 input psi filename:./perm/s2-bomd.movecs Warning - Gram-Schmidt being performed on psi: 12.0000000000000 7.95390810283379 12.0000000000000 4.04609189716621 number of processors used: 2 processor grid : 2 x 1 parallel mapping : hilbert parallel mapping : balanced options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = LDA (Vosko et al) parameterization elements involved in the cluster: 1: S core charge: 6.0000 lmax= 2 comment : Hamann pseudopotential pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.843 0.955 0.955 total charge: 0.000 atomic composition: S : 2 number of electrons: spin up= 7 ( 7 per task) down= 5 ( 5 per task) (fourier space) number of orbitals : spin up= 7 ( 7 per task) down= 5 ( 5 per task) (fourier space) supercell: cell_name: cell_default lattice: a1=< 20.000 0.000 0.000 > a2=< 0.000 20.000 0.000 > a3=< 0.000 0.000 20.000 > reciprocal: b1=< 0.314 0.000 0.000 > b2=< 0.000 0.314 0.000 > b3=< 0.000 0.000 0.314 > lattice: a= 20.000 b= 20.000 c= 20.000 alpha= 90.000 beta= 90.000 gamma= 90.000 omega= 8000.0 density cutoff= 30.000 fft= 50x 50x 50( 31335 waves 15668 per task) wavefnc cutoff= 15.000 fft= 50x 50x 50( 11060 waves 5530 per task) ewald summation: cut radius= 6.37 and 1 madelung= 1.76011888 technical parameters: time step= 5.80 ficticious mass= 400000.0 tolerance=.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) == Energy Calculation == ============ Grassmann lmbfgs iteration ============ >>> ITERATION STARTED AT Tue Aug 10 15:41:36 2010 <<< iter. Energy DeltaE DeltaRho ------------------------------------------------------ - 15 steepest descent iterations performed 10 -0.1993134867E+02 -0.22934E+00 0.10437E+01 - 10 steepest descent iterations performed 20 -0.2028238962E+02 -0.31716E-02 0.51553E-02 30 -0.2034037491E+02 -0.12827E-01 0.98166E-02 - 10 steepest descent iterations performed 40 -0.2040952281E+02 -0.86197E-03 0.59873E-03 50 -0.2041178393E+02 -0.50792E-04 0.77800E-04 60 -0.2041189515E+02 -0.85957E-06 0.26750E-05 70 -0.2041189672E+02 -0.75520E-07 0.18532E-07 *** tolerance ok. iteration terminated >>> ITERATION ENDED AT Tue Aug 10 15:41:58 2010 <<< == Summary Of Results == number of electrons: spin up= 7.00000 down= 5.00000 (real space) total energy : -0.2041189672E+02 ( -0.10206E+02/ion) total orbital energy: -0.4907546453E+01 ( -0.40896E+00/electron) hartree energy : 0.1651320391E+02 ( 0.13761E+01/electron) exc-corr energy : -0.4294211808E+01 ( -0.35785E+00/electron) ion-ion energy : -0.3103701257E+00 ( -0.15519E+00/ion) kinetic (planewave) : 0.7453123441E+01 ( 0.62109E+00/electron) V_local (planewave) : -0.4438319979E+02 ( -0.36986E+01/electron) V_nl (planewave) : 0.4609557651E+01 ( 0.38413E+00/electron) V_Coul (planewave) : 0.3302640782E+02 ( 0.27522E+01/electron) V_xc. (planewave) : -0.5613435576E+01 ( -0.46779E+00/electron) Virial Coefficient : -0.1658455008E+01 orbital energies: -0.2060464E+00 ( -5.607eV) -0.2060468E+00 ( -5.607eV) -0.3223469E+00 ( -8.772eV) -0.2915946E+00 ( -7.935eV) -0.3223475E+00 ( -8.772eV) -0.2915951E+00 ( -7.935eV) -0.3562787E+00 ( -9.695eV) -0.3341330E+00 ( -9.092eV) -0.5699889E+00 ( -15.510eV) -0.5316549E+00 ( -14.467eV) -0.7497389E+00 ( -20.402eV) -0.7257747E+00 ( -19.749eV) Total PSPW energy : -0.2041189672E+02 === Spin Contamination === = 2.00000000000000 = 2.00304460940712 == Center of Charge == spin up ( 0.0000, 0.0000, -0.0001 ) spin down ( 0.0000, 0.0000, 0.0001 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) crystal ( 0.0000, 0.0000, 0.0000 ) == Crystal Dipole == mu = ( 0.0002, -0.0002, 0.0000 ) au |mu| = 0.0003 au, 0.0007 Debye == Molecular Dipole wrt Center of Mass == mu = ( 0.0002, -0.0002, 0.0000 ) au |mu| = 0.0003 au, 0.0007 Debye output psi filename:./perm/s2-bomd.movecs == Timing == cputime in seconds prologue : 0.237376E+01 main loop : 0.217059E+02 epilogue : 0.960310E-01 total : 0.241757E+02 cputime/step: 0.119263E+00 ( 182 evalulations, 66 linesearches) Time spent doing total step FFTs : 0.106724E+01 0.586393E-02 dot products : 0.685108E+00 0.376433E-02 geodesic : 0.195585E+01 0.107464E-01 ffm_dgemm : 0.238275E+00 0.130920E-02 fmf_dgemm : 0.166293E+01 0.913695E-02 m_diagonalize : 0.399656E-02 0.219591E-04 - m_tredq : 0.000000E+00 0.000000E+00 - m_getdiags : 0.000000E+00 0.000000E+00 - m_tqliq : 0.000000E+00 0.000000E+00 - m_eigsrt : 0.000000E+00 0.000000E+00 exchange correlation : 0.334427E+01 0.183751E-01 local pseudopotentials : 0.589132E-03 0.323699E-05 non-local pseudopotentials : 0.760957E+00 0.418108E-02 hartree potentials : 0.592434E-01 0.325513E-03 ion-ion interaction : 0.373149E-02 0.205027E-04 structure factors : 0.230284E-01 0.126530E-03 phase factors : 0.109673E-04 0.602596E-07 masking and packing : 0.143174E+01 0.786668E-02 queue fft : 0.116734E+02 0.641398E-01 queue fft (serial) : 0.552551E+01 0.303600E-01 queue fft (message passing): 0.576138E+01 0.316560E-01 HFX potential : 0.000000E+00 0.000000E+00 qmmm LJ : 0.000000E+00 0.000000E+00 qmmm residual Q : 0.000000E+00 0.000000E+00 >>> JOB COMPLETED AT Tue Aug 10 15:41:58 2010 <<< Task times cpu: 23.4s wall: 24.2s NWChem Input Module ------------------- >>>> PSPW Parallel Module - Born-Oppenheimer <<<< **************************************************** * * * NWPW PSPW Calculation * * * * [ Born-Oppenheimer molecular ] * * [ dynamics simulation ] * * * * [ NorthWest Chemistry implementation ] * * * * version #1.00 07/18/06 * * * * This code was developed by Eric J. Bylaska. * * * **************************************************** >>> JOB STARTED AT Tue Aug 10 15:41:58 2010 <<< ================ input data ======================== input psi filename:./perm/s2-bomd.movecs number of processors used: 2 processor grid : 2 x 1 parallel mapping : hilbert parallel mapping : balanced options: boundary conditions = periodic (version3) electron spin = unrestricted exchange-correlation = LDA (Vosko et al) parameterization elements involved in the cluster: 1: S core charge: 6.0 lmax= 2 comment : Hamann pseudopotential pseudpotential type : 0 highest angular component : 2 local potential used : 2 number of non-local projections: 4 cutoff = 0.843 0.955 0.955 total charge: 0.000 atomic composition: S : 2 initial position of ions (au): 1 S ( 0.00000 0.00000 -1.84248 ) - atomic mass= 31.972 2 S ( 0.00000 0.00000 1.84248 ) - atomic mass= 31.972 G.C. ( 0.00000 0.00000 0.00000 ) C.O.M. ( 0.00000 0.00000 0.00000 ) initial velocity of ions after scaling (au): 1 S ( 0.00000 0.00000 0.00000 ) 2 S ( 0.00000 0.00000 0.00000 ) G.C. ( 0.00000 0.00000 0.00000 ) C.O.M. ( 0.00000 0.00000 0.00000 ) number of active electrons: spin up= 7 down= 5 (fourier space) number of active orbitals : spin up= 7 down= 5 (fourier space) supercell: cell_name: cell_default lattice: a1=< 20.000 0.000 0.000 > a2=< 0.000 20.000 0.000 > a3=< 0.000 0.000 20.000 > reciprocal: b1=< 0.314 0.000 0.000 > b2=< 0.000 0.314 0.000 > b3=< 0.000 0.000 0.314 > lattice: a= 20.000 b= 20.000 c= 20.000 alpha= 90.000 beta= 90.000 gamma= 90.000 omega= 8000.0 density cutoff= 30.000 fft= 50x 50x 50( 31335 waves 15668 per task) wavefnc cutoff= 15.000 fft= 50x 50x 50( 11060 waves 5530 per task) ewald summation: cut radius= 6.37 and 1 madelung= 1.76011888 technical parameters for minimizer: time step= 5.80 ficticious mass= 400000.0 tolerance=.100E-06 (energy) 0.100E-06 (density) maximum iterations = 1000 ( 10 inner 100 outer ) molecular dynamics parameters: translation contrained time step= 10.00 iterations= 500 ( 1 inner 500 outer ) integration agorithm= Verlet cooling/heatting rate= 0.10000E+01 (ion) initial kinetic energy= 0.00000E+00 (ion) 0.00000E+00 (c.o.m.) after scaling= 0.00000E+00 (ion) increased energy= 0.00000E+00 (ion) Constant Energy Simulation >>> ITERATION STARTED AT Tue Aug 10 15:41:58 2010 <<< iter. KE+Energy Energy KE_ion Temperature ------------------------------------------------------------------------------------ 1 -0.2041189678E+02 -0.2041189830E+02 0.15256E-05 0.48 2 -0.2041189681E+02 -0.2041190287E+02 0.60598E-05 1.60 3 -0.2041189692E+02 -0.2041191053E+02 0.13610E-04 3.35 4 -0.2041189694E+02 -0.2041192112E+02 0.24180E-04 5.73 5 -0.2041189704E+02 -0.2041193474E+02 0.37703E-04 8.74 6 -0.2041189706E+02 -0.2041195125E+02 0.54193E-04 12.38 7 -0.2041189712E+02 -0.2041197071E+02 0.73583E-04 16.65 8 -0.2041189721E+02 -0.2041199304E+02 0.95828E-04 21.52 9 -0.2041189722E+02 -0.2041201806E+02 0.12084E-03 27.00 10 -0.2041189722E+02 -0.2041204576E+02 0.14854E-03 33.07 11 -0.2041189723E+02 -0.2041207605E+02 0.17882E-03 39.73 12 -0.2041189723E+02 -0.2041210884E+02 0.21161E-03 46.95 13 -0.2041189724E+02 -0.2041214402E+02 0.24678E-03 54.73 14 -0.2041189725E+02 -0.2041218147E+02 0.28422E-03 63.05 15 -0.2041189726E+02 -0.2041222107E+02 0.32382E-03 71.89 16 -0.2041189727E+02 -0.2041226270E+02 0.36543E-03 81.24 17 -0.2041189728E+02 -0.2041230620E+02 0.40892E-03 91.07 18 -0.2041189729E+02 -0.2041235144E+02 0.45416E-03 101.37 19 -0.2041189730E+02 -0.2041239826E+02 0.50097E-03 112.13 20 -0.2041189731E+02 -0.2041244651E+02 0.54921E-03 123.30 21 -0.2041189732E+02 -0.2041249602E+02 0.59870E-03 134.89 22 -0.2041189733E+02 -0.2041254661E+02 0.64929E-03 146.85 23 -0.2041189734E+02 -0.2041259812E+02 0.70078E-03 159.17 24 -0.2041189735E+02 -0.2041265036E+02 0.75300E-03 171.83 25 -0.2041189736E+02 -0.2041270314E+02 0.80577E-03 184.79 26 -0.2041189738E+02 -0.2041275627E+02 0.85890E-03 198.04 27 -0.2041189739E+02 -0.2041280957E+02 0.91218E-03 211.54 28 -0.2041189740E+02 -0.2041286284E+02 0.96544E-03 225.27 29 -0.2041189742E+02 -0.2041291587E+02 0.10185E-02 239.20 30 -0.2041189743E+02 -0.2041296848E+02 0.10711E-02 253.31 31 -0.2041189744E+02 -0.2041302046E+02 0.11230E-02 267.55 32 -0.2041189746E+02 -0.2041307160E+02 0.11741E-02 281.92 33 -0.2041189747E+02 -0.2041312172E+02 0.12243E-02 296.37 34 -0.2041189748E+02 -0.2041317061E+02 0.12731E-02 310.87 35 -0.2041189749E+02 -0.2041321807E+02 0.13206E-02 325.40 36 -0.2041189751E+02 -0.2041326391E+02 0.13664E-02 339.93 37 -0.2041189752E+02 -0.2041330795E+02 0.14104E-02 354.43 38 -0.2041189753E+02 -0.2041334998E+02 0.14525E-02 368.86 39 -0.2041189754E+02 -0.2041338984E+02 0.14923E-02 383.20 40 -0.2041189756E+02 -0.2041342735E+02 0.15298E-02 397.42 41 -0.2041189757E+02 -0.2041346233E+02 0.15648E-02 411.48 42 -0.2041189758E+02 -0.2041349463E+02 0.15971E-02 425.37 43 -0.2041189759E+02 -0.2041352409E+02 0.16265E-02 439.05 44 -0.2041189760E+02 -0.2041355057E+02 0.16530E-02 452.49 45 -0.2041189761E+02 -0.2041357393E+02 0.16763E-02 465.67 46 -0.2041189762E+02 -0.2041359405E+02 0.16964E-02 478.56 47 -0.2041189762E+02 -0.2041361081E+02 0.17132E-02 491.13 48 -0.2041189763E+02 -0.2041362413E+02 0.17265E-02 503.36 49 -0.2041189764E+02 -0.2041363390E+02 0.17363E-02 515.22 50 -0.2041189764E+02 -0.2041364006E+02 0.17424E-02 526.70 51 -0.2041189764E+02 -0.2041364255E+02 0.17449E-02 537.76 52 -0.2041189765E+02 -0.2041364132E+02 0.17437E-02 548.39 53 -0.2041189765E+02 -0.2041363634E+02 0.17387E-02 558.57 54 -0.2041189765E+02 -0.2041362760E+02 0.17300E-02 568.28 55 -0.2041189765E+02 -0.2041361511E+02 0.17175E-02 577.50 56 -0.2041189765E+02 -0.2041359889E+02 0.17012E-02 586.22 57 -0.2041189765E+02 -0.2041357897E+02 0.16813E-02 594.42 58 -0.2041189764E+02 -0.2041355541E+02 0.16578E-02 602.09 59 -0.2041189764E+02 -0.2041352828E+02 0.16306E-02 609.22 60 -0.2041189763E+02 -0.2041349768E+02 0.16000E-02 615.80 61 -0.2041189763E+02 -0.2041346370E+02 0.15661E-02 621.82 62 -0.2041189762E+02 -0.2041342649E+02 0.15289E-02 627.27 63 -0.2041189761E+02 -0.2041338617E+02 0.14886E-02 632.16 64 -0.2041189760E+02 -0.2041334292E+02 0.14453E-02 636.48 65 -0.2041189759E+02 -0.2041329691E+02 0.13993E-02 640.23 66 -0.2041189758E+02 -0.2041324833E+02 0.13508E-02 643.40 67 -0.2041189757E+02 -0.2041319739E+02 0.12998E-02 646.01 68 -0.2041189755E+02 -0.2041314431E+02 0.12468E-02 648.06 69 -0.2041189754E+02 -0.2041308933E+02 0.11918E-02 649.56 70 -0.2041189752E+02 -0.2041303270E+02 0.11352E-02 650.51 71 -0.2041189751E+02 -0.2041297468E+02 0.10772E-02 650.92 72 -0.2041189749E+02 -0.2041291555E+02 0.10181E-02 650.81 73 -0.2041189747E+02 -0.2041285557E+02 0.95810E-03 650.19 74 -0.2041189745E+02 -0.2041279505E+02 0.89759E-03 649.08 75 -0.2041189744E+02 -0.2041273427E+02 0.83683E-03 647.50 76 -0.2041189742E+02 -0.2041267353E+02 0.77611E-03 645.45 77 -0.2041189740E+02 -0.2041261314E+02 0.71574E-03 642.97 78 -0.2041189738E+02 -0.2041255340E+02 0.65602E-03 640.08 79 -0.2041189736E+02 -0.2041249462E+02 0.59726E-03 636.79 80 -0.2041189734E+02 -0.2041243709E+02 0.53974E-03 633.14 81 -0.2041189733E+02 -0.2041238111E+02 0.48378E-03 629.14 82 -0.2041189731E+02 -0.2041232698E+02 0.42967E-03 624.83 83 -0.2041189729E+02 -0.2041227497E+02 0.37768E-03 620.23 84 -0.2041189728E+02 -0.2041222537E+02 0.32809E-03 615.37 85 -0.2041189726E+02 -0.2041217843E+02 0.28117E-03 610.28 86 -0.2041189725E+02 -0.2041213441E+02 0.23717E-03 604.99 87 -0.2041189723E+02 -0.2041209355E+02 0.19632E-03 599.53 88 -0.2041189722E+02 -0.2041205606E+02 0.15884E-03 593.92 89 -0.2041189721E+02 -0.2041202216E+02 0.12495E-03 588.19 90 -0.2041189720E+02 -0.2041199202E+02 0.94817E-04 582.39 91 -0.2041189719E+02 -0.2041196580E+02 0.68614E-04 576.53 92 -0.2041189718E+02 -0.2041194362E+02 0.46440E-04 570.65 93 -0.2041189722E+02 -0.2041192568E+02 0.28466E-04 564.77 94 -0.2041189724E+02 -0.2041191207E+02 0.14837E-04 558.92 95 -0.2041189723E+02 -0.2041190276E+02 0.55340E-05 553.13 96 -0.2041189728E+02 -0.2041189797E+02 0.69815E-06 547.44 97 -0.2041189732E+02 -0.2041189779E+02 0.46674E-06 541.85 98 -0.2041189735E+02 -0.2041190214E+02 0.47989E-05 536.41 99 -0.2041189737E+02 -0.2041191090E+02 0.13524E-04 531.13 100 -0.2041189746E+02 -0.2041192401E+02 0.26552E-04 526.04 101 -0.2041189749E+02 -0.2041194138E+02 0.43889E-04 521.15 102 -0.2041189751E+02 -0.2041196293E+02 0.65413E-04 516.49 103 -0.2041189763E+02 -0.2041198867E+02 0.91035E-04 512.08 104 -0.2041189764E+02 -0.2041201828E+02 0.12064E-03 507.93 105 -0.2041189766E+02 -0.2041205169E+02 0.15403E-03 504.06 106 -0.2041189767E+02 -0.2041208869E+02 0.19102E-03 500.47 107 -0.2041189768E+02 -0.2041212910E+02 0.23142E-03 497.19 108 -0.2041189770E+02 -0.2041217270E+02 0.27500E-03 494.22 109 -0.2041189771E+02 -0.2041221924E+02 0.32153E-03 491.58 110 -0.2041189773E+02 -0.2041226849E+02 0.37076E-03 489.26 111 -0.2041189775E+02 -0.2041232017E+02 0.42243E-03 487.27 112 -0.2041189776E+02 -0.2041237402E+02 0.47626E-03 485.62 113 -0.2041189778E+02 -0.2041242975E+02 0.53197E-03 484.31 114 -0.2041189780E+02 -0.2041248708E+02 0.58928E-03 483.33 115 -0.2041189782E+02 -0.2041254570E+02 0.64788E-03 482.69 116 -0.2041189784E+02 -0.2041260531E+02 0.70747E-03 482.39 117 -0.2041189786E+02 -0.2041266561E+02 0.76775E-03 482.41 118 -0.2041189788E+02 -0.2041272630E+02 0.82843E-03 482.75 119 -0.2041189790E+02 -0.2041278708E+02 0.88918E-03 483.41 120 -0.2041189791E+02 -0.2041284764E+02 0.94972E-03 484.37 121 -0.2041189793E+02 -0.2041290769E+02 0.10098E-02 485.63 122 -0.2041189795E+02 -0.2041296694E+02 0.10690E-02 487.17 123 -0.2041189797E+02 -0.2041302511E+02 0.11271E-02 488.98 124 -0.2041189798E+02 -0.2041308193E+02 0.11839E-02 491.05 125 -0.2041189800E+02 -0.2041313713E+02 0.12391E-02 493.36 126 -0.2041189801E+02 -0.2041319046E+02 0.12924E-02 495.91 127 -0.2041189803E+02 -0.2041324169E+02 0.13437E-02 498.66 128 -0.2041189804E+02 -0.2041329059E+02 0.13925E-02 501.61 129 -0.2041189805E+02 -0.2041333695E+02 0.14389E-02 504.74 130 -0.2041189807E+02 -0.2041338057E+02 0.14825E-02 508.04 131 -0.2041189808E+02 -0.2041342128E+02 0.15232E-02 511.48 132 -0.2041189809E+02 -0.2041345891E+02 0.15608E-02 515.04 133 -0.2041189809E+02 -0.2041349332E+02 0.15952E-02 518.72 134 -0.2041189810E+02 -0.2041352438E+02 0.16263E-02 522.48 135 -0.2041189811E+02 -0.2041355198E+02 0.16539E-02 526.32 136 -0.2041189811E+02 -0.2041357602E+02 0.16779E-02 530.22 137 -0.2041189811E+02 -0.2041359643E+02 0.16983E-02 534.15 138 -0.2041189812E+02 -0.2041361315E+02 0.17150E-02 538.10 139 -0.2041189812E+02 -0.2041362614E+02 0.17280E-02 542.05 140 -0.2041189812E+02 -0.2041363538E+02 0.17373E-02 545.98 141 -0.2041189812E+02 -0.2041364086E+02 0.17427E-02 549.89 142 -0.2041189812E+02 -0.2041364260E+02 0.17445E-02 553.75 143 -0.2041189811E+02 -0.2041364061E+02 0.17425E-02 557.55 144 -0.2041189811E+02 -0.2041363495E+02 0.17368E-02 561.27 145 -0.2041189811E+02 -0.2041362566E+02 0.17276E-02 564.89 146 -0.2041189810E+02 -0.2041361282E+02 0.17147E-02 568.42 147 -0.2041189810E+02 -0.2041359652E+02 0.16984E-02 571.83 148 -0.2041189809E+02 -0.2041357685E+02 0.16788E-02 575.10 149 -0.2041189808E+02 -0.2041355393E+02 0.16558E-02 578.24 150 -0.2041189807E+02 -0.2041352787E+02 0.16298E-02 581.23 151 -0.2041189806E+02 -0.2041349882E+02 0.16008E-02 584.06 152 -0.2041189806E+02 -0.2041346690E+02 0.15688E-02 586.71 153 -0.2041189805E+02 -0.2041343228E+02 0.15342E-02 589.20 154 -0.2041189804E+02 -0.2041339512E+02 0.14971E-02 591.49 155 -0.2041189802E+02 -0.2041335558E+02 0.14576E-02 593.60 156 -0.2041189801E+02 -0.2041331383E+02 0.14158E-02 595.52 157 -0.2041189800E+02 -0.2041327007E+02 0.13721E-02 597.23 158 -0.2041189799E+02 -0.2041322447E+02 0.13265E-02 598.75 159 -0.2041189798E+02 -0.2041317723E+02 0.12793E-02 600.05 160 -0.2041189797E+02 -0.2041312854E+02 0.12306E-02 601.15 161 -0.2041189795E+02 -0.2041307859E+02 0.11806E-02 602.05 162 -0.2041189794E+02 -0.2041302758E+02 0.11296E-02 602.73 163 -0.2041189793E+02 -0.2041297572E+02 0.10778E-02 603.20 164 -0.2041189792E+02 -0.2041292320E+02 0.10253E-02 603.47 165 -0.2041189791E+02 -0.2041287023E+02 0.97232E-03 603.54 166 -0.2041189789E+02 -0.2041281700E+02 0.91910E-03 603.40 167 -0.2041189788E+02 -0.2041276370E+02 0.86582E-03 603.06 168 -0.2041189787E+02 -0.2041271055E+02 0.81268E-03 602.53 169 -0.2041189786E+02 -0.2041265772E+02 0.75987E-03 601.81 170 -0.2041189785E+02 -0.2041260542E+02 0.70757E-03 600.90 171 -0.2041189783E+02 -0.2041255382E+02 0.65598E-03 599.82 172 -0.2041189782E+02 -0.2041250310E+02 0.60528E-03 598.56 173 -0.2041189781E+02 -0.2041245345E+02 0.55564E-03 597.14 174 -0.2041189780E+02 -0.2041240504E+02 0.50724E-03 595.56 175 -0.2041189779E+02 -0.2041235803E+02 0.46024E-03 593.82 176 -0.2041189778E+02 -0.2041231258E+02 0.41480E-03 591.95 177 -0.2041189778E+02 -0.2041226884E+02 0.37107E-03 589.94 178 -0.2041189777E+02 -0.2041222697E+02 0.32921E-03 587.80 179 -0.2041189776E+02 -0.2041218710E+02 0.28934E-03 585.55 180 -0.2041189775E+02 -0.2041214936E+02 0.25161E-03 583.20 181 -0.2041189774E+02 -0.2041211388E+02 0.21614E-03 580.74 182 -0.2041189774E+02 -0.2041208077E+02 0.18304E-03 578.20 183 -0.2041189773E+02 -0.2041205015E+02 0.15242E-03 575.58 184 -0.2041189773E+02 -0.2041202210E+02 0.12438E-03 572.90 185 -0.2041189772E+02 -0.2041199673E+02 0.99006E-04 570.15 186 -0.2041189772E+02 -0.2041197410E+02 0.76388E-04 567.36 187 -0.2041189771E+02 -0.2041195425E+02 0.56542E-04 564.53 188 -0.2041189776E+02 -0.2041193734E+02 0.39583E-04 561.68 189 -0.2041189778E+02 -0.2041192340E+02 0.25620E-04 558.81 190 -0.2041189780E+02 -0.2041191246E+02 0.14667E-04 555.93 191 -0.2041189779E+02 -0.2041190446E+02 0.66727E-05 553.06 192 -0.2041189784E+02 -0.2041189958E+02 0.17381E-05 550.20 193 -0.2041189789E+02 -0.2041189789E+02 0.24633E-08 547.36 194 -0.2041189792E+02 -0.2041189941E+02 0.14827E-05 544.56 195 -0.2041189795E+02 -0.2041190404E+02 0.60879E-05 541.80 196 -0.2041189798E+02 -0.2041191165E+02 0.13672E-04 539.09 197 -0.2041189808E+02 -0.2041192231E+02 0.24223E-04 536.45 198 -0.2041189811E+02 -0.2041193589E+02 0.37785E-04 533.87 199 -0.2041189812E+02 -0.2041195240E+02 0.54278E-04 531.37 200 -0.2041189817E+02 -0.2041197184E+02 0.73664E-04 528.96 201 -0.2041189824E+02 -0.2041199414E+02 0.95905E-04 526.64 202 -0.2041189825E+02 -0.2041201916E+02 0.12092E-03 524.42 203 -0.2041189825E+02 -0.2041204686E+02 0.14861E-03 522.31 204 -0.2041189826E+02 -0.2041207715E+02 0.17889E-03 520.32 205 -0.2041189826E+02 -0.2041210994E+02 0.21167E-03 518.44 206 -0.2041189827E+02 -0.2041214511E+02 0.24684E-03 516.69 207 -0.2041189828E+02 -0.2041218255E+02 0.28427E-03 515.07 208 -0.2041189829E+02 -0.2041222214E+02 0.32386E-03 513.58 209 -0.2041189830E+02 -0.2041226375E+02 0.36546E-03 512.24 210 -0.2041189831E+02 -0.2041230724E+02 0.40894E-03 511.03 211 -0.2041189832E+02 -0.2041235247E+02 0.45415E-03 509.97 212 -0.2041189833E+02 -0.2041239928E+02 0.50095E-03 509.07 213 -0.2041189834E+02 -0.2041244750E+02 0.54917E-03 508.31 214 -0.2041189835E+02 -0.2041249699E+02 0.59864E-03 507.70 215 -0.2041189836E+02 -0.2041254756E+02 0.64920E-03 507.25 216 -0.2041189837E+02 -0.2041259904E+02 0.70067E-03 506.95 217 -0.2041189838E+02 -0.2041265125E+02 0.75287E-03 506.81 218 -0.2041189839E+02 -0.2041270401E+02 0.80561E-03 506.82 219 -0.2041189841E+02 -0.2041275711E+02 0.85871E-03 506.98 220 -0.2041189842E+02 -0.2041281038E+02 0.91196E-03 507.29 221 -0.2041189843E+02 -0.2041286362E+02 0.96518E-03 507.75 222 -0.2041189845E+02 -0.2041291662E+02 0.10182E-02 508.36 223 -0.2041189846E+02 -0.2041296919E+02 0.10707E-02 509.11 224 -0.2041189847E+02 -0.2041302114E+02 0.11227E-02 509.99 225 -0.2041189849E+02 -0.2041307225E+02 0.11738E-02 511.02 226 -0.2041189850E+02 -0.2041312233E+02 0.12238E-02 512.17 227 -0.2041189851E+02 -0.2041317118E+02 0.12727E-02 513.45 228 -0.2041189852E+02 -0.2041321861E+02 0.13201E-02 514.85 229 -0.2041189854E+02 -0.2041326441E+02 0.13659E-02 516.36 230 -0.2041189855E+02 -0.2041330841E+02 0.14099E-02 517.98 231 -0.2041189856E+02 -0.2041335041E+02 0.14518E-02 519.70 232 -0.2041189857E+02 -0.2041339023E+02 0.14917E-02 521.51 233 -0.2041189859E+02 -0.2041342770E+02 0.15291E-02 523.41 234 -0.2041189860E+02 -0.2041346265E+02 0.15640E-02 525.39 235 -0.2041189861E+02 -0.2041349491E+02 0.15963E-02 527.43 236 -0.2041189862E+02 -0.2041352433E+02 0.16257E-02 529.54 237 -0.2041189863E+02 -0.2041355078E+02 0.16521E-02 531.70 238 -0.2041189864E+02 -0.2041357410E+02 0.16755E-02 533.90 239 -0.2041189865E+02 -0.2041359419E+02 0.16955E-02 536.14 240 -0.2041189865E+02 -0.2041361093E+02 0.17123E-02 538.40 241 -0.2041189866E+02 -0.2041362421E+02 0.17256E-02 540.68 242 -0.2041189867E+02 -0.2041363396E+02 0.17353E-02 542.96 243 -0.2041189867E+02 -0.2041364009E+02 0.17414E-02 545.25 244 -0.2041189867E+02 -0.2041364255E+02 0.17439E-02 547.52 245 -0.2041189868E+02 -0.2041364130E+02 0.17426E-02 549.76 246 -0.2041189868E+02 -0.2041363630E+02 0.17376E-02 551.98 247 -0.2041189868E+02 -0.2041362755E+02 0.17289E-02 554.16 248 -0.2041189868E+02 -0.2041361504E+02 0.17164E-02 556.29 249 -0.2041189868E+02 -0.2041359880E+02 0.17001E-02 558.36 250 -0.2041189868E+02 -0.2041357887E+02 0.16802E-02 560.36 251 -0.2041189867E+02 -0.2041355530E+02 0.16566E-02 562.29 252 -0.2041189867E+02 -0.2041352817E+02 0.16295E-02 564.13 253 -0.2041189866E+02 -0.2041349756E+02 0.15989E-02 565.89 254 -0.2041189866E+02 -0.2041346359E+02 0.15649E-02 567.54 255 -0.2041189865E+02 -0.2041342638E+02 0.15277E-02 569.10 256 -0.2041189864E+02 -0.2041338607E+02 0.14874E-02 570.54 257 -0.2041189863E+02 -0.2041334282E+02 0.14442E-02 571.86 258 -0.2041189862E+02 -0.2041329682E+02 0.13982E-02 573.06 259 -0.2041189861E+02 -0.2041324826E+02 0.13496E-02 574.14 260 -0.2041189860E+02 -0.2041319733E+02 0.12987E-02 575.08 261 -0.2041189858E+02 -0.2041314427E+02 0.12457E-02 575.89 262 -0.2041189857E+02 -0.2041308932E+02 0.11907E-02 576.56 263 -0.2041189855E+02 -0.2041303271E+02 0.11342E-02 577.09 264 -0.2041189854E+02 -0.2041297472E+02 0.10762E-02 577.47 265 -0.2041189852E+02 -0.2041291561E+02 0.10171E-02 577.72 266 -0.2041189850E+02 -0.2041285567E+02 0.95717E-03 577.82 267 -0.2041189848E+02 -0.2041279518E+02 0.89670E-03 577.77 268 -0.2041189847E+02 -0.2041273444E+02 0.83597E-03 577.59 269 -0.2041189845E+02 -0.2041267374E+02 0.77529E-03 577.26 270 -0.2041189843E+02 -0.2041261339E+02 0.71496E-03 576.80 271 -0.2041189841E+02 -0.2041255369E+02 0.65528E-03 576.20 272 -0.2041189839E+02 -0.2041249495E+02 0.59656E-03 575.47 273 -0.2041189837E+02 -0.2041243746E+02 0.53909E-03 574.61 274 -0.2041189836E+02 -0.2041238153E+02 0.48317E-03 573.63 275 -0.2041189834E+02 -0.2041232744E+02 0.42910E-03 572.54 276 -0.2041189832E+02 -0.2041227548E+02 0.37716E-03 571.33 277 -0.2041189831E+02 -0.2041222592E+02 0.32762E-03 570.02 278 -0.2041189829E+02 -0.2041217903E+02 0.28074E-03 568.61 279 -0.2041189828E+02 -0.2041213506E+02 0.23678E-03 567.11 280 -0.2041189826E+02 -0.2041209424E+02 0.19598E-03 565.54 281 -0.2041189825E+02 -0.2041205680E+02 0.15855E-03 563.89 282 -0.2041189824E+02 -0.2041202293E+02 0.12469E-03 562.17 283 -0.2041189823E+02 -0.2041199283E+02 0.94601E-04 560.40 284 -0.2041189822E+02 -0.2041196666E+02 0.68437E-04 558.59 285 -0.2041189821E+02 -0.2041194451E+02 0.46300E-04 556.74 286 -0.2041189825E+02 -0.2041192661E+02 0.28361E-04 554.86 287 -0.2041189827E+02 -0.2041191303E+02 0.14765E-04 552.96 288 -0.2041189826E+02 -0.2041190375E+02 0.54921E-05 551.06 289 -0.2041189831E+02 -0.2041189899E+02 0.68413E-06 549.16 290 -0.2041189835E+02 -0.2041189883E+02 0.47735E-06 547.28 291 -0.2041189837E+02 -0.2041190320E+02 0.48301E-05 545.41 292 -0.2041189840E+02 -0.2041191197E+02 0.13572E-04 543.58 293 -0.2041189849E+02 -0.2041192511E+02 0.26613E-04 541.79 294 -0.2041189852E+02 -0.2041194248E+02 0.43961E-04 540.05 295 -0.2041189854E+02 -0.2041196404E+02 0.65492E-04 538.36 296 -0.2041189866E+02 -0.2041198978E+02 0.91119E-04 536.74 297 -0.2041189868E+02 -0.2041201940E+02 0.12072E-03 535.20 298 -0.2041189869E+02 -0.2041205280E+02 0.15411E-03 533.73 299 -0.2041189870E+02 -0.2041208980E+02 0.19110E-03 532.36 300 -0.2041189871E+02 -0.2041213020E+02 0.23149E-03 531.08 301 -0.2041189873E+02 -0.2041217378E+02 0.27506E-03 529.89 302 -0.2041189874E+02 -0.2041222031E+02 0.32157E-03 528.81 303 -0.2041189876E+02 -0.2041226954E+02 0.37078E-03 527.84 304 -0.2041189878E+02 -0.2041232121E+02 0.42243E-03 526.99 305 -0.2041189879E+02 -0.2041237503E+02 0.47624E-03 526.25 306 -0.2041189881E+02 -0.2041243074E+02 0.53193E-03 525.63 307 -0.2041189883E+02 -0.2041248803E+02 0.58920E-03 525.13 308 -0.2041189885E+02 -0.2041254662E+02 0.64777E-03 524.75 309 -0.2041189887E+02 -0.2041260620E+02 0.70733E-03 524.50 310 -0.2041189889E+02 -0.2041266647E+02 0.76758E-03 524.38 311 -0.2041189891E+02 -0.2041272712E+02 0.82822E-03 524.37 312 -0.2041189893E+02 -0.2041278786E+02 0.88894E-03 524.49 313 -0.2041189894E+02 -0.2041284838E+02 0.94944E-03 524.73 314 -0.2041189896E+02 -0.2041290840E+02 0.10094E-02 525.09 315 -0.2041189898E+02 -0.2041296761E+02 0.10686E-02 525.56 316 -0.2041189900E+02 -0.2041302574E+02 0.11267E-02 526.15 317 -0.2041189901E+02 -0.2041308251E+02 0.11835E-02 526.84 318 -0.2041189903E+02 -0.2041313767E+02 0.12386E-02 527.64 319 -0.2041189904E+02 -0.2041319097E+02 0.12919E-02 528.55 320 -0.2041189906E+02 -0.2041324215E+02 0.13431E-02 529.54 321 -0.2041189907E+02 -0.2041329101E+02 0.13919E-02 530.63 322 -0.2041189908E+02 -0.2041333733E+02 0.14382E-02 531.80 323 -0.2041189910E+02 -0.2041338092E+02 0.14818E-02 533.04 324 -0.2041189911E+02 -0.2041342159E+02 0.15225E-02 534.36 325 -0.2041189912E+02 -0.2041345918E+02 0.15601E-02 535.74 326 -0.2041189912E+02 -0.2041349356E+02 0.15944E-02 537.19 327 -0.2041189913E+02 -0.2041352458E+02 0.16254E-02 538.68 328 -0.2041189914E+02 -0.2041355214E+02 0.16530E-02 540.21 329 -0.2041189914E+02 -0.2041357615E+02 0.16770E-02 541.79 330 -0.2041189914E+02 -0.2041359654E+02 0.16974E-02 543.39 331 -0.2041189915E+02 -0.2041361323E+02 0.17141E-02 545.01 332 -0.2041189915E+02 -0.2041362620E+02 0.17271E-02 546.65 333 -0.2041189915E+02 -0.2041363542E+02 0.17363E-02 548.30 334 -0.2041189915E+02 -0.2041364088E+02 0.17417E-02 549.94 335 -0.2041189915E+02 -0.2041364260E+02 0.17435E-02 551.58 336 -0.2041189914E+02 -0.2041364060E+02 0.17415E-02 553.21 337 -0.2041189914E+02 -0.2041363492E+02 0.17358E-02 554.82 338 -0.2041189914E+02 -0.2041362562E+02 0.17265E-02 556.40 339 -0.2041189913E+02 -0.2041361278E+02 0.17136E-02 557.94 340 -0.2041189913E+02 -0.2041359647E+02 0.16973E-02 559.45 341 -0.2041189912E+02 -0.2041357680E+02 0.16777E-02 560.91 342 -0.2041189911E+02 -0.2041355387E+02 0.16548E-02 562.33 343 -0.2041189910E+02 -0.2041352782E+02 0.16287E-02 563.68 344 -0.2041189909E+02 -0.2041349876E+02 0.15997E-02 564.98 345 -0.2041189909E+02 -0.2041346686E+02 0.15678E-02 566.20 346 -0.2041189908E+02 -0.2041343225E+02 0.15332E-02 567.36 347 -0.2041189907E+02 -0.2041339509E+02 0.14960E-02 568.45 348 -0.2041189905E+02 -0.2041335556E+02 0.14565E-02 569.45 349 -0.2041189904E+02 -0.2041331384E+02 0.14148E-02 570.38 350 -0.2041189903E+02 -0.2041327009E+02 0.13711E-02 571.22 351 -0.2041189902E+02 -0.2041322451E+02 0.13255E-02 571.98 352 -0.2041189901E+02 -0.2041317729E+02 0.12783E-02 572.64 353 -0.2041189900E+02 -0.2041312862E+02 0.12296E-02 573.22 354 -0.2041189898E+02 -0.2041307869E+02 0.11797E-02 573.70 355 -0.2041189897E+02 -0.2041302771E+02 0.11287E-02 574.10 356 -0.2041189896E+02 -0.2041297587E+02 0.10769E-02 574.39 357 -0.2041189895E+02 -0.2041292338E+02 0.10244E-02 574.60 358 -0.2041189894E+02 -0.2041287044E+02 0.97150E-03 574.70 359 -0.2041189892E+02 -0.2041281724E+02 0.91831E-03 574.72 360 -0.2041189891E+02 -0.2041276398E+02 0.86506E-03 574.64 361 -0.2041189890E+02 -0.2041271085E+02 0.81195E-03 574.47 362 -0.2041189889E+02 -0.2041265806E+02 0.75917E-03 574.21 363 -0.2041189888E+02 -0.2041260579E+02 0.70691E-03 573.86 364 -0.2041189886E+02 -0.2041255422E+02 0.65535E-03 573.42 365 -0.2041189885E+02 -0.2041250354E+02 0.60468E-03 572.89 366 -0.2041189884E+02 -0.2041245392E+02 0.55508E-03 572.29 367 -0.2041189883E+02 -0.2041240554E+02 0.50671E-03 571.60 368 -0.2041189882E+02 -0.2041235856E+02 0.45974E-03 570.84 369 -0.2041189881E+02 -0.2041231315E+02 0.41433E-03 570.01 370 -0.2041189881E+02 -0.2041226945E+02 0.37064E-03 569.10 371 -0.2041189880E+02 -0.2041222761E+02 0.32881E-03 568.13 372 -0.2041189879E+02 -0.2041218777E+02 0.28898E-03 567.09 373 -0.2041189878E+02 -0.2041215006E+02 0.25128E-03 566.00 374 -0.2041189877E+02 -0.2041211461E+02 0.21584E-03 564.86 375 -0.2041189877E+02 -0.2041208154E+02 0.18277E-03 563.66 376 -0.2041189876E+02 -0.2041205094E+02 0.15218E-03 562.42 377 -0.2041189876E+02 -0.2041202292E+02 0.12417E-03 561.14 378 -0.2041189875E+02 -0.2041199758E+02 0.98825E-04 559.82 379 -0.2041189875E+02 -0.2041197498E+02 0.76234E-04 558.48 380 -0.2041189874E+02 -0.2041195515E+02 0.56414E-04 557.11 381 -0.2041189879E+02 -0.2041193827E+02 0.39479E-04 555.71 382 -0.2041189881E+02 -0.2041192435E+02 0.25540E-04 554.30 383 -0.2041189883E+02 -0.2041191344E+02 0.14609E-04 552.88 384 -0.2041189882E+02 -0.2041190545E+02 0.66351E-05 551.46 385 -0.2041189887E+02 -0.2041190059E+02 0.17199E-05 550.03 386 -0.2041189892E+02 -0.2041189892E+02 0.18108E-08 548.61 387 -0.2041189895E+02 -0.2041190045E+02 0.14973E-05 547.20 388 -0.2041189898E+02 -0.2041190509E+02 0.61150E-05 545.80 389 -0.2041189901E+02 -0.2041191272E+02 0.13709E-04 544.43 390 -0.2041189911E+02 -0.2041192338E+02 0.24268E-04 543.07 391 -0.2041189914E+02 -0.2041193697E+02 0.37837E-04 541.75 392 -0.2041189915E+02 -0.2041195349E+02 0.54336E-04 540.46 393 -0.2041189920E+02 -0.2041197293E+02 0.73725E-04 539.20 394 -0.2041189927E+02 -0.2041199524E+02 0.95968E-04 537.99 395 -0.2041189928E+02 -0.2041202026E+02 0.12098E-03 536.83 396 -0.2041189928E+02 -0.2041204795E+02 0.14867E-03 535.71 397 -0.2041189929E+02 -0.2041207824E+02 0.17895E-03 534.65 398 -0.2041189929E+02 -0.2041211102E+02 0.21172E-03 533.64 399 -0.2041189930E+02 -0.2041214618E+02 0.24688E-03 532.70 400 -0.2041189931E+02 -0.2041218361E+02 0.28430E-03 531.82 401 -0.2041189932E+02 -0.2041222319E+02 0.32388E-03 531.01 402 -0.2041189933E+02 -0.2041226479E+02 0.36547E-03 530.26 403 -0.2041189934E+02 -0.2041230827E+02 0.40893E-03 529.59 404 -0.2041189935E+02 -0.2041235348E+02 0.45413E-03 528.99 405 -0.2041189936E+02 -0.2041240026E+02 0.50091E-03 528.46 406 -0.2041189937E+02 -0.2041244847E+02 0.54910E-03 528.02 407 -0.2041189938E+02 -0.2041249793E+02 0.59856E-03 527.65 408 -0.2041189939E+02 -0.2041254848E+02 0.64909E-03 527.36 409 -0.2041189940E+02 -0.2041259994E+02 0.70054E-03 527.16 410 -0.2041189941E+02 -0.2041265212E+02 0.75271E-03 527.03 411 -0.2041189942E+02 -0.2041270485E+02 0.80543E-03 526.98 412 -0.2041189944E+02 -0.2041275793E+02 0.85849E-03 527.02 413 -0.2041189945E+02 -0.2041281117E+02 0.91172E-03 527.14 414 -0.2041189946E+02 -0.2041286437E+02 0.96491E-03 527.34 415 -0.2041189948E+02 -0.2041291734E+02 0.10179E-02 527.61 416 -0.2041189949E+02 -0.2041296988E+02 0.10704E-02 527.97 417 -0.2041189950E+02 -0.2041302179E+02 0.11223E-02 528.40 418 -0.2041189952E+02 -0.2041307287E+02 0.11734E-02 528.91 419 -0.2041189953E+02 -0.2041312291E+02 0.12234E-02 529.49 420 -0.2041189954E+02 -0.2041317173E+02 0.12722E-02 530.14 421 -0.2041189955E+02 -0.2041321912E+02 0.13196E-02 530.86 422 -0.2041189957E+02 -0.2041326489E+02 0.13653E-02 531.64 423 -0.2041189958E+02 -0.2041330885E+02 0.14093E-02 532.49 424 -0.2041189959E+02 -0.2041335082E+02 0.14512E-02 533.39 425 -0.2041189960E+02 -0.2041339060E+02 0.14910E-02 534.35 426 -0.2041189962E+02 -0.2041342803E+02 0.15284E-02 535.36 427 -0.2041189963E+02 -0.2041346294E+02 0.15633E-02 536.42 428 -0.2041189964E+02 -0.2041349517E+02 0.15955E-02 537.51 429 -0.2041189965E+02 -0.2041352456E+02 0.16249E-02 538.65 430 -0.2041189966E+02 -0.2041355097E+02 0.16513E-02 539.82 431 -0.2041189967E+02 -0.2041357427E+02 0.16746E-02 541.02 432 -0.2041189968E+02 -0.2041359433E+02 0.16946E-02 542.24 433 -0.2041189968E+02 -0.2041361103E+02 0.17113E-02 543.48 434 -0.2041189969E+02 -0.2041362429E+02 0.17246E-02 544.74 435 -0.2041189970E+02 -0.2041363401E+02 0.17343E-02 546.00 436 -0.2041189970E+02 -0.2041364012E+02 0.17404E-02 547.27 437 -0.2041189970E+02 -0.2041364255E+02 0.17428E-02 548.53 438 -0.2041189971E+02 -0.2041364128E+02 0.17416E-02 549.79 439 -0.2041189971E+02 -0.2041363626E+02 0.17366E-02 551.03 440 -0.2041189971E+02 -0.2041362749E+02 0.17278E-02 552.25 441 -0.2041189971E+02 -0.2041361497E+02 0.17153E-02 553.46 442 -0.2041189971E+02 -0.2041359872E+02 0.16990E-02 554.63 443 -0.2041189971E+02 -0.2041357878E+02 0.16791E-02 555.77 444 -0.2041189970E+02 -0.2041355521E+02 0.16555E-02 556.87 445 -0.2041189970E+02 -0.2041352807E+02 0.16284E-02 557.92 446 -0.2041189969E+02 -0.2041349746E+02 0.15978E-02 558.93 447 -0.2041189969E+02 -0.2041346349E+02 0.15638E-02 559.89 448 -0.2041189968E+02 -0.2041342628E+02 0.15266E-02 560.79 449 -0.2041189967E+02 -0.2041338598E+02 0.14863E-02 561.63 450 -0.2041189966E+02 -0.2041334275E+02 0.14431E-02 562.41 451 -0.2041189965E+02 -0.2041329676E+02 0.13971E-02 563.11 452 -0.2041189964E+02 -0.2041324821E+02 0.13486E-02 563.75 453 -0.2041189963E+02 -0.2041319730E+02 0.12977E-02 564.32 454 -0.2041189961E+02 -0.2041314426E+02 0.12447E-02 564.80 455 -0.2041189960E+02 -0.2041308933E+02 0.11897E-02 565.21 456 -0.2041189958E+02 -0.2041303275E+02 0.11332E-02 565.54 457 -0.2041189957E+02 -0.2041297479E+02 0.10752E-02 565.79 458 -0.2041189955E+02 -0.2041291571E+02 0.10162E-02 565.95 459 -0.2041189953E+02 -0.2041285580E+02 0.95627E-03 566.04 460 -0.2041189951E+02 -0.2041279535E+02 0.89583E-03 566.04 461 -0.2041189950E+02 -0.2041273464E+02 0.83514E-03 565.95 462 -0.2041189948E+02 -0.2041267398E+02 0.77450E-03 565.79 463 -0.2041189946E+02 -0.2041261367E+02 0.71421E-03 565.54 464 -0.2041189944E+02 -0.2041255401E+02 0.65457E-03 565.21 465 -0.2041189942E+02 -0.2041249531E+02 0.59589E-03 564.81 466 -0.2041189940E+02 -0.2041243787E+02 0.53846E-03 564.33 467 -0.2041189939E+02 -0.2041238198E+02 0.48259E-03 563.77 468 -0.2041189937E+02 -0.2041232793E+02 0.42856E-03 563.15 469 -0.2041189935E+02 -0.2041227602E+02 0.37666E-03 562.45 470 -0.2041189934E+02 -0.2041222650E+02 0.32717E-03 561.70 471 -0.2041189932E+02 -0.2041217965E+02 0.28033E-03 560.88 472 -0.2041189931E+02 -0.2041213572E+02 0.23642E-03 560.01 473 -0.2041189929E+02 -0.2041209495E+02 0.19566E-03 559.09 474 -0.2041189928E+02 -0.2041205755E+02 0.15827E-03 558.13 475 -0.2041189927E+02 -0.2041202372E+02 0.12445E-03 557.12 476 -0.2041189926E+02 -0.2041199366E+02 0.94399E-04 556.08 477 -0.2041189925E+02 -0.2041196752E+02 0.68272E-04 555.00 478 -0.2041189924E+02 -0.2041194541E+02 0.46169E-04 553.90 479 -0.2041189928E+02 -0.2041192754E+02 0.28263E-04 552.79 480 -0.2041189930E+02 -0.2041191400E+02 0.14698E-04 551.66 481 -0.2041189929E+02 -0.2041190474E+02 0.54537E-05 550.52 482 -0.2041189934E+02 -0.2041190001E+02 0.67146E-06 549.38 483 -0.2041189938E+02 -0.2041189987E+02 0.48710E-06 548.25 484 -0.2041189940E+02 -0.2041190426E+02 0.48584E-05 547.12 485 -0.2041189943E+02 -0.2041191305E+02 0.13614E-04 546.02 486 -0.2041189952E+02 -0.2041192619E+02 0.26668E-04 544.93 487 -0.2041189955E+02 -0.2041194357E+02 0.44024E-04 543.87 488 -0.2041189957E+02 -0.2041196513E+02 0.65560E-04 542.84 489 -0.2041189969E+02 -0.2041199088E+02 0.91189E-04 541.85 490 -0.2041189971E+02 -0.2041202050E+02 0.12079E-03 540.90 491 -0.2041189972E+02 -0.2041205389E+02 0.15418E-03 540.00 492 -0.2041189973E+02 -0.2041209089E+02 0.19116E-03 539.15 493 -0.2041189974E+02 -0.2041213128E+02 0.23154E-03 538.36 494 -0.2041189976E+02 -0.2041217485E+02 0.27509E-03 537.62 495 -0.2041189977E+02 -0.2041222136E+02 0.32159E-03 536.94 496 -0.2041189979E+02 -0.2041227057E+02 0.37078E-03 536.34 497 -0.2041189981E+02 -0.2041232221E+02 0.42241E-03 535.79 498 -0.2041189982E+02 -0.2041237602E+02 0.47619E-03 535.32 499 -0.2041189984E+02 -0.2041243170E+02 0.53185E-03 534.92 500 -0.2041189986E+02 -0.2041248896E+02 0.58910E-03 534.60 *** arived at the Maximum iteration. terminated. Elapsed time of simulation was 120.944 fs ============= summary of results ================= final position of ions: 1 S ( -0.00003 0.00022 -1.74537 ) - atomic mass= 31.972 2 S ( 0.00003 -0.00022 1.74537 ) - atomic mass= 31.972 G.C. ( 0.00000 0.00000 0.00000 ) C.O.M. ( 0.00000 0.00000 0.00000 ) final velocity of ions: 1 S ( 0.00000 0.00000 -0.00010 ) 2 S ( 0.00000 0.00000 0.00010 ) G.C. ( 0.00000 0.00000 0.00000 ) C.O.M. ( 0.00000 0.00000 0.00000 ) == Summary Of Results == number of electrons: spin up= 7.00000 down= 5.00000 (real space) total energy : -0.2041248896E+02 ( -0.10206E+02/ion) total orbital energy: -0.4969399250E+01 ( -0.41412E+00/electron) hartree energy : 0.1699511472E+02 ( 0.14163E+01/electron) exc-corr energy : -0.4338365042E+01 ( -0.36153E+00/electron) ion-ion energy : 0.2188170914E+00 ( 0.10941E+00/ion) kinetic (planewave) : 0.7584364385E+01 ( 0.63203E+00/electron) V_local (planewave) : -0.4550760392E+02 ( -0.37923E+01/electron) V_nl (planewave) : 0.4635183805E+01 ( 0.38627E+00/electron) V_coul (planewave) : 0.3399022944E+02 ( 0.28325E+01/electron) V_xc (planewave) : -0.5671572962E+01 ( -0.47263E+00/electron) Virial Coefficient : -0.1655216310E+01 orbital energies: -0.1962285E+00 ( -5.340eV) -0.1962285E+00 ( -5.340eV) -0.3342998E+00 ( -9.097eV) -0.3042598E+00 ( -8.279eV) -0.3342998E+00 ( -9.097eV) -0.3042598E+00 ( -8.279eV) -0.3592474E+00 ( -9.776eV) -0.3357308E+00 ( -9.136eV) -0.5588872E+00 ( -15.208eV) -0.5201014E+00 ( -14.153eV) -0.7740222E+00 ( -21.062eV) -0.7518341E+00 ( -20.459eV) Kinetic energy (ion) : 0.5890999095E-03 ( 0.29455E-03/ion) Temperature : 534.6 K (ion) : 0.0 K (c.o.m.) Eaverage : -0.2041274664E+02 Evariance : 0.3821111818E-06 Cv - f*kb/(2*nion) : 0.2110980482E-06 Translation force removed: ( 0.00000 0.00000 0.00000) ============= Ion Gradients ================= Ion Forces: 1 S ( 0.00000 0.00000 -0.02882 ) 2 S ( 0.00000 0.00000 0.02882 ) C.O.M. ( 0.00000 0.00000 0.00000 ) =============================================== |F| = 0.407564E-01 |F|/nion = 0.203782E-01 max|Fatom|= 0.288191E-01 ( 1.482eV/Angstrom) Total PSPW energy : -0.2041248896E+02 === Spin Contamination === = 2.00000000000000 = 2.00289928697927 == Center of Charge == spin up ( 0.0000, 0.0000, 0.0000 ) spin down ( 0.0000, 0.0000, 0.0000 ) total ( 0.0000, 0.0000, 0.0000 ) ionic ( 0.0000, 0.0000, 0.0000 ) crystal ( 0.0000, 0.0000, 0.0000 ) == Crystal Dipole == mu = ( 0.0000, 0.0000, 0.0000 ) au |mu| = 0.0000 au, 0.0000 Debye == Molecular Dipole wrt Center of Mass == mu = ( 0.0000, 0.0000, 0.0000 ) au |mu| = 0.0000 au, 0.0000 Debye output psi filename:./perm/s2-bomd.movecs ===== Calculating pair distribution function ==== frames used: 500 natoms: 2 delr: 1.000000000000000E-002 volume: 8000.00000000000 maxbin: 999 MOTION filename: ./perm/s2-bomd.ion_motion g(r) filename: ./perm/PAIR_DISTRIBUTION ===== Calculating radial distribution functions ==== frames used: 500 natoms: 2 nkatoms: 1 atoms: S delr: 1.000000000000000E-002 volume: 8000.00000000000 a1: 20.000 0.000 0.000 a2: 0.000 20.000 0.000 a3: 0.000 0.000 20.000 maxbin: 999 MOTION filename: ./perm/s2-bomd.ion_motion creating gr filename:./perm/SS.gr creating hist filename:./perm/SS.hist ===== Calculating velocity auto correlation ==== frames used: 500 natoms: 2 volume: 8000.00000000000 time interval (au) : 10.0000000000000 Max time (au) : 4990.00000000000 frequency interval (cm-1): 275.976477774143 Max frequency (cm-1) : 68994.1194435358 MOTION filename: ./perm/s2-bomd.ion_motion filename: ./perm/VEL_CORRELATION s(k) filename: ./perm/POWER_SPECTRUM == Timing == cputime in seconds prologue : 0.294209E+00 main loop : 0.139011E+04 epilogue : 0.215233E+00 total : 0.139062E+04 cputime/step: 0.119057E+00 ( 11676 evalulations, 4901 linesearches) Time spent doing total step FFTs : 0.683384E+02 0.585289E-02 dot products : 0.391666E+02 0.335445E-02 geodesic : 0.127394E+03 0.109107E-01 ffm_dgemm : 0.119895E+02 0.102685E-02 fmf_dgemm : 0.115739E+03 0.991256E-02 m_diagonalize : 0.311176E+00 0.266509E-04 - m_tredq : 0.000000E+00 0.000000E+00 - m_getdiags : 0.000000E+00 0.000000E+00 - m_tqliq : 0.000000E+00 0.000000E+00 - m_eigsrt : 0.000000E+00 0.000000E+00 exchange correlation : 0.214557E+03 0.183759E-01 local pseudopotentials : 0.687694E+00 0.588980E-04 non-local pseudopotentials : 0.556925E+02 0.476983E-02 hartree potentials : 0.408362E+01 0.349745E-03 ion-ion interaction : 0.102148E+01 0.874852E-04 structure factors : 0.235273E+01 0.201502E-03 phase factors : 0.196814E-02 0.168563E-06 masking and packing : 0.905270E+02 0.775326E-02 queue fft : 0.738079E+03 0.632134E-01 queue fft (serial) : 0.352001E+03 0.301474E-01 queue fft (message passing): 0.361980E+03 0.310020E-01 HFX potential : 0.000000E+00 0.000000E+00 qmmm LJ : 0.000000E+00 0.000000E+00 qmmm residual Q : 0.000000E+00 0.000000E+00 >>> JOB COMPLETED AT Tue Aug 10 16:05:09 2010 <<< Task times cpu: 1370.4s wall: 1390.6s NWChem Input Module ------------------- Summary of allocated global arrays ----------------------------------- No active global arrays GA Statistics for process 0 ------------------------------ create destroy get put acc scatter gather read&inc calls: 0 0 0 0 0 0 0 0 number of processes/call 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 bytes total: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 bytes remote: 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 0.00e+00 Max memory consumed for GA by this process: 0 bytes ACKNOWLEDGEMENT --------------- Please use the following acknowledgement where appropriate for results obtained with NWChem: High Performance Computational Chemistry Group, "NWChem, A Computational Chemistry Package for Parallel Computers, Version 5.1.1" (2008), Pacific Northwest National Laboratory, Richland, Washington 99352-0999, USA. CITATION -------- Please use the following citation when publishing results obtained with NWChem: E. J. Bylaska, W. A. de Jong, N. Govind, K. Kowalski, T. P. Straatsma, M. Valiev, H. J. J. Van Dam, D. Wang, E. Apra, T. L. Windus, J. Hammond, J. Autschbach, P. Nichols, S. Hirata, M. T. Hackler, Y. Zhao, P.-D. Fan, R. J. Harrison, M. Dupuis, D. M. A. Smith, J. Nieplocha, V. Tipparaju, M. Krishnan, A. Vazquez-Mayagoitia, L. Jensen, M. Swart, Q. Wu, T. Van Voorhis, A. A. Auer, M. Nooijen, L. D. Crosby, E. Brown, G. Cisneros, G. I. Fann, H. Fruchtl, J. Garza, K. Hirao, R. Kendall, J. A. Nichols, K. Tsemekhman, K. Wolinski, J. Anchell, D. Bernholdt, P. Borowski, T. Clark, D. Clerc, H. Dachsel, M. Deegan, K. Dyall, D. Elwood, E. Glendening, M. Gutowski, A. Hess, J. Jaffe, B. Johnson, J. Ju, R. Kobayashi, R. Kutteh, Z. Lin, R. Littlefield, X. Long, B. Meng, T. Nakajima, S. Niu, L. Pollack, M. Rosing, G. Sandrone, M. Stave, H. Taylor, G. Thomas, J. van Lenthe, A. Wong, and Z. Zhang, "NWChem, A Computational Chemistry Package for Parallel Computers, Version 5.1.1" (2008), Pacific Northwest National Laboratory, Richland, Washington 99352-0999, USA. Total times cpu: 1393.8s wall: 1414.9s MA_summarize_allocated_blocks: starting scan ... MA_summarize_allocated_blocks: scan completed: 0 heap blocks, 0 stack blocks MA usage statistics: allocation statistics: heap stack ---- ----- current number of blocks 0 0 maximum number of blocks 235 11 current total bytes 0 0 maximum total bytes 49804024 4205376 maximum total K-bytes 49805 4206 maximum total M-bytes 50 5