NRG Ljubljana (c) Rok Zitko, rok.zitko@ijs.si, 2005-2018 Mathematica version: 11.3.0 for Linux x86 (64-bit) (March 7, 2018) sneg version: 1.250 Loading module initialparse.m Options: {} Loading module models.m "models started" Loading module custommodels.m models $Id: custommodels.m,v 1.1 2015/11/09 12:23:47 rokzitko Exp rokzitko $ custommodels.m done Loading module kondo.m def1ch, NRDOTS=0 COEFCHANNELS:1 H0=coefzeta[1, 0]*(-1 + nc[f[0, 0, 0], f[1, 0, 0]] + nc[f[0, 0, 1], f[1, 0, 1]]) adddots, nrdots=0 params={gammaPol -> Sqrt[gammaA*theta0]/Sqrt[Pi], gammaPolCh[ch_] :> Sqrt[1/Pi*theta0Ch[ch]*gammaA], hybV[i_, j_] :> Sqrt[1/Pi]*V[i, j], coefzeta[ch_, j__] :> N[bandrescale*zeta[ch][j]], coefxi[ch_, j__] :> N[bandrescale*xi[ch][j]], coefrung[ch_, j__] :> N[bandrescale*zetaR[ch][j]], coefdelta[ch_, j__] :> N[bandrescale*scdelta[ch][j]], coefkappa[ch_, j__] :> N[bandrescale*sckappa[ch][j]], U -> 0.1, delta -> 0., t -> 0., gammaPol2 -> Sqrt[extraGamma2*gammaA*thetaCh[1]]/Sqrt[Pi], gammaPol2to2 -> Sqrt[extraGamma2to2*gammaA*thetaCh[2]]/Sqrt[Pi], gammaPolch1 -> Sqrt[extraGamma1*gammaA*thetaCh[1]]/Sqrt[Pi], gammaPolch2 -> Sqrt[extraGamma2*gammaA*thetaCh[2]]/Sqrt[Pi], gammaPolch3 -> Sqrt[extraGamma3*gammaA*thetaCh[3]]/Sqrt[Pi], Jspin -> extraJspin*gammaA, Jcharge -> extraJcharge*gammaA, Jcharge1 -> extraJcharge1*gammaA, Jcharge2 -> extraJcharge2*gammaA, Jkondo -> 0.2*gammaA, Jkondo1 -> extraJkondo1*gammaA, Jkondo2 -> extraJkondo2*gammaA, Jkondo3 -> extraJkondo3*gammaA, Jkondo1P -> extraJkondo1P*gammaA, Jkondo2P -> extraJkondo2P*gammaA, Jkondo1Z -> extraJkondo1Z*gammaA, Jkondo2Z -> extraJkondo2Z*gammaA, JkondoP -> extraJkondoP*gammaA, JkondoZ -> extraJkondoZ*gammaA, Jkondo1ch2 -> extraJkondo1ch2*gammaA, Jkondo2ch2 -> extraJkondo2ch2*gammaA, gep -> extrag, dd -> extrad, hybV11 -> Sqrt[extraGamma11*gammaA*thetaCh[1]]/Sqrt[Pi], hybV12 -> Sqrt[extraGamma12*gammaA*thetaCh[2]]/Sqrt[Pi], hybV21 -> Sqrt[extraGamma21*gammaA*thetaCh[1]]/Sqrt[Pi], hybV22 -> Sqrt[extraGamma22*gammaA*thetaCh[2]]/Sqrt[Pi], spin -> 1/2, Jkondo -> 0.2, B -> 0.0001} NRDOTS:0 CHANNELS:1 basis:{f[0]} lrchain:{} lrextrarule:{} NROPS:1 Hamiltonian generated. -coefzeta[1, 0] + coefzeta[1, 0]*nc[f[0, 0, 0], f[1, 0, 0]] + coefzeta[1, 0]*nc[f[0, 0, 1], f[1, 0, 1]] - (B*nc[ket[-1/2], bra[-1/2]])/2 + (B*nc[ket[1/2], bra[1/2]])/2 + (Jkondo*nc[f[0, 0, 0], f[1, 0, 0], ket[-1/2], bra[-1/2]])/4 - (Jkondo*nc[f[0, 0, 0], f[1, 0, 0], ket[1/2], bra[1/2]])/4 + (Jkondo*nc[f[0, 0, 0], f[1, 0, 1], ket[1/2], bra[-1/2]])/2 + (Jkondo*nc[f[0, 0, 1], f[1, 0, 0], ket[-1/2], bra[1/2]])/2 - (Jkondo*nc[f[0, 0, 1], f[1, 0, 1], ket[-1/2], bra[-1/2]])/4 + (Jkondo*nc[f[0, 0, 1], f[1, 0, 1], ket[1/2], bra[1/2]])/4 H-conj[H]=0 SCALE[0]=1.0510537250399226 faktor=1.6479184330021646 Generating basis Basis states generated. BASIS NR=4 MAKESPINKET=1/2 Basis: basis.model..QSZ PREC=1000 Tmin=1.*^-10 Tmin=1.*^-10 ==> Nmax=42 DISCNMAX=42 mMAX=84 Diagonalisation. Discretization checksum [-1] (channel 1): 2.784154609541307270170357727625724949`10.*^-41 BAND="flat" thetaCh={"2."} Discretization (channel 1) "xitable" (channel 1) 0.5049098776 0.3180107496 0.195230373 0.1153712823 0.06714955356 0.0388746386 0.02246477381 0.01297399467 0.007491300293 0.004325250724 0.002497212862 0.001441771944 0.0008324084643 0.0004805914519 0.0002774696428 0.0001601971804 0.00009248988666 0.00005339906124 0.00003082996243 0.00001779968712 0.00001027665415 5.933229041e-6 3.425551384e-6 1.977743014e-6 1.141850461e-6 6.592476713e-7 3.806168205e-7 2.197492238e-7 1.268722735e-7 7.324974125e-8 4.229075783e-8 2.441658042e-8 1.409691928e-8 8.138860139e-9 4.698973092e-9 2.71295338e-9 1.566324364e-9 9.043177933e-10 5.221081214e-10 3.014392644e-10 1.740360405e-10 1.004797548e-10 5.801201349e-11 "zetatable" (channel 1) 0.e-999 0.e-998 0.e-998 0.e-997 0.e-996 0.e-995 0.e-994 0.e-993 0.e-992 0.e-991 0.e-991 0.e-990 0.e-989 0.e-988 0.e-987 0.e-986 0.e-985 0.e-984 0.e-984 0.e-983 0.e-982 0.e-981 0.e-980 0.e-979 0.e-978 0.e-977 0.e-976 0.e-976 0.e-975 0.e-974 0.e-973 0.e-972 0.e-971 0.e-970 0.e-969 0.e-968 0.e-968 0.e-967 0.e-966 0.e-965 0.e-964 0.e-963 0.e-962 Precision last xi:952.2651260013099 Precision last zeta: 0. Discretization done. --EOF-- {{# Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, 2005-2015}, {# symtype , QSZ}, {# Using sneg version , 1.250}, {#!8}, {# Number of channels, impurities, chain sites, subspaces: }, {1, 0, 42, 7}} maketable[] exnames={B, d, g, Gamma1, Gamma11, Gamma12, Gamma2, Gamma21, Gamma22, Gamma2to2, Gamma3, Jcharge, Jcharge1, Jcharge2, Jkondo, Jkondo1, Jkondo1ch2, Jkondo1P, Jkondo1Z, Jkondo2, Jkondo2ch2, Jkondo2P, Jkondo2Z, Jkondo3, JkondoP, JkondoZ, Jspin, spin} thetaCh={"2."} theta0Ch={"0.2"} gammaPolCh={"0.252313252202016"} checkdefinitions[] -> 0.2 calcgsenergy[] diagvc[{-1, 0}] Generating matrix: ham.model..QSZ_-1.0 hamil={{-B/2 - coefzeta[1, 0]}} dim={1, 1} det[vec]=1. 1-abs=0. orthogonality check=0. diagvc[{-1, 2}] Generating matrix: ham.model..QSZ_-1.2 hamil={{(B - 2*coefzeta[1, 0])/2}} dim={1, 1} det[vec]=1. 1-abs=0. orthogonality check=0. diagvc[{0, -1}] Generating matrix: ham.model..QSZ_0.-1 hamil={{(-2*B + Jkondo)/4}} dim={1, 1} det[vec]=1. 1-abs=0. orthogonality check=0. diagvc[{0, 1}] Generating matrix: ham.model..QSZ_0.1 hamil={{(2*B - Jkondo)/4, Jkondo/2}, {Jkondo/2, (-2*B - Jkondo)/4}} dim={2, 2} det[vec]=-0.9999999999999999 1-abs=1.1102230246251565*^-16 orthogonality check=0. diagvc[{0, 3}] Generating matrix: ham.model..QSZ_0.3 hamil={{(2*B + Jkondo)/4}} dim={1, 1} det[vec]=1. 1-abs=0. orthogonality check=0. diagvc[{1, 0}] Generating matrix: ham.model..QSZ_1.0 hamil={{-B/2 + coefzeta[1, 0]}} dim={1, 1} det[vec]=1. 1-abs=0. orthogonality check=0. diagvc[{1, 2}] Generating matrix: ham.model..QSZ_1.2 hamil={{B/2 + coefzeta[1, 0]}} dim={1, 1} det[vec]=1. 1-abs=0. orthogonality check=0. Lowest energies (absolute):{-0.15000001249999922, -0.00005, -0.00005, 0.00005, 0.00005, 0.04995, 0.050000012499999226, 0.050050000000000004} Lowest energies (GS shifted):{0., 0.14995001249999923, 0.14995001249999923, 0.15005001249999922, 0.15005001249999922, 0.19995001249999922, 0.20000002499999844, 0.20005001249999924} Scale factor SCALE(Ninit):1.0510537250399226 Lowest energies (shifted and scaled):{0., 0.14266636321973325, 0.14266636321973325, 0.14276150583482286, 0.14276150583482286, 0.19023767076455053, 0.19028525396492152, 0.1903328133796402} makeireducf GENERAL ireducTable: f[0]{} Loading module operators.m "operators.m started" s: SZ op.model..QSZ.SZ (-nc[ket[-1/2], bra[-1/2]] + nc[ket[1/2], bra[1/2]])/2 operators.m done Loading module customoperators.m "customoperators $Id: customoperators.m,v 1.1 2015/11/09 12:23:54 rokzitko Exp rokzitko $" Customoperators done. Loading module modeloperators.m Can't load modeloperators.m. Continuing. -- maketable[] done -- Timing report {basis, 0.002976`3.925177920369814} {ham, 0.0267819999999999999`4.034289959187249} {maketable, 0.500097`6.150599242789984} {xi, 0.502004`6.152252171141006} {_, 0} data gammaPol=0.252313252202016 "Success!"