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 ../model.m
def1ch, NRDOTS=1
COEFCHANNELS:2
H0=coefzeta[2, 0]*(-1/2 + nc[f[0, 0, 0], f[1, 0, 0]]) + coefzeta[1, 0]*(-1/2 + nc[f[0, 0, 1], f[1, 0, 1]])
adddots, nrdots=1
"selfopd[CR,UP]="-nc[d[0, 1], epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]])] + nc[epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]]), d[0, 1]] - 0.05*nc[d[0, 0], d[0, 1], d[1, 0]]
"selfopd[CR,DO]="-nc[d[0, 0], epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]])] + nc[epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]]), d[0, 0]] + 0.05*nc[d[0, 0], d[0, 1], d[1, 1]]
"selfopd[AN,UP]="-nc[d[1, 1], epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]])] + nc[epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]]), d[1, 1]] - 0.05*nc[d[0, 0], d[1, 0], d[1, 1]]
"selfopd[AN,DO]="-nc[d[1, 0], epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]])] + nc[epsilon*(nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]]), d[1, 0]] + 0.05*nc[d[0, 1], d[1, 0], d[1, 1]]
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.05, 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 -> extraJkondo*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], 
   U -> 0.05, epsilon -> -0.025, GammaU -> 0.003, GammaD -> 0.05}
NRDOTS:1
CHANNELS:1
basis:{d[], f[0]}
lrchain:{}
lrextrarule:{}
NROPS:2
Hamiltonian generated. -coefzeta[1, 0]/2 - coefzeta[2, 0]/2 + epsilon*nc[d[0, 0], d[1, 0]] + hybV[2, 2]*nc[d[0, 0], f[1, 0, 0]] + epsilon*nc[d[0, 1], d[1, 1]] + hybV[1, 1]*nc[d[0, 1], f[1, 0, 1]] + hybV[2, 2]*nc[f[0, 0, 0], d[1, 0]] + 
   coefzeta[2, 0]*nc[f[0, 0, 0], f[1, 0, 0]] + hybV[1, 1]*nc[f[0, 0, 1], d[1, 1]] + coefzeta[1, 0]*nc[f[0, 0, 1], f[1, 0, 1]] - U*nc[d[0, 0], d[0, 1], d[1, 0], d[1, 1]]
H-conj[H]=0
SCALE[0]=1.715663207404297
faktor=0.8242955588659627
Generating basis
Basis states generated.
BASIS NR=16
Basis: basis.model..QSZ
PREC=30
DISCNMAX=30
mMAX=80
"band=manual_V, importing V, COEFCHANNELS="2
"V[1,1]="-0.07741600514660126886`18.888830756913833
"V[1,2]="0
"V[2,1]="0
"V[2,2]="-0.3160495175564158599`18.49975513172201
Diagonalisation.
Loading discretization data from files.
"nrch="1
"xi="{0.5693244560603186732`18.75535983988729, 0.5412410214981017154`18.733390705026693, 0.5904425846578158499`18.77117767268855, 0.5169751794174348714`18.71346969260882, 0.3416212279788689576`18.533544849451523, 
   0.2217640163312762369`18.345891078438186, 0.1530018487268470506`18.184696678445164, 0.1075851280929436687`18.031752241278692, 0.07593498761202688296`18.880441926431878, 0.05365252349944660964`18.729590153442903, 
   0.03792436582842798748`18.57891832724342, 0.0268118811505536625`18.42832728562199, 0.0189572229469709766`18.27777471764655, 0.01340420216690603981`18.127240969395224, 0.009477997362884266197`18.97671658352724, 
   0.006701883098144649151`18.826196848107806, 0.004738921272614505728`18.675679493853846, 0.003350913917608767251`18.525163271288662, 0.002369450879951075922`18.374647709959955, 0.001675453397051370873`18.224132352383222, 
   0.00118472425043976144`18.07361727823791, 0.0008377261770810412456`18.923102086391697, 0.0005923620180968751587`18.772587204035933, 0.0004188629609934379724`18.622071958481634, 0.0002961810425688772043`18.47155725752345, 
   0.0002094313978023195082`18.321041791335656, 0.0001480905705438362698`18.17052740631896, 0.0001047156244037063141`18.02001148670325, 0.0000740453376059605758`18.869497717607373, 0.00005235773789499322298`18.71898087468269, 
   0.0000370227211208864259`18.56846833565825}
"zeta="{-0.06601070754445456656`18.819614387819083, 0.008624389719958163933`18.93572837323728, -0.008749899663157218593`18.942003072909515, 0.0124899490977765041`18.096560668430193, 0.02591507756878054167`18.41355251310389, 
   0.01276297478892132221`18.105951911363192, 0.005336854111126844956`18.727285331040694, 0.002592821447369992108`18.413772610447758, 0.001325124328335814707`18.122256627387614, 0.0006828556598936225`18.83432891343782, 
   0.0003523244805971346223`18.546942820356595, 0.000181744630905970601`18.259461589806993, 0.00009369145276153558735`18.971699973088263, 0.00004826599692506304952`18.68364128088938, 0.00002484518955543801012`18.395242314527586, 
   0.00001277857918831260114`18.106482568612805, 6.567578600430061392`18.817405279102385*^-6, 3.372943389809536583`18.528009052218064*^-6, 1.731004410828119904`18.23829817451657*^-6, 8.878175727887964662`18.948323736870556*^-7, 
   4.550697493772038891`18.65807796687267*^-7, 2.331057701659095965`18.36755302392898*^-7, 1.193406381699212878`18.07678835588896*^-7, 6.106434375763226941`18.785787694212114*^-8, 3.12280587491866844`18.4945449877294*^-8, 
   1.596278244531300245`18.20310859474168*^-8, 8.155742974209679301`18.911463530661603*^-9, 4.164952525675344739`18.619610055453553*^-9, 2.126005823009399054`18.32756444969687*^-9, 1.084808557536724891`18.0353531024783*^-9, 
   5.532631159419261048`18.742931718331217*^-10}
