0
ERROR: Bottom of the band has a negative slope
[-1.78250019 -1.7825425 -1.78262535 -1.78274518 -1.78289668 -1.7830728
-1.78326476 -1.78346212 -1.78365276 -1.78397332 -1.78452752 -1.78506972
-1.78558128 -1.78604213 -1.78643085 -1.78672474 -1.7868999 -1.7869313
-1.78679289 -1.78645766]
[ -4.23128288e-05 -8.28468952e-05 -1.19829629e-04 -1.51500824e-04
-1.76118758e-04 -1.91966253e-04 -1.97356635e-04 -1.90639590e-04
-3.20564685e-04 -5.54193916e-04 -5.42200868e-04 -5.11560512e-04
-4.60851662e-04 -3.88721742e-04 -2.93891329e-04 -1.75158436e-04
-3.14025103e-05 1.38411861e-04 3.35231550e-04]
[False False True True True True True True True True True True
True True True True False False False]
ERROR: Chemical potential exceeds the evaporation threshold
mu0 = -0.597
T = 0.240
Eth = -0.455
ERROR: Chemical potential exceeds the bottom of the band along 100
mu0 = -0.597
T = 0.240
E100 = -0.455
ERROR: Radial density profile along 111 has a positive slope
radius check start = 0.001
[ 1.01039381e+00 1.01042412e+00 1.01048476e+00 1.01057582e+00
1.01069736e+00 1.01084950e+00 1.01103235e+00 1.01124602e+00
1.01149062e+00 1.01177669e+00 1.01211331e+00 1.01248498e+00
1.01289184e+00 1.01333394e+00 1.01381113e+00 1.01432312e+00
1.01486935e+00 1.01544895e+00 1.01606075e+00 1.01670312e+00
1.01740547e+00 1.01820037e+00 1.01903726e+00 1.01990552e+00
1.02080131e+00 1.02171991e+00 1.02265570e+00 1.02363860e+00
1.02473970e+00 1.02585024e+00 1.02696081e+00 1.02806062e+00
1.02913750e+00 1.03029839e+00 1.03156743e+00 1.03279356e+00
1.03395906e+00 1.03504464e+00 1.03606002e+00 1.03724693e+00
1.03831273e+00 1.03923188e+00 1.04002114e+00 1.04071099e+00
1.04152001e+00 1.04209776e+00 1.04241043e+00 1.04242282e+00
1.04239719e+00 1.04226318e+00 1.04176776e+00 1.04087017e+00
1.03951789e+00 1.03830525e+00 1.03655424e+00 1.03420793e+00
1.03120483e+00 1.02794440e+00 1.02428703e+00 1.01980780e+00
1.01441387e+00 1.00812222e+00 1.00164715e+00 9.94019942e-01
9.85125487e-01 9.74841321e-01 9.63812193e-01 9.51716700e-01
9.37975745e-01 9.22488561e-01 9.05188515e-01 8.87696156e-01
8.68463887e-01 8.47548824e-01 8.25063158e-01 8.01644744e-01
7.78674486e-01 7.54999253e-01 7.30963718e-01 7.06912770e-01
6.83676098e-01 6.62445105e-01 6.41732782e-01 6.21413967e-01
6.01122504e-01 5.80519315e-01 5.61439652e-01 5.41347983e-01
5.19692556e-01 4.96151555e-01 4.70582327e-01 4.43709754e-01
4.17612468e-01 3.90208041e-01 3.61858699e-01 3.33107820e-01
3.04459200e-01 2.76385278e-01 2.51391883e-01 2.27725224e-01
2.05251969e-01 1.84149965e-01 1.64536031e-01 1.46471231e-01
1.29968576e-01 1.15707971e-01 1.03180481e-01 9.18166110e-02
8.15620294e-02 7.23481365e-02 6.41004541e-02 5.67418523e-02
5.01950930e-02 4.43847547e-02 3.95686860e-02 3.52907068e-02
3.14742816e-02 2.80743252e-02 2.50488412e-02 2.23590451e-02
1.99693795e-02 1.78474482e-02 1.59638923e-02 1.42922273e-02
1.28086397e-02 1.15493293e-02 1.04442361e-02 9.45509621e-03
8.56955096e-03 7.77649764e-03 7.06598212e-03 6.42909562e-03
5.85787682e-03 5.34522004e-03 4.88478993e-03 4.47094277e-03
4.09865443e-03 3.76345462e-03 3.46136743e-03 3.18885747e-03
2.94278149e-03 2.72034494e-03 2.52515205e-03 2.35144631e-03
2.19320137e-03 2.04890328e-03 1.91719486e-03 1.79685913e-03
1.68680456e-03 1.58605176e-03 1.49372158e-03 1.40902452e-03
1.33125116e-03 1.25976373e-03 1.19398847e-03 1.13340890e-03
1.07755977e-03 1.02602169e-03 9.78416291e-04 9.34401985e-04]
[False False False True True True True True True True True True
True True True True True True True True True True True True
True True True True True True True True True True True True
True True True True True True True True True True False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False False False False False False False False False False
False False False]
0.95511238578
2708644.2959
1
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.200093376909
23637.5956254
2
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.732341759488
2217546.76447
3
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.70144345962
299586.152208
4
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.0435736677116
300268.631553
5
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.456879877276
300072.944547
6
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.429467084639
300112.686314
7
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.704155561003
299850.856949
8
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.0429467084639
300134.643761
9
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.450210098046
300176.99891
10
OK: Bottom of the band has positive slope up to r111 = 10 um
OK: Chemical potential is below evaporation threshold.
OK: Chemical potential is below the bottom of the band along 100
OK: Radial density profile along 111 decreases monotonically.
0.423197492163
300784.859313
scubic.py:271: RuntimeWarning: invalid value encountered in sqrt
nu.append( (1./2./np.pi) * np.sqrt(c2) * factor )