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=0.01 Gamma=0.0006 delta=0 t=0.
epsilon=-0.005
QS ../dqd.m 
Options:
Options: {}
Loading module models.m
"models $Id: models.m 877 2009-12-16 16:34:11Z zitko $"
Loading module custommodels.m
models $Id: custommodels.m,v 1.2 2011/10/21 16:59:38 rokzitko Exp rokzitko $
custommodels.m done
Loading module ../dqd.m
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
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 + gammaPol*nc[a[0, 0], f[1, 0, 0]] + 
   delta*nc[a[0, 1], a[1, 1]] - (U*nc[a[0, 1], a[1, 1]])/2 + 
   gammaPol*nc[a[0, 1], f[1, 0, 1]] + 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]] + 
   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]] + gammaPol*nc[f[0, 0, 0], a[1, 0]] + 
   gammaPolCh[1]*nc[f[0, 0, 0], d[1, 0]] + 
   coefzeta[1, 0]*nc[f[0, 0, 0], f[1, 0, 0]] + 
   gammaPol*nc[f[0, 0, 1], a[1, 1]] + 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.model..QS
PREC=1000
Tmin=1.*^-10
Tmin=1.*^-10 ==> Nmax=66
DISCNMAX=66
mMAX=132
Diagonalisation.
Discretization checksum [-1] (channel 1): \
9.183549615799121156005754197048794358`10.*^-41
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.
"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
Precision last xi:933.5204393616123
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: , ../dqd.m,  Variant: , ,  Channels: , 1}, 
   {# U=, 0.01,  Gamma=, 0.0006,  Delta=, 0,  t=, 0.}, {#}, 
   {# Gamma^(1/2)=, Sqrt[theta0]/Sqrt[Pi]}, {# Lambda=, 2.,  BAND=, flat}, 
   {# DISCRETIZATION=, Z,  z=, 1.}, {# ops=, n_d^2 n_a^2 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.0012"}
gammaPolCh={"0.019544100476116797"}
checkdefinitions[] -> 0.12635280380893438
calcgsenergy[]
diagvc[{-3, 1}]
Generating matrix and saving to ham.model..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.model..QS_-2.2
hamil={{U, gammaPolCh[1], gammaPol}, 
   {gammaPolCh[1], delta + U/2 - coefzeta[1, 0], 0}, 
   {gammaPol, 0, delta + U/2 - coefzeta[1, 0]}}
dim={3, 3}
det[vec]=-0.9999999999999999 1-abs=1.1102230246251565*^-16
orthogonality check=3*Abs[-1 + 
      MyDot[{{0.6745063324698264, -0.5220350598657642, 
         -0.5220350598657643}, 
        {0., -0.7071067811865476, 0.7071067811865475}, 
        {0.7382690616964144, 0.47694800164268214, 0.47694800164268225}}, 
       {{0.6745063324698264, 0., 0.7382690616964144}, 
        {-0.5220350598657642, -0.7071067811865476, 0.47694800164268214}, 
        {-0.5220350598657643, 0.7071067811865475, 0.47694800164268225}}]] + 
   6*Abs[MyDot[{{0.6745063324698264, -0.5220350598657642, 
        -0.5220350598657643}, {0., -0.7071067811865476, 0.7071067811865475}, 
       {0.7382690616964144, 0.47694800164268214, 0.47694800164268225}}, 
      {{0.6745063324698264, 0., 0.7382690616964144}, 
       {-0.5220350598657642, -0.7071067811865476, 0.47694800164268214}, 
       {-0.5220350598657643, 0.7071067811865475, 0.47694800164268225}}]]
diagvc[{-1, 1}]
Generating matrix and saving to ham.model..QS_-1.1
hamil={{U + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], Sqrt[2]*gammaPol, 0, 0, 
    0}, {Sqrt[2]*gammaPolCh[1], delta + U/2, 0, Sqrt[2]*gammaPolCh[1], 
    gammaPol, 0}, {Sqrt[2]*gammaPol, 0, delta + U/2, 0, gammaPolCh[1], 
    Sqrt[2]*gammaPol}, {0, Sqrt[2]*gammaPolCh[1], 0, 
    2*delta + U - coefzeta[1, 0], 0, 0}, 
   {0, gammaPol, gammaPolCh[1], 0, 2*delta - coefzeta[1, 0], 0}, 
   {0, 0, Sqrt[2]*gammaPol, 0, 0, 2*delta + U - coefzeta[1, 0]}}
dim={6, 6}
det[vec]=-1.0000000000000002 1-abs=-2.220446049250313*^-16
orthogonality check=6*Abs[-1 + 
      MyDot[{{-0.4704027095450446, 0.5039440838066555, 0.5039440838066556, 
         -0.2352013547725223, -0.40020326969124287, -0.23520135477252235}, 
        {1.5087419797821271*^-16, -0.5220350598657644, 0.5220350598657642, 
         0.47694800164268236, 1.6184470817437574*^-16, 
         -0.47694800164268214}, 
        {-0.409480948076921, 0.05544303014664071, 0.05544303014664078, 
         -0.2047404740384601, 0.8615916495413338, -0.2047404740384602}, 
        {-0.5773502691896261, -1.5325028301493154*^-16, 
         -1.799521496441761*^-17, 0.5773502691896257, 
         -3.354491875516997*^-16, 0.5773502691896254}, 
        {-9.926578512486655*^-17, -0.4769480016426821, 
         0.47694800164268203, -0.5220350598657642, 
         1.1618819014753942*^-17, 0.5220350598657645}, 
        {-0.5269851142882245, -0.49291625130891076, -0.492916251308911, 
         -0.26349255714411235, -0.31224537205390523, -0.2634925571441123}}, 
       {{-0.4704027095450446, 1.5087419797821271*^-16, -0.409480948076921, 
         -0.5773502691896261, -9.926578512486655*^-17, -0.5269851142882245}, 
        {0.5039440838066555, -0.5220350598657644, 0.05544303014664071, 
         -1.5325028301493154*^-16, -0.4769480016426821, 
         -0.49291625130891076}, 
        {0.5039440838066556, 0.5220350598657642, 0.05544303014664078, 
         -1.799521496441761*^-17, 0.47694800164268203, -0.492916251308911}, 
        {-0.2352013547725223, 0.47694800164268236, -0.2047404740384601, 
         0.5773502691896257, -0.5220350598657642, -0.26349255714411235}, 
        {-0.40020326969124287, 1.6184470817437574*^-16, 0.8615916495413338, 
         -3.354491875516997*^-16, 1.1618819014753942*^-17, 
         -0.31224537205390523}, 
        {-0.23520135477252235, -0.47694800164268214, -0.2047404740384602, 
         0.5773502691896254, 0.5220350598657645, -0.2634925571441123}}]] + 
   30*Abs[MyDot[{{-0.4704027095450446, 0.5039440838066555, 
        0.5039440838066556, -0.2352013547725223, -0.40020326969124287, 
        -0.23520135477252235}, 
       {1.5087419797821271*^-16, -0.5220350598657644, 0.5220350598657642, 
        0.47694800164268236, 1.6184470817437574*^-16, -0.47694800164268214}, 
       {-0.409480948076921, 0.05544303014664071, 0.05544303014664078, 
        -0.2047404740384601, 0.8615916495413338, -0.2047404740384602}, 
       {-0.5773502691896261, -1.5325028301493154*^-16, 
        -1.799521496441761*^-17, 0.5773502691896257, 
        -3.354491875516997*^-16, 0.5773502691896254}, 
       {-9.926578512486655*^-17, -0.4769480016426821, 0.47694800164268203, 
        -0.5220350598657642, 1.1618819014753942*^-17, 0.5220350598657645}, 
       {-0.5269851142882245, -0.49291625130891076, -0.492916251308911, 
        -0.26349255714411235, -0.31224537205390523, -0.2634925571441123}}, 
      {{-0.4704027095450446, 1.5087419797821271*^-16, -0.409480948076921, 
        -0.5773502691896261, -9.926578512486655*^-17, -0.5269851142882245}, 
       {0.5039440838066555, -0.5220350598657644, 0.05544303014664071, 
        -1.5325028301493154*^-16, -0.4769480016426821, -0.49291625130891076}, 
       {0.5039440838066556, 0.5220350598657642, 0.05544303014664078, 
        -1.799521496441761*^-17, 0.47694800164268203, -0.492916251308911}, 
       {-0.2352013547725223, 0.47694800164268236, -0.2047404740384601, 
        0.5773502691896257, -0.5220350598657642, -0.26349255714411235}, 
       {-0.40020326969124287, 1.6184470817437574*^-16, 0.8615916495413338, 
        -3.354491875516997*^-16, 1.1618819014753942*^-17, 
        -0.31224537205390523}, {-0.23520135477252235, -0.47694800164268214, 
        -0.2047404740384602, 0.5773502691896254, 0.5220350598657645, 
        -0.2634925571441123}}]]
diagvc[{-1, 3}]
Generating matrix and saving to ham.model..QS_-1.3
hamil={{delta + U/2, 0, -gammaPol}, {0, delta + U/2, gammaPolCh[1]}, 
   {-gammaPol, gammaPolCh[1], 2*delta - coefzeta[1, 0]}}
dim={3, 3}
det[vec]=-1. 1-abs=0.
orthogonality check=3*Abs[-1 + 
      MyDot[{{0.47694800164268236, -0.47694800164268225, 
         0.7382690616964142}, {-0.7071067811865475, -0.7071067811865477, 
         -1.4484941024406337*^-16}, 
        {-0.5220350598657643, 0.5220350598657644, 0.6745063324698263}}, 
       {{0.47694800164268236, -0.7071067811865475, -0.5220350598657643}, 
        {-0.47694800164268225, -0.7071067811865477, 0.5220350598657644}, 
        {0.7382690616964142, -1.4484941024406337*^-16, 0.6745063324698263}}]] \
+ 6*Abs[MyDot[{{0.47694800164268236, -0.47694800164268225, 
        0.7382690616964142}, {-0.7071067811865475, -0.7071067811865477, 
        -1.4484941024406337*^-16}, 
       {-0.5220350598657643, 0.5220350598657644, 0.6745063324698263}}, 
      {{0.47694800164268236, -0.7071067811865475, -0.5220350598657643}, 
       {-0.47694800164268225, -0.7071067811865477, 0.5220350598657644}, 
       {0.7382690616964142, -1.4484941024406337*^-16, 0.6745063324698263}}]]
diagvc[{0, 2}]
Generating matrix and saving to ham.model..QS_0.2
hamil={{delta + U/2 + coefzeta[1, 0], 0, -gammaPolCh[1], 
    -(gammaPol/Sqrt[2]), 0, -(Sqrt[3/2]*gammaPol), 0, 0}, 
   {0, delta + U/2 + coefzeta[1, 0], 0, -(gammaPolCh[1]/Sqrt[2]), -gammaPol, 
    Sqrt[3/2]*gammaPolCh[1], 0, 0}, 
   {-gammaPolCh[1], 0, 2*delta + U, 0, 0, 0, gammaPol, 0}, 
   {-(gammaPol/Sqrt[2]), -(gammaPolCh[1]/Sqrt[2]), 0, 2*delta, 0, 0, 
    -(gammaPolCh[1]/Sqrt[2]), -(gammaPol/Sqrt[2])}, 
   {0, -gammaPol, 0, 0, 2*delta + U, 0, 0, gammaPolCh[1]}, 
   {-(Sqrt[3/2]*gammaPol), Sqrt[3/2]*gammaPolCh[1], 0, 0, 0, 2*delta, 
    Sqrt[3/2]*gammaPolCh[1], -(Sqrt[3/2]*gammaPol)}, 
   {0, 0, gammaPol, -(gammaPolCh[1]/Sqrt[2]), 0, Sqrt[3/2]*gammaPolCh[1], 
    (6*delta + U - 2*coefzeta[1, 0])/2, 0}, 
   {0, 0, 0, -(gammaPol/Sqrt[2]), gammaPolCh[1], -(Sqrt[3/2]*gammaPol), 0, 
    (6*delta + U - 2*coefzeta[1, 0])/2}}
dim={8, 8}
det[vec]=1.0000000000000002 1-abs=-2.220446049250313*^-16
orthogonality check=8*Abs[-1 + 
      MyDot[{{-0.3485444238575831, 0.34854442385758366, 
         -0.2207908199734328, 1.7184579567826543*^-16, 
         0.22079081997343297, -0.6454223160242133, 0.3485444238575832, 
         -0.3485444238575833}, 
        {0.33725316623491225, 0.33725316623491286, 
         -9.747953736116397*^-16, 0.738269061696414, 
         -8.058984862367899*^-16, 2.1638303670558363*^-17, 
         0.33725316623491397, 0.337253166234914}, 
        {-0.3691345308482078, -0.369134530848208, -0.47694800164268175, 
         -1.1654089152046685*^-15, -0.4769480016426823, 
         1.4813047012023047*^-18, 0.3691345308482067, 0.3691345308482066}, 
        {0.5000000000000002, -0.5000000000000001, 
         -1.0740639822002585*^-16, 8.695132016744295*^-17, 
         1.3359829807546412*^-16, -9.632964910710786*^-19, 
         0.5000000000000001, -0.49999999999999983}, 
        {-0.0392041425862201, 0.039204142586219785, -0.6092372980043806, 
         1.8607946320307096*^-16, 0.6092372980043804, 0.5015096910897742, 
         0.03920414258621973, -0.03920414258621967}, 
        {0.36913453084820624, 0.36913453084820635, 
         1.5906872173621088*^-15, -0.6745063324698266, 
         1.922615712468546*^-15, 2.9916187474620113*^-16, 
         0.369134530848208, 0.369134530848208}, 
        {0.33725316623491464, 0.3372531662349146, -0.5220350598657638, 
         -2.5343496832391764*^-15, -0.5220350598657646, 
         -2.1109755111471655*^-17, -0.33725316623491175, 
         -0.3372531662349118}, 
        {0.35634227899852805, -0.35634227899852855, -0.2829864458517064, 
         -5.0977500499571324*^-17, 0.28298644585170635, 
         -0.5761233060040005, -0.3563422789985279, 0.35634227899852816}}, 
       {{-0.3485444238575831, 0.33725316623491225, -0.3691345308482078, 
         0.5000000000000002, -0.0392041425862201, 0.36913453084820624, 
         0.33725316623491464, 0.35634227899852805}, 
        {0.34854442385758366, 0.33725316623491286, -0.369134530848208, 
         -0.5000000000000001, 0.039204142586219785, 0.36913453084820635, 
         0.3372531662349146, -0.35634227899852855}, 
        {-0.2207908199734328, -9.747953736116397*^-16, 
         -0.47694800164268175, -1.0740639822002585*^-16, 
         -0.6092372980043806, 1.5906872173621088*^-15, -0.5220350598657638, 
         -0.2829864458517064}, 
        {1.7184579567826543*^-16, 0.738269061696414, 
         -1.1654089152046685*^-15, 8.695132016744295*^-17, 
         1.8607946320307096*^-16, -0.6745063324698266, 
         -2.5343496832391764*^-15, -5.0977500499571324*^-17}, 
        {0.22079081997343297, -8.058984862367899*^-16, -0.4769480016426823, 
         1.3359829807546412*^-16, 0.6092372980043804, 
         1.922615712468546*^-15, -0.5220350598657646, 0.28298644585170635}, 
        {-0.6454223160242133, 2.1638303670558363*^-17, 
         1.4813047012023047*^-18, -9.632964910710786*^-19, 
         0.5015096910897742, 2.9916187474620113*^-16, 
         -2.1109755111471655*^-17, -0.5761233060040005}, 
        {0.3485444238575832, 0.33725316623491397, 0.3691345308482067, 
         0.5000000000000001, 0.03920414258621973, 0.369134530848208, 
         -0.33725316623491175, -0.3563422789985279}, 
        {-0.3485444238575833, 0.337253166234914, 0.3691345308482066, 
         -0.49999999999999983, -0.03920414258621967, 0.369134530848208, 
         -0.3372531662349118, 0.35634227899852816}}]] + 
   56*Abs[MyDot[{{-0.3485444238575831, 0.34854442385758366, 
        -0.2207908199734328, 1.7184579567826543*^-16, 0.22079081997343297, 
        -0.6454223160242133, 0.3485444238575832, -0.3485444238575833}, 
       {0.33725316623491225, 0.33725316623491286, -9.747953736116397*^-16, 
        0.738269061696414, -8.058984862367899*^-16, 
        2.1638303670558363*^-17, 0.33725316623491397, 0.337253166234914}, 
       {-0.3691345308482078, -0.369134530848208, -0.47694800164268175, 
        -1.1654089152046685*^-15, -0.4769480016426823, 
        1.4813047012023047*^-18, 0.3691345308482067, 0.3691345308482066}, 
       {0.5000000000000002, -0.5000000000000001, -1.0740639822002585*^-16, 
        8.695132016744295*^-17, 1.3359829807546412*^-16, 
        -9.632964910710786*^-19, 0.5000000000000001, -0.49999999999999983}, 
       {-0.0392041425862201, 0.039204142586219785, -0.6092372980043806, 
        1.8607946320307096*^-16, 0.6092372980043804, 0.5015096910897742, 
        0.03920414258621973, -0.03920414258621967}, 
       {0.36913453084820624, 0.36913453084820635, 1.5906872173621088*^-15, 
        -0.6745063324698266, 1.922615712468546*^-15, 
        2.9916187474620113*^-16, 0.369134530848208, 0.369134530848208}, 
       {0.33725316623491464, 0.3372531662349146, -0.5220350598657638, 
        -2.5343496832391764*^-15, -0.5220350598657646, 
        -2.1109755111471655*^-17, -0.33725316623491175, -0.3372531662349118}\
, {0.35634227899852805, -0.35634227899852855, -0.2829864458517064, 
        -5.0977500499571324*^-17, 0.28298644585170635, -0.5761233060040005, 
        -0.3563422789985279, 0.35634227899852816}}, 
      {{-0.3485444238575831, 0.33725316623491225, -0.3691345308482078, 
        0.5000000000000002, -0.0392041425862201, 0.36913453084820624, 
        0.33725316623491464, 0.35634227899852805}, 
       {0.34854442385758366, 0.33725316623491286, -0.369134530848208, 
        -0.5000000000000001, 0.039204142586219785, 0.36913453084820635, 
        0.3372531662349146, -0.35634227899852855}, 
       {-0.2207908199734328, -9.747953736116397*^-16, -0.47694800164268175, 
        -1.0740639822002585*^-16, -0.6092372980043806, 
        1.5906872173621088*^-15, -0.5220350598657638, -0.2829864458517064}, 
       {1.7184579567826543*^-16, 0.738269061696414, 
        -1.1654089152046685*^-15, 8.695132016744295*^-17, 
        1.8607946320307096*^-16, -0.6745063324698266, 
        -2.5343496832391764*^-15, -5.0977500499571324*^-17}, 
       {0.22079081997343297, -8.058984862367899*^-16, -0.4769480016426823, 
        1.3359829807546412*^-16, 0.6092372980043804, 1.922615712468546*^-15, 
        -0.5220350598657646, 0.28298644585170635}, 
       {-0.6454223160242133, 2.1638303670558363*^-17, 
        1.4813047012023047*^-18, -9.632964910710786*^-19, 
        0.5015096910897742, 2.9916187474620113*^-16, 
        -2.1109755111471655*^-17, -0.5761233060040005}, 
       {0.3485444238575832, 0.33725316623491397, 0.3691345308482067, 
        0.5000000000000001, 0.03920414258621973, 0.369134530848208, 
        -0.33725316623491175, -0.3563422789985279}, 
       {-0.3485444238575833, 0.337253166234914, 0.3691345308482066, 
        -0.49999999999999983, -0.03920414258621967, 0.369134530848208, 
        -0.3372531662349118, 0.35634227899852816}}]]
diagvc[{0, 4}]
Generating matrix and saving to ham.model..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.model..QS_1.1
hamil={{2*delta + U + coefzeta[1, 0], 0, 0, Sqrt[2]*gammaPol, 0, 0}, 
   {0, 2*delta + coefzeta[1, 0], 0, -gammaPolCh[1], -gammaPol, 0}, 
   {0, 0, 2*delta + U + coefzeta[1, 0], 0, Sqrt[2]*gammaPolCh[1], 0}, 
   {Sqrt[2]*gammaPol, -gammaPolCh[1], 0, (6*delta + U)/2, 0, 
    Sqrt[2]*gammaPol}, {0, -gammaPol, Sqrt[2]*gammaPolCh[1], 0, 
    (6*delta + U)/2, Sqrt[2]*gammaPolCh[1]}, 
   {0, 0, 0, Sqrt[2]*gammaPol, Sqrt[2]*gammaPolCh[1], 
    4*delta + U - coefzeta[1, 0]}}
dim={6, 6}
det[vec]=1.0000000000000004 1-abs=-4.440892098500626*^-16
orthogonality check=6*Abs[-1 + 
      MyDot[{{0.2352013547725224, -0.40020326969124304, 
         0.2352013547725224, -0.5039440838066552, -0.5039440838066556, 
         0.47040270954504476}, 
        {0.4769480016426824, 5.20617123369026*^-17, -0.47694800164268203, 
         -0.5220350598657643, 0.5220350598657643, 1.2729865179820583*^-16}, 
        {0.20474047403846005, 0.8615916495413339, 0.2047404740384602, 
         -0.0554430301466407, -0.05544303014664088, 0.409480948076921}, 
        {0.5773502691896255, 2.1261585767267003*^-18, 0.5773502691896262, 
         4.2934406030425963*^-17, -1.432231069853351*^-16, 
         -0.5773502691896255}, 
        {-0.5220350598657646, 1.1384289196882621*^-18, 0.5220350598657642, 
         -0.476948001642682, 0.4769480016426822, -1.9641566181699297*^-16}, 
        {0.26349255714411257, -0.31224537205390546, 0.26349255714411257, 
         0.49291625130891076, 0.4929162513089105, 0.5269851142882247}}, 
       {{0.2352013547725224, 0.4769480016426824, 0.20474047403846005, 
         0.5773502691896255, -0.5220350598657646, 0.26349255714411257}, 
        {-0.40020326969124304, 5.20617123369026*^-17, 0.8615916495413339, 
         2.1261585767267003*^-18, 1.1384289196882621*^-18, 
         -0.31224537205390546}, 
        {0.2352013547725224, -0.47694800164268203, 0.2047404740384602, 
         0.5773502691896262, 0.5220350598657642, 0.26349255714411257}, 
        {-0.5039440838066552, -0.5220350598657643, -0.0554430301466407, 
         4.2934406030425963*^-17, -0.476948001642682, 0.49291625130891076}, 
        {-0.5039440838066556, 0.5220350598657643, -0.05544303014664088, 
         -1.432231069853351*^-16, 0.4769480016426822, 0.4929162513089105}, 
        {0.47040270954504476, 1.2729865179820583*^-16, 0.409480948076921, 
         -0.5773502691896255, -1.9641566181699297*^-16, 0.5269851142882247}}]\
] + 30*Abs[MyDot[{{0.2352013547725224, -0.40020326969124304, 
        0.2352013547725224, -0.5039440838066552, -0.5039440838066556, 
        0.47040270954504476}, {0.4769480016426824, 5.20617123369026*^-17, 
        -0.47694800164268203, -0.5220350598657643, 0.5220350598657643, 
        1.2729865179820583*^-16}, 
       {0.20474047403846005, 0.8615916495413339, 0.2047404740384602, 
        -0.0554430301466407, -0.05544303014664088, 0.409480948076921}, 
       {0.5773502691896255, 2.1261585767267003*^-18, 0.5773502691896262, 
        4.2934406030425963*^-17, -1.432231069853351*^-16, 
        -0.5773502691896255}, {-0.5220350598657646, 
        1.1384289196882621*^-18, 0.5220350598657642, -0.476948001642682, 
        0.4769480016426822, -1.9641566181699297*^-16}, 
       {0.26349255714411257, -0.31224537205390546, 0.26349255714411257, 
        0.49291625130891076, 0.4929162513089105, 0.5269851142882247}}, 
      {{0.2352013547725224, 0.4769480016426824, 0.20474047403846005, 
        0.5773502691896255, -0.5220350598657646, 0.26349255714411257}, 
       {-0.40020326969124304, 5.20617123369026*^-17, 0.8615916495413339, 
        2.1261585767267003*^-18, 1.1384289196882621*^-18, 
        -0.31224537205390546}, {0.2352013547725224, -0.47694800164268203, 
        0.2047404740384602, 0.5773502691896262, 0.5220350598657642, 
        0.26349255714411257}, {-0.5039440838066552, -0.5220350598657643, 
        -0.0554430301466407, 4.2934406030425963*^-17, -0.476948001642682, 
        0.49291625130891076}, {-0.5039440838066556, 0.5220350598657643, 
        -0.05544303014664088, -1.432231069853351*^-16, 0.4769480016426822, 
        0.4929162513089105}, {0.47040270954504476, 1.2729865179820583*^-16, 
        0.409480948076921, -0.5773502691896255, -1.9641566181699297*^-16, 
        0.5269851142882247}}]]
diagvc[{1, 3}]
Generating matrix and saving to ham.model..QS_1.3
hamil={{2*delta + coefzeta[1, 0], -gammaPolCh[1], gammaPol}, 
   {-gammaPolCh[1], (6*delta + U)/2, 0}, {gammaPol, 0, (6*delta + U)/2}}
dim={3, 3}
det[vec]=0.9999999999999999 1-abs=1.1102230246251565*^-16
orthogonality check=3*Abs[-1 + 
      MyDot[{{-0.7382690616964144, -0.47694800164268214, 
         0.47694800164268225}, 
        {0., 0.7071067811865476, 0.7071067811865475}, 
        {-0.6745063324698264, 0.5220350598657642, -0.5220350598657643}}, 
       {{-0.7382690616964144, 0., -0.6745063324698264}, 
        {-0.47694800164268214, 0.7071067811865476, 0.5220350598657642}, 
        {0.47694800164268225, 0.7071067811865475, -0.5220350598657643}}]] + 
   6*Abs[MyDot[{{-0.7382690616964144, -0.47694800164268214, 
        0.47694800164268225}, {0., 0.7071067811865476, 0.7071067811865475}, 
       {-0.6745063324698264, 0.5220350598657642, -0.5220350598657643}}, 
      {{-0.7382690616964144, 0., -0.6745063324698264}, 
       {-0.47694800164268214, 0.7071067811865476, 0.5220350598657642}, 
       {0.47694800164268225, 0.7071067811865475, -0.5220350598657643}}]]
diagvc[{2, 2}]
Generating matrix and saving to ham.model..QS_2.2
hamil={{3*delta + U/2 + coefzeta[1, 0], 0, -gammaPol}, 
   {0, 3*delta + U/2 + coefzeta[1, 0], -gammaPolCh[1]}, 
   {-gammaPol, -gammaPolCh[1], 4*delta + U}}
dim={3, 3}
det[vec]=0.9999999999999999 1-abs=1.1102230246251565*^-16
orthogonality check=3*Abs[-1 + 
      MyDot[{{0.5220350598657642, 0.5220350598657642, 0.6745063324698264}, 
        {-0.7071067811865477, 0.7071067811865474, 
         -1.7954387976359953*^-16}, 
        {-0.4769480016426824, -0.47694800164268236, 0.7382690616964143}}, 
       {{0.5220350598657642, -0.7071067811865477, -0.4769480016426824}, 
        {0.5220350598657642, 0.7071067811865474, -0.47694800164268236}, 
        {0.6745063324698264, -1.7954387976359953*^-16, 0.7382690616964143}}]] \
+ 6*Abs[MyDot[{{0.5220350598657642, 0.5220350598657642, 0.6745063324698264}, 
       {-0.7071067811865477, 0.7071067811865474, -1.7954387976359953*^-16}, 
       {-0.4769480016426824, -0.47694800164268236, 0.7382690616964143}}, 
      {{0.5220350598657642, -0.7071067811865477, -0.4769480016426824}, 
       {0.5220350598657642, 0.7071067811865474, -0.47694800164268236}, 
       {0.6745063324698264, -1.7954387976359953*^-16, 0.7382690616964143}}]]
diagvc[{3, 1}]
Generating matrix and saving to ham.model..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.051705348447752555, -0.04922065637225061, 
   -0.04922065637225057, -0.025252364346864172, -0.02525236434686417, 
   -0.025252364346864144, -0.020252364346864192, -0.02025236434686418, 
   -0.02025236434686417, -0.020252364346864168, -0.020252364346864164, 0., 
   0.002515307924498015, 0.0025153079244980216, 0.004999999999999996, 
   0.004999999999999998, 0.005, 0.005, 0.005000000000000003, 
   0.007484692075501967}
Lowest energies (GS shifted):{0., 0.0024846920755019486, 
   0.0024846920755019833, 0.026452984100888383, 0.026452984100888387, 
   0.02645298410088841, 0.03145298410088836, 0.03145298410088837, 
   0.031452984100888384, 0.031452984100888384, 0.03145298410088839, 
   0.051705348447752555, 0.05422065637225057, 0.05422065637225058, 
   0.05670534844775255, 0.05670534844775255, 0.05670534844775255, 
   0.05670534844775255, 0.05670534844775256, 0.05919004052325452}
Scale factor SCALE(Ninit):1.0201394465967895
Lowest energies (shifted and scaled):
  {0., 0.0024356396410225513, 0.0024356396410225856, 0.025930753083939842, 
   0.025930753083939845, 0.02593075308393987, 0.03083204380128256, 
   0.030832043801282568, 0.030832043801282582, 0.030832043801282582, 
   0.03083204380128259, 0.05068458887678825, 0.0531502399531084, 
   0.05315023995310841, 0.055585879594130985, 0.055585879594130985, 
   0.055585879594130985, 0.055585879594130985, 0.05558587959413099, 
   0.05802151923515356}
ireducTable: f[0]
Loading module operators.m
"operators.m $Id: operators.m 877 2009-12-16 16:34:11Z zitko $"
n_d^2 op.model..QS.n_d^2 nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]] - 
   2*nc[d[0, 0], d[0, 1], d[1, 0], d[1, 1]]
n_a^2 op.model..QS.n_a^2 nc[a[0, 0], a[1, 0]] + nc[a[0, 1], a[1, 1]] - 
   2*nc[a[0, 0], a[0, 1], a[1, 0], a[1, 1]]
SaSd op.model..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.108788`5.488125986089707}
{ham, 0.304861`4.935646863357802}
{maketable, 2.295935`6.812504582125688}
{xi, 3.423227`6.985980692393903}
{_, 0}
data
gammaPol=0.019544100476116797

"Success!"