"nrch="2
"xi="{0.5721845308761048798`18.75753611238927, 0.5408066988610586989`18.73304206248309, 0.5907412525709155871`18.771397299515904, 0.516637841469106962`18.71318621317941, 0.3410615059841057284`18.532832705387815, 
   0.2215637648057046938`18.345498736057838, 0.1529384769595524274`18.184516760945396, 0.1075623157965448773`18.031660143930903, 0.07592645862306388049`18.880393143900978, 0.05364932222800904843`18.7295642397317, 
   0.03792316623041345341`18.578904589715407, 0.02681143292452412008`18.428320025269116, 0.01895705578890664531`18.277770888175052, 0.0134041399145835155`18.127238952423042, 0.009477974227689846265`18.97671552344048, 
   0.006701874517726179661`18.826196292080382, 0.004738918096552468227`18.675679202786192, 0.003350912744489300231`18.525163119246734, 0.002369450447197457392`18.374647630640933, 0.001675453237541592658`18.22413231103667, 
   0.001184724191724331353`18.073617256714094, 0.0008377261554967741029`18.92310207520197, 0.000592362010172679896`18.77258719822625, 0.0004188629580884751895`18.622071955469647, 0.0002961810415048384817`18.471557255963237, 
   0.0002094313974127634992`18.32104179052784, 0.0001480905704013458395`18.170527405901087, 0.000104715624351634591`18.020011486487288, 0.00007404533758694982254`18.86949771749587, 0.00005235773788805991272`18.718980874625178, 
   0.00003702272111835944879`18.56846833562861}
"zeta="{0.03300548326266206089`18.51858609602628, -0.004597739346966279485`18.662544346797784, 0.004393236305534689579`18.642784563936765, -0.006213491104262195037`18.793335680901397, -0.01283082805312386193`18.108254685002983, 
   -0.00631893677241592705`18.800644009814338, -0.002652075390376710824`18.423585865569546, -0.00129167080749406811`18.111151844382494, -0.0006609444096815082336`18.820164933640218, -0.0003408134905320127188`18.532516777255616, 
   -0.0001759077974132691995`18.24528509072808, -0.00009075948804598057149`18.957892037461548, -0.00004679322980610927936`18.670183022507015, -0.00002410780779446445858`18.38215772028815, -0.00001241028395874306986`18.093781718670716, 
   -6.383210662114918721`18.805039177483913*^-6, -3.280770319091711653`18.51597582726425*^-6, -1.684967673734687271`18.22659157330628*^-6, -8.647516514053711398`18.936891400053913*^-7, -4.435338949224263893`18.646926814281148*^-7, 
   -2.273475310990564514`18.356690242306783*^-7, -1.164592401379888615`18.066173952138527*^-7, -5.962350381203530823`18.775417494030062*^-8, -3.050876605848379709`18.48442464275127*^-8, -1.560232699961318192`18.193189375772718*^-8, 
   -7.975540133935943633`18.901760104800626*^-9, -4.074945996127746119`18.610121857557193*^-9, -2.081013513775810823`18.318274900461127*^-9, -1.062271550432422904`18.026235550436017*^-9, -5.420385909475623237`18.734030207643457*^-10, 
   -2.764487191396183192`18.44161458207767*^-10}
