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.2131570881878404
faktor=1.1657299587521546
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.07741600515104565816`18.888830756938766
"V[1,2]="0
"V[2,1]="0
"V[2,2]="-0.3160495175745600127`18.499755131746944
Diagonalisation.
Loading discretization data from files.
"nrch="1
"xi="{0.5829792498644170706`18.765653097075038, 0.5063102955121788762`18.704416758591, 0.4067960944949384205`18.609376774763394, 0.2954306426607566372`18.47045553920982, 0.2104017863188205606`18.323049422673307, 
   0.1499300661081554575`18.175888732484935, 0.1065879770469266852`18.02770821973133, 0.07559172745766201684`18.878474270156367, 0.05353249833925766837`18.728617512197825, 0.0378821502549768277`18.57843462232087, 
   0.02679699422148569196`18.428086082602228, 0.01895196593452728312`18.277654267086803, 0.01340234489606898297`18.127180789881965, 0.00947734060004963054`18.9766864887354, 0.006701651239453967381`18.826181822973773, 
   0.00473883903937700416`18.67567195759209, 0.003350885143300713805`18.525159541984607, 0.002369440406831054089`18.374645790347135, 0.001675449985907770216`18.22413146817937, 0.001184722759582000195`18.07361673172111, 
   0.0008377259499169357939`18.923101968625136, 0.0005923616491485490581`18.772586933538722, 0.0004188631251166735907`18.622072128651375, 0.0002961808121419566586`18.471556919645046, 0.0002094316137352320636`18.321042239111993, 
   0.0001480903591245977626`18.17052678630464, 0.000104715845841514061`18.020012405086888, 0.0000740451273362124565`18.8694964843208, 0.00005235796001866340849`18.718982717139667, 0.00003702251158379867181`18.56846587767959, 
   0.00002617901719575167334`18.417953338376936}
"zeta="{-0.06786225246857814641`18.831628270552827, 0.01482927473344288795`18.171119911179023, 0.01461109917120870239`18.164682888489704, 0.01412981638591657985`18.15013651831697, 0.008560263672416412728`18.932487141979944, 
   0.00463081236084496306`18.665657183865715, 0.002472097501770702168`18.393065595723385, 0.00130398212868691808`18.115271639351548, 0.0006811129958038703459`18.833219166807442, 0.000353619489766962016`18.548536193176403, 
   0.0001829444382788862351`18.262319210972095, 0.00009444107063288228695`18.97516090182174, 0.0000486764752536322784`18.68731912269444, 0.00002505913018017294478`18.39896599225641, 0.0000128884984426033513`18.110202323398614, 
   6.622912043241889214`18.82104898736029*^-6, 3.400691141165698547`18.53156719006901*^-6, 1.745039723345499249`18.241805317503452*^-6, 8.948371381287679921`18.951744000178383*^-7, 4.585700936184351952`18.661405728108704*^-7, 
   2.348642777069945818`18.370816966652047*^-7, 1.202164449725248673`18.079963880913887*^-7, 6.150035057371447394`18.788877591414373*^-8, 3.144822800927558449`18.497596179591326*^-8, 1.607272964678714951`18.20608963966731*^-8, 
   8.210601568439954655`18.914374977860355*^-9, 4.192435978990855298`18.62246643939183*^-9, 2.139738480591479371`18.33036069700507*^-9, 1.091603235662045037`18.03806481430831*^-9, 5.567204987468703523`18.74563721255672*^-10, 
   2.838012360377617847`18.453014282605636*^-10}
