NRG Ljubljana  (c) Rok Zitko, rok.zitko@ijs.si, 2005-2012
initial $Id: initial.m,v 1.1 2012/05/15 09:58:19 rokzitko Exp rokzitko $
Mathematica version: 8.0 for Mac OS X x86 (64-bit) (November 6, 2010)
sneg version: 1.233
Loading module initialparse.m
Discretization: Z
U=1 Gamma=0.03 delta=0 t=0.0001
epsilon=-1/2
QS SIDE 
Options:
Options: {}
Loading module models.m
"models $Id: models.m 877 2009-12-16 16:34:11Z zitko $"
def1ch, NRDOTS=2
COEFCHANNELS:1
Hc=gammaPolCh[1]*(nc[d[0, 0], f[1, 0, 0]] + nc[d[0, 1], f[1, 0, 1]] + 
     nc[f[0, 0, 0], d[1, 0]] + nc[f[0, 0, 1], d[1, 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=2
Loading module custommodels.m
models $Id: custommodels.m,v 1.2 2011/10/21 16:59:38 rokzitko Exp rokzitko $
custommodels.m done
NRDOTS:2
CHANNELS:1
basis:{a[], d[], f[0]}
lrchain:{}
lrextrarule:{}
NROPS:3
Hamiltonian generated. U - coefzeta[1, 0] + delta*nc[a[0, 0], a[1, 0]] - 
   (U*nc[a[0, 0], a[1, 0]])/2 + t*nc[a[0, 0], d[1, 0]] + 
   delta*nc[a[0, 1], a[1, 1]] - (U*nc[a[0, 1], a[1, 1]])/2 + 
   t*nc[a[0, 1], d[1, 1]] + t*nc[d[0, 0], a[1, 0]] + 
   delta*nc[d[0, 0], d[1, 0]] - (U*nc[d[0, 0], d[1, 0]])/2 + 
   gammaPolCh[1]*nc[d[0, 0], f[1, 0, 0]] + t*nc[d[0, 1], a[1, 1]] + 
   delta*nc[d[0, 1], d[1, 1]] - (U*nc[d[0, 1], d[1, 1]])/2 + 
   gammaPolCh[1]*nc[d[0, 1], f[1, 0, 1]] + 
   gammaPolCh[1]*nc[f[0, 0, 0], d[1, 0]] + 
   coefzeta[1, 0]*nc[f[0, 0, 0], f[1, 0, 0]] + 
   gammaPolCh[1]*nc[f[0, 0, 1], d[1, 1]] + 
   coefzeta[1, 0]*nc[f[0, 0, 1], f[1, 0, 1]] - 
   U*nc[a[0, 0], a[0, 1], a[1, 0], a[1, 1]] - 
   U*nc[d[0, 0], d[0, 1], d[1, 0], d[1, 1]]
H-conj[H]=0
SCALE[0]=1.0201394465967895
faktor=1.3862943611198906
Generating basis
Basis states generated.
BASIS NR=64
Saving basis to basis.SIDE..QS
PREC=1000
Tmin=1.*^-16
Tmin=1.*^-16 ==> Nmax=106
DISCNMAX=106
mMAX=212
Diagonalisation.
Discretization checksum [-1] (channel 1): \
7.5964541966078389979785938156496`10.*^-65
BAND="flat" thetaCh={"2."}
Discretization (channel 1)
"xitable" (channel 1)
0.755928946
0.81461308
0.8924940297
0.942104181
0.9699342623
0.9846755156
0.9922632047
0.9961127494
0.9980516341
0.9990246284
0.9995120166
0.9997559339
0.9998779483
0.9999389695
0.9999694836
0.9999847415
0.9999923707
0.9999961853
0.9999980927
0.9999990463
0.9999995232
0.9999997616
0.9999998808
0.9999999404
0.9999999702
0.9999999851
0.9999999925
0.9999999963
0.9999999981
0.9999999991
0.9999999995
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
"zetatable" (channel 1)
0.e-999
0.e-998
0.e-997
0.e-996
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-986
0.e-985
0.e-984
0.e-983
0.e-982
0.e-981
0.e-980
0.e-980
0.e-979
0.e-978
0.e-977
0.e-976
0.e-975
0.e-975
0.e-974
0.e-973
0.e-972
0.e-971
0.e-970
0.e-969
0.e-969
0.e-968
0.e-967
0.e-966
0.e-965
0.e-964
0.e-964
0.e-963
0.e-962
0.e-961
0.e-960
0.e-959
0.e-958
0.e-958
0.e-957
0.e-956
0.e-955
0.e-954
0.e-953
0.e-953
0.e-952
0.e-951
0.e-950
0.e-949
0.e-948
0.e-947
0.e-947
0.e-946
0.e-945
0.e-944
0.e-943
0.e-942
0.e-942
0.e-941
0.e-940
0.e-939
0.e-938
0.e-937
0.e-937
0.e-936
0.e-935
0.e-934
0.e-933
0.e-932
0.e-931
0.e-931
0.e-930
0.e-929
0.e-928
0.e-919
0.e-893
0.e-867
0.e-841
0.e-814
0.e-788
0.e-761
0.e-733
0.e-706
0.e-678
0.e-649
0.e-621
0.e-592
0.e-563
0.e-533
0.e-503
0.e-473
0.e-443
0.e-412
0.e-381
0.e-350
0.e-318
Precision last xi:286.8284471949013
Precision last zeta: 0.
Discretization done.
--EOF--

           {{# Input file for NRG Ljubljana, Rok Zitko, rok.zitko@ijs.si, \

maketable[]

2005-2012}, {# symtype , QS}, {# $Id: initial.m,v 1.1 2012/05/15 09:58:19 \
rokzitko Exp rokzitko $}, {# Using sneg version , 1.233}, 
   {# Model: , SIDE,  Variant: , ,  Channels: , 1}, 
   {# U=, 1,  Gamma=, 0.03,  Delta=, 0,  t=, 0.0001}, {#}, 
   {# Gamma^(1/2)=, Sqrt[theta0]/Sqrt[Pi]}, {# Lambda=, 2.,  BAND=, flat}, 
   {# DISCRETIZATION=, Z,  z=, 1.}, {# ops=, SaSd}, {#!7}, 
   {# Number of channels, impurities, subspaces: }, {1, 2, 10}}
exnames={an1x, an1y, an1z, an2x, an2y, an2z, anB2, anB4, anB6, anD, anE, 
   anodelta, B, beta, Bf, Borb, Bx, By, Bz, delta1, delta2, delta3, deltaeps, 
   eps, eps0, eps1, eps2, eps3, epsilon, F0, F2, F4, g, g1, g2, Gamma1, 
   Gamma2, Gamma2to2, Gamma3, gammaa, gch, ggg, gl, gs, J0, J1, J12, J2, Jad, 
   Jcharge, Jcharge1, Jcharge2, JH, jhop, Jkondo, Jkondo1, Jkondo1ch2, 
   Jkondo1P, Jkondo1Z, Jkondo2, Jkondo2ch2, Jkondo2P, Jkondo2Z, Jkondo3, 
   JkondoDM, JkondoP, JkondoZ, Jspin, k, kappa, kFR, lr, lslambda, maj1do, 
   maj1up, maj2do, maj2up, majdelta, majdeltado, majdeltaup, n0, omega, 
   pairdelta, phi, Pl, Pr, Ra, t0, t12, ta, tab, tad, Tad, tb, tpp, tpp12, 
   tpp13, tpp23, U1, U12, U2, U3, Uad, Uf, Uprime, V, V1, V2, V3}
thetaCh={"2."}
theta0Ch={"0.06"}
gammaPolCh={"0.13819765978853418"}
checkdefinitions[] -> -2.4468093608458634
calcgsenergy[]
diagvc[{-3, 1}]
Generating matrix and saving to ham.SIDE..QS_-3.1
hamil={{U - coefzeta[1, 0]}}
dim={1, 1}
det[vec]=1. 1-abs=0.
orthogonality check=Abs[-1 + MyDot[{{1.}}, {{1.}}]]
diagvc[{-2, 2}]
Generating matrix and saving to ham.SIDE..QS_-2.2
hamil={{U, gammaPolCh[1], 0}, {gammaPolCh[1], delta + U/2 - coefzeta[1, 0], 
    t}, {0, t, delta + U/2 - coefzeta[1, 0]}}
dim={3, 3}
det[vec]=1.0000000000000004 1-abs=-4.440892098500626*^-16
orthogonality check=3*Abs[-1 + 
      MyDot[{{-0.24981630019504283, 0.9682894407169533, 
         -0.0027157232757237717}, 
        {-0.000723593626038112, 0.0026179649072038957, 0.9999963113292014}, 
        {0.9682929786842023, 0.24981734378500364, 0.000046637760304660004}}, 
       {{-0.24981630019504283, -0.000723593626038112, 0.9682929786842023}, 
        {0.9682894407169533, 0.0026179649072038957, 0.24981734378500364}, 
        {-0.0027157232757237717, 0.9999963113292014, 
         0.000046637760304660004}}]] + 
   6*Abs[MyDot[{{-0.24981630019504283, 0.9682894407169533, 
        -0.0027157232757237717}, 
       {-0.000723593626038112, 0.0026179649072038957, 0.9999963113292014}, 
       {0.9682929786842023, 0.24981734378500364, 0.000046637760304660004}}, 
      {{-0.24981630019504283, -0.000723593626038112, 0.9682929786842023}, 
       {0.9682894407169533, 0.0026179649072038957, 0.24981734378500364}, 
       {-0.0027157232757237717, 0.9999963113292014, 0.000046637760304660004}}]\
]
diagvc[{-1, 1}]
Generating matrix and saving to ham.SIDE..QS_-1.1
hamil={{U + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], 0, 0, 0, 0}, 
   {Sqrt[2]*gammaPolCh[1], delta + U/2, t, Sqrt[2]*gammaPolCh[1], 0, 0}, 
   {0, t, delta + U/2, 0, gammaPolCh[1], 0}, 
   {0, Sqrt[2]*gammaPolCh[1], 0, 2*delta + U - coefzeta[1, 0], Sqrt[2]*t, 0}, 
   {0, 0, gammaPolCh[1], Sqrt[2]*t, 2*delta - coefzeta[1, 0], Sqrt[2]*t}, 
   {0, 0, 0, 0, Sqrt[2]*t, 2*delta + U - coefzeta[1, 0]}}
dim={6, 6}
det[vec]=-0.9999999999999996 1-abs=4.440892098500626*^-16
orthogonality check=6*Abs[-1 + 
      MyDot[{{-0.000020764739017871614, 0.00011003371542703169, 
         -0.24981733070924536, -0.00015298767088339904, 
         0.9682929555927706, -0.0001322229318654701}, 
        {-0.2868766054104093, 0.9140040448743016, -0.0004416600500244338, 
         -0.2868765442506637, -0.0002692894524697976, 
         6.115974557024613*^-8}, 
        {-0.0001723304827126727, 0.00040943736089502346, 
         0.9682928773959391, -0.0002484150020636783, 0.2498172106734797, 
         -0.00007608451935047234}, 
        {-0.5773502679277233, 3.613752299039022*^-16, 
         -1.5323785280643059*^-13, 0.5773502679277241, 
         -5.568594037627531*^-13, -0.5773502717134298}, 
        {0.40824838119964435, 6.515442850245887*^-8, 
         -0.00004977428916745137, -0.40824818904628707, 
         -0.00018008374880979915, -0.8164965648921346}, 
        {0.6462984685970459, 0.40570485110119153, 0.00008556004226248933, 
         0.6462985745828072, 0.0000919444011115157, 1.0598576177402735*^-7}}\
, {{-0.000020764739017871614, -0.2868766054104093, -0.0001723304827126727, 
         -0.5773502679277233, 0.40824838119964435, 0.6462984685970459}, 
        {0.00011003371542703169, 0.9140040448743016, 
         0.00040943736089502346, 3.613752299039022*^-16, 
         6.515442850245887*^-8, 0.40570485110119153}, 
        {-0.24981733070924536, -0.0004416600500244338, 0.9682928773959391, 
         -1.5323785280643059*^-13, -0.00004977428916745137, 
         0.00008556004226248933}, 
        {-0.00015298767088339904, -0.2868765442506637, 
         -0.0002484150020636783, 0.5773502679277241, -0.40824818904628707, 
         0.6462985745828072}, {0.9682929555927706, -0.0002692894524697976, 
         0.2498172106734797, -5.568594037627531*^-13, 
         -0.00018008374880979915, 0.0000919444011115157}, 
        {-0.0001322229318654701, 6.115974557024613*^-8, 
         -0.00007608451935047234, -0.5773502717134298, -0.8164965648921346, 
         1.0598576177402735*^-7}}]] + 
   30*Abs[MyDot[{{-0.000020764739017871614, 0.00011003371542703169, 
        -0.24981733070924536, -0.00015298767088339904, 0.9682929555927706, 
        -0.0001322229318654701}, 
       {-0.2868766054104093, 0.9140040448743016, -0.0004416600500244338, 
        -0.2868765442506637, -0.0002692894524697976, 6.115974557024613*^-8}, 
       {-0.0001723304827126727, 0.00040943736089502346, 0.9682928773959391, 
        -0.0002484150020636783, 0.2498172106734797, -0.00007608451935047234}\
, {-0.5773502679277233, 3.613752299039022*^-16, -1.5323785280643059*^-13, 
        0.5773502679277241, -5.568594037627531*^-13, -0.5773502717134298}, 
       {0.40824838119964435, 6.515442850245887*^-8, 
        -0.00004977428916745137, -0.40824818904628707, 
        -0.00018008374880979915, -0.8164965648921346}, 
       {0.6462984685970459, 0.40570485110119153, 0.00008556004226248933, 
        0.6462985745828072, 0.0000919444011115157, 1.0598576177402735*^-7}}, 
      {{-0.000020764739017871614, -0.2868766054104093, 
        -0.0001723304827126727, -0.5773502679277233, 0.40824838119964435, 
        0.6462984685970459}, {0.00011003371542703169, 0.9140040448743016, 
        0.00040943736089502346, 3.613752299039022*^-16, 
        6.515442850245887*^-8, 0.40570485110119153}, 
       {-0.24981733070924536, -0.0004416600500244338, 0.9682928773959391, 
        -1.5323785280643059*^-13, -0.00004977428916745137, 
        0.00008556004226248933}, 
       {-0.00015298767088339904, -0.2868765442506637, 
        -0.0002484150020636783, 0.5773502679277241, -0.40824818904628707, 
        0.6462985745828072}, {0.9682929555927706, -0.0002692894524697976, 
        0.2498172106734797, -5.568594037627531*^-13, 
        -0.00018008374880979915, 0.0000919444011115157}, 
       {-0.0001322229318654701, 6.115974557024613*^-8, 
        -0.00007608451935047234, -0.5773502717134298, -0.8164965648921346, 
        1.0598576177402735*^-7}}]]
diagvc[{-1, 3}]
Generating matrix and saving to ham.SIDE..QS_-1.3
hamil={{delta + U/2, t, 0}, {t, delta + U/2, gammaPolCh[1]}, 
   {0, gammaPolCh[1], 2*delta - coefzeta[1, 0]}}
dim={3, 3}
det[vec]=1.0000000000000002 1-abs=-2.220446049250313*^-16
orthogonality check=3*Abs[-1 + 
      MyDot[{{0.00004663776030466001, -0.24981734378500361, 
         0.968292978684202}, {0.9999963113292014, -0.002617964907203922, 
         -0.0007235936260381192}, 
        {0.0027157232757238007, 0.9682894407169532, 0.2498163001950428}}, 
       {{0.00004663776030466001, 0.9999963113292014, 0.0027157232757238007}, 
        {-0.24981734378500361, -0.002617964907203922, 0.9682894407169532}, 
        {0.968292978684202, -0.0007235936260381192, 0.2498163001950428}}]] + 
   6*Abs[MyDot[{{0.00004663776030466001, -0.24981734378500361, 
        0.968292978684202}, {0.9999963113292014, -0.002617964907203922, 
        -0.0007235936260381192}, 
       {0.0027157232757238007, 0.9682894407169532, 0.2498163001950428}}, 
      {{0.00004663776030466001, 0.9999963113292014, 0.0027157232757238007}, 
       {-0.24981734378500361, -0.002617964907203922, 0.9682894407169532}, 
       {0.968292978684202, -0.0007235936260381192, 0.2498163001950428}}]]
diagvc[{0, 2}]
Generating matrix and saving to ham.SIDE..QS_0.2
hamil={{delta + U/2 + coefzeta[1, 0], t, -gammaPolCh[1], 0, 0, 0, 0, 0}, 
   {t, delta + U/2 + coefzeta[1, 0], 0, -(gammaPolCh[1]/Sqrt[2]), 0, 
    Sqrt[3/2]*gammaPolCh[1], 0, 0}, 
   {-gammaPolCh[1], 0, 2*delta + U, Sqrt[2]*t, 0, 0, 0, 0}, 
   {0, -(gammaPolCh[1]/Sqrt[2]), Sqrt[2]*t, 2*delta, Sqrt[2]*t, 0, 
    -(gammaPolCh[1]/Sqrt[2]), 0}, 
   {0, 0, 0, Sqrt[2]*t, 2*delta + U, 0, 0, gammaPolCh[1]}, 
   {0, Sqrt[3/2]*gammaPolCh[1], 0, 0, 0, 2*delta, Sqrt[3/2]*gammaPolCh[1], 
    0}, {0, 0, 0, -(gammaPolCh[1]/Sqrt[2]), 0, Sqrt[3/2]*gammaPolCh[1], 
    (6*delta + U - 2*coefzeta[1, 0])/2, -t}, 
   {0, 0, 0, 0, gammaPolCh[1], 0, -t, (6*delta + U - 2*coefzeta[1, 0])/2}}
dim={8, 8}
det[vec]=1.0000000000000004 1-abs=-4.440892098500626*^-16
orthogonality check=8*Abs[-1 + 
      MyDot[{{-0.00006050008608251001, 0.28687665137393115, 
         -0.00006501450951817033, 0.4570021797019532, 
         -0.00006501450951809992, -0.79155073480247, 0.2868766513739313, 
         0.00006050008608245874}, 
        {-0.000035195737398820544, -4.607113813191728*^-8, 
         -0.00012733843996512264, 0.8660253207325078, 
         -0.00012733843996507363, 0.5000001089425044, 
         -4.6071138373043845*^-8, 0.00003519573739875454}, 
        {-0.6846840291099598, 0.0019203063438494629, -0.17664679977118994, 
         1.917487721991862*^-13, 0.17664680006650768, 
         -1.7636625513918332*^-13, -0.0019203063443331506, 
         -0.6846840302546146}, 
        {0.6846864603536281, -0.0002895159359657518, 0.17664744387198453, 
         0.00022945581208363893, 0.17664744357666784, 
         -0.00021106087488727827, -0.00028951593275473023, 
         -0.6846864592089774}, 
        {0.0018511807387937985, 0.7071041729024125, 0.0005116579597952318, 
         -4.603905068302046*^-16, -0.0005116579597949823, 
         2.323990368177095*^-16, -0.7071041729024122, 
         0.0018511807387929053}, 
        {0.0003123008163518747, 0.6462984581625518, 
         0.00019041639794336451, -0.2028523065565367, 
         0.00019041639794311824, 0.3513506511986275, 0.6462984581625525, 
         -0.00031230081635170585}, 
        {0.17664753045151063, 0.000032977873312455134, 
         -0.6846865027325991, 8.444782186359761*^-12, 0.6846865600732407, 
         -1.0648322874221963*^-12, -0.00003297787982830297, 
         0.176647545245263}, {-0.1766475359993067, -0.00007780558771853452, 
         0.6846865437451344, 0.0002016743512698568, 0.6846864864044897, 
         -0.000025431507618109802, -0.00007780558495615368, 
         0.17664752120555277}}, 
       {{-0.00006050008608251001, -0.000035195737398820544, 
         -0.6846840291099598, 0.6846864603536281, 0.0018511807387937985, 
         0.0003123008163518747, 0.17664753045151063, -0.1766475359993067}, 
        {0.28687665137393115, -4.607113813191728*^-8, 
         0.0019203063438494629, -0.0002895159359657518, 0.7071041729024125, 
         0.6462984581625518, 0.000032977873312455134, 
         -0.00007780558771853452}, 
        {-0.00006501450951817033, -0.00012733843996512264, 
         -0.17664679977118994, 0.17664744387198453, 0.0005116579597952318, 
         0.00019041639794336451, -0.6846865027325991, 0.6846865437451344}, 
        {0.4570021797019532, 0.8660253207325078, 1.917487721991862*^-13, 
         0.00022945581208363893, -4.603905068302046*^-16, 
         -0.2028523065565367, 8.444782186359761*^-12, 0.0002016743512698568}\
, {-0.00006501450951809992, -0.00012733843996507363, 0.17664680006650768, 
         0.17664744357666784, -0.0005116579597949823, 
         0.00019041639794311824, 0.6846865600732407, 0.6846864864044897}, 
        {-0.79155073480247, 0.5000001089425044, -1.7636625513918332*^-13, 
         -0.00021106087488727827, 2.323990368177095*^-16, 
         0.3513506511986275, -1.0648322874221963*^-12, 
         -0.000025431507618109802}, 
        {0.2868766513739313, -4.6071138373043845*^-8, 
         -0.0019203063443331506, -0.00028951593275473023, 
         -0.7071041729024122, 0.6462984581625525, -0.00003297787982830297, 
         -0.00007780558495615368}, 
        {0.00006050008608245874, 0.00003519573739875454, 
         -0.6846840302546146, -0.6846864592089774, 0.0018511807387929053, 
         -0.00031230081635170585, 0.176647545245263, 0.17664752120555277}}]] \
+ 56*Abs[MyDot[{{-0.00006050008608251001, 0.28687665137393115, 
        -0.00006501450951817033, 0.4570021797019532, 
        -0.00006501450951809992, -0.79155073480247, 0.2868766513739313, 
        0.00006050008608245874}, 
       {-0.000035195737398820544, -4.607113813191728*^-8, 
        -0.00012733843996512264, 0.8660253207325078, 
        -0.00012733843996507363, 0.5000001089425044, 
        -4.6071138373043845*^-8, 0.00003519573739875454}, 
       {-0.6846840291099598, 0.0019203063438494629, -0.17664679977118994, 
        1.917487721991862*^-13, 0.17664680006650768, 
        -1.7636625513918332*^-13, -0.0019203063443331506, 
        -0.6846840302546146}, {0.6846864603536281, -0.0002895159359657518, 
        0.17664744387198453, 0.00022945581208363893, 0.17664744357666784, 
        -0.00021106087488727827, -0.00028951593275473023, 
        -0.6846864592089774}, {0.0018511807387937985, 0.7071041729024125, 
        0.0005116579597952318, -4.603905068302046*^-16, 
        -0.0005116579597949823, 2.323990368177095*^-16, 
        -0.7071041729024122, 0.0018511807387929053}, 
       {0.0003123008163518747, 0.6462984581625518, 0.00019041639794336451, 
        -0.2028523065565367, 0.00019041639794311824, 0.3513506511986275, 
        0.6462984581625525, -0.00031230081635170585}, 
       {0.17664753045151063, 0.000032977873312455134, -0.6846865027325991, 
        8.444782186359761*^-12, 0.6846865600732407, 
        -1.0648322874221963*^-12, -0.00003297787982830297, 
        0.176647545245263}, {-0.1766475359993067, -0.00007780558771853452, 
        0.6846865437451344, 0.0002016743512698568, 0.6846864864044897, 
        -0.000025431507618109802, -0.00007780558495615368, 
        0.17664752120555277}}, {{-0.00006050008608251001, 
        -0.000035195737398820544, -0.6846840291099598, 0.6846864603536281, 
        0.0018511807387937985, 0.0003123008163518747, 0.17664753045151063, 
        -0.1766475359993067}, {0.28687665137393115, -4.607113813191728*^-8, 
        0.0019203063438494629, -0.0002895159359657518, 0.7071041729024125, 
        0.6462984581625518, 0.000032977873312455134, -0.00007780558771853452}\
, {-0.00006501450951817033, -0.00012733843996512264, -0.17664679977118994, 
        0.17664744387198453, 0.0005116579597952318, 0.00019041639794336451, 
        -0.6846865027325991, 0.6846865437451344}, 
       {0.4570021797019532, 0.8660253207325078, 1.917487721991862*^-13, 
        0.00022945581208363893, -4.603905068302046*^-16, 
        -0.2028523065565367, 8.444782186359761*^-12, 0.0002016743512698568}, 
       {-0.00006501450951809992, -0.00012733843996507363, 
        0.17664680006650768, 0.17664744357666784, -0.0005116579597949823, 
        0.00019041639794311824, 0.6846865600732407, 0.6846864864044897}, 
       {-0.79155073480247, 0.5000001089425044, -1.7636625513918332*^-13, 
        -0.00021106087488727827, 2.323990368177095*^-16, 0.3513506511986275, 
        -1.0648322874221963*^-12, -0.000025431507618109802}, 
       {0.2868766513739313, -4.6071138373043845*^-8, -0.0019203063443331506, 
        -0.00028951593275473023, -0.7071041729024122, 0.6462984581625525, 
        -0.00003297787982830297, -0.00007780558495615368}, 
       {0.00006050008608245874, 0.00003519573739875454, -0.6846840302546146, 
        -0.6846864592089774, 0.0018511807387929053, -0.00031230081635170585, 
        0.176647545245263, 0.17664752120555277}}]]
diagvc[{0, 4}]
Generating matrix and saving to ham.SIDE..QS_0.4
hamil={{2*delta}}
dim={1, 1}
det[vec]=1. 1-abs=0.
orthogonality check=Abs[-1 + MyDot[{{1.}}, {{1.}}]]
diagvc[{1, 1}]
Generating matrix and saving to ham.SIDE..QS_1.1
hamil={{2*delta + U + coefzeta[1, 0], Sqrt[2]*t, 0, 0, 0, 0}, 
   {Sqrt[2]*t, 2*delta + coefzeta[1, 0], Sqrt[2]*t, -gammaPolCh[1], 0, 0}, 
   {0, Sqrt[2]*t, 2*delta + U + coefzeta[1, 0], 0, Sqrt[2]*gammaPolCh[1], 0}, 
   {0, -gammaPolCh[1], 0, (6*delta + U)/2, -t, 0}, 
   {0, 0, Sqrt[2]*gammaPolCh[1], -t, (6*delta + U)/2, Sqrt[2]*gammaPolCh[1]}, 
   {0, 0, 0, 0, Sqrt[2]*gammaPolCh[1], 4*delta + U - coefzeta[1, 0]}}
dim={6, 6}
det[vec]=-1.0000000000000002 1-abs=-2.220446049250313*^-16
orthogonality check=6*Abs[-1 + 
      MyDot[{{0.00013222293186547653, -0.9682929555927707, 
         0.00015298767088338768, -0.24981733070924522, 
         -0.00011003371542705252, 0.000020764739017938753}, 
        {-6.115974574268144*^-8, 0.0002692894524696393, 
         0.2868765442506641, -0.0004416600500242631, -0.9140040448743015, 
         0.2868766054104094}, {-0.000076084519350465, 0.24981721067348003, 
         -0.0002484150020636044, -0.968292877395939, 
         0.00040943736089466806, -0.00017233048271281766}, 
        {0.5773502579028013, -2.4895277837174417*^-12, -0.577350274833038, 
         6.883890927975091*^-13, -2.675963322159202*^-16, 
         0.5773502748330378}, {-0.8164965746577241, 
         -0.00018008374880982083, -0.40824817928069745, 
         0.000049774289167236666, 6.515442915498355*^-8, 
         0.4082483714340557}, {-1.0598576365083804*^-7, 
         -0.00009194440111157459, -0.6462985745828063, 
         0.00008556004226283188, -0.405704851101192, -0.6462984685970466}}, 
       {{0.00013222293186547653, -6.115974574268144*^-8, 
         -0.000076084519350465, 0.5773502579028013, -0.8164965746577241, 
         -1.0598576365083804*^-7}, 
        {-0.9682929555927707, 0.0002692894524696393, 0.24981721067348003, 
         -2.4895277837174417*^-12, -0.00018008374880982083, 
         -0.00009194440111157459}, 
        {0.00015298767088338768, 0.2868765442506641, 
         -0.0002484150020636044, -0.577350274833038, -0.40824817928069745, 
         -0.6462985745828063}, 
        {-0.24981733070924522, -0.0004416600500242631, -0.968292877395939, 
         6.883890927975091*^-13, 0.000049774289167236666, 
         0.00008556004226283188}, 
        {-0.00011003371542705252, -0.9140040448743015, 
         0.00040943736089466806, -2.675963322159202*^-16, 
         6.515442915498355*^-8, -0.405704851101192}, 
        {0.000020764739017938753, 0.2868766054104094, 
         -0.00017233048271281766, 0.5773502748330378, 0.4082483714340557, 
         -0.6462984685970466}}]] + 
   30*Abs[MyDot[{{0.00013222293186547653, -0.9682929555927707, 
        0.00015298767088338768, -0.24981733070924522, 
        -0.00011003371542705252, 0.000020764739017938753}, 
       {-6.115974574268144*^-8, 0.0002692894524696393, 0.2868765442506641, 
        -0.0004416600500242631, -0.9140040448743015, 0.2868766054104094}, 
       {-0.000076084519350465, 0.24981721067348003, -0.0002484150020636044, 
        -0.968292877395939, 0.00040943736089466806, -0.00017233048271281766}\
, {0.5773502579028013, -2.4895277837174417*^-12, -0.577350274833038, 
        6.883890927975091*^-13, -2.675963322159202*^-16, 0.5773502748330378}\
, {-0.8164965746577241, -0.00018008374880982083, -0.40824817928069745, 
        0.000049774289167236666, 6.515442915498355*^-8, 0.4082483714340557}, 
       {-1.0598576365083804*^-7, -0.00009194440111157459, 
        -0.6462985745828063, 0.00008556004226283188, -0.405704851101192, 
        -0.6462984685970466}}, {{0.00013222293186547653, 
        -6.115974574268144*^-8, -0.000076084519350465, 0.5773502579028013, 
        -0.8164965746577241, -1.0598576365083804*^-7}, 
       {-0.9682929555927707, 0.0002692894524696393, 0.24981721067348003, 
        -2.4895277837174417*^-12, -0.00018008374880982083, 
        -0.00009194440111157459}, 
       {0.00015298767088338768, 0.2868765442506641, -0.0002484150020636044, 
        -0.577350274833038, -0.40824817928069745, -0.6462985745828063}, 
       {-0.24981733070924522, -0.0004416600500242631, -0.968292877395939, 
        6.883890927975091*^-13, 0.000049774289167236666, 
        0.00008556004226283188}, 
       {-0.00011003371542705252, -0.9140040448743015, 
        0.00040943736089466806, -2.675963322159202*^-16, 
        6.515442915498355*^-8, -0.405704851101192}, 
       {0.000020764739017938753, 0.2868766054104094, -0.00017233048271281766, 
        0.5773502748330378, 0.4082483714340557, -0.6462984685970466}}]]
diagvc[{1, 3}]
Generating matrix and saving to ham.SIDE..QS_1.3
hamil={{2*delta + coefzeta[1, 0], -gammaPolCh[1], 0}, 
   {-gammaPolCh[1], (6*delta + U)/2, -t}, {0, -t, (6*delta + U)/2}}
dim={3, 3}
det[vec]=-0.9999999999999998 1-abs=2.220446049250313*^-16
orthogonality check=3*Abs[-1 + 
      MyDot[{{-0.9682929786842019, -0.24981734378500364, 
         -0.00004663776030466}, 
        {-0.0007235936260381085, 0.00261796490720388, 0.9999963113292013}, 
        {0.24981630019504286, -0.968289440716953, 0.002715723275723755}}, 
       {{-0.9682929786842019, -0.0007235936260381085, 0.24981630019504286}, 
        {-0.24981734378500364, 0.00261796490720388, -0.968289440716953}, 
        {-0.00004663776030466, 0.9999963113292013, 0.002715723275723755}}]] + 
   6*Abs[MyDot[{{-0.9682929786842019, -0.24981734378500364, 
        -0.00004663776030466}, 
       {-0.0007235936260381085, 0.00261796490720388, 0.9999963113292013}, 
       {0.24981630019504286, -0.968289440716953, 0.002715723275723755}}, 
      {{-0.9682929786842019, -0.0007235936260381085, 0.24981630019504286}, 
       {-0.24981734378500364, 0.00261796490720388, -0.968289440716953}, 
       {-0.00004663776030466, 0.9999963113292013, 0.002715723275723755}}]]
diagvc[{2, 2}]
Generating matrix and saving to ham.SIDE..QS_2.2
hamil={{3*delta + U/2 + coefzeta[1, 0], -t, 0}, 
   {-t, 3*delta + U/2 + coefzeta[1, 0], -gammaPolCh[1]}, 
   {0, -gammaPolCh[1], 4*delta + U}}
dim={3, 3}
det[vec]=-1. 1-abs=0.
orthogonality check=3*Abs[-1 + 
      MyDot[{{-0.002715723275723763, -0.9682894407169532, 
         -0.2498163001950426}, 
        {-0.9999963113292014, 0.0026179649072038884, 
         0.0007235936260381096}, 
        {0.00004663776030465997, -0.2498173437850035, 0.968292978684202}}, 
       {{-0.002715723275723763, -0.9999963113292014, 
         0.00004663776030465997}, 
        {-0.9682894407169532, 0.0026179649072038884, -0.2498173437850035}, 
        {-0.2498163001950426, 0.0007235936260381096, 0.968292978684202}}]] + 
   6*Abs[MyDot[{{-0.002715723275723763, -0.9682894407169532, 
        -0.2498163001950426}, {-0.9999963113292014, 0.0026179649072038884, 
        0.0007235936260381096}, 
       {0.00004663776030465997, -0.2498173437850035, 0.968292978684202}}, 
      {{-0.002715723275723763, -0.9999963113292014, 0.00004663776030465997}, 
       {-0.9682894407169532, 0.0026179649072038884, -0.2498173437850035}, 
       {-0.2498163001950426, 0.0007235936260381096, 0.968292978684202}}]]
diagvc[{3, 1}]
Generating matrix and saving to ham.SIDE..QS_3.1
hamil={{4*delta + U + coefzeta[1, 0]}}
dim={1, 1}
det[vec]=1. 1-abs=0.
orthogonality check=Abs[-1 + MyDot[{{1.}}, {{1.}}]]
Lowest energies (absolute):{-0.1226853653356458, -0.03565471649266649, 
   -0.03565471649266638, -0.03565467585295268, -0.03565467585295265, 
   -3.1191420868996954*^-8, 0., 0.3773146170966112, 0.37731461709661174, 
   0.4643450623495911, 0.46434506234959116, 0.4643450623495913, 
   0.46434529713101147, 0.4999997382025437, 0.4999997382025439, 
   0.5000002617974564, 0.5000002617974566, 0.5000002617974572, 
   0.5356547028689883, 0.5356547028689884}
Lowest energies (GS shifted):{0., 0.08703064884297931, 0.08703064884297941, 
   0.08703068948269312, 0.08703068948269314, 0.12268533414422493, 
   0.1226853653356458, 0.49999998243225696, 0.4999999824322575, 
   0.5870304276852369, 0.5870304276852369, 0.587030427685237, 
   0.5870306624666572, 0.6226851035381895, 0.6226851035381897, 
   0.6226856271331022, 0.6226856271331024, 0.6226856271331029, 
   0.6583400682046341, 0.6583400682046342}
Scale factor SCALE(Ninit):1.0201394465967895
Lowest energies (shifted and scaled):
  {0., 0.08531250225968197, 0.08531250225968207, 0.08531254209709237, 
   0.0853125420970924, 0.1202632978790363, 0.1202633284546806, 
   0.4901290545133504, 0.490129054513351, 0.5754413572022776, 
   0.5754413572022776, 0.5754413572022777, 0.5754415873486768, 
   0.6103921435598657, 0.6103921435598659, 0.6103926568180428, 
   0.610392656818043, 0.6103926568180436, 0.6453432130292314, 
   0.6453432130292315}
ireducTable: f[0]
Loading module operators.m
"operators.m $Id: operators.m 877 2009-12-16 16:34:11Z zitko $"
SaSd op.SIDE..QS.SaSd (-nc[a[0, 0], d[0, 0], a[1, 0], d[1, 0]] + 
     nc[a[0, 0], d[0, 1], a[1, 0], d[1, 1]] - 
     2*nc[a[0, 0], d[0, 1], a[1, 1], d[1, 0]] - 
     2*nc[a[0, 1], d[0, 0], a[1, 0], d[1, 1]] + 
     nc[a[0, 1], d[0, 0], a[1, 1], d[1, 0]] - 
     nc[a[0, 1], d[0, 1], a[1, 1], d[1, 1]])/4
$Id: operators.m 877 2009-12-16 16:34:11Z zitko $ DONE
Loading module customoperators.m
"customoperators $Id: customoperators.m,v 1.3 2010/04/07 09:12:36 rok Exp rok \
$"
Customoperators done.
Loading module modeloperators.m
Can't load modeloperators.m. Continuing.
-- maketable[] done --
Timing report
{basis, 0.083974`5.3756898342672725}
{ham, 0.22926`4.811873281568025}
{maketable, 1.559546`6.644543182625387}
{xi, 5.33842`7.178957732394243}
{_, 0}
data
gammaPol=0.13819765978853418

"Success!"