BAND="manual_V" thetaCh={"0.005993237853", "0.09988729755"}
Discretization (channel 1)
"xitable" (channel 1)
0.5693244561
0.5412410215
0.5904425847
0.5169751794
0.341621228
0.2217640163
0.1530018487
0.1075851281
0.07593498761
0.0536525235
0.03792436583
0.02681188115
0.01895722295
0.01340420217
0.009477997363
0.006701883098
0.004738921273
0.003350913918
0.00236945088
0.001675453397
0.00118472425
0.0008377261771
0.0005923620181
0.000418862961
0.0002961810426
0.0002094313978
0.0001480905705
0.0001047156244
0.00007404533761
0.00005235773789
0.00003702272112
"zetatable" (channel 1)
-0.06601070754
0.00862438972
-0.008749899663
0.0124899491
0.02591507757
0.01276297479
0.005336854111
0.002592821447
0.001325124328
0.0006828556599
0.0003523244806
0.0001817446309
0.00009369145276
0.00004826599693
0.00002484518956
0.00001277857919
6.5675786e-6
3.37294339e-6
1.731004411e-6
8.878175728e-7
4.550697494e-7
2.331057702e-7
1.193406382e-7
6.106434376e-8
3.122805875e-8
1.596278245e-8
8.155742974e-9
4.164952526e-9
2.126005823e-9
1.084808558e-9
5.532631159e-10
Precision last xi:18.56846833565825
Precision last zeta: 18.742931718331217
Discretization (channel 2)
"xitable" (channel 2)
0.5721845309
0.5408066989
0.5907412526
0.5166378415
0.341061506
0.2215637648
0.152938477
0.1075623158
0.07592645862
0.05364932223
0.03792316623
0.02681143292
0.01895705579
0.01340413991
0.009477974228
0.006701874518
0.004738918097
0.003350912744
0.002369450447
0.001675453238
0.001184724192
0.0008377261555
0.0005923620102
0.0004188629581
0.0002961810415
0.0002094313974
0.0001480905704
0.0001047156244
0.00007404533759
0.00005235773789
0.00003702272112
"zetatable" (channel 2)
0.03300548326
-0.004597739347
0.004393236306
-0.006213491104
-0.01283082805
-0.006318936772
-0.00265207539
-0.001291670807
-0.0006609444097
-0.0003408134905
-0.0001759077974
-0.00009075948805
-0.00004679322981
-0.00002410780779
-0.00001241028396
-6.383210662e-6
-3.280770319e-6
-1.684967674e-6
-8.647516514e-7
-4.435338949e-7
-2.273475311e-7
-1.164592401e-7
-5.962350381e-8
-3.050876606e-8
-1.5602327e-8
-7.975540134e-9
-4.074945996e-9
-2.081013514e-9
-1.06227155e-9
-5.420385909e-10
-2.764487191e-10
Precision last xi:18.56846833562861
Precision last zeta: 18.44161458207767
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, 1, 30, 9}}

maketable[]

exnames={d, epsilon, g, Gamma1, Gamma11, Gamma12, Gamma2, Gamma21, Gamma22, Gamma2to2, Gamma3, GammaD, GammaU, Jcharge, Jcharge1, Jcharge2, Jkondo, Jkondo1, Jkondo1ch2, Jkondo1P, Jkondo1Z, Jkondo2, Jkondo2ch2, Jkondo2P, Jkondo2Z, Jkondo3, 
   JkondoP, JkondoZ, Jspin, U}