"nrch="2
"xi="{0.585938731203994867`18.767852206306245, 0.5058315886558011343`18.704005947093272, 0.4064588691635462325`18.6090166045449, 0.2951524431277023908`18.470046382449603, 0.210267783544837078`18.32277273678282, 
   0.1498749235748607633`18.17572897464886, 0.1065661596105750286`18.027619315125502, 0.07558328974372030551`18.878425790550722, 0.05352928493824177236`18.728591441973943, 0.03788093798650661592`18.578420724220372, 
   0.02679653970666521015`18.42807871629275, 0.018951796291269047`18.277650379602992, 0.01340228174838656479`18.127178743616565, 0.009477317137162211136`18.976685413558766, 0.006701642539845489334`18.826181259203103, 
   0.004738835820033975889`18.675671662552865, 0.003350883954212654451`18.525159387871753, 0.00236943996844386686`18.374645709995185, 0.001675449824470473092`18.224131426333102, 0.001184722700182529393`18.07361670994651, 
   0.0008377259280876383097`18.923101957308376, 0.0005923616411361322647`18.77258692766436, 0.0004188631221794167619`18.622072125605904, 0.0002961808110665226459`18.47155691806812, 0.0002094316133417721407`18.32104223829608, 
   0.000148090358980724297`18.17052678588271, 0.0001047158457889483483`18.020012404868883, 0.00007404512731702294541`18.86949648420825, 0.00005235796001166491057`18.718982717081616, 0.0000370225115812486283`18.56846587764968, 
   0.00002617901719482303385`18.41795333836153}
"zeta="{0.03393085530865209976`18.530594807217952, -0.007682925763276970445`18.885526636812337, -0.007204019835183774444`18.857574899913597, -0.006998879477190387416`18.84502851489621, -0.004245723321876953768`18.62795168938316, 
   -0.002302175639183133993`18.36213845405407, -0.001231221171378656717`18.090336074768413, -0.0006501560895576558204`18.813017634636005, -0.0003398123450392781194`18.531239152299506, -0.0001764876359975103186`18.24671428589303, 
   -0.00009132517246120590369`18.960590501019507, -0.00004715087273336133534`18.673489735668795, -0.00002430445860982540827`18.385685951459354, -0.00001251295952018517919`18.097360039611555, -6.436022693527665028`18.80861756675957*^-6, 
   -3.307363203194255905`18.519481890251377*^-6, -1.698304624409916882`18.230015592170925*^-6, -8.715008434655855587`18.940267811769903*^-7, -4.469094466983836132`18.650219534802716*^-7, -2.290305875973325763`18.359893487307748*^-7, 
   -1.173049363067027918`18.069316288037566*^-7, -6.00446282681636628`18.778474160484212*^-8, -3.07183801143594133`18.487398310124323*^-8, -1.570821687472151807`18.19612688870512*^-8, -8.028416418627888567`18.90462989042507*^-9, 
   -4.101326591268056144`18.612924353812755*^-9, -2.094230913431284391`18.32102456602479*^-9, -1.068875697722485578`18.028927202933712*^-9, -5.453048553060268265`18.736639364591696*^-10, -2.781118627603717535`18.44421951393245*^-10, 
   -1.417764279933518783`18.151604030258934*^-10}
