Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Last revisionBoth sides next revision
forum:start [2020/02/22 19:08] – Maurits W. Haverkortforum:start [2020/12/02 16:00] – BODRY TEGOMO CHIOGO
Line 1: Line 1:
-~~NOTOC~~ 
-====== Forum ====== 
  
-### +====== Ground state  in ligand field calculations for Cerium  ======
-If you can not find your answer in the documentation feel free to ask your question here. If you want to be notified of an answer by e-mail you can register yourself to get e-mail notifications with the [[forum:start?do=subscribe|E-mail]] button on the right top of this page. +
-###+
  
-<form> +<WRAP center box 100%> 
-action   template forum:data:template "forum:data:%Y:@@Title@@"  +Dear developers, 
-thanks   "Thanks for asking" +I wrote a simple code to calculate the ground state of cerium compounds using Ligand field calculation without crytal elcetric field. 
-fieldset "Create new question" +Without spin-orbit ,I think  the code works correctly and gives me the energies and states with expected weights for  configurations $L^{14} f^0$,$L^{13}  f^{1}$, and $L^{12}  f^{2}$
-textbox  "Title"/^[^\/:]+$/ +
-textarea "Question" +
-submit +
-</form>+
  
-===== Questions =====+ —The $L^{12}  f^{2}$ configuration has an energy $2*\Delta+U_{ff}$ 
 + 
 + —The $L^{13}  f^{1}$ configuration has an energy $\Delta$ 
 + 
 + —The $L^{14} f^0$   configuration has an energy 0 
 + 
 +**Please I would like to known if it possible to calculate in Quanty the weight of the configuration $L^{13}  f^{1}_{5/2}$ and $L^{13}  f^{1}_{7/2}$ when I add spin-orbit (zeta_4f = 0.087 for $Ce^{3+}$ )** 
 + 
 + —The $L^{12}  f^{2}$ configuration has an energy $2*\Delta+U_{ff}$ 
 + 
 + —The $L^{13}  f^{1}_{7/2}$ configuration has an energy $\Delta$ +  $\Delta$so 
 + 
 + 
 + —The $L^{13}  f^{1}_{5/2}$ configuration has an energy $\Delta$ 
 + 
 + —The $L^{14} f^0$   configuration has an energy 0 
 + 
 + 
 +**Secondly, I don't understand why Quanty gime in output 322 states. I expected to have 106 states with  6 $f^1_{5/2}$ states and  8 $f^1_{7/2}$ states. ** 
 + 
 +All the best, 
 + 
 +Bodry 
 + 
 +here is the scrip and the output: 
 + 
 +<code> 
 +Verbosity(0*0000) 
 + 
 +NF = 38 
 +NB = 0 
 +IndexDn_3d = {0,2,4,6,8} 
 +IndexUp_3d = {1,3,5,7,9} 
 +IndexDn_4f = {10,12,14,16,18,20,22} 
 +IndexUp_4f = {11,13,15,17,19,21,23} 
 +IndexDn_Lf = {24,26,28,30,32,34,36} 
 +IndexUp_Lf = {25,27,29,31,33,35,37} 
 + 
 + 
 + 
 +OppSx_4f = NewOperator("Sx", NF, IndexUp_4f,IndexDn_4f) 
 +OppSy_4f = NewOperator("Sy", NF, IndexUp_4f,IndexDn_4f) 
 +OppSz_4f = NewOperator("Sz", NF, IndexUp_4f,IndexDn_4f) 
 +OppSsqr_4f = NewOperator("Ssqr", NF, IndexUp_4f,IndexDn_4f) 
 +OppSplus_4f = NewOperator("Splus", NF, IndexUp_4f,IndexDn_4f) 
 +OppSmin_4f = NewOperator("Smin", NF, IndexUp_4f,IndexDn_4f) 
 + 
 +OppLx_4f = NewOperator("Lx", NF, IndexUp_4f,IndexDn_4f) 
 +OppLy_4f = NewOperator("Ly", NF, IndexUp_4f,IndexDn_4f) 
 +OppLz_4f = NewOperator("Lz", NF, IndexUp_4f,IndexDn_4f) 
 +OppLsqr_4f = NewOperator("Lsqr", NF, IndexUp_4f,IndexDn_4f) 
 +OppLplus_4f = NewOperator("Lplus", NF, IndexUp_4f,IndexDn_4f) 
 +OppLmin_4f = NewOperator("Lmin", NF, IndexUp_4f,IndexDn_4f) 
 + 
 +OppJx_4f = NewOperator("Jx", NF, IndexUp_4f,IndexDn_4f) 
 +OppJy_4f = NewOperator("Jy", NF, IndexUp_4f,IndexDn_4f) 
 +OppJz_4f = NewOperator("Jz", NF, IndexUp_4f,IndexDn_4f) 
 +OppJsqr_4f = NewOperator("Jsqr", NF, IndexUp_4f,IndexDn_4f) 
 +OppJplus_4f = NewOperator("Jplus", NF, IndexUp_4f,IndexDn_4f) 
 +OppJmin_4f = NewOperator("Jmin", NF, IndexUp_4f,IndexDn_4f) 
 + 
 +Oppldots_4f = NewOperator("ldots", NF, IndexUp_4f,IndexDn_4f) 
 +Oppldots_3d = NewOperator('ldots', NF, IndexUp_3d, IndexDn_3d) 
 + 
 +-- Angular momentum operator for the figand f-shell 
 + 
 +OppSx_Lf = NewOperator("Sx", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppSy_Lf = NewOperator("Sy", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppSz_Lf = NewOperator("Sz", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppSsqr_Lf = NewOperator("Ssqr", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppSplus_Lf = NewOperator("Splus", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppSmin_Lf = NewOperator("Smin", NF, IndexUp_Lf,IndexDn_Lf) 
 + 
 +OppLx_Lf = NewOperator("Lx", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppLy_Lf = NewOperator("Ly", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppLz_Lf = NewOperator("Lz", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppLsqr_Lf = NewOperator("Lsqr", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppLplus_Lf = NewOperator("Lplus", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppLmin_Lf = NewOperator("Lmin", NF, IndexUp_Lf,IndexDn_Lf) 
 + 
 +OppJx_Lf = NewOperator("Jx", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppJy_Lf = NewOperator("Jy", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppJz_Lf = NewOperator("Jz", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppJsqr_Lf = NewOperator("Jsqr", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppJplus_Lf = NewOperator("Jplus", NF, IndexUp_Lf,IndexDn_Lf) 
 +OppJmin_Lf = NewOperator("Jmin", NF, IndexUp_Lf,IndexDn_Lf) 
 + 
 +-- SUM OF THE OPERATOR 
 + 
 +OppSx = OppSx_4f + OppSx_Lf 
 +OppSy = OppSy_4f + OppSy_Lf 
 +OppSz = OppSz_4f + OppSz_Lf 
 +OppSsqr = OppSx*OppSx + OppSy*OppSy + OppSz*OppSz 
 + 
 +OppLx = OppLx_4f + OppLx_Lf 
 +OppLy = OppLy_4f + OppLy_Lf 
 +OppLz = OppLz_4f + OppLz_Lf 
 +OppLsqr = OppLx*OppLx  + OppLy*OppLy + OppLz*OppLz 
 + 
 +OppJx = OppJx_4f + OppJx_Lf 
 +OppJy = OppJy_4f + OppJy_Lf 
 +OppJz = OppJz_4f + OppJz_Lf 
 +OppJsqr = OppJx*OppJx + OppJy*OppJy + OppJz*OppJz 
 + 
 + 
 + 
 + 
 + OppF0_4f = NewOperator("U", NF, IndexUp_4f, IndexDn_4f, {1, 0, 0, 0}) 
 + OppF2_4f = NewOperator("U", NF, IndexUp_4f, IndexDn_4f, {0, 1, 0, 0}) 
 + OppF4_4f = NewOperator("U", NF, IndexUp_4f, IndexDn_4f, {0, 0, 1, 0}) 
 + OppF6_4f = NewOperator("U", NF, IndexUp_4f, IndexDn_4f, {0, 0, 0, 1}) 
 +  
 + OppF0_3d_4f = NewOperator('U', NF, IndexUp_3d, IndexDn_3d, IndexUp_4f, IndexDn_4f, {1, 0, 0}, {0, 0, 0}); 
 + OppF2_3d_4f = NewOperator('U', NF, IndexUp_3d, IndexDn_3d, IndexUp_4f, IndexDn_4f, {0, 1, 0}, {0, 0, 0}); 
 + OppF4_3d_4f = NewOperator('U', NF, IndexUp_3d, IndexDn_3d, IndexUp_4f, IndexDn_4f, {0, 0, 1}, {0, 0, 0}); 
 +  
 + OppG1_3d_4f = NewOperator("U", NF, IndexUp_3d, IndexDn_3d, IndexUp_4f, IndexDn_4f, {0, 0, 0}, {1, 0, 0}) 
 + OppG3_3d_4f = NewOperator("U", NF, IndexUp_3d, IndexDn_3d, IndexUp_4f, IndexDn_4f, {0, 0, 0}, {0, 1, 0}) 
 + OppG5_3d_4f = NewOperator("U", NF, IndexUp_3d, IndexDn_3d, IndexUp_4f, IndexDn_4f, {0, 0, 0}, {0, 0, 1}) 
 +  
 +  
 +OppNUp_4f = NewOperator("Number", NF, IndexUp_4f,IndexUp_4f,{1,1,1,1,1,1,1}) 
 +OppNDn_4f = NewOperator("Number", NF, IndexDn_4f,IndexDn_4f,{1,1,1,1,1,1,1}) 
 + 
 +OppN_4f = OppNUp_4f + OppNDn_4f 
 + 
 +OppNUp_Lf = NewOperator("Number", NF, IndexUp_Lf,IndexUp_Lf,{1,1,1,1,1,1,1}) 
 +OppNDn_Lf = NewOperator("Number", NF, IndexDn_Lf,IndexDn_Lf,{1,1,1,1,1,1,1}) 
 +OppN_Lf = OppNUp_Lf + OppNDn_Lf 
 + 
 +OppNUp_3d = NewOperator("Number", NF, IndexUp_3d, IndexUp_3d, {1, 1, 1, 1, 1}) 
 +OppNDn_3d = NewOperator("Number", NF, IndexDn_3d, IndexDn_3d, {1, 1, 1, 1, 1}) 
 +OppN_3d = OppNUp_3d  + OppNDn_3d 
 +      
 + zeta_4f_i = 0.087* 1.0 
 + U_4f_4f_i = 8.5 
 + F2_4f_4f_i = 0 * 0.8 
 + F4_4f_4f_i = 0 * 0.8 
 + F6_4f_4f_i = 0 * 0.8 
 + F0_4f_4f_i = U_4f_4f_i + 4 / 195 * F2_4f_4f_i + 2 / 143 * F4_4f_4f_i + 100 / 5577 * F6_4f_4f_i 
 +  
 +  
 + zeta_4f_f = 0.119 * 1.0 
 + zeta_3d_f = 19.393 * 1.0 
 +  
 + U_4f_4f_f = 8.5 
 + U_3d_4f_f =- 10 
 + F2_4f_4f_f = 10.10 * 0.8 
 + F4_4f_4f_f = 6.35 * 0.8 
 + F6_4f_4f_f = 4.57 * 0.8 
 + F0_4f_4f_f = U_4f_4f_f + 4 / 195 * F2_4f_4f_f + 2 / 143 * F4_4f_4f_f + 100 / 5577 * F6_4f_4f_f 
 +  
 + F2_3d_4f_f = 9.8 * 1 
 + F4_3d_4f_f = 4.6* 1 
 + G1_3d_4f_f = 7.3 * 1 
 + G3_3d_4f_f = 4.2 * 1 
 + G5_3d_4f_f = 3 * 1 
 + F0_3d_4f_f = U_3d_4f_f + 3 / 70 * G1_3d_4f_f + 2 / 105 * G3_3d_4f_f + 5 / 231 * G5_3d_4f_f 
 + 
 + Bz = 0.000001 
 +  
 +H_i = Chop(F0_4f_4f_i *OppF0_4f + F2_4f_4f_i *OppF2_4f+ F4_4f_4f_i *OppF4_4f + F6_4f_4f_i *OppF6_4f +  zeta_4f_i * Oppldots_4f +  Bz*(2*OppSz_4f + OppLz_4f)) 
 + 
 +H_f = Chop(F0_4f_4f_f *OppF0_4f + F2_4f_4f_f *OppF2_4f+ F4_4f_4f_f *OppF4_4f + F6_4f_4f_f *OppF6_4f +  G1_3d_4f_f *OppG1_3d_4f + G3_3d_4f_f *OppG3_3d_4f + G5_3d_4f_f *OppG5_3d_4f + 
 + F0_3d_4f_f *OppF0_3d_4f+ F2_3d_4f_f *OppF2_3d_4f + F4_3d_4f_f *OppF4_3d_4f+  zeta_4f_f * Oppldots_4f +  zeta_3d_f * Oppldots_3d +  Bz*(2*OppSz_4f + OppLz_4f)) 
 +  
 +  
 + 
 +  
 +Delta_4f_Lf_i = -2.3 
 + 
 +NElectrons_4f = 0 
 + 
 +e_4f_i = (28 * Delta_4f_Lf_i - 27 * U_4f_4f_i * NElectrons_4f - U_4f_4f_i * NElectrons_4f^2) / (2 * (14 + NElectrons_4f)) 
 +e_Lf_i = NElectrons_4f * (-2 * Delta_4f_Lf_i + U_4f_4f_i * NElectrons_4f + U_4f_4f_i) / (2 * (NElectrons_4f + 14)) 
 + 
 + 
 + 
 +H_i = H_i + Chop( e_4f_i * OppN_4f + e_Lf_i * OppN_Lf) 
 + 
 +Delta_4f_L1_f = -2.3 
 +e_4f_f = (28 * Delta_4f_L1_f - 460 * U_3d_4f_f - U_4f_4f_f * NElectrons_4f^2 - 47 * U_4f_4f_f * NElectrons_4f) / (2 * (NElectrons_4f + 24)) 
 +e_3d_f = (28 * Delta_4f_L1_f - 2 * U_3d_4f_f * NElectrons_4f^2 - 30 * U_3d_4f_f * NElectrons_4f - 28 * U_3d_4f_f + U_4f_4f_f * NElectrons_4f^2 + U_4f_4f_f * NElectrons_4f) / (2 * (NElectrons_4f + 24)) 
 +e_Lf_f = (-2 * Delta_4f_L1_f * NElectrons_4f - 20 * Delta_4f_L1_f + 20 * U_3d_4f_f * NElectrons_4f + 20 * U_3d_4f_f + U_4f_4f_f * NElectrons_4f^2 + U_4f_4f_f * NElectrons_4f) / (2 * (NElectrons_4f + 24)) 
 + 
 +H_f = H_f + Chop( e_4f_f * OppN_4f + e_3d_f * OppN_3d + e_Lf_f * OppN_Lf) 
 + 
 + 
 +NElectrons_3d = 10 
 + 
 + 
 + 
 +V = 0.23 
 + 
 + 
 + 
 +OppV = NewOperator("Number", NF, IndexUp_4f,IndexUp_Lf, {1,1,1,1,1,1,1}) + NewOperator("Number", NF, IndexDn_4f,IndexDn_Lf, {1,1,1,1,1,1,1})+ NewOperator("Number", NF, IndexUp_Lf,IndexUp_4f, {1,1,1,1,1,1,1})+ NewOperator("Number", NF, IndexDn_Lf,IndexDn_4f, {1,1,1,1,1,1,1}) 
 + 
 +H_i  = H_i +Chop(V*OppV) 
 + 
 + 
 + 
 +H_f  = H_f +Chop(V*OppV) 
 + 
 + 
 +NConfigurations = 3 
 + 
 +InitialRestrictions = {NF, NB,{"1111111111 00000000000000 00000000000000",NElectrons_3d,NElectrons_3d}, {"0000000000 11111111111111 00000000000000",NElectrons_4f,NElectrons_4f}, {"0000000000 00000000000000 11111111111111",14,14}} 
 + 
 + 
 +FinalRestrictions = {NF, NB,{"1111111111 00000000000000 00000000000000",NElectrons_3d -1,NElectrons_3d -1}, {"0000000000 11111111111111 00000000000000",NElectrons_4f,NElectrons_4f}, {"0000000000 00000000000000 11111111111111",14,14}} 
 + 
 +CalculationRestrictions = {NF, NB, {'0000000000 00000000000000 11111111111111', 14 - (NConfigurations - 1), 14}} 
 + 
 +Npsis = 322 
 + 
 +psiList = Eigensystem(H_i, InitialRestrictions, Npsis, {{'restrictions', CalculationRestrictions}}) 
 + 
 + 
 + 
 +print(psiList) 
 + 
 + 
 +Oppf1    =  NewOperator ( "Identity"    , NF )  
 +Oppf1.Restrictions =  { NF , NB , { "0000000000 11111111111111 00000000000000" ,  1 ,  1 } } 
 +Oppf0    =  NewOperator ( "Identity"    , NF )  
 +Oppf0.Restrictions =  { NF , NB , { "0000000000 11111111111111 00000000000000" ,  0 ,  0 } }  
 +Oppf2    =  NewOperator ( "Identity"    , NF )  
 +Oppf2.Restrictions =  { NF , NB , { "0000000000 11111111111111 00000000000000" ,  2 , 2 } } 
 + 
 +oppList={H_i, OppSsqr, OppLsqr, OppJsqr, OppSz_4f, OppLz_4f, Oppldots_4f, OppF2_4f, OppF4_4f, OppF6_4f, OppN_4f, OppN_Lf,Oppf0,Oppf1,Oppf2} 
 + 
 + 
 + 
 +print("  #    <E>      <S^2>    <L^2>    <J^2>  <S_z^4f>  <L_z^4f>  <l.s>    <F[2]>   <F[4]>   <F[4]>   <N^4f> <N^Lf>  <W_4f^0>  <W_4f^1>  <W_4f^2>  "); 
 +for i = 1,#psiList do 
 +  io.write(string.format("%3i ",i)) 
 +  for j = 1,#oppList do 
 +    expectationvalue = Chop(psiList[i]*oppList[j]*psiList[i]) 
 +    io.write(string.format("%8.3f ",expectationvalue)) 
 +  end 
 +  io.write("\n") 
 +end 
 + 
 +</code> 
 + 
 +<code> 
 + #    <E>      <S^2>    <L^2>    <J^2>  <S_z^4f>  <L_z^4f>  <l.s>    <F[2]>   <F[4]>   <F[4]>   <N^4f> <N^Lf>  <W_4f^0>  <W_4f^1>  <W_4f^2> 
 +  1   -2.828    0.251    0.251    0.000    0.000   -0.000   -1.170   -0.001   -0.000   -0.000    0.954   13.046    0.072    0.902    0.026 
 +  2   -2.575    1.755   11.755   12.000    0.357   -2.850   -2.001   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 +  3   -2.575    1.592   19.592   20.000    0.214   -1.710   -2.001   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 +  4   -2.575    1.510   23.510   24.000    0.071   -0.570   -2.001   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 +  5   -2.575    1.510   23.510   24.000   -0.071    0.570   -2.001   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 +  6   -2.575    1.592   19.592   20.000   -0.214    1.710   -2.001   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 +  7   -2.575    1.755   11.755   12.000   -0.357    2.850   -2.001   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 +  8   -2.574    1.330   13.110    9.163    0.308   -2.466   -1.996   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 +  9   -2.574    1.478   22.610   16.436    0.153   -1.224   -1.996   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 10   -2.574    1.527   25.702   18.803   -0.000    0.000   -1.996   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 11   -2.574    1.478   22.610   16.436   -0.153    1.224   -1.996   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 12   -2.574    1.330   13.110    9.163   -0.308    2.466   -1.996   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 13   -2.390    1.749    1.749    0.000   -0.000   -0.000    0.650   -0.001   -0.000   -0.000    1.003   12.997    0.014    0.969    0.017 
 + 14   -2.271    1.110    9.944   12.368   -0.455   -2.727    1.495   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 15   -2.271    1.349   18.554   24.087   -0.301   -1.807    1.495   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 16   -2.271    1.489   23.590   30.941   -0.150   -0.901    1.495   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 17   -2.271    1.535   25.254   33.207    0.000    0.000    1.495   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 18   -2.271    1.489   23.590   30.941    0.150    0.901    1.495   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 19   -2.271    1.349   18.554   24.087    0.301    1.807    1.495   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 20   -2.271    1.110    9.944   12.368    0.455    2.727    1.495   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 21   -2.270    1.755   11.755   12.000   -0.357   -2.136    1.499   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 22   -2.270    1.592   19.592   20.000   -0.214   -1.282    1.499   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 23   -2.270    1.510   23.510   24.000   -0.071   -0.427    1.499   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 24   -2.270    1.510   23.510   24.000    0.071    0.427    1.499   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 25   -2.270    1.592   19.592   20.000    0.214    1.282    1.499   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 26   -2.270    1.755   11.755   12.000    0.357    2.136    1.499   -0.000   -0.000   -0.000    1.016   12.984    0.000    0.984    0.016 
 + 27    0.253    0.001    0.001    0.000    0.000    0.000    0.016   -0.000   -0.000   -0.000    0.096   13.904    0.912    0.079    0.008 
 + 28    3.552    3.612   16.612   15.000    0.571   -4.571   -4.000   -0.114   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 29    3.552    2.755   31.719   26.750    0.521   -4.165   -4.000   -0.110   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 30    3.552    2.650   20.992   20.266    0.485   -3.881   -4.000   -0.095   -0.019   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 31    3.552    3.469   35.469   36.000    0.429   -3.429   -4.000   -0.114   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 32    3.552    2.561   19.933   13.511    0.384   -3.070   -4.000   -0.084   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 33    3.552    2.804   36.778   33.313    0.372   -2.977   -4.000   -0.086   -0.022   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 34    3.552    3.256   37.771   31.068    0.363   -2.907   -4.000   -0.089   -0.021   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 35    3.552    2.886   40.691   29.404    0.310   -2.479   -4.000   -0.069   -0.025   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 36    3.552    2.755   21.756   21.000    0.286   -2.286   -4.000   -0.026   -0.036   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 37    3.552    3.061   35.061   36.000    0.286   -2.286   -4.000   -0.065   -0.027   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 38    3.552    2.840   37.264   33.663    0.233   -1.862   -4.000   -0.063   -0.025   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 39    3.552    2.844   43.111   37.319    0.225   -1.800   -4.000   -0.042   -0.030   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 40    3.552    3.229   34.760   28.746    0.208   -1.664   -4.000   -0.035   -0.031   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 41    3.552    2.548   26.150   25.246    0.161   -1.287   -4.000   -0.018   -0.020   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 42    3.552    2.841   37.818   27.205    0.157   -1.255   -4.000   -0.035   -0.022   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 43    3.552    2.980   42.980   44.000    0.143   -1.143   -4.000   -0.075   -0.025   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 44    3.552    2.735   30.735   32.000    0.143   -1.143   -4.000   -0.016   -0.039   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 45    3.552    2.923   33.166   29.522    0.091   -0.730   -4.000   -0.014   -0.024   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 46    3.552    2.861   45.146   39.777    0.072   -0.574   -4.000   -0.016   -0.010   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 47    3.552    3.060   42.780   36.893    0.067   -0.535   -4.000   -0.031   -0.003   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 48    3.552    3.086   25.220   21.437    0.051   -0.406   -4.000    0.051   -0.008   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 49    3.552    2.608   22.896   15.780    0.000   -0.000   -4.000   -0.048   -0.025   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 50    3.552    2.988   47.218   34.401    0.000   -0.000   -4.000   -0.014   -0.006   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 51    3.552    2.685   27.857   19.578    0.000   -0.000   -4.000    0.027   -0.013   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 52    3.552    2.469   23.469   23.000   -0.000    0.000   -4.000    0.052    0.009   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 53    3.552    2.816   38.816   40.000   -0.000    0.000   -4.000    0.003    0.020   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 54    3.552    2.490   22.490   24.000   -0.000    0.000   -4.000    0.082    0.002   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 55    3.552    3.086   25.219   21.437   -0.051    0.406   -4.000    0.051   -0.008   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 56    3.552    3.060   42.780   36.893   -0.067    0.535   -4.000   -0.031   -0.003   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 57    3.552    2.861   45.146   39.777   -0.072    0.574   -4.000   -0.016   -0.010   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 58    3.552    2.923   33.166   29.522   -0.091    0.730   -4.000   -0.014   -0.024   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 59    3.552    2.980   42.979   44.000   -0.143    1.143   -4.000   -0.075   -0.025   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 60    3.552    2.735   30.735   32.000   -0.143    1.143   -4.000   -0.016   -0.039   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 61    3.552    2.841   37.819   27.205   -0.157    1.255   -4.000   -0.035   -0.022   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 62    3.552    2.548   26.150   25.246   -0.161    1.287   -4.000   -0.018   -0.020   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 63    3.552    3.229   34.759   28.746   -0.208    1.664   -4.000   -0.035   -0.031   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 64    3.552    2.844   43.111   37.319   -0.225    1.800   -4.000   -0.042   -0.030   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 65    3.552    2.840   37.265   33.663   -0.233    1.862   -4.000   -0.063   -0.025   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 66    3.552    2.755   21.757   21.002   -0.286    2.286   -4.000   -0.026   -0.036   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 67    3.552    3.061   35.060   35.998   -0.286    2.286   -4.000   -0.065   -0.027   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 68    3.552    2.886   40.692   29.405   -0.310    2.479   -4.000   -0.069   -0.025   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 69    3.552    3.256   37.771   31.068   -0.363    2.907   -4.000   -0.089   -0.021   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 70    3.552    2.804   36.778   33.313   -0.372    2.977   -4.000   -0.086   -0.022   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 71    3.552    2.561   19.933   13.511   -0.384    3.069   -4.000   -0.084   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 72    3.552    3.469   35.469   36.000   -0.429    3.429   -4.000   -0.114   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 73    3.552    2.650   20.992   20.266   -0.485    3.881   -4.000   -0.095   -0.019   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 74    3.552    2.755   31.719   26.750   -0.521    4.165   -4.000   -0.110   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 75    3.552    3.612   16.612   15.000   -0.571    4.571   -4.000   -0.114   -0.016   -0.000    2.000   12.000    0.000    0.000    1.000 
 + 76    3.580    2.682   14.463    9.163    0.158   -1.265   -3.923   -0.029   -0.019   -0.001    1.996   12.004    0.000    0.004    0.996 
 + 77    3.580    2.831   23.963   16.436    0.078   -0.628   -3.923   -0.049   -0.015   -0.001    1.996   12.004    0.000    0.004    0.996 
 + 78    3.580    2.879   27.055   18.803   -0.000    0.000   -3.923   -0.068   -0.026   -0.001    1.996   12.004    0.000    0.004    0.996 
 + 79    3.580    2.831   23.963   16.436   -0.079    0.628   -3.923   -0.049   -0.015   -0.001    1.996   12.004    0.000    0.004    0.996 
 + 80    3.580    2.682   14.463    9.163   -0.158    1.265   -3.923   -0.029   -0.019   -0.001    1.996   12.004    0.000    0.004    0.996 
 + 81    3.588    3.026   13.026   12.000    0.289   -2.295   -3.907   -0.045   -0.019   -0.001    1.995   12.005    0.000    0.005    0.995 
 + 82    3.588    2.903   20.903   20.000    0.173   -1.377   -3.907   -0.045   -0.019   -0.001    1.995   12.005    0.000    0.005    0.995 
 + 83    3.588    2.842   24.842   24.000    0.058   -0.459   -3.907   -0.045   -0.019   -0.001    1.995   12.005    0.000    0.005    0.995 
 + 84    3.588    2.842   24.842   24.000   -0.058    0.459   -3.907   -0.045   -0.019   -0.001    1.995   12.005    0.000    0.005    0.995 
 + 85    3.588    2.903   20.903   20.000   -0.173    1.377   -3.907   -0.045   -0.019   -0.001    1.995   12.005    0.000    0.005    0.995 
 + 86    3.588    3.026   13.026   12.000   -0.289    2.295   -3.907   -0.045   -0.019   -0.001    1.995   12.005    0.000    0.005    0.995 
 + 87    3.622    2.811    2.811    0.000   -0.000    0.000   -3.806   -0.043   -0.018   -0.001    1.991   12.009    0.000    0.009    0.991 
 + 88    3.856    1.911   12.134   12.313   -0.069   -5.412   -0.500    0.057   -0.003   -0.002    2.000   12.000    0.000    0.000    1.000 
 + 89    3.856    1.315   21.930   21.530   -0.147   -5.207   -0.500    0.038   -0.008   -0.002    2.000   12.000    0.000    0.000    1.000 
 + 90    3.856    2.561   29.582   32.630   -0.234   -4.404   -0.500   -0.024   -0.019   -0.003    2.000   12.000    0.000    0.000    1.000 
 + 91    3.856    2.378   32.344   28.804   -0.302   -3.961   -0.500   -0.033   -0.018   -0.007    2.000   12.000    0.000    0.000    1.000 
 + 92    3.856    2.166   31.151   33.249    0.007   -4.284   -0.500   -0.010   -0.015   -0.005    2.000   12.000    0.000    0.000    1.000 
 + 93    3.856    2.082   30.082   32.000   -0.143   -3.857   -0.500   -0.003    0.005   -0.005    2.000   12.000    0.000    0.000    1.000 
 + 94    3.856    2.618   33.901   36.996   -0.380   -3.279   -0.500   -0.061   -0.017   -0.008    2.000   12.000    0.000    0.000    1.000 
 + 95    3.856    2.082   30.082   32.000    0.143   -4.143   -0.500   -0.048   -0.021   -0.003    2.000   12.000    0.000    0.000    1.000 
 + 96    3.856    2.755   31.719   26.750   -0.194   -3.451   -0.500   -0.008   -0.002   -0.009    2.000   12.000    0.000    0.000    1.000 
 + 97    3.856    2.579   35.589   31.171   -0.456   -2.734   -0.500   -0.053   -0.015   -0.016    2.000   12.000    0.000    0.000    1.000 
 + 98    3.856    3.133   32.839   39.588    0.196   -3.822   -0.500   -0.050   -0.019   -0.006    2.000   12.000    0.000    0.000    1.000 
 + 99    3.856    2.650   20.992   20.266    0.057   -3.453   -0.500   -0.032   -0.017   -0.010    2.000   12.000    0.000    0.000    1.000 
 +100    3.856    2.162   40.499   40.523   -0.148   -3.038   -0.500   -0.006   -0.008   -0.014    2.000   12.000    0.000    0.000    1.000 
 +101    3.856    2.108   21.593   22.079   -0.039   -3.234   -0.500   -0.004   -0.005   -0.013    2.000   12.000    0.000    0.000    1.000 
 +102    3.856    3.061   35.061   36.000   -0.286   -2.714   -0.500   -0.025   -0.025   -0.015    2.000   12.000    0.000    0.000    1.000 
 +103    3.856    2.625   33.784   36.824   -0.524   -2.144   -0.500   -0.059   -0.014   -0.021    2.000   12.000    0.000    0.000    1.000 
 +104    3.856    2.608   36.489   40.104    0.159   -3.376   -0.500   -0.036   -0.016   -0.010    2.000   12.000    0.000    0.000    1.000 
 +105    3.856    3.256   37.771   31.068   -0.351   -2.193   -0.500   -0.017   -0.019   -0.024    2.000   12.000    0.000    0.000    1.000 
 +106    3.856    2.692   41.499   47.530   -0.221   -2.327   -0.500   -0.006   -0.016   -0.021    2.000   12.000    0.000    0.000    1.000 
 +107    3.856    2.563   32.529   28.804   -0.609   -1.507   -0.500   -0.026   -0.012   -0.038    2.000   12.000    0.000    0.000    1.000 
 +108    3.856    3.061   35.061   36.000    0.286   -3.286   -0.500   -0.070   -0.017   -0.008    2.000   12.000    0.000    0.000    1.000 
 +109    3.856    3.469   35.469   36.000   -0.429   -1.571   -0.500   -0.011   -0.016   -0.034    2.000   12.000    0.000    0.000    1.000 
 +110    3.856    2.695   44.076   42.890   -0.302   -1.811   -0.500   -0.006   -0.015   -0.028    2.000   12.000    0.000    0.000    1.000 
 +111    3.856    2.804   36.778   33.313    0.229   -2.834   -0.500   -0.046   -0.017   -0.013    2.000   12.000    0.000    0.000    1.000 
 +112    3.856    2.706   29.790   32.717   -0.668   -1.005   -0.500   -0.015   -0.011   -0.047    2.000   12.000    0.000    0.000    1.000 
 +113    3.856    3.289   37.322   45.646    0.351   -2.894   -0.500   -0.063   -0.015   -0.014    2.000   12.000    0.000    0.000    1.000 
 +114    3.856    2.449   42.449   44.000   -0.143   -1.857   -0.500    0.030   -0.004   -0.029    2.000   12.000    0.000    0.000    1.000 
 +115    3.856    2.383   45.615   47.377    0.003   -2.130   -0.500    0.009   -0.009   -0.024    2.000   12.000    0.000    0.000    1.000 
 +116    3.856    3.229   34.760   28.746   -0.506   -0.950   -0.500   -0.006   -0.008   -0.046    2.000   12.000    0.000    0.000    1.000 
 +117    3.856    2.860   42.571   48.796   -0.368   -1.207   -0.500   -0.002   -0.015   -0.037    2.000   12.000    0.000    0.000    1.000 
 +118    3.856    2.449   42.449   44.000    0.143   -2.143   -0.500    0.003   -0.019   -0.024    2.000   12.000    0.000    0.000    1.000 
 +119    3.856    2.783   38.282   42.370    0.309   -2.473   -0.500   -0.044   -0.015   -0.019    2.000   12.000    0.000    0.000    1.000 
 +120    3.856    2.427   23.042   21.531   -0.765   -0.262   -0.500    0.051   -0.002   -0.082    2.000   12.000    0.000    0.000    1.000 
 +121    3.856    2.844   43.111   37.319   -0.204   -1.371   -0.500    0.024   -0.008   -0.035    2.000   12.000    0.000    0.000    1.000 
 +122    3.856    3.959   31.959   32.000   -0.571   -0.429   -0.500   -0.003   -0.001   -0.057    2.000   12.000    0.000    0.000    1.000 
 +123    3.856    3.469   35.469   36.000    0.429   -2.429   -0.500   -0.068   -0.014   -0.017    2.000   12.000    0.000    0.000    1.000 
 +124    3.856    2.777   41.114   40.523   -0.455   -0.584   -0.500   -0.002   -0.008   -0.045    2.000   12.000    0.000    0.000    1.000 
 +125    3.856    2.869   21.904   24.650   -0.811    0.135   -0.500    0.073    0.001   -0.095    2.000   12.000    0.000    0.000    1.000 
 +126    3.856    2.987   42.720   50.425    0.204   -1.777   -0.500    0.004   -0.017   -0.030    2.000   12.000    0.000    0.000    1.000 
 +127    3.856    3.224   43.224   44.000   -0.286   -0.714   -0.500    0.022   -0.015   -0.046    2.000   12.000    0.000    0.000    1.000 
 +128    3.856    2.840   37.264   33.663    0.376   -2.005   -0.500   -0.045   -0.015   -0.023    2.000   12.000    0.000    0.000    1.000 
 +129    3.856    2.550   48.826   49.744   -0.151   -0.903   -0.500    0.025   -0.007   -0.037    2.000   12.000    0.000    0.000    1.000 
 +130    3.856    2.548   26.150   25.246    0.018   -1.144   -0.500    0.005   -0.014   -0.030    2.000   12.000    0.000    0.000    1.000 
 +131    3.856    2.174   26.044   26.617   -0.013   -1.077   -0.500    0.020   -0.008   -0.033    2.000   12.000    0.000    0.000    1.000 
 +132    3.856    2.946   38.501   44.577   -0.512   -0.074   -0.500    0.000   -0.005   -0.053    2.000   12.000    0.000    0.000    1.000 
 +133    3.856    3.086   25.220   21.437   -0.663    0.308   -0.500    0.009   -0.007   -0.060    2.000   12.000    0.000    0.000    1.000 
 +134    3.856    3.247   35.644   43.356    0.503   -1.979   -0.500   -0.047   -0.013   -0.028    2.000   12.000    0.000    0.000    1.000 
 +135    3.856    2.678   47.530   49.643    0.153   -1.227   -0.500    0.014   -0.011   -0.034    2.000   12.000    0.000    0.000    1.000 
 +136    3.856    3.060   42.780   36.893   -0.362   -0.107   -0.500    0.011   -0.010   -0.047    2.000   12.000    0.000    0.000    1.000 
 +137    3.856    4.531   24.531   24.000   -0.714    0.714   -0.500    0.014   -0.008   -0.063    2.000   12.000    0.000    0.000    1.000 
 +138    3.856    3.224   43.224   44.000    0.286   -1.286   -0.500    0.008   -0.018   -0.038    2.000   12.000    0.000    0.000    1.000 
 +139    3.856    2.838   45.722   53.237   -0.212   -0.270   -0.500    0.029   -0.009   -0.045    2.000   12.000    0.000    0.000    1.000 
 +140    3.856    2.781   36.663   40.104    0.460   -1.569   -0.500   -0.032   -0.012   -0.032    2.000   12.000    0.000    0.000    1.000 
 +141    3.856    2.658   31.643   33.250   -0.611    0.661   -0.500   -0.001   -0.008   -0.053    2.000   12.000    0.000    0.000    1.000 
 +142    3.856    3.551   39.551   40.000   -0.429    0.429   -0.500    0.008   -0.007   -0.051    2.000   12.000    0.000    0.000    1.000 
 +143    3.856    3.959   31.959   32.000    0.571   -1.571   -0.500   -0.041   -0.012   -0.034    2.000   12.000    0.000    0.000    1.000 
 +144    3.856    2.861   45.146   39.777    0.215   -0.717   -0.500    0.020   -0.011   -0.041    2.000   12.000    0.000    0.000    1.000 
 +145    3.856    2.835   46.068   47.377   -0.304    0.324   -0.500    0.011   -0.010   -0.042    2.000   12.000    0.000    0.000    1.000 
 +146    3.856    3.108   30.570   36.447   -0.656    1.064   -0.500   -0.003   -0.009   -0.054    2.000   12.000    0.000    0.000    1.000 
 +147    3.856    2.571   46.571   48.000   -0.143    0.143   -0.500    0.046   -0.003   -0.048    2.000   12.000    0.000    0.000    1.000 
 +148    3.856    3.082   43.633   51.571    0.359   -0.847   -0.500    0.006   -0.012   -0.044    2.000   12.000    0.000    0.000    1.000 
 +149    3.856    2.923   33.166   29.522    0.520   -1.159   -0.500   -0.026   -0.011   -0.038    2.000   12.000    0.000    0.000    1.000 
 +150    3.856    2.466   50.361   52.010   -0.000   -0.000   -0.500    0.035   -0.003   -0.040    2.000   12.000    0.000    0.000    1.000 
 +151    3.857    2.923   33.166   29.522   -0.520    1.159   -0.500   -0.026   -0.011   -0.038    2.000   12.000    0.000    0.000    1.000 
 +152    3.857    3.082   43.633   51.571   -0.359    0.847   -0.500    0.006   -0.012   -0.044    2.000   12.000    0.000    0.000    1.000 
 +153    3.857    2.571   46.571   48.000    0.143   -0.143   -0.500    0.046   -0.003   -0.048    2.000   12.000    0.000    0.000    1.000 
 +154    3.857    3.108   30.571   36.448    0.656   -1.064   -0.500   -0.003   -0.009   -0.054    2.000   12.000    0.000    0.000    1.000 
 +155    3.857    2.835   46.068   47.377    0.304   -0.324   -0.500    0.011   -0.010   -0.042    2.000   12.000    0.000    0.000    1.000 
 +156    3.857    2.861   45.146   39.777   -0.215    0.717   -0.500    0.020   -0.011   -0.041    2.000   12.000    0.000    0.000    1.000 
 +157    3.857    3.551   39.551   40.000    0.429   -0.429   -0.500    0.008   -0.007   -0.051    2.000   12.000    0.000    0.000    1.000 
 +158    3.857    3.959   31.959   32.000   -0.571    1.571   -0.500   -0.041   -0.012   -0.034    2.000   12.000    0.000    0.000    1.000 
 +159    3.857    2.658   31.643   33.250    0.611   -0.661   -0.500   -0.001   -0.008   -0.053    2.000   12.000    0.000    0.000    1.000 
 +160    3.857    2.781   36.663   40.104   -0.460    1.569   -0.500   -0.032   -0.012   -0.032    2.000   12.000    0.000    0.000    1.000 
 +161    3.857    2.838   45.722   53.237    0.212    0.270   -0.500    0.029   -0.009   -0.045    2.000   12.000    0.000    0.000    1.000 
 +162    3.857    4.531   24.531   24.000    0.714   -0.714   -0.500    0.014   -0.008   -0.063    2.000   12.000    0.000    0.000    1.000 
 +163    3.857    3.224   43.224   44.000   -0.286    1.286   -0.500    0.008   -0.018   -0.038    2.000   12.000    0.000    0.000    1.000 
 +164    3.857    3.060   42.780   36.893    0.362    0.107   -0.500    0.011   -0.010   -0.047    2.000   12.000    0.000    0.000    1.000 
 +165    3.857    2.678   47.530   49.643   -0.153    1.227   -0.500    0.014   -0.011   -0.034    2.000   12.000    0.000    0.000    1.000 
 +166    3.857    3.247   35.644   43.356   -0.503    1.979   -0.500   -0.047   -0.013   -0.028    2.000   12.000    0.000    0.000    1.000 
 +167    3.857    3.086   25.219   21.437    0.663   -0.308   -0.500    0.009   -0.007   -0.060    2.000   12.000    0.000    0.000    1.000 
 +168    3.857    2.946   38.501   44.578    0.512    0.074   -0.500    0.000   -0.005   -0.053    2.000   12.000    0.000    0.000    1.000 
 +169    3.857    2.174   26.044   26.617    0.013    1.076   -0.500    0.020   -0.008   -0.033    2.000   12.000    0.000    0.000    1.000 
 +170    3.857    2.548   26.150   25.246   -0.018    1.144   -0.500    0.005   -0.014   -0.030    2.000   12.000    0.000    0.000    1.000 
 +171    3.857    2.550   48.826   49.744    0.151    0.903   -0.500    0.025   -0.007   -0.037    2.000   12.000    0.000    0.000    1.000 
 +172    3.857    2.840   37.265   33.663   -0.376    2.005   -0.500   -0.045   -0.015   -0.023    2.000   12.000    0.000    0.000    1.000 
 +173    3.857    3.224   43.224   44.000    0.286    0.714   -0.500    0.022   -0.015   -0.046    2.000   12.000    0.000    0.000    1.000 
 +174    3.857    2.987   42.720   50.426   -0.204    1.777   -0.500    0.004   -0.017   -0.030    2.000   12.000    0.000    0.000    1.000 
 +175    3.857    2.869   21.904   24.650    0.811   -0.135   -0.500    0.073    0.001   -0.095    2.000   12.000    0.000    0.000    1.000 
 +176    3.857    2.777   41.114   40.523    0.455    0.584   -0.500   -0.002   -0.008   -0.045    2.000   12.000    0.000    0.000    1.000 
 +177    3.857    3.959   31.959   32.000    0.571    0.429   -0.500   -0.003   -0.001   -0.057    2.000   12.000    0.000    0.000    1.000 
 +178    3.857    3.469   35.469   36.000   -0.429    2.429   -0.500   -0.068   -0.014   -0.017    2.000   12.000    0.000    0.000    1.000 
 +179    3.857    2.844   43.111   37.319    0.204    1.371   -0.500    0.024   -0.008   -0.035    2.000   12.000    0.000    0.000    1.000 
 +180    3.857    2.427   23.042   21.531    0.765    0.262   -0.500    0.051   -0.002   -0.082    2.000   12.000    0.000    0.000    1.000 
 +181    3.857    2.783   38.283   42.370   -0.309    2.472   -0.500   -0.044   -0.015   -0.019    2.000   12.000    0.000    0.000    1.000 
 +182    3.857    2.449   42.449   44.000   -0.143    2.143   -0.500    0.003   -0.019   -0.024    2.000   12.000    0.000    0.000    1.000 
 +183    3.857    2.860   42.571   48.796    0.368    1.207   -0.500   -0.002   -0.015   -0.037    2.000   12.000    0.000    0.000    1.000 
 +184    3.857    3.229   34.759   28.746    0.506    0.950   -0.500   -0.006   -0.008   -0.046    2.000   12.000    0.000    0.000    1.000 
 +185    3.857    2.383   45.616   47.377   -0.003    2.130   -0.500    0.009   -0.009   -0.024    2.000   12.000    0.000    0.000    1.000 
 +186    3.857    2.449   42.449   44.000    0.143    1.857   -0.500    0.030   -0.004   -0.029    2.000   12.000    0.000    0.000    1.000 
 +187    3.857    3.289   37.322   45.646   -0.351    2.894   -0.500   -0.063   -0.015   -0.014    2.000   12.000    0.000    0.000    1.000 
 +188    3.857    2.706   29.790   32.718    0.668    1.005   -0.500   -0.015   -0.011   -0.047    2.000   12.000    0.000    0.000    1.000 
 +189    3.857    2.804   36.778   33.313   -0.229    2.834   -0.500   -0.046   -0.017   -0.013    2.000   12.000    0.000    0.000    1.000 
 +190    3.857    2.695   44.076   42.890    0.302    1.811   -0.500   -0.006   -0.015   -0.028    2.000   12.000    0.000    0.000    1.000 
 +191    3.857    3.469   35.469   36.000    0.429    1.571   -0.500   -0.011   -0.016   -0.034    2.000   12.000    0.000    0.000    1.000 
 +192    3.857    3.061   35.061   36.000   -0.286    3.286   -0.500   -0.070   -0.017   -0.008    2.000   12.000    0.000    0.000    1.000 
 +193    3.857    2.563   32.529   28.804    0.609    1.507   -0.500   -0.026   -0.012   -0.038    2.000   12.000    0.000    0.000    1.000 
 +194    3.857    2.692   41.499   47.530    0.221    2.327   -0.500   -0.006   -0.016   -0.021    2.000   12.000    0.000    0.000    1.000 
 +195    3.857    3.256   37.771   31.068    0.351    2.193   -0.500   -0.017   -0.019   -0.024    2.000   12.000    0.000    0.000    1.000 
 +196    3.857    2.608   36.490   40.104   -0.159    3.376   -0.500   -0.036   -0.016   -0.010    2.000   12.000    0.000    0.000    1.000 
 +197    3.857    2.625   33.785   36.824    0.524    2.144   -0.500   -0.059   -0.014   -0.021    2.000   12.000    0.000    0.000    1.000 
 +198    3.857    3.061   35.061   36.000    0.286    2.714   -0.500   -0.025   -0.025   -0.015    2.000   12.000    0.000    0.000    1.000 
 +199    3.857    2.108   21.593   22.079    0.039    3.234   -0.500   -0.004   -0.005   -0.013    2.000   12.000    0.000    0.000    1.000 
 +200    3.857    2.162   40.499   40.523    0.148    3.038   -0.500   -0.006   -0.008   -0.014    2.000   12.000    0.000    0.000    1.000 
 +201    3.857    2.650   20.992   20.266   -0.057    3.452   -0.500   -0.032   -0.017   -0.010    2.000   12.000    0.000    0.000    1.000 
 +202    3.857    3.133   32.839   39.589   -0.196    3.821   -0.500   -0.050   -0.019   -0.006    2.000   12.000    0.000    0.000    1.000 
 +203    3.857    2.579   35.589   31.171    0.456    2.734   -0.500   -0.053   -0.015   -0.016    2.000   12.000    0.000    0.000    1.000 
 +204    3.857    2.755   31.719   26.750    0.194    3.451   -0.500   -0.008   -0.002   -0.009    2.000   12.000    0.000    0.000    1.000 
 +205    3.857    2.082   30.082   32.000   -0.143    4.143   -0.500   -0.048   -0.021   -0.003    2.000   12.000    0.000    0.000    1.000 
 +206    3.857    2.618   33.901   36.997    0.380    3.279   -0.500   -0.061   -0.017   -0.008    2.000   12.000    0.000    0.000    1.000 
 +207    3.857    2.082   30.082   32.000    0.143    3.857   -0.500   -0.003    0.005   -0.005    2.000   12.000    0.000    0.000    1.000 
 +208    3.857    2.166   31.151   33.250   -0.007    4.284   -0.500   -0.010   -0.015   -0.005    2.000   12.000    0.000    0.000    1.000 
 +209    3.857    2.378   32.344   28.804    0.302    3.961   -0.500   -0.033   -0.018   -0.007    2.000   12.000    0.000    0.000    1.000 
 +210    3.857    2.561   29.582   32.630    0.234    4.404   -0.500   -0.024   -0.019   -0.003    2.000   12.000    0.000    0.000    1.000 
 +211    3.857    1.315   21.930   21.531    0.147    5.207   -0.500    0.038   -0.008   -0.002    2.000   12.000    0.000    0.000    1.000 
 +212    3.857    1.911   12.134   12.313    0.069    5.412   -0.500    0.057   -0.003   -0.002    2.000   12.000    0.000    0.000    1.000 
 +213    3.892    3.026   13.026   12.000   -0.425   -1.580   -0.407   -0.006   -0.009   -0.034    1.995   12.005    0.000    0.005    0.995 
 +214    3.892    2.903   20.903   20.000   -0.255   -0.948   -0.407   -0.006   -0.012   -0.032    1.995   12.005    0.000    0.005    0.995 
 +215    3.892    2.842   24.842   24.000   -0.085   -0.316   -0.407   -0.017   -0.014   -0.026    1.995   12.005    0.000    0.005    0.995 
 +216    3.892    2.842   24.842   24.000    0.085    0.316   -0.407   -0.017   -0.014   -0.026    1.995   12.005    0.000    0.005    0.995 
 +217    3.892    2.903   20.903   20.000    0.255    0.948   -0.407   -0.006   -0.012   -0.032    1.995   12.005    0.000    0.005    0.995 
 +218    3.892    3.026   13.026   12.000    0.425    1.580   -0.407   -0.006   -0.009   -0.034    1.995   12.005    0.000    0.005    0.995 
 +219    3.900    2.239   11.073   12.368   -0.452   -2.709   -0.394   -0.005   -0.011   -0.028    1.994   12.006    0.000    0.006    0.994 
 +220    3.900    2.478   19.683   24.087   -0.299   -1.795   -0.394   -0.005   -0.011   -0.028    1.994   12.006    0.000    0.006    0.994 
 +221    3.900    2.618   24.719   30.941   -0.149   -0.895   -0.394   -0.006   -0.011   -0.028    1.994   12.006    0.000    0.006    0.994 
 +222    3.900    2.664   26.383   33.207   -0.000    0.000   -0.394   -0.006   -0.011   -0.028    1.994   12.006    0.000    0.006    0.994 
 +223    3.900    2.618   24.719   30.941    0.149    0.895   -0.394   -0.006   -0.011   -0.028    1.994   12.006    0.000    0.006    0.994 
 +224    3.900    2.478   19.683   24.087    0.299    1.795   -0.394   -0.005   -0.011   -0.028    1.994   12.006    0.000    0.006    0.994 
 +225    3.900    2.239   11.073   12.368    0.452    2.710   -0.394   -0.005   -0.011   -0.028    1.994   12.006    0.000    0.006    0.994 
 +226    3.922    2.233   12.233   12.000    0.405   -2.549   -0.592   -0.006   -0.010   -0.032    1.989   12.011    0.000    0.011    0.989 
 +227    3.922    2.086   20.086   20.000    0.243   -1.529   -0.592   -0.006   -0.011   -0.031    1.989   12.011    0.000    0.011    0.989 
 +228    3.922    2.012   24.012   24.000    0.081   -0.510   -0.592   -0.014   -0.013   -0.027    1.989   12.011    0.000    0.011    0.989 
 +229    3.922    2.012   24.012   24.000   -0.081    0.510   -0.592   -0.014   -0.013   -0.027    1.989   12.011    0.000    0.011    0.989 
 +230    3.922    2.085   20.085   20.000   -0.243    1.529   -0.592   -0.006   -0.011   -0.031    1.989   12.011    0.000    0.011    0.989 
 +231    3.922    2.233   12.233   12.000   -0.405    2.549   -0.592   -0.006   -0.010   -0.032    1.989   12.011    0.000    0.011    0.989 
 +232    3.929    1.692   13.472    9.162    0.306   -2.450   -0.581   -0.005   -0.011   -0.028    1.988   12.012    0.000    0.012    0.988 
 +233    3.929    1.840   22.972   16.436    0.152   -1.216   -0.581   -0.006   -0.011   -0.028    1.988   12.012    0.000    0.012    0.988 
 +234    3.929    1.888   26.064   18.803    0.000    0.000   -0.581   -0.006   -0.011   -0.028    1.988   12.012    0.000    0.012    0.988 
 +235    3.929    1.840   22.972   16.436   -0.152    1.216   -0.581   -0.006   -0.011   -0.028    1.988   12.012    0.000    0.012    0.988 
 +236    3.929    1.692   13.472    9.163   -0.306    2.450   -0.581   -0.005   -0.011   -0.028    1.988   12.012    0.000    0.012    0.988 
 +237    3.973    1.822    1.822    0.000    0.000    0.000   -0.454   -0.008   -0.011   -0.026    1.981   12.019    0.001    0.017    0.982 
 +238    4.161    1.911   12.134   12.313   -0.783   -4.698    3.000   -0.076   -0.034   -0.003    2.000   12.000    0.000    0.000    1.000 
 +239    4.161    2.116   16.167   19.671   -0.704   -4.226    3.000   -0.069   -0.033   -0.004    2.000   12.000    0.000    0.000    1.000 
 +240    4.161    2.561   29.582   32.630   -0.663   -3.976    3.000   -0.077   -0.030   -0.005    2.000   12.000    0.000    0.000    1.000 
 +241    4.161    2.627   34.586   44.743   -0.607   -3.642    3.000   -0.076   -0.024   -0.007    2.000   12.000    0.000    0.000    1.000 
 +242    4.161    3.612   16.612   15.000   -0.571   -3.429    3.000   -0.079   -0.004   -0.014    2.000   12.000    0.000    0.000    1.000 
 +243    4.161    2.618   33.901   36.996   -0.523   -3.136    3.000   -0.075   -0.022   -0.008    2.000   12.000    0.000    0.000    1.000 
 +244    4.161    3.133   32.839   39.588   -0.518   -3.107    3.000   -0.075   -0.010   -0.012    2.000   12.000    0.000    0.000    1.000 
 +245    4.161    2.108   21.593   22.079   -0.468   -2.806    3.000   -0.056   -0.012   -0.012    2.000   12.000    0.000    0.000    1.000 
 +246    4.161    2.680   36.494   47.339   -0.458   -2.745    3.000   -0.067   -0.018   -0.009    2.000   12.000    0.000    0.000    1.000 
 +247    4.161    3.469   35.469   36.000   -0.429   -2.571    3.000   -0.079   -0.020   -0.008    2.000   12.000    0.000    0.000    1.000 
 +248    4.161    2.625   33.784   36.824   -0.381   -2.287    3.000   -0.064   -0.012   -0.012    2.000   12.000    0.000    0.000    1.000 
 +249    4.161    2.692   41.499   47.530   -0.364   -2.184    3.000   -0.043   -0.021   -0.011    2.000   12.000    0.000    0.000    1.000 
 +250    4.161    3.289   37.322   45.646   -0.363   -2.180    3.000   -0.049   -0.026   -0.008    2.000   12.000    0.000    0.000    1.000 
 +251    4.161    2.246   20.866   26.068   -0.350   -2.102    3.000   -0.049   -0.013   -0.009    2.000   12.000    0.000    0.000    1.000 
 +252    4.161    2.628   34.616   44.783   -0.307   -1.840    3.000   -0.045   -0.011   -0.013    2.000   12.000    0.000    0.000    1.000 
 +253    4.161    2.916   44.969   58.875   -0.302   -1.810    3.000   -0.030   -0.023   -0.011    2.000   12.000    0.000    0.000    1.000 
 +254    4.161    3.061   35.061   36.000   -0.286   -1.714    3.000   -0.011   -0.039   -0.008    2.000   12.000    0.000    0.000    1.000 
 +255    4.161    2.755   21.755   21.000   -0.286   -1.714    3.000    0.011   -0.023   -0.016    2.000   12.000    0.000    0.000    1.000 
 +256    4.161    2.706   29.790   32.717   -0.239   -1.434    3.000   -0.038   -0.008   -0.015    2.000   12.000    0.000    0.000    1.000 
 +257    4.161    2.860   42.571   48.796   -0.225   -1.350    3.000   -0.027   -0.029   -0.008    2.000   12.000    0.000    0.000    1.000 
 +258    4.161    2.987   42.720   50.425   -0.225   -1.349    3.000   -0.013   -0.021   -0.012    2.000   12.000    0.000    0.000    1.000 
 +259    4.161    3.247   35.644   43.356   -0.211   -1.265    3.000   -0.013   -0.029   -0.010    2.000   12.000    0.000    0.000    1.000 
 +260    4.161    2.174   26.044   26.617   -0.156   -0.934    3.000    0.004   -0.012   -0.008    2.000   12.000    0.000    0.000    1.000 
 +261    4.161    2.484   29.419   37.709   -0.155   -0.929    3.000   -0.005   -0.008   -0.013    2.000   12.000    0.000    0.000    1.000 
 +262    4.161    2.883   43.770   57.243   -0.152   -0.911    3.000   -0.016   -0.021   -0.007    2.000   12.000    0.000    0.000    1.000 
 +263    4.161    2.735   30.735   32.000   -0.143   -0.857    3.000   -0.007   -0.029   -0.012    2.000   12.000    0.000    0.000    1.000 
 +264    4.161    2.980   42.980   44.000   -0.143   -0.857    3.000   -0.057   -0.021   -0.010    2.000   12.000    0.000    0.000    1.000 
 +265    4.161    2.869   21.904   24.650   -0.097   -0.579    3.000    0.011   -0.008   -0.014    2.000   12.000    0.000    0.000    1.000 
 +266    4.161    2.946   38.501   44.577   -0.084   -0.502    3.000   -0.003   -0.018   -0.005    2.000   12.000    0.000    0.000    1.000 
 +267    4.161    3.082   43.633   51.571   -0.070   -0.419    3.000   -0.016   -0.011   -0.001    2.000   12.000    0.000    0.000    1.000 
 +268    4.161    2.838   45.722   53.237   -0.069   -0.413    3.000    0.001   -0.005   -0.004    2.000   12.000    0.000    0.000    1.000 
 +269    4.161    3.108   30.571   36.448   -0.058   -0.350    3.000   -0.002    0.006   -0.003    2.000   12.000    0.000    0.000    1.000 
 +270    4.161    2.287   22.344   28.080   -0.000   -0.000    3.000   -0.024   -0.019   -0.010    2.000   12.000    0.000    0.000    1.000 
 +271    4.161    2.732   38.335   49.845   -0.000   -0.000    3.000   -0.004   -0.005   -0.005    2.000   12.000    0.000    0.000    1.000 
 +272    4.161    3.009   48.333   63.453   -0.000   -0.000    3.000    0.001   -0.005   -0.002    2.000   12.000    0.000    0.000    1.000 
 +273    4.161    2.240   20.641   25.762   -0.000   -0.000    3.000    0.063    0.000   -0.006    2.000   12.000    0.000    0.000    1.000 
 +274    4.161    2.490   22.490   24.000    0.000   -0.000    3.000    0.002    0.029    0.003    2.000   12.000    0.000    0.000    1.000 
 +275    4.161    2.816   38.816   40.000    0.000   -0.000    3.000    0.020    0.002    0.011    2.000   12.000    0.000    0.000    1.000 
 +276    4.161    2.469   23.469   23.000    0.000   -0.000    3.000    0.057    0.014    0.003    2.000   12.000    0.000    0.000    1.000 
 +277    4.161    3.108   30.570   36.447    0.058    0.350    3.000   -0.002    0.006   -0.003    2.000   12.000    0.000    0.000    1.000 
 +278    4.161    2.838   45.722   53.237    0.069    0.413    3.000    0.001   -0.005   -0.004    2.000   12.000    0.000    0.000    1.000 
 +279    4.161    3.082   43.633   51.571    0.070    0.418    3.000   -0.016   -0.011   -0.001    2.000   12.000    0.000    0.000    1.000 
 +280    4.161    2.946   38.501   44.578    0.084    0.502    3.000   -0.003   -0.018   -0.005    2.000   12.000    0.000    0.000    1.000 
 +281    4.161    2.869   21.904   24.650    0.097    0.579    3.000    0.011   -0.008   -0.014    2.000   12.000    0.000    0.000    1.000 
 +282    4.161    2.735   30.735   32.000    0.143    0.857    3.000   -0.007   -0.029   -0.012    2.000   12.000    0.000    0.000    1.000 
 +283    4.161    2.980   42.979   44.000    0.143    0.857    3.000   -0.057   -0.021   -0.010    2.000   12.000    0.000    0.000    1.000 
 +284    4.161    2.883   43.770   57.243    0.152    0.911    3.000   -0.016   -0.021   -0.007    2.000   12.000    0.000    0.000    1.000 
 +285    4.161    2.484   29.419   37.710    0.155    0.929    3.000   -0.005   -0.008   -0.013    2.000   12.000    0.000    0.000    1.000 
 +286    4.161    2.174   26.044   26.617    0.156    0.934    3.000    0.004   -0.012   -0.008    2.000   12.000    0.000    0.000    1.000 
 +287    4.161    3.247   35.644   43.356    0.211    1.265    3.000   -0.013   -0.029   -0.010    2.000   12.000    0.000    0.000    1.000 
 +288    4.161    2.987   42.720   50.426    0.225    1.349    3.000   -0.013   -0.021   -0.012    2.000   12.000    0.000    0.000    1.000 
 +289    4.161    2.860   42.571   48.796    0.225    1.350    3.000   -0.027   -0.029   -0.008    2.000   12.000    0.000    0.000    1.000 
 +290    4.161    2.706   29.790   32.718    0.239    1.434    3.000   -0.038   -0.008   -0.015    2.000   12.000    0.000    0.000    1.000 
 +291    4.161    3.061   35.061   36.000    0.286    1.714    3.000   -0.011   -0.039   -0.008    2.000   12.000    0.000    0.000    1.000 
 +292    4.161    2.755   21.755   21.000    0.286    1.714    3.000    0.011   -0.023   -0.016    2.000   12.000    0.000    0.000    1.000 
 +293    4.161    2.916   44.969   58.875    0.302    1.810    3.000   -0.030   -0.023   -0.011    2.000   12.000    0.000    0.000    1.000 
 +294    4.161    2.628   34.616   44.783    0.307    1.840    3.000   -0.045   -0.011   -0.013    2.000   12.000    0.000    0.000    1.000 
 +295    4.161    2.246   20.866   26.068    0.350    2.101    3.000   -0.049   -0.013   -0.009    2.000   12.000    0.000    0.000    1.000 
 +296    4.161    3.289   37.322   45.646    0.363    2.180    3.000   -0.049   -0.026   -0.008    2.000   12.000    0.000    0.000    1.000 
 +297    4.161    2.692   41.499   47.530    0.364    2.184    3.000   -0.043   -0.021   -0.011    2.000   12.000    0.000    0.000    1.000 
 +298    4.161    2.625   33.785   36.824    0.381    2.287    3.000   -0.064   -0.012   -0.012    2.000   12.000    0.000    0.000    1.000 
 +299    4.161    3.469   35.469   36.000    0.429    2.571    3.000   -0.079   -0.020   -0.008    2.000   12.000    0.000    0.000    1.000 
 +300    4.161    2.680   36.494   47.339    0.458    2.745    3.000   -0.067   -0.018   -0.009    2.000   12.000    0.000    0.000    1.000 
 +301    4.161    2.108   21.593   22.079    0.468    2.805    3.000   -0.056   -0.012   -0.012    2.000   12.000    0.000    0.000    1.000 
 +302    4.161    3.133   32.839   39.589    0.518    3.107    3.000   -0.075   -0.010   -0.012    2.000   12.000    0.000    0.000    1.000 
 +303    4.161    2.618   33.901   36.997    0.523    3.136    3.000   -0.075   -0.022   -0.008    2.000   12.000    0.000    0.000    1.000 
 +304    4.161    3.612   16.612   15.000    0.571    3.429    3.000   -0.079   -0.004   -0.014    2.000   12.000    0.000    0.000    1.000 
 +305    4.161    2.627   34.587   44.743    0.607    3.642    3.000   -0.076   -0.024   -0.007    2.000   12.000    0.000    0.000    1.000 
 +306    4.161    2.561   29.582   32.630    0.663    3.976    3.000   -0.077   -0.030   -0.005    2.000   12.000    0.000    0.000    1.000 
 +307    4.161    2.116   16.167   19.672    0.704    4.226    3.000   -0.069   -0.033   -0.004    2.000   12.000    0.000    0.000    1.000 
 +308    4.161    1.911   12.134   12.313    0.783    4.698    3.000   -0.076   -0.034   -0.003    2.000   12.000    0.000    0.000    1.000 
 +309    4.219    1.695   10.529   12.368   -0.309   -1.854    2.899   -0.026   -0.013   -0.010    1.990   12.010    0.000    0.010    0.990 
 +310    4.219    1.934   19.139   24.087   -0.205   -1.228    2.899   -0.026   -0.018   -0.008    1.990   12.010    0.000    0.010    0.990 
 +311    4.219    2.074   24.175   30.941   -0.102   -0.612    2.899   -0.040   -0.015   -0.008    1.990   12.010    0.000    0.010    0.990 
 +312    4.219    2.120   25.840   33.207    0.000    0.000    2.899   -0.047   -0.021   -0.011    1.990   12.010    0.000    0.010    0.990 
 +313    4.219    2.074   24.175   30.941    0.102    0.613    2.899   -0.040   -0.015   -0.008    1.990   12.010    0.000    0.010    0.990 
 +314    4.219    1.934   19.139   24.087    0.205    1.228    2.899   -0.026   -0.018   -0.008    1.990   12.010    0.000    0.010    0.990 
 +315    4.219    1.695   10.529   12.368    0.309    1.854    2.899   -0.026   -0.013   -0.010    1.990   12.010    0.000    0.010    0.990 
 +316    4.226    2.233   12.233   12.000   -0.309   -1.834    2.908   -0.033   -0.016   -0.009    1.989   12.011    0.000    0.011    0.989 
 +317    4.226    2.086   20.086   20.000   -0.186   -1.100    2.908   -0.033   -0.016   -0.009    1.989   12.011    0.000    0.011    0.989 
 +318    4.226    2.012   24.012   24.000   -0.062   -0.367    2.908   -0.033   -0.016   -0.009    1.989   12.011    0.000    0.011    0.989 
 +319    4.226    2.012   24.012   24.000    0.062    0.367    2.908   -0.033   -0.016   -0.009    1.989   12.011    0.000    0.011    0.989 
 +320    4.226    2.085   20.085   20.000    0.186    1.101    2.908   -0.033   -0.016   -0.009    1.989   12.011    0.000    0.011    0.989 
 +321    4.226    2.233   12.233   12.000    0.309    1.835    2.908   -0.033   -0.016   -0.009    1.989   12.011    0.000    0.011    0.989 
 +322    4.296    1.082    1.082    0.000    0.000    0.000    2.763   -0.032   -0.016   -0.009    1.975   12.025    0.001    0.023    0.976 
 +</code> 
 + 
 + 
 +</WRAP>
  
-{{namespace>forum:data&link&noeditbtn&firstsectiononly&readmore&order=modified&depth=2}} 

You are here: Quanty - a quantum many body script language » Forum

Print/export