thetaCh={"0.005993237853", "0.09988729755"}
theta0Ch={"0.0005993237852858595", "0.009988729754764322"}
gammaPolCh={"0.013811976175825837", "0.056387156617085016"}
checkdefinitions[] -> -0.5604809109293273
calcgsenergy[]
diagvc[{-2, 1}]
Generating matrix: ham.model..QSZ_-2.1
hamil={{(-coefzeta[1, 0] - coefzeta[2, 0])/2}}
dim={1, 1}
det[vec]=1. 1-abs=0.
orthogonality check=0.
diagvc[{-1, 0}]
Generating matrix: ham.model..QSZ_-1.0
hamil={{(-coefzeta[1, 0] + coefzeta[2, 0])/2, hybV[2, 2]}, {hybV[2, 2], epsilon - coefzeta[1, 0]/2 - coefzeta[2, 0]/2}}
dim={2, 2}
det[vec]=-1. 1-abs=0.
orthogonality check=0.
diagvc[{-1, 2}]
Generating matrix: ham.model..QSZ_-1.2
hamil={{(coefzeta[1, 0] - coefzeta[2, 0])/2, hybV[1, 1]}, {hybV[1, 1], epsilon - coefzeta[1, 0]/2 - coefzeta[2, 0]/2}}
dim={2, 2}
det[vec]=1.0000000000000002 1-abs=-2.220446049250313*^-16
orthogonality check=8.881784197001252*^-16
diagvc[{0, -1}]
Generating matrix: ham.model..QSZ_0.-1
hamil={{(2*epsilon - coefzeta[1, 0] + coefzeta[2, 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={{(coefzeta[1, 0] + coefzeta[2, 0])/2, -hybV[2, 2], hybV[1, 1], 0}, {-hybV[2, 2], (2*epsilon + coefzeta[1, 0] - coefzeta[2, 0])/2, 0, -hybV[1, 1]}, {hybV[1, 1], 0, (2*epsilon - coefzeta[1, 0] + coefzeta[2, 0])/2, hybV[2, 2]}, 
   {0, -hybV[1, 1], hybV[2, 2], 2*epsilon + U - coefzeta[1, 0]/2 - coefzeta[2, 0]/2}}
dim={4, 4}
det[vec]=1. 1-abs=0.
orthogonality check=1.609823385706477*^-15
diagvc[{0, 3}]
Generating matrix: ham.model..QSZ_0.3
hamil={{(2*epsilon + coefzeta[1, 0] - coefzeta[2, 0])/2}}
dim={1, 1}
det[vec]=1. 1-abs=0.
orthogonality check=0.
diagvc[{1, 0}]
Generating matrix: ham.model..QSZ_1.0
hamil={{(2*epsilon + coefzeta[1, 0] + coefzeta[2, 0])/2, -hybV[1, 1]}, {-hybV[1, 1], (4*epsilon + 2*U - coefzeta[1, 0] + coefzeta[2, 0])/2}}
dim={2, 2}
det[vec]=1.0000000000000002 1-abs=-2.220446049250313*^-16
orthogonality check=4.440892098500626*^-16
diagvc[{1, 2}]
Generating matrix: ham.model..QSZ_1.2
hamil={{(2*epsilon + coefzeta[1, 0] + coefzeta[2, 0])/2, -hybV[2, 2]}, {-hybV[2, 2], (4*epsilon + 2*U + coefzeta[1, 0] - coefzeta[2, 0])/2}}
dim={2, 2}
det[vec]=-1. 1-abs=0.
orthogonality check=0.
diagvc[{2, 1}]
Generating matrix: ham.model..QSZ_2.1
hamil={{(4*epsilon + 2*U + coefzeta[1, 0] + coefzeta[2, 0])/2}}
dim={1, 1}
det[vec]=1. 1-abs=0.
orthogonality check=0.
Lowest energies (absolute):{-0.24738856940681223, -0.22386212054647653, -0.16014976886682353, -0.1366796373244448, -0.07725392184632682, -0.07450809540355832, -0.05907216685652461, -0.016502612140896253, 0.016502612140896253, 
   0.01924843858366476, 0.024508095403558312, 0.06707765011918668, 0.11306721898751425, 0.13285141300202197, 0.2011604764112781, 0.22100098774374283}
Lowest energies (GS shifted):{0., 0.023526448860335697, 0.0872388005399887, 0.11070893208236743, 0.1701346475604854, 0.1728804740032539, 0.18831640255028761, 0.230885957265916, 0.26389118154770846, 0.26663700799047696, 0.27189666481037056, 
   0.3144662195259989, 0.36045578839432646, 0.38023998240883417, 0.4485490458180903, 0.4683895571505551}
Scale factor SCALE(Ninit):1.715663207404297
Lowest energies (shifted and scaled):{0., 0.013712743129771901, 0.050848441677534224, 0.06452835941493662, 0.09916552784149849, 0.10076597391443304, 0.10976303608865044, 0.13457533872002395, 0.15381292809033373, 0.15541337416326828, 
   0.1584790439271208, 0.1832913465584943, 0.21009705566844672, 0.22162856950468507, 0.2614435303399203, 0.2730078695685282}
makeireducf GENERAL
ireducTable: f[0]{}
Loading module operators.m
"operators.m started"
s: n_d op.model..QSZ.n_d nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]]
d: A_d d
ireducTable: d{}
ireducTable: Chop[Expand[komutator[Hselfd /. params, d[#1, #2]]]] & {}
s: SZd op.model..QSZ.SZd (-nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]])/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.005284`4.1745078024384625}
{ham, 0.0477619999999999999`4.176384988333275}
{maketable, 0.548544`6.190756462582266}
{xi, 0.052258`5.169697778037972}
{_, 0}
data
gammaPol=0.013811976175825838
"Success!"