Page 1 of 1

fatal internal in brackting...

Posted: Mon Oct 22, 2007 4:16 pm
by Danny
I'm running some tests to time my vasp executables(different ways of compiling version 4.6.31) and see how they respond to different combinations of processors and threads.
All goes well for all combinations I try except 1 cpu 4threads and 1 cpu 8threads(1 cpu 2 threads, works nicely and N cpu 4&8 threads, N>1 also does alright) in this case the job crashes with the error:

Code: Select all

ZBRENT: fatal internal in brackting
             system-shutdown; contact gK immediatly
With all the versions I tried already only for the combinations stated above the job goes bad.
(all vasp parameters are exactly the same)

Are there any suggestions I could implement in my compilation to solve this problem or are they not related to the vasp compilation completely?

thx
Danny

fatal internal in brackting...

Posted: Wed Oct 24, 2007 9:39 am
by admin
the error messsage says that the geometry optimization cannot find a reasonable next geometry by BRENT's algorithm (root bisectioning). This has nothing to do with the hardware or the code itself, but is due to the ionic optimization "history" of your system.
please check
1)( best to do in XDATCAR ): did the system converge into an unreasonable geometry? if yes: set POTIM to a smaller value (default=0.5) and/or change the optimization algorithm (IBRION)
2) check EDIFFG: a calculation can be regarded as converged if the maximum remaining forces are not larger that 0.01eV/Ã…. (EDIFFG = -0.01). if a much higher precision is afforded, you may have to set IBRION=1; ADDGRID=.True.
3) please in any case also check if each ionic step is converged electronically (ie the forces are calculated after the EDIFF-criterium is matched, not because the electronic convergence is not reached after 60 el. steps and vasp therefore breaks the electronic scf iterations for this geometry. If the latter is the case, the calculated forces (and therefore the next ionic configuration) are unreasonable.

fatal internal in brackting...

Posted: Mon Oct 29, 2007 3:16 pm
by Danny
[quote="admin"]the error messsage says that the geometry optimization cannot find a reasonable next geometry by BRENT's algorithm (root bisectioning). This has nothing to do with the hardware or the code itself, but is due to the ionic optimization "history" of your system.
please check
[/quote]
if this has nothing to do with the hardware nor the code then doing exactly the same calculation(i.e. same input files) should give the same error? Am I correct?
If so, then why does the same calculation run on 16 CPU's run perfectly, and when I run it on 16CPU's each running 4 threads or 8CPU's running 2 threads each, the calculations gives the above error?
All of them should fail if it were only a geometry problem.

Danny

fatal internal in brackting...

Posted: Wed Oct 31, 2007 2:41 pm
by Danny
output when the job is run on 1cpu, 1thread/cpu

Code: Select all

step_1
job started on Wednesday 24 October 2007, at 15:28:10
 running on    1 nodes
 distr:  one band on    1 nodes,    1 groups
 vasp.4.6.31 08Feb07 complex 
 POSCAR found :  1 types and    2 ions
 LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro
 POSCAR, INCAR and KPOINTS ok, starting setup
 FFT: planning ... 1
 reading WAVECAR
 entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.541995558798E+01    0.54200E+01   -0.32487E+03  7280   0.801E+02
DAV:   2    -0.108298552793E+02   -0.16250E+02   -0.15628E+02 11120   0.778E+01
DAV:   3    -0.110322582370E+02   -0.20240E+00   -0.20240E+00  9760   0.111E+01
DAV:   4    -0.110330766749E+02   -0.81844E-03   -0.81844E-03 11392   0.747E-01
DAV:   5    -0.110330776903E+02   -0.10154E-05   -0.10154E-05  9992   0.213E-02    0.301E+00
DAV:   6    -0.109131780101E+02    0.11990E+00   -0.84721E-02  8776   0.151E+00    0.179E+00
DAV:   7    -0.108573305462E+02    0.55847E-01   -0.16141E-01  9248   0.218E+00    0.158E-01
DAV:   8    -0.108582608010E+02   -0.93025E-03   -0.42091E-03  8280   0.480E-01    0.272E-02
DAV:   9    -0.108582986809E+02   -0.37880E-04   -0.10941E-04 10544   0.838E-02
   1 F= -.10858299E+02 E0= -.10858299E+02  d E =-.108583E+02
 curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
 trial: gam= 0.00000 g(F)=  0.320E-45 g(S)=  0.161E-01 ort = 0.000E+00 (trialstep = 0.100E+01)
 search vector abs. value=  0.161E-01
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108674017214E+02   -0.10867E+02   -0.24938E-01  9040   0.283E+00    0.316E-01
DAV:   2    -0.108630875918E+02    0.43141E-02   -0.38499E-03 10152   0.427E-01    0.180E-01
DAV:   3    -0.108608719497E+02    0.22156E-02   -0.47378E-03  8664   0.347E-01    0.236E-02
DAV:   4    -0.108608612291E+02    0.10721E-04   -0.55128E-05  6192   0.677E-02
   2 F= -.10860861E+02 E0= -.10860861E+02  d E =-.256255E-02
 trial-energy change:   -0.002563  1 .order   -0.002850   -0.016107    0.010406
 step:   0.6075(harm=  0.6075)  dis= 0.00000  next Energy=   -10.863191 (dE=-0.489E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108640013077E+02   -0.10864E+02   -0.39184E-02  9000   0.113E+00    0.126E-01
DAV:   2    -0.108633291668E+02    0.67214E-03   -0.59649E-04 10176   0.168E-01    0.718E-02
DAV:   3    -0.108629924546E+02    0.33671E-03   -0.73965E-04  8616   0.137E-01    0.939E-03
DAV:   4    -0.108629944707E+02   -0.20161E-05   -0.81055E-06  3808   0.267E-02
   3 F= -.10862994E+02 E0= -.10862994E+02  d E =-.469579E-02
 curvature:  -0.30 expect dE=-0.306E-05 dE for cont linesearch -0.306E-05
 ZBRENT: interpolating
 opt :   0.5927  next Energy=   -10.862997 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629988933E+02   -0.10863E+02   -0.51208E-05  3752   0.421E-02    0.441E-03
DAV:   2    -0.108629982717E+02    0.62153E-06   -0.25033E-06  3664   0.125E-02
   4 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469959E-02
 curvature:  -0.30 expect dE=-0.168E-06 dE for cont linesearch -0.168E-06
 ZBRENT: interpolating
 opt :   0.5881  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629978226E+02   -0.10863E+02   -0.34251E-06  3648   0.124E-02    0.181E-03
DAV:   2    -0.108629977824E+02    0.40198E-07   -0.35836E-07  3640   0.453E-03
   5 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469910E-02
 curvature:  -0.39 expect dE=-0.226E-06 dE for cont linesearch -0.226E-06
 ZBRENT: interpolating
 opt :   0.5904  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629981381E+02   -0.10863E+02   -0.19219E-06  3640   0.791E-03    0.636E-04
DAV:   2    -0.108629981287E+02    0.93158E-08   -0.80961E-08  3640   0.245E-03
   6 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469945E-02
 curvature:  -0.19 expect dE=-0.117E-06 dE for cont linesearch -0.117E-06
 ZBRENT: interpolating
 opt :   0.5916  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629982012E+02   -0.10863E+02   -0.26870E-07  3640   0.328E-03    0.426E-04
DAV:   2    -0.108629981785E+02    0.22707E-07   -0.19972E-08  3640   0.111E-03
   7 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469950E-02
 curvature:  -0.09 expect dE=-0.229E-07 dE for cont linesearch -0.229E-07
 trial: gam=-0.00406 g(F)=  0.160E-45 g(S)=  0.268E-06 ort = 0.657E-04 (trialstep = 0.918E+00)
 search vector abs. value=  0.104E-09
 reached required accuracy - stopping structural energy minimisation
Total resource usage:
 Tasks  : 1
 Memory : 23.00 MB
 CPU    : 1.000000

Average per task :
 Memory : 23.00 MB
 CPU    : 1.00
job completed on Wednesday 24 October 2007, at 15:46:23

fatal internal in brackting...

Posted: Wed Oct 31, 2007 2:43 pm
by Danny
the job when it is run on 4cpu's, 1thread/cpu

Code: Select all

step_2
job started on Wednesday 24 October 2007, at 19:02:47
 running on    4 nodes
 distr:  one band on    1 nodes,    4 groups
 vasp.4.6.31 08Feb07 complex 
 POSCAR found :  1 types and    2 ions
 LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro
 POSCAR, INCAR and KPOINTS ok, starting setup
 FFT: planning ... 1
 reading WAVECAR
 entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.541995558798E+01    0.54200E+01   -0.32487E+03  7280   0.801E+02
DAV:   2    -0.108298552793E+02   -0.16250E+02   -0.15628E+02 11120   0.778E+01
DAV:   3    -0.110322582371E+02   -0.20240E+00   -0.20240E+00  9760   0.111E+01
DAV:   4    -0.110330766749E+02   -0.81844E-03   -0.81844E-03 11392   0.747E-01
DAV:   5    -0.110330776903E+02   -0.10154E-05   -0.10154E-05  9992   0.213E-02    0.301E+00
DAV:   6    -0.109131780101E+02    0.11990E+00   -0.84721E-02  8776   0.151E+00    0.179E+00
DAV:   7    -0.108573305462E+02    0.55847E-01   -0.16141E-01  9248   0.218E+00    0.158E-01
DAV:   8    -0.108582608010E+02   -0.93025E-03   -0.42091E-03  8280   0.480E-01    0.272E-02
DAV:   9    -0.108582986809E+02   -0.37880E-04   -0.10941E-04 10544   0.838E-02
   1 F= -.10858299E+02 E0= -.10858299E+02  d E =-.108583E+02
 curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
 trial: gam= 0.00000 g(F)=  0.820E-45 g(S)=  0.161E-01 ort = 0.000E+00 (trialstep = 0.100E+01)
 search vector abs. value=  0.161E-01
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108674017165E+02   -0.10867E+02   -0.24938E-01  9040   0.283E+00    0.316E-01
DAV:   2    -0.108630876032E+02    0.43141E-02   -0.38499E-03 10152   0.427E-01    0.180E-01
DAV:   3    -0.108608719695E+02    0.22156E-02   -0.47378E-03  8664   0.347E-01    0.236E-02
DAV:   4    -0.108608612489E+02    0.10721E-04   -0.55127E-05  6192   0.677E-02
   2 F= -.10860861E+02 E0= -.10860861E+02  d E =-.256257E-02
 trial-energy change:   -0.002563  1 .order   -0.002850   -0.016107    0.010406
 step:   0.6075(harm=  0.6075)  dis= 0.00000  next Energy=   -10.863191 (dE=-0.489E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108640013018E+02   -0.10864E+02   -0.39184E-02  9000   0.113E+00    0.126E-01
DAV:   2    -0.108633291650E+02    0.67214E-03   -0.59649E-04 10176   0.168E-01    0.718E-02
DAV:   3    -0.108629924549E+02    0.33671E-03   -0.73964E-04  8616   0.137E-01    0.939E-03
DAV:   4    -0.108629944710E+02   -0.20161E-05   -0.81055E-06  3808   0.267E-02
   3 F= -.10862994E+02 E0= -.10862994E+02  d E =-.469579E-02
 curvature:  -0.30 expect dE=-0.306E-05 dE for cont linesearch -0.306E-05
 ZBRENT: interpolating
 opt :   0.5927  next Energy=   -10.862997 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629988933E+02   -0.10863E+02   -0.51211E-05  3752   0.421E-02    0.441E-03
DAV:   2    -0.108629982718E+02    0.62157E-06   -0.25034E-06  3664   0.125E-02
   4 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469959E-02
 curvature:  -0.30 expect dE=-0.168E-06 dE for cont linesearch -0.168E-06
 ZBRENT: interpolating
 opt :   0.5881  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629978226E+02   -0.10863E+02   -0.34226E-06  3648   0.124E-02    0.181E-03
DAV:   2    -0.108629977825E+02    0.40149E-07   -0.35818E-07  3640   0.453E-03
   5 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469910E-02
 curvature:  -0.39 expect dE=-0.226E-06 dE for cont linesearch -0.226E-06
 ZBRENT: interpolating
 opt :   0.5904  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629981381E+02   -0.10863E+02   -0.19216E-06  3640   0.791E-03    0.636E-04
DAV:   2    -0.108629981288E+02    0.93201E-08   -0.80950E-08  3640   0.245E-03
   6 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469945E-02
 curvature:  -0.19 expect dE=-0.116E-06 dE for cont linesearch -0.116E-06
 ZBRENT: interpolating
 opt :   0.5916  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629982012E+02   -0.10863E+02   -0.26843E-07  3640   0.328E-03    0.425E-04
DAV:   2    -0.108629981785E+02    0.22696E-07   -0.19955E-08  3640   0.111E-03
   7 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469950E-02
 curvature:  -0.09 expect dE=-0.229E-07 dE for cont linesearch -0.229E-07
 trial: gam=-0.00406 g(F)=  0.128E-45 g(S)=  0.268E-06 ort = 0.657E-04 (trialstep = 0.918E+00)
 search vector abs. value=  0.104E-09
 reached required accuracy - stopping structural energy minimisation
Total resource usage:
 Tasks  : 4
 Memory : 4107.00 MB
 CPU    : 3.990000

Average per task :
 Memory : 1026.75 MB
 CPU    : .99
job completed on Wednesday 24 October 2007, at 19:11:02
(also runs correctly)
<span class='smallblacktext'>[ Edited Wed Oct 31 2007, 03:45PM ]</span>

fatal internal in brackting...

Posted: Wed Oct 31, 2007 2:46 pm
by Danny
and now multithreaded: 1cpu and 2 threads/cpu:

Code: Select all

step_6
job started on Thursday 25 October 2007, at 01:40:48
 running on    1 nodes
 distr:  one band on    1 nodes,    1 groups
 vasp.4.6.31 08Feb07 complex 
 POSCAR found :  1 types and    2 ions
 LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro
 POSCAR, INCAR and KPOINTS ok, starting setup
 FFT: planning ... 1
 reading WAVECAR
 entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.541995558798E+01    0.54200E+01   -0.32487E+03  7280   0.801E+02
DAV:   2    -0.108298552793E+02   -0.16250E+02   -0.15628E+02 11120   0.778E+01
DAV:   3    -0.110322582370E+02   -0.20240E+00   -0.20240E+00  9760   0.111E+01
DAV:   4    -0.110330766749E+02   -0.81844E-03   -0.81844E-03 11392   0.747E-01
DAV:   5    -0.110330776903E+02   -0.10154E-05   -0.10154E-05  9992   0.213E-02    0.301E+00
DAV:   6    -0.109131780101E+02    0.11990E+00   -0.84721E-02  8776   0.151E+00    0.179E+00
DAV:   7    -0.108573305462E+02    0.55847E-01   -0.16141E-01  9248   0.218E+00    0.158E-01
DAV:   8    -0.108582608010E+02   -0.93025E-03   -0.42091E-03  8280   0.480E-01    0.272E-02
DAV:   9    -0.108582986809E+02   -0.37880E-04   -0.10941E-04 10544   0.838E-02
   1 F= -.10858299E+02 E0= -.10858299E+02  d E =-.108583E+02
 curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
 trial: gam= 0.00000 g(F)=  0.539E-46 g(S)=  0.161E-01 ort = 0.000E+00 (trialstep = 0.100E+01)
 search vector abs. value=  0.161E-01
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108674017216E+02   -0.10867E+02   -0.24938E-01  9040   0.283E+00    0.316E-01
DAV:   2    -0.108630875912E+02    0.43141E-02   -0.38499E-03 10152   0.427E-01    0.180E-01
DAV:   3    -0.108608719486E+02    0.22156E-02   -0.47378E-03  8664   0.347E-01    0.236E-02
DAV:   4    -0.108608612280E+02    0.10721E-04   -0.55128E-05  6192   0.677E-02
   2 F= -.10860861E+02 E0= -.10860861E+02  d E =-.256255E-02
 trial-energy change:   -0.002563  1 .order   -0.002850   -0.016107    0.010406
 step:   0.6075(harm=  0.6075)  dis= 0.00000  next Energy=   -10.863191 (dE=-0.489E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108640013075E+02   -0.10864E+02   -0.39184E-02  9000   0.113E+00    0.126E-01
DAV:   2    -0.108633291667E+02    0.67214E-03   -0.59649E-04 10176   0.168E-01    0.718E-02
DAV:   3    -0.108629924545E+02    0.33671E-03   -0.73965E-04  8616   0.137E-01    0.939E-03
DAV:   4    -0.108629944707E+02   -0.20161E-05   -0.81055E-06  3808   0.267E-02
   3 F= -.10862994E+02 E0= -.10862994E+02  d E =-.469579E-02
 curvature:  -0.30 expect dE=-0.306E-05 dE for cont linesearch -0.306E-05
 ZBRENT: interpolating
 opt :   0.5927  next Energy=   -10.862997 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629988933E+02   -0.10863E+02   -0.51212E-05  3752   0.421E-02    0.441E-03
DAV:   2    -0.108629982718E+02    0.62158E-06   -0.25035E-06  3664   0.125E-02
   4 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469959E-02
 curvature:  -0.30 expect dE=-0.168E-06 dE for cont linesearch -0.168E-06
 ZBRENT: interpolating
 opt :   0.5881  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629978226E+02   -0.10863E+02   -0.34249E-06  3648   0.124E-02    0.181E-03
DAV:   2    -0.108629977824E+02    0.40188E-07   -0.35835E-07  3640   0.453E-03
   5 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469910E-02
 curvature:  -0.39 expect dE=-0.226E-06 dE for cont linesearch -0.226E-06
 ZBRENT: interpolating
 opt :   0.5904  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629981381E+02   -0.10863E+02   -0.19221E-06  3640   0.791E-03    0.636E-04
DAV:   2    -0.108629981287E+02    0.93258E-08   -0.80969E-08  3640   0.245E-03
   6 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469945E-02
 curvature:  -0.19 expect dE=-0.116E-06 dE for cont linesearch -0.116E-06
 ZBRENT: interpolating
 opt :   0.5916  next Energy=   -10.862998 (dE=-0.470E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.108629982012E+02   -0.10863E+02   -0.26865E-07  3640   0.328E-03    0.425E-04
DAV:   2    -0.108629981785E+02    0.22704E-07   -0.19969E-08  3640   0.111E-03
   7 F= -.10862998E+02 E0= -.10862998E+02  d E =-.469950E-02
 curvature:  -0.09 expect dE=-0.230E-07 dE for cont linesearch -0.230E-07
 trial: gam=-0.00406 g(F)=  0.215E-46 g(S)=  0.268E-06 ort = 0.657E-04 (trialstep = 0.918E+00)
 search vector abs. value=  0.104E-09
 reached required accuracy - stopping structural energy minimisation
Total resource usage:
 Tasks  : 1
 Memory : 10317.00 MB
 CPU    : 1.100000

Average per task :
 Memory : 10317.00 MB
 CPU    : 1.10
job completed on Thursday 25 October 2007, at 01:58:15

still working fine

fatal internal in brackting...

Posted: Wed Oct 31, 2007 2:57 pm
by Danny
but when I run the same job on 1 CPU, 4threads/cpu I get:

Code: Select all

step_7
job started on Thursday 25 October 2007, at 01:58:26
 running on    1 nodes
 distr:  one band on    1 nodes,    1 groups
 vasp.4.6.31 08Feb07 complex 
 POSCAR found :  1 types and    2 ions
 LDA part: xc-table for Ceperly-Alder, Vosko type interpolation para-ferro
 POSCAR, INCAR and KPOINTS ok, starting setup
 FFT: planning ... 1
 reading WAVECAR
 entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.527107267689E+01    0.52711E+01   -0.32438E+03  7280   0.801E+02
DAV:   2    -0.108411462071E+02   -0.16112E+02   -0.15532E+02 11072   0.776E+01
DAV:   3    -0.110418128896E+02   -0.20067E+00   -0.20067E+00  9688   0.112E+01
DAV:   4    -0.110426010265E+02   -0.78814E-03   -0.78814E-03 11368   0.743E-01
DAV:   5    -0.110426020903E+02   -0.10638E-05   -0.10456E-05  9936   0.212E-02BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00114 new    8.00000
    0.312E+00
DAV:   6    -0.105930682893E+02    0.44953E+00   -0.65238E-01  7528   0.878E+00    0.184E+00
DAV:   7    -0.107330651923E+02   -0.14000E+00   -0.13668E+00  7688   0.971E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00209 new    8.00000
    0.109E+00
DAV:   8    -0.105508792177E+02    0.18219E+00   -0.66775E-01  7456   0.802E+00    0.287E-01
DAV:   9    -0.107105430165E+02   -0.15966E+00   -0.58651E-01  7464   0.779E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.106E+00
DAV:  10    -0.105266830017E+02    0.18386E+00   -0.60392E-01  7456   0.778E+00    0.198E-01
DAV:  11    -0.107138464510E+02   -0.18716E+00   -0.56045E-01  7456   0.740E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  12    -0.105279104520E+02    0.18594E+00   -0.58286E-01  7456   0.777E+00    0.183E-01
DAV:  13    -0.107153801812E+02   -0.18747E+00   -0.55949E-01  7448   0.743E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  14    -0.105287152099E+02    0.18666E+00   -0.58338E-01  7448   0.778E+00    0.182E-01
DAV:  15    -0.107164543969E+02   -0.18774E+00   -0.55812E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  16    -0.105290047349E+02    0.18745E+00   -0.58388E-01  7448   0.779E+00    0.182E-01
DAV:  17    -0.107167642591E+02   -0.18776E+00   -0.55821E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  18    -0.105290423363E+02    0.18772E+00   -0.58407E-01  7448   0.779E+00    0.182E-01
DAV:  19    -0.107168176894E+02   -0.18778E+00   -0.55828E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  20    -0.105290532327E+02    0.18776E+00   -0.58407E-01  7448   0.779E+00    0.182E-01
DAV:  21    -0.107168170020E+02   -0.18776E+00   -0.55826E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  22    -0.105290591325E+02    0.18776E+00   -0.58403E-01  7448   0.779E+00    0.182E-01
DAV:  23    -0.107168080961E+02   -0.18775E+00   -0.55823E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  24    -0.105290621651E+02    0.18775E+00   -0.58399E-01  7448   0.779E+00    0.182E-01
DAV:  25    -0.107167999073E+02   -0.18774E+00   -0.55819E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  26    -0.105290637166E+02    0.18774E+00   -0.58396E-01  7448   0.779E+00    0.182E-01
DAV:  27    -0.107167934737E+02   -0.18773E+00   -0.55817E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  28    -0.105290646100E+02    0.18773E+00   -0.58394E-01  7448   0.779E+00    0.182E-01
DAV:  29    -0.107167884424E+02   -0.18772E+00   -0.55815E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  30    -0.105290652170E+02    0.18772E+00   -0.58393E-01  7448   0.779E+00    0.182E-01
DAV:  31    -0.107167843973E+02   -0.18772E+00   -0.55813E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  32    -0.104878493280E+02    0.22894E+00   -0.62364E-01  7448   0.803E+00    0.182E-01
DAV:  33    -0.107170040754E+02   -0.22915E+00   -0.59508E-01  7448   0.762E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.105E+00
DAV:  34    -0.105291755982E+02    0.18783E+00   -0.58274E-01  7448   0.778E+00    0.182E-01
DAV:  35    -0.107168309649E+02   -0.18766E+00   -0.55794E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  36    -0.105291026726E+02    0.18773E+00   -0.58394E-01  7448   0.779E+00    0.182E-01
DAV:  37    -0.107168081816E+02   -0.18771E+00   -0.55806E-01  7448   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  38    -0.105290920798E+02    0.18772E+00   -0.58389E-01  7448   0.779E+00    0.182E-01
DAV:  39    -0.107814089451E+02   -0.25232E+00   -0.53641E-01  7448   0.730E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  40    -0.105636956921E+02    0.21771E+00   -0.68046E-01  7368   0.877E+00    0.187E-01
DAV:  41    -0.108562991494E+02   -0.29260E+00   -0.19305E-01  7752   0.444E+00    0.603E-02
DAV:  42    -0.108627995033E+02   -0.65004E-02   -0.51979E-03 12408   0.796E-01    0.490E-02
DAV:  43    -0.107896392114E+02    0.73160E-01   -0.72666E-01  7560   0.875E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00211 new    8.00000
    0.105E+00
DAV:  44    -0.105640892162E+02    0.22555E+00   -0.70443E-01  7368   0.884E+00    0.186E-01
DAV:  45    -0.108563641011E+02   -0.29227E+00   -0.19144E-01  7616   0.437E+00    0.498E-02
DAV:  46    -0.108626789148E+02   -0.63148E-02   -0.48171E-03 11000   0.768E-01    0.381E-02
DAV:  47    -0.107898656308E+02    0.72813E-01   -0.72584E-01  7552   0.875E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.105E+00
DAV:  48    -0.105638856477E+02    0.22598E+00   -0.70477E-01  7376   0.884E+00    0.188E-01
DAV:  49    -0.108563339038E+02   -0.29245E+00   -0.19156E-01  7608   0.437E+00    0.439E-02
DAV:  50    -0.104879086702E+02    0.36843E+00   -0.19460E-01  7408   0.443E+00    0.513E-02
DAV:  51    -0.108617241346E+02   -0.37382E+00   -0.18669E-01  7344   0.425E+00    0.201E-02
DAV:  52    -0.108611853194E+02    0.53882E-03   -0.25988E-03 10968   0.545E-01    0.167E-02
DAV:  53    -0.107894728757E+02    0.71712E-01   -0.68635E-01  7544   0.863E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.105E+00
DAV:  54    -0.105635503134E+02    0.22592E+00   -0.70605E-01  7376   0.885E+00    0.193E-01
DAV:  55    -0.108562954977E+02   -0.29275E+00   -0.19229E-01  7624   0.438E+00    0.384E-02
DAV:  56    -0.104876183483E+02    0.36868E+00   -0.19493E-01  7416   0.444E+00    0.617E-02
DAV:  57    -0.108617330848E+02   -0.37411E+00   -0.18709E-01  7352   0.425E+00    0.999E-03
DAV:  58    -0.108611896736E+02    0.54341E-03   -0.26000E-03 10840   0.544E-01    0.812E-03
DAV:  59    -0.107889304236E+02    0.72259E-01   -0.68607E-01  7544   0.863E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.105E+00
DAV:  60    -0.105634185250E+02    0.22551E+00   -0.70620E-01  7376   0.886E+00
   1 F= -.10563419E+02 E0= -.10563419E+02  d E =-.105634E+02
 curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
 trial: gam= 0.00000 g(F)=  0.455E-38 g(S)=  0.957E-02 ort = 0.000E+00 (trialstep = 0.100E+01)
 search vector abs. value=  0.957E-02
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.104816064924E+02   -0.10482E+02   -0.19011E-01 10072   0.289E+00    0.249E-01
DAV:   2    -0.107802669300E+02   -0.29866E+00   -0.75286E-01  7664   0.872E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.103E+00
DAV:   3    -0.109457948706E+02   -0.16553E+00   -0.82736E-01  7464   0.909E+00    0.294E-01
DAV:   4    -0.108533647835E+02    0.92430E-01   -0.24942E-01  7416   0.538E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00201 new    8.00000
    0.102E+00
DAV:   5    -0.105060208594E+02    0.34734E+00   -0.18130E-01  7464   0.542E+00    0.214E-01
DAV:   6    -0.108604977199E+02   -0.35448E+00   -0.13291E-01  7392   0.418E+00    0.679E-02
DAV:   7    -0.108663190994E+02   -0.58214E-02   -0.30123E-03 11792   0.613E-01    0.577E-02
DAV:   8    -0.108258096888E+02    0.40509E-01   -0.11801E-01  7520   0.387E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.103E+00
DAV:   9    -0.109167586138E+02   -0.90949E-01   -0.20761E-01  7432   0.545E+00    0.180E-01
DAV:  10    -0.107793610324E+02    0.13740E+00   -0.71861E-01  7464   0.861E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00205 new    8.00000
    0.104E+00
DAV:  11    -0.105708952002E+02    0.20847E+00   -0.69002E-01  7368   0.879E+00    0.208E-01
DAV:  12    -0.108602857371E+02   -0.28939E+00   -0.19009E-01  7624   0.432E+00    0.420E-02
DAV:  13    -0.108666006255E+02   -0.63149E-02   -0.47623E-03 11208   0.752E-01    0.254E-02
DAV:  14    -0.108253969358E+02    0.41204E-01   -0.11847E-01  7512   0.388E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.103E+00
DAV:  15    -0.109171383862E+02   -0.91741E-01   -0.20908E-01  7424   0.545E+00    0.181E-01
DAV:  16    -0.107792061041E+02    0.13793E+00   -0.72383E-01  7448   0.862E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00205 new    8.00000
    0.104E+00
DAV:  17    -0.105711423416E+02    0.20806E+00   -0.69107E-01  7376   0.879E+00    0.203E-01
DAV:  18    -0.108602984144E+02   -0.28916E+00   -0.18941E-01  7624   0.431E+00    0.388E-02
DAV:  19    -0.108665492551E+02   -0.62508E-02   -0.46410E-03 11312   0.750E-01    0.155E-02
DAV:  20    -0.107879255033E+02    0.78624E-01   -0.70752E-01  7536   0.857E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.103E+00
DAV:  21    -0.105721737584E+02    0.21575E+00   -0.69330E-01  7384   0.869E+00    0.197E-01
DAV:  22    -0.108603766667E+02   -0.28820E+00   -0.18565E-01  7592   0.427E+00    0.365E-02
DAV:  23    -0.105011780785E+02    0.35920E+00   -0.18760E-01  7416   0.433E+00    0.709E-02
DAV:  24    -0.107196012929E+02   -0.21842E+00   -0.77577E-01  7568   0.885E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.103E+00
DAV:  25    -0.105370748438E+02    0.18253E+00   -0.57904E-01  7504   0.755E+00    0.181E-01
DAV:  26    -0.107110520362E+02   -0.17398E+00   -0.54171E-01  7496   0.733E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.104E+00
DAV:  27    -0.105367353069E+02    0.17432E+00   -0.56199E-01  7496   0.755E+00    0.183E-01
DAV:  28    -0.107113141606E+02   -0.17458E+00   -0.53635E-01  7496   0.726E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  29    -0.105369923733E+02    0.17432E+00   -0.56188E-01  7496   0.757E+00    0.182E-01
DAV:  30    -0.107115851174E+02   -0.17459E+00   -0.53707E-01  7496   0.725E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  31    -0.105371946896E+02    0.17439E+00   -0.56227E-01  7496   0.758E+00    0.182E-01
DAV:  32    -0.107117996794E+02   -0.17460E+00   -0.53759E-01  7496   0.724E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  33    -0.105373461618E+02    0.17445E+00   -0.56259E-01  7496   0.759E+00    0.182E-01
DAV:  34    -0.107119771272E+02   -0.17463E+00   -0.53794E-01  7496   0.724E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  35    -0.105374634090E+02    0.17451E+00   -0.56283E-01  7496   0.759E+00    0.182E-01
DAV:  36    -0.107121298527E+02   -0.17467E+00   -0.53819E-01  7496   0.724E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  37    -0.105375569756E+02    0.17457E+00   -0.56300E-01  7496   0.759E+00    0.182E-01
DAV:  38    -0.107122654792E+02   -0.17471E+00   -0.53835E-01  7496   0.724E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  39    -0.105376337806E+02    0.17463E+00   -0.56312E-01  7496   0.760E+00    0.181E-01
DAV:  40    -0.107123891965E+02   -0.17476E+00   -0.53846E-01  7496   0.724E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  41    -0.105376987461E+02    0.17469E+00   -0.56319E-01  7496   0.760E+00    0.181E-01
DAV:  42    -0.107125047404E+02   -0.17481E+00   -0.53853E-01  7496   0.724E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  43    -0.105720734519E+02    0.14043E+00   -0.69478E-01  7416   0.877E+00    0.191E-01
DAV:  44    -0.108603354600E+02   -0.28826E+00   -0.18765E-01  7712   0.435E+00    0.410E-02
DAV:  45    -0.108666162824E+02   -0.62808E-02   -0.48343E-03 12408   0.763E-01    0.243E-02
DAV:  46    -0.108257048153E+02    0.40911E-01   -0.11781E-01  7512   0.388E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.102E+00
DAV:  47    -0.109177012156E+02   -0.91996E-01   -0.20955E-01  7424   0.545E+00    0.185E-01
DAV:  48    -0.107802349097E+02    0.13747E+00   -0.72456E-01  7448   0.862E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.103E+00
DAV:  49    -0.105715296089E+02    0.20871E+00   -0.68936E-01  7360   0.879E+00    0.193E-01
DAV:  50    -0.108603269592E+02   -0.28880E+00   -0.18926E-01  7624   0.431E+00    0.382E-02
DAV:  51    -0.108665487234E+02   -0.62218E-02   -0.46288E-03 11208   0.748E-01    0.192E-02
DAV:  52    -0.107892064318E+02    0.77342E-01   -0.70623E-01  7536   0.857E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.103E+00
DAV:  53    -0.105723878579E+02    0.21682E+00   -0.69180E-01  7384   0.868E+00    0.191E-01
DAV:  54    -0.104417986610E+02    0.13059E+00   -0.83371E-02  9400   0.321E+00    0.650E-02
DAV:  55    -0.108663239277E+02   -0.42453E+00   -0.24937E-01  7568   0.510E+00    0.210E-02
DAV:  56    -0.107891767263E+02    0.77147E-01   -0.70044E-01  7520   0.850E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.103E+00
DAV:  57    -0.105725912008E+02    0.21659E+00   -0.69425E-01  7384   0.871E+00    0.189E-01
DAV:  58    -0.108604124204E+02   -0.28782E+00   -0.18506E-01  7608   0.428E+00    0.390E-02
DAV:  59    -0.108665707897E+02   -0.61584E-02   -0.45289E-03 11792   0.741E-01    0.233E-02
DAV:  60    -0.108258720640E+02    0.40699E-01   -0.11758E-01  7512   0.388E+00
BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
   2 F= -.10825872E+02 E0= -.10825872E+02  d E =-.262454E+00
 trial-energy change:   -0.262454  1 .order    0.003944   -0.009567    0.017455
 step:   0.3540(harm=  0.3540)  dis= 0.00000  next Energy=   -10.565112 (dE=-0.169E-02)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.113321302861E+02   -0.11332E+02   -0.79906E-01  8528   0.876E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00171 new    8.00000
    0.110E+00
DAV:   2    -0.108761451457E+02    0.45599E+00   -0.59537E-01  7576   0.840E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00170 new    8.00000
    0.104E+00
DAV:   3    -0.105419016653E+02    0.33424E+00   -0.78511E-01  7880   0.941E+00    0.333E-01
DAV:   4    -0.108598287362E+02   -0.31793E+00   -0.18531E-01  7712   0.468E+00    0.976E-02
DAV:   5    -0.108648894627E+02   -0.50607E-02   -0.44196E-03 12160   0.714E-01    0.555E-02
DAV:   6    -0.107896394314E+02    0.75250E-01   -0.71040E-01  7536   0.866E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.105E+00
DAV:   7    -0.105689547746E+02    0.22068E+00   -0.69993E-01  7368   0.879E+00    0.193E-01
DAV:   8    -0.104335585401E+02    0.13540E+00   -0.86404E-02  9384   0.327E+00    0.544E-02
DAV:   9    -0.108655178088E+02   -0.43196E+00   -0.25534E-01  7568   0.518E+00    0.383E-02
DAV:  10    -0.107915664130E+02    0.73951E-01   -0.71291E-01  7512   0.862E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00211 new    8.00000
    0.104E+00
DAV:  11    -0.105689892039E+02    0.22258E+00   -0.70206E-01  7368   0.881E+00    0.186E-01
DAV:  12    -0.108596158001E+02   -0.29063E+00   -0.18954E-01  7600   0.435E+00    0.545E-02
DAV:  13    -0.108659139344E+02   -0.62981E-02   -0.48785E-03 11656   0.771E-01    0.466E-02
DAV:  14    -0.107921493589E+02    0.73765E-01   -0.71879E-01  7536   0.868E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00211 new    8.00000
    0.104E+00
DAV:  15    -0.105692183700E+02    0.22293E+00   -0.69918E-01  7376   0.878E+00    0.183E-01
DAV:  16    -0.108596472141E+02   -0.29043E+00   -0.18942E-01  7592   0.433E+00    0.583E-02
DAV:  17    -0.104940640798E+02    0.36558E+00   -0.19200E-01  7408   0.439E+00    0.278E-02
DAV:  18    -0.108649896966E+02   -0.37093E+00   -0.18442E-01  7352   0.422E+00    0.435E-02
DAV:  19    -0.108644610873E+02    0.52861E-03   -0.25858E-03 10824   0.541E-01    0.444E-02
DAV:  20    -0.107924721038E+02    0.71989E-01   -0.67904E-01  7536   0.856E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00211 new    8.00000
    0.104E+00
DAV:  21    -0.105693300461E+02    0.22314E+00   -0.69941E-01  7384   0.879E+00    0.182E-01
DAV:  22    -0.108596739549E+02   -0.29034E+00   -0.18937E-01  7592   0.434E+00    0.614E-02
DAV:  23    -0.108660049651E+02   -0.63310E-02   -0.48168E-03 11192   0.764E-01    0.564E-02
DAV:  24    -0.107923032524E+02    0.73702E-01   -0.71777E-01  7536   0.868E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00211 new    8.00000
    0.104E+00
DAV:  25    -0.105693532545E+02    0.22295E+00   -0.69854E-01  7376   0.878E+00    0.181E-01
DAV:  26    -0.108596937151E+02   -0.29034E+00   -0.18915E-01  7592   0.433E+00    0.648E-02
DAV:  27    -0.108660409096E+02   -0.63472E-02   -0.47878E-03 11120   0.760E-01    0.607E-02
DAV:  28    -0.107923729621E+02    0.73668E-01   -0.71745E-01  7536   0.868E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.104E+00
DAV:  29    -0.105694098074E+02    0.22296E+00   -0.69834E-01  7384   0.877E+00    0.180E-01
DAV:  30    -0.108597177200E+02   -0.29031E+00   -0.18905E-01  7592   0.433E+00    0.682E-02
DAV:  31    -0.108660799642E+02   -0.63622E-02   -0.47709E-03 11232   0.759E-01    0.649E-02
DAV:  32    -0.107924140919E+02    0.73666E-01   -0.71708E-01  7536   0.867E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.104E+00
DAV:  33    -0.105694644298E+02    0.22295E+00   -0.69813E-01  7384   0.877E+00    0.180E-01
DAV:  34    -0.108597418491E+02   -0.29028E+00   -0.18895E-01  7592   0.433E+00    0.716E-02
DAV:  35    -0.108661185847E+02   -0.63767E-02   -0.47562E-03 11136   0.758E-01    0.690E-02
DAV:  36    -0.107924515706E+02    0.73667E-01   -0.71676E-01  7536   0.867E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.104E+00
DAV:  37    -0.105695165986E+02    0.22293E+00   -0.69792E-01  7384   0.877E+00    0.179E-01
DAV:  38    -0.108597655965E+02   -0.29025E+00   -0.18885E-01  7592   0.433E+00    0.749E-02
DAV:  39    -0.108661555843E+02   -0.63900E-02   -0.47425E-03 11256   0.756E-01    0.729E-02
DAV:  40    -0.107924967310E+02    0.73659E-01   -0.71643E-01  7536   0.867E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.104E+00
DAV:  41    -0.105695657691E+02    0.22293E+00   -0.69771E-01  7384   0.877E+00    0.179E-01
DAV:  42    -0.108597885694E+02   -0.29022E+00   -0.18877E-01  7592   0.433E+00    0.781E-02
DAV:  43    -0.104943002570E+02    0.36549E+00   -0.19182E-01  7416   0.439E+00    0.155E-02
DAV:  44    -0.108650391235E+02   -0.37074E+00   -0.18403E-01  7352   0.421E+00    0.672E-02
DAV:  45    -0.108645148482E+02    0.52428E-03   -0.25905E-03 10792   0.542E-01    0.691E-02
DAV:  46    -0.108260530841E+02    0.38462E-01   -0.10288E-01  7512   0.378E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.103E+00
DAV:  47    -0.109200216834E+02   -0.93969E-01   -0.21767E-01  7432   0.554E+00    0.211E-01
DAV:  48    -0.107839760805E+02    0.13605E+00   -0.74300E-01  7432   0.874E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00202 new    8.00000
    0.104E+00
DAV:  49    -0.105691489808E+02    0.21483E+00   -0.69455E-01  7368   0.887E+00    0.180E-01
DAV:  50    -0.108597534805E+02   -0.29060E+00   -0.19204E-01  7608   0.437E+00    0.845E-02
DAV:  51    -0.108662896992E+02   -0.65362E-02   -0.48051E-03 11296   0.764E-01    0.828E-02
DAV:  52    -0.107926905055E+02    0.73599E-01   -0.71651E-01  7536   0.867E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.104E+00
DAV:  53    -0.105696617782E+02    0.22303E+00   -0.69725E-01  7392   0.877E+00    0.178E-01
DAV:  54    -0.108598238806E+02   -0.29016E+00   -0.18847E-01  7592   0.433E+00    0.831E-02
DAV:  55    -0.108662361280E+02   -0.64122E-02   -0.46996E-03 11088   0.753E-01    0.823E-02
DAV:  56    -0.107927757993E+02    0.73460E-01   -0.71569E-01  7536   0.867E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00209 new    8.00000
    0.104E+00
DAV:  57    -0.105696532265E+02    0.22312E+00   -0.69721E-01  7392   0.877E+00    0.178E-01
DAV:  58    -0.108598325797E+02   -0.29018E+00   -0.18847E-01  7592   0.433E+00    0.844E-02
DAV:  59    -0.108662550898E+02   -0.64225E-02   -0.47054E-03 11312   0.753E-01    0.843E-02
DAV:  60    -0.108259111440E+02    0.40344E-01   -0.11941E-01  7512   0.389E+00
BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00209 new    8.00000
   3 F= -.10825911E+02 E0= -.10825911E+02  d E =-.262493E+00
 curvature:  -0.18 expect dE=-0.557E-03 dE for cont linesearch -0.557E-03
 ZBRENT: interpolating
 opt :   0.2250  next Energy=   -10.826265 (dE=-0.263E+00)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.113261096322E+02   -0.11326E+02   -0.68873E-01  7984   0.854E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00169 new    8.00000
    0.108E+00
DAV:   2    -0.108643916511E+02    0.46172E+00   -0.72668E-01  7400   0.885E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00169 new    8.00000
    0.102E+00
DAV:   3    -0.109532980618E+02   -0.88906E-01   -0.33559E-01  7352   0.600E+00    0.355E-01
DAV:   4    -0.108172651549E+02    0.13603E+00   -0.93612E-01  7384   0.907E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00206 new    8.00000
    0.104E+00
DAV:   5    -0.105721871366E+02    0.24508E+00   -0.71264E-01  7352   0.880E+00    0.244E-01
DAV:   6    -0.108580180616E+02   -0.28583E+00   -0.19593E-01  7608   0.438E+00    0.716E-02
DAV:   7    -0.108644424013E+02   -0.64243E-02   -0.46727E-03 11872   0.751E-01    0.601E-02
DAV:   8    -0.108247089685E+02    0.39733E-01   -0.12106E-01  7520   0.396E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00215 new    8.00000
    0.105E+00
DAV:   9    -0.109159607549E+02   -0.91252E-01   -0.21505E-01  7424   0.555E+00    0.182E-01
DAV:  10    -0.107790629558E+02    0.13690E+00   -0.74045E-01  7432   0.876E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00207 new    8.00000
    0.106E+00
DAV:  11    -0.105653905365E+02    0.21367E+00   -0.70411E-01  7352   0.893E+00    0.215E-01
DAV:  12    -0.108585015082E+02   -0.29311E+00   -0.19590E-01  7616   0.440E+00    0.383E-02
DAV:  13    -0.108648984049E+02   -0.63969E-02   -0.50793E-03 11064   0.789E-01    0.170E-02
DAV:  14    -0.107890759439E+02    0.75822E-01   -0.72506E-01  7552   0.872E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.105E+00
DAV:  15    -0.105668500218E+02    0.22223E+00   -0.70440E-01  7384   0.882E+00    0.197E-01
DAV:  16    -0.108586427983E+02   -0.29179E+00   -0.19142E-01  7592   0.435E+00    0.446E-02
DAV:  17    -0.108649522823E+02   -0.63095E-02   -0.48652E-03 11184   0.771E-01    0.256E-02
DAV:  18    -0.107892919886E+02    0.75660E-01   -0.72330E-01  7536   0.872E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.104E+00
DAV:  19    -0.105669610800E+02    0.22233E+00   -0.70346E-01  7376   0.881E+00    0.194E-01
DAV:  20    -0.108586685932E+02   -0.29171E+00   -0.19098E-01  7592   0.435E+00    0.450E-02
DAV:  21    -0.104919797458E+02    0.36669E+00   -0.19324E-01  7408   0.441E+00    0.602E-02
DAV:  22    -0.108640446483E+02   -0.37206E+00   -0.18533E-01  7352   0.423E+00    0.184E-02
DAV:  23    -0.108635146520E+02    0.53000E-03   -0.25709E-03 10864   0.541E-01    0.176E-02
DAV:  24    -0.107891306027E+02    0.74384E-01   -0.68305E-01  7536   0.859E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.104E+00
DAV:  25    -0.105670746963E+02    0.22206E+00   -0.70433E-01  7384   0.882E+00    0.192E-01
DAV:  26    -0.104328235601E+02    0.13425E+00   -0.86063E-02  9416   0.329E+00    0.684E-02
DAV:  27    -0.108648482097E+02   -0.43202E+00   -0.25656E-01  7560   0.520E+00    0.389E-02
DAV:  28    -0.108241232491E+02    0.40725E-01   -0.12267E-01  7528   0.396E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00213 new    8.00000
    0.104E+00
DAV:  29    -0.109182491645E+02   -0.94126E-01   -0.21872E-01  7424   0.558E+00    0.195E-01
DAV:  30    -0.107809521423E+02    0.13730E+00   -0.74675E-01  7432   0.877E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.104E+00
DAV:  31    -0.105665628884E+02    0.21439E+00   -0.69994E-01  7360   0.890E+00    0.191E-01
DAV:  32    -0.108586600707E+02   -0.29210E+00   -0.19361E-01  7616   0.439E+00    0.544E-02
DAV:  33    -0.108651803563E+02   -0.65203E-02   -0.48768E-03 11552   0.772E-01    0.471E-02
DAV:  34    -0.107894689553E+02    0.75711E-01   -0.72196E-01  7544   0.871E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.104E+00
DAV:  35    -0.105674089245E+02    0.22206E+00   -0.70227E-01  7376   0.881E+00    0.187E-01
DAV:  36    -0.108587641260E+02   -0.29136E+00   -0.19024E-01  7592   0.435E+00    0.585E-02
DAV:  37    -0.108652113736E+02   -0.64472E-02   -0.47832E-03 11080   0.761E-01    0.535E-02
DAV:  38    -0.107897395629E+02    0.75472E-01   -0.72080E-01  7536   0.870E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00212 new    8.00000
    0.104E+00
DAV:  39    -0.105674972267E+02    0.22224E+00   -0.70171E-01  7384   0.880E+00    0.185E-01
DAV:  40    -0.108587990755E+02   -0.29130E+00   -0.19008E-01  7592   0.435E+00    0.637E-02
DAV:  41    -0.108652559748E+02   -0.64569E-02   -0.47653E-03 11264   0.759E-01    0.598E-02
DAV:  42    -0.108243333034E+02    0.40923E-01   -0.12026E-01  7512   0.391E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00211 new    8.00000
    0.103E+00
DAV:  43    -0.109189704363E+02   -0.94637E-01   -0.21861E-01  7432   0.555E+00    0.206E-01
DAV:  44    -0.107816495405E+02    0.13732E+00   -0.74794E-01  7432   0.877E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  45    -0.105669838170E+02    0.21467E+00   -0.69839E-01  7368   0.889E+00    0.184E-01
DAV:  46    -0.108587713849E+02   -0.29179E+00   -0.19309E-01  7616   0.438E+00    0.718E-02
DAV:  47    -0.108653556702E+02   -0.65843E-02   -0.48309E-03 11224   0.766E-01    0.701E-02
DAV:  48    -0.108244079790E+02    0.40948E-01   -0.12024E-01  7520   0.391E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00211 new    8.00000
    0.103E+00
DAV:  49    -0.109192356595E+02   -0.94828E-01   -0.21909E-01  7432   0.554E+00    0.211E-01
DAV:  50    -0.107817898548E+02    0.13745E+00   -0.74863E-01  7432   0.877E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  51    -0.105671534907E+02    0.21464E+00   -0.69786E-01  7368   0.889E+00    0.183E-01
DAV:  52    -0.108588304487E+02   -0.29168E+00   -0.19290E-01  7616   0.438E+00    0.797E-02
DAV:  53    -0.108654407946E+02   -0.66103E-02   -0.48073E-03 11264   0.764E-01    0.789E-02
DAV:  54    -0.107906293396E+02    0.74811E-01   -0.71930E-01  7536   0.870E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.104E+00
DAV:  55    -0.105677946083E+02    0.22283E+00   -0.70017E-01  7392   0.879E+00    0.180E-01
DAV:  56    -0.108589405409E+02   -0.29115E+00   -0.18947E-01  7592   0.434E+00    0.826E-02
DAV:  57    -0.108654603089E+02   -0.65198E-02   -0.47082E-03 11136   0.753E-01    0.830E-02
DAV:  58    -0.107907987371E+02    0.74662E-01   -0.71833E-01  7536   0.869E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00210 new    8.00000
    0.104E+00
DAV:  59    -0.105678807004E+02    0.22292E+00   -0.69985E-01  7392   0.879E+00    0.179E-01
DAV:  60    -0.104345206230E+02    0.13336E+00   -0.84999E-02  9376   0.328E+00
   4 F= -.10434521E+02 E0= -.10434521E+02  d E =0.128898E+00
 curvature:  -0.11 expect dE=-0.508E-02 dE for cont linesearch -0.508E-02
 ZBRENT: interpolating
 opt :   0.3271  next Energy=   -10.825985 (dE=-0.263E+00)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.107258123863E+02   -0.10726E+02   -0.95505E-01  7680   0.967E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00209 new    8.00000
    0.104E+00
DAV:   2    -0.108631782598E+02   -0.13737E+00   -0.39402E-02 11264   0.234E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00167 new    8.00000
    0.102E+00
DAV:   3    -0.105331698320E+02    0.33001E+00   -0.89159E-01  7544   0.938E+00    0.245E-01
DAV:   4    -0.107075736298E+02   -0.17440E+00   -0.74288E-01  7488   0.827E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.105E+00
DAV:   5    -0.105476644920E+02    0.15991E+00   -0.55613E-01  7464   0.765E+00    0.197E-01
DAV:   6    -0.107174337876E+02   -0.16977E+00   -0.54358E-01  7456   0.737E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:   7    -0.104944584849E+02    0.22298E+00   -0.61409E-01  7456   0.798E+00    0.182E-01
DAV:   8    -0.107184552202E+02   -0.22400E+00   -0.58721E-01  7456   0.755E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:   9    -0.105348721994E+02    0.18358E+00   -0.57617E-01  7456   0.772E+00    0.181E-01
DAV:  10    -0.107183020563E+02   -0.18343E+00   -0.55214E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  11    -0.105344582665E+02    0.18384E+00   -0.57681E-01  7456   0.773E+00    0.181E-01
DAV:  12    -0.107181158920E+02   -0.18366E+00   -0.55142E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  13    -0.104936203991E+02    0.22450E+00   -0.61606E-01  7456   0.796E+00    0.181E-01
DAV:  14    -0.107181376558E+02   -0.22452E+00   -0.58777E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  15    -0.105343940708E+02    0.18374E+00   -0.57564E-01  7456   0.771E+00    0.181E-01
DAV:  16    -0.107181176155E+02   -0.18372E+00   -0.55130E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  17    -0.104936010841E+02    0.22452E+00   -0.61596E-01  7456   0.796E+00    0.181E-01
DAV:  18    -0.107181331537E+02   -0.22453E+00   -0.58769E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  19    -0.105343599707E+02    0.18377E+00   -0.57556E-01  7456   0.771E+00    0.181E-01
DAV:  20    -0.107181070136E+02   -0.18375E+00   -0.55113E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  21    -0.104935957848E+02    0.22451E+00   -0.61588E-01  7456   0.796E+00    0.181E-01
DAV:  22    -0.107181214621E+02   -0.22453E+00   -0.58766E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  23    -0.105343425768E+02    0.18378E+00   -0.57554E-01  7456   0.771E+00    0.181E-01
DAV:  24    -0.107180984524E+02   -0.18376E+00   -0.55105E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  25    -0.105342600337E+02    0.18384E+00   -0.57655E-01  7456   0.772E+00    0.181E-01
DAV:  26    -0.107180391553E+02   -0.18378E+00   -0.55107E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  27    -0.104935741351E+02    0.22447E+00   -0.61584E-01  7456   0.796E+00    0.181E-01
DAV:  28    -0.107180784181E+02   -0.22450E+00   -0.58770E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  29    -0.105343129906E+02    0.18377E+00   -0.57554E-01  7456   0.771E+00    0.181E-01
DAV:  30    -0.107180750632E+02   -0.18376E+00   -0.55103E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  31    -0.105342442958E+02    0.18383E+00   -0.57653E-01  7456   0.772E+00    0.181E-01
DAV:  32    -0.107180263518E+02   -0.18378E+00   -0.55107E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  33    -0.105342186110E+02    0.18381E+00   -0.57652E-01  7456   0.772E+00    0.181E-01
DAV:  34    -0.107180004004E+02   -0.18378E+00   -0.55109E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  35    -0.104935597585E+02    0.22444E+00   -0.61583E-01  7456   0.796E+00    0.181E-01
DAV:  36    -0.107180466417E+02   -0.22449E+00   -0.58775E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  37    -0.105342883190E+02    0.18376E+00   -0.57555E-01  7456   0.771E+00    0.181E-01
DAV:  38    -0.107180559861E+02   -0.18377E+00   -0.55101E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  39    -0.104935785168E+02    0.22448E+00   -0.61582E-01  7456   0.796E+00    0.181E-01
DAV:  40    -0.107180722724E+02   -0.22449E+00   -0.58773E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  41    -0.105342960453E+02    0.18378E+00   -0.57555E-01  7456   0.771E+00    0.181E-01
DAV:  42    -0.107180688953E+02   -0.18377E+00   -0.55097E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  43    -0.104935839454E+02    0.22448E+00   -0.61580E-01  7456   0.796E+00    0.181E-01
DAV:  44    -0.107180780921E+02   -0.22449E+00   -0.58773E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  45    -0.105342966113E+02    0.18378E+00   -0.57555E-01  7456   0.771E+00    0.181E-01
DAV:  46    -0.107180729447E+02   -0.18378E+00   -0.55094E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  47    -0.105342359498E+02    0.18384E+00   -0.57651E-01  7456   0.772E+00    0.181E-01
DAV:  48    -0.107180273591E+02   -0.18379E+00   -0.55095E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  49    -0.105342347233E+02    0.18379E+00   -0.57635E-01  7456   0.772E+00    0.181E-01
DAV:  50    -0.107180049530E+02   -0.18377E+00   -0.55097E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  51    -0.105342214092E+02    0.18378E+00   -0.57637E-01  7456   0.772E+00    0.181E-01
DAV:  52    -0.107179853262E+02   -0.18376E+00   -0.55099E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  53    -0.105342143743E+02    0.18377E+00   -0.57640E-01  7456   0.772E+00    0.181E-01
DAV:  54    -0.107179722567E+02   -0.18376E+00   -0.55101E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  55    -0.104935679667E+02    0.22440E+00   -0.61570E-01  7456   0.796E+00    0.181E-01
DAV:  56    -0.107180108086E+02   -0.22444E+00   -0.58798E-01  7456   0.756E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  57    -0.105342212038E+02    0.18379E+00   -0.57571E-01  7456   0.771E+00    0.181E-01
DAV:  58    -0.107179899875E+02   -0.18377E+00   -0.55080E-01  7456   0.736E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
    0.104E+00
DAV:  59    -0.104935729797E+02    0.22442E+00   -0.61573E-01  7456   0.796E+00    0.181E-01
DAV:  60    -0.107180210471E+02   -0.22445E+00   -0.58798E-01  7456   0.756E+00
BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00203 new    8.00000
   5 F= -.10718021E+02 E0= -.10718021E+02  d E =-.154603E+00
 curvature:   0.00 expect dE=-0.366E-03 dE for cont linesearch -0.366E-03
 ZBRENT: interpolating
 opt :   0.2566  next Energy=   -10.434850 (dE= 0.129E+00)
 bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.106799631364E+02   -0.10680E+02   -0.99529E-03 11640   0.116E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00170 new    8.00000
    0.106E+00
DAV:   2    -0.108256057407E+02   -0.14564E+00   -0.27423E-02  7368   0.222E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00168 new    8.00000
    0.102E+00


..... I removed 2000lines off the same story....


DAV:  48    -0.105290862108E+02    0.18764E+00   -0.58376E-01  7456   0.779E+00    0.182E-01
DAV:  49    -0.107167379298E+02   -0.18765E+00   -0.55798E-01  7456   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  50    -0.105290856052E+02    0.18765E+00   -0.58376E-01  7456   0.779E+00    0.182E-01
DAV:  51    -0.107167399889E+02   -0.18765E+00   -0.55798E-01  7456   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  52    -0.105290793399E+02    0.18766E+00   -0.58374E-01  7456   0.779E+00    0.182E-01
DAV:  53    -0.107167214272E+02   -0.18764E+00   -0.55795E-01  7456   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  54    -0.105290755848E+02    0.18765E+00   -0.58371E-01  7456   0.779E+00    0.182E-01
DAV:  55    -0.107167147982E+02   -0.18764E+00   -0.55793E-01  7456   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  56    -0.105290744331E+02    0.18764E+00   -0.58371E-01  7456   0.779E+00    0.182E-01
DAV:  57    -0.107167117298E+02   -0.18764E+00   -0.55793E-01  7456   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  58    -0.105290740417E+02    0.18764E+00   -0.58371E-01  7456   0.779E+00    0.182E-01
DAV:  59    -0.107167144539E+02   -0.18764E+00   -0.55794E-01  7456   0.742E+00BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00204 new    8.00000
    0.105E+00
DAV:  60    -0.105290743818E+02    0.18764E+00   -0.58372E-01  7456   0.779E+00
BRMIX: very serious problems
 the old and the new charge density differ
 old charge density:     8.00000 new    8.00212
  20 F= -.10529074E+02 E0= -.10529074E+02  d E =0.343441E-01
 curvature:   0.00 expect dE= 0.000E+00 dE for cont linesearch  0.000E+00
 ZBRENT: fatal internal in brackting
         system-shutdown; contact gK immediately
Total resource usage:
 Tasks  : 1
 Memory : 10344.00 MB
 CPU    : 1.470000

Average per task :
 Memory : 10344.00 MB
 CPU    : 1.47
job completed on Thursday 25 October 2007, at 11:59:33

(same goes for more than 4 threads...but when I run the job on say 4 cpu's, 4 threads or 8 threads /cpu it works perfectly again.)

As can be seen from the last posts, geometry can't really be the cause, because then the job should crash in all cases given here.

any thoughts what could be the cause and how it could be resolved?

thx
Danny

fatal internal in brackting...

Posted: Mon Nov 05, 2007 2:52 pm
by Danny
Hi,

I'm still interested in what might be causing this problem with VASP.

thanx
Danny

fatal internal in brackting...

Posted: Mon Nov 19, 2007 4:28 pm
by Danny
Hi,

I'm still interested in what could cause this problem, since it is clear from the above examples that it has no connection to the actual geometry of the system.

Or should I direct my question elsewhere?
I really want to understand the cause of this problem.

thanx
Danny

fatal internal in brackting...

Posted: Fri Nov 23, 2007 2:08 pm
by admin
maybe the reason is the following: some compilers automatically support OPEN_MP if several processes are running. This should NOT be done, because all parallelization is done explicitely in the vasp code itself. (the errors at the first mixing step after the pre-conditioning electronic steps may be due to some parallelization-related problem)

Therefore please add a line
export OMP_NUM_THREADS=1
in your .bashrc file to suppress any additional parallelization which is not encoded in vasp itself.

fatal internal in brackting...

Posted: Tue Nov 27, 2007 3:30 pm
by Danny
Thank you for the confirmation. I already suspected something like this, but I was told MPI and the Open_MP(done by the compiler) would not clash ...it seems they do, and thus it is a code-related problem.

Danny

fatal internal in brackting...

Posted: Wed Nov 28, 2007 10:15 am
by admin
well, in VASP all reasonable parallelization is hard-coded, all that has to be done is to tell the compiler to do no additional Open-MPI parallelization (this behaviour is compiler-dependent), or to set the environment-variables such that any additional parallelization is suppressed.

fatal internal in brackting...

Posted: Fri May 29, 2009 9:15 pm
by cylee
Does the default MKL_DYNAMIC in MKL 10.1 detect the optimal threading and solve the parallelization inefficiency for VASP?
http://software.intel.com/en-us/article ... threading/

I have Intel E5410 8 cores per node (2 quads). When it's using 4 cores on the node, VASP will be threaded to 120-150% CPU for all 4 VASP processes. And with 8 processes, it comes back to 100% CPU per process.

I optimize a H2 molecule for VASP benchmarking. Within the same node, ppn=1 as 100% walltime:
ppn=2: 0.70
ppn=3: 0.54
ppn=4: 0.52
ppn=5: 0.45
ppn=6: 0.62
ppn=7: 0.81
ppn=8: 0.60

Is this normal for an 8-core node?
Should I use a bigger system as benchmark?
Or use OMP_NUM_THREADS = 1?

Thanks in advance for your help.

CY

fatal internal in brackting...

Posted: Fri May 29, 2009 9:21 pm
by cylee
Here is a ppt presentation by Intel:
MKL PPT

On page 11 they state that if the thread count is larger than physical core count, it will be automatically reduced to the latter.

Thoughts ?!