Differences

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

Link to this comparison view

Next revision		Previous revision
Next revisionBoth sides next revision
forum:data:2020:stranger_things [2020/01/07 18:57] – Created from the form at forum:start Riccardo Piomboforum:data:2020:stranger_things [2020/01/07 19:25] – Riccardo Piombo
Line 9: Line 9:
  
 ---------------------------------------------------------------------------------------------- ----------------------------------------------------------------------------------------------
-"Verbosity(0)+<code> 
 + 
 +Verbosity(0)
  
 NF=20 NF=20
Line 30: Line 32:
 TdLb1g = 2.6     TdLb1g = 2.6    
  
--- Racah parameters [1]+-- Racah parameters 
 B = 0.09 B = 0.09
 C = 0.54 C = 0.54
  
 -- Slater-Koster integrals for monopole and  -- Slater-Koster integrals for monopole and 
--- multipole part of Coulomb interaction [2]+-- multipole part of Coulomb interaction 
 F2dd = 49.0*B + 7*C F2dd = 49.0*B + 7*C
 F4dd = 441.0*C/35.0 F4dd = 441.0*C/35.0
Line 113: Line 115:
 print("") print("")
 print("N-body states energies computation") print("N-body states energies computation")
-Npsi_i = 2+Npsi_i = 1
 StartRestrictions1 = {NF, NB, {"1111111111 0000000000",n4d,n4d}, {"0000000000 1111111111",nL,nL}} StartRestrictions1 = {NF, NB, {"1111111111 0000000000",n4d,n4d}, {"0000000000 1111111111",nL,nL}}
 psiList_N = Eigensystem(H_tot, StartRestrictions1, Npsi_i) psiList_N = Eigensystem(H_tot, StartRestrictions1, Npsi_i)
---psiList_N = {psiList_N}+psiList_N = {psiList_N}
 print("Done") print("Done")
  
Line 134: Line 136:
 print("") print("")
 print("(N-1)-body states energies computation") print("(N-1)-body states energies computation")
-StartRestrictions_rem = {NF, NB, {"1111111111 0000000000",n4d-1,n4d-1}, {"0000000000 1111111111",nL,nL}}+StartRestrictions_rem = {NF, NB, {"1111111111 0000000000",8,8}, {"0000000000 1111111111",10,10}}
 psiList_Nminus1 = Eigensystem(H_tot, StartRestrictions_rem, 1) psiList_Nminus1 = Eigensystem(H_tot, StartRestrictions_rem, 1)
 psiList_Nminus1 = {psiList_Nminus1} psiList_Nminus1 = {psiList_Nminus1}
Line 147: Line 149:
 -- chemical potential μ+ to add a particle -- chemical potential μ+ to add a particle
 print("(N+1)-body states energies computation") print("(N+1)-body states energies computation")
-StartRestrictions_add = {NF, NB, {"1111111111 0000000000",n4d+1,n4d+1}, {"0000000000 1111111111",nL,nL}}+StartRestrictions_add = {NF, NB, {"1111111111 0000000000",10,10}, {"0000000000 1111111111",10,10}}
 psiList_Nplus1 = Eigensystem(H_tot, StartRestrictions_add, 1) psiList_Nplus1 = Eigensystem(H_tot, StartRestrictions_add, 1)
 psiList_Nplus1 = {psiList_Nplus1} psiList_Nplus1 = {psiList_Nplus1}
Line 162: Line 164:
 mean_mu = (add_mu + remove_mu)/2.0 mean_mu = (add_mu + remove_mu)/2.0
 print("mean μ = [(μ+) + (μ-)]/2 = ", mean_mu) print("mean μ = [(μ+) + (μ-)]/2 = ", mean_mu)
-"+ 
 +</code> 
 </WRAP> </WRAP>
  
 ~~DISCUSSION|Answers~~ ~~DISCUSSION|Answers~~
  

You are here: Quanty - a quantum many body script language » Forum » data » 2020 » Wrong (missing) ground state Part 2

Print/export