BAND="manual_V" thetaCh={"0.005993237854", "0.09988729756"}
Discretization (channel 1)
"xitable" (channel 1)
0.5829792499
0.5063102955
0.4067960945
0.2954306427
0.2104017863
0.1499300661
0.106587977
0.07559172746
0.05353249834
0.03788215025
0.02679699422
0.01895196593
0.0134023449
0.0094773406
0.006701651239
0.004738839039
0.003350885143
0.002369440407
0.001675449986
0.00118472276
0.0008377259499
0.0005923616491
0.0004188631251
0.0002961808121
0.0002094316137
0.0001480903591
0.0001047158458
0.00007404512734
0.00005235796002
0.00003702251158
0.0000261790172
"zetatable" (channel 1)
-0.06786225247
0.01482927473
0.01461109917
0.01412981639
0.008560263672
0.004630812361
0.002472097502
0.001303982129
0.0006811129958
0.0003536194898
0.0001829444383
0.00009444107063
0.00004867647525
0.00002505913018
0.00001288849844
6.622912043e-6
3.400691141e-6
1.745039723e-6
8.948371381e-7
4.585700936e-7
2.348642777e-7
1.20216445e-7
6.150035057e-8
3.144822801e-8
1.607272965e-8
8.210601568e-9
4.192435979e-9
2.139738481e-9
1.091603236e-9
5.567204987e-10
2.83801236e-10
Precision last xi:18.417953338376936
Precision last zeta: 18.453014282605636
Discretization (channel 2)
"xitable" (channel 2)
0.5859387312
0.5058315887
0.4064588692
0.2951524431
0.2102677835
0.1498749236
0.1065661596
0.07558328974
0.05352928494
0.03788093799
0.02679653971
0.01895179629
0.01340228175
0.009477317137
0.00670164254
0.00473883582
0.003350883954
0.002369439968
0.001675449824
0.0011847227
0.0008377259281
0.0005923616411
0.0004188631222
0.0002961808111
0.0002094316133
0.000148090359
0.0001047158458
0.00007404512732
0.00005235796001
0.00003702251158
0.00002617901719
"zetatable" (channel 2)
0.03393085531
-0.007682925763
-0.007204019835
-0.006998879477
-0.004245723322
-0.002302175639
-0.001231221171
-0.0006501560896
-0.000339812345
-0.000176487636
-0.00009132517246
-0.00004715087273
-0.00002430445861
-0.00001251295952
-6.436022694e-6
-3.307363203e-6
-1.698304624e-6
-8.715008435e-7
-4.469094467e-7
-2.290305876e-7
-1.173049363e-7
-6.004462827e-8
-3.071838011e-8
-1.570821687e-8
-8.028416419e-9
-4.101326591e-9
-2.094230913e-9
-1.068875698e-9
-5.453048553e-10
-2.781118628e-10
-1.41776428e-10
Precision last xi:18.41795333836153
Precision last zeta: 18.151604030258934
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.005993237854", "0.09988729756"}
theta0Ch={"0.0005993237853546729", "0.009988729755911212"}
gammaPolCh={"0.013811976176618772", "0.05638715662032216"}
checkdefinitions[] -> -0.5609439973938823
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. 1-abs=0.
orthogonality check=0.
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.0000000000000002 1-abs=-2.220446049250313*^-16
orthogonality check=1.9984014443252818*^-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. 1-abs=0.
orthogonality check=2.220446049250313*^-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.24802740509691323, -0.22479887657437553, -0.15929885396849777, -0.13612669213130027, -0.07811725219056392, -0.07589655388861513, -0.05928060235868207, -0.016965698579963023, 0.016965698579963023, 
   0.019186396881911837, 0.02589655388861512, 0.06821145766733416, 0.11257786228855032, 0.13193662410579743, 0.20216110643707594, 0.22157623493966322}
Lowest energies (GS shifted):{0., 0.023228528522537706, 0.08872855112841546, 0.11190071296561296, 0.1699101529063493, 0.1721308512082981, 0.18874680273823116, 0.2310617065169502, 0.26499310367687623, 0.26721380197882505, 
   0.27392395898552835, 0.3162388627642474, 0.36060526738546356, 0.37996402920271066, 0.4501885115339892, 0.46960364003657645}
Scale factor SCALE(Ninit):1.2131570881878404
Lowest energies (shifted and scaled):{0., 0.019147172900119174, 0.07313855063976438, 0.09223926073148962, 0.14005618444694037, 0.14188669619481797, 0.15558315124727387, 0.19046313850591254, 0.2184326384909567, 0.2202631502388343, 
   0.22579430285875315, 0.2606742901173918, 0.2972453204095106, 0.31320266180061135, 0.37108839070994554, 0.3870921949094402}
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.00522`4.1692154964982375}
{ham, 0.0426800000000000001`4.12752689480908}
{maketable, 0.563148`6.202167539566682}
{xi, 0.052111`5.168474400769715}
{_, 0}
data
gammaPol=0.013811976176618772
"Success!"