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 Gamma=0.1 delta=0 t=0.
epsilon=0
{{omega, 0.1}, {g, 0}}
{extrag, extraomega}
QS model.m 
Options:Nph=10
Options: {Nph=10}
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 model.m
def1ch, NRDOTS=1
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=1
NRDOTS:1
CHANNELS:1
basis:{d[], f[0]}
lrchain:{}
lrextrarule:{}
NROPS:2
Hamiltonian generated. U/2 - coefzeta[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]] + 
   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]] - gep*nc[ket[0], bra[1]] - 
   gep*nc[ket[1], bra[0]] + omega*nc[ket[1], bra[1]] - 
   Sqrt[2]*gep*nc[ket[1], bra[2]] - Sqrt[2]*gep*nc[ket[2], bra[1]] + 
   2*omega*nc[ket[2], bra[2]] - Sqrt[3]*gep*nc[ket[2], bra[3]] - 
   Sqrt[3]*gep*nc[ket[3], bra[2]] + 3*omega*nc[ket[3], bra[3]] - 
   2*gep*nc[ket[3], bra[4]] - 2*gep*nc[ket[4], bra[3]] + 
   4*omega*nc[ket[4], bra[4]] - Sqrt[5]*gep*nc[ket[4], bra[5]] - 
   Sqrt[5]*gep*nc[ket[5], bra[4]] + 5*omega*nc[ket[5], bra[5]] - 
   Sqrt[6]*gep*nc[ket[5], bra[6]] - Sqrt[6]*gep*nc[ket[6], bra[5]] + 
   6*omega*nc[ket[6], bra[6]] - Sqrt[7]*gep*nc[ket[6], bra[7]] - 
   Sqrt[7]*gep*nc[ket[7], bra[6]] + 7*omega*nc[ket[7], bra[7]] - 
   2*Sqrt[2]*gep*nc[ket[7], bra[8]] - 2*Sqrt[2]*gep*nc[ket[8], bra[7]] + 
   8*omega*nc[ket[8], bra[8]] - 3*gep*nc[ket[8], bra[9]] - 
   3*gep*nc[ket[9], bra[8]] + 9*omega*nc[ket[9], bra[9]] - 
   Sqrt[10]*gep*nc[ket[9], bra[10]] - Sqrt[10]*gep*nc[ket[10], bra[9]] + 
   10*omega*nc[ket[10], bra[10]] - U*nc[d[0, 0], d[0, 1], d[1, 0], d[1, 1]] + 
   gep*nc[d[0, 0], d[1, 0], ket[0], bra[1]] + 
   gep*nc[d[0, 0], d[1, 0], ket[1], bra[0]] + 
   Sqrt[2]*gep*nc[d[0, 0], d[1, 0], ket[1], bra[2]] + 
   Sqrt[2]*gep*nc[d[0, 0], d[1, 0], ket[2], bra[1]] + 
   Sqrt[3]*gep*nc[d[0, 0], d[1, 0], ket[2], bra[3]] + 
   Sqrt[3]*gep*nc[d[0, 0], d[1, 0], ket[3], bra[2]] + 
   2*gep*nc[d[0, 0], d[1, 0], ket[3], bra[4]] + 
   2*gep*nc[d[0, 0], d[1, 0], ket[4], bra[3]] + 
   Sqrt[5]*gep*nc[d[0, 0], d[1, 0], ket[4], bra[5]] + 
   Sqrt[5]*gep*nc[d[0, 0], d[1, 0], ket[5], bra[4]] + 
   Sqrt[6]*gep*nc[d[0, 0], d[1, 0], ket[5], bra[6]] + 
   Sqrt[6]*gep*nc[d[0, 0], d[1, 0], ket[6], bra[5]] + 
   Sqrt[7]*gep*nc[d[0, 0], d[1, 0], ket[6], bra[7]] + 
   Sqrt[7]*gep*nc[d[0, 0], d[1, 0], ket[7], bra[6]] + 
   2*Sqrt[2]*gep*nc[d[0, 0], d[1, 0], ket[7], bra[8]] + 
   2*Sqrt[2]*gep*nc[d[0, 0], d[1, 0], ket[8], bra[7]] + 
   3*gep*nc[d[0, 0], d[1, 0], ket[8], bra[9]] + 
   3*gep*nc[d[0, 0], d[1, 0], ket[9], bra[8]] + 
   Sqrt[10]*gep*nc[d[0, 0], d[1, 0], ket[9], bra[10]] + 
   Sqrt[10]*gep*nc[d[0, 0], d[1, 0], ket[10], bra[9]] + 
   gep*nc[d[0, 1], d[1, 1], ket[0], bra[1]] + 
   gep*nc[d[0, 1], d[1, 1], ket[1], bra[0]] + 
   Sqrt[2]*gep*nc[d[0, 1], d[1, 1], ket[1], bra[2]] + 
   Sqrt[2]*gep*nc[d[0, 1], d[1, 1], ket[2], bra[1]] + 
   Sqrt[3]*gep*nc[d[0, 1], d[1, 1], ket[2], bra[3]] + 
   Sqrt[3]*gep*nc[d[0, 1], d[1, 1], ket[3], bra[2]] + 
   2*gep*nc[d[0, 1], d[1, 1], ket[3], bra[4]] + 
   2*gep*nc[d[0, 1], d[1, 1], ket[4], bra[3]] + 
   Sqrt[5]*gep*nc[d[0, 1], d[1, 1], ket[4], bra[5]] + 
   Sqrt[5]*gep*nc[d[0, 1], d[1, 1], ket[5], bra[4]] + 
   Sqrt[6]*gep*nc[d[0, 1], d[1, 1], ket[5], bra[6]] + 
   Sqrt[6]*gep*nc[d[0, 1], d[1, 1], ket[6], bra[5]] + 
   Sqrt[7]*gep*nc[d[0, 1], d[1, 1], ket[6], bra[7]] + 
   Sqrt[7]*gep*nc[d[0, 1], d[1, 1], ket[7], bra[6]] + 
   2*Sqrt[2]*gep*nc[d[0, 1], d[1, 1], ket[7], bra[8]] + 
   2*Sqrt[2]*gep*nc[d[0, 1], d[1, 1], ket[8], bra[7]] + 
   3*gep*nc[d[0, 1], d[1, 1], ket[8], bra[9]] + 
   3*gep*nc[d[0, 1], d[1, 1], ket[9], bra[8]] + 
   Sqrt[10]*gep*nc[d[0, 1], d[1, 1], ket[9], bra[10]] + 
   Sqrt[10]*gep*nc[d[0, 1], d[1, 1], ket[10], bra[9]]
H-conj[H]=0
SCALE[0]=1.0201394465967895
faktor=1.3862943611198906
Generating basis
Basis states generated.
BASIS NR=16
Adding phonons
Saving basis to basis.model..QS
PREC=1000
Tmin=1.*^-6
Tmin=1.*^-6 ==> Nmax=39
DISCNMAX=39
mMAX=80
Diagonalisation.
Discretization checksum [-1] (channel 1): \
4.135903062765138374357043460349814`10.*^-25
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.
"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
Precision last xi:960.4026965462326
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: , model.m,  Variant: , ,  Channels: , 1}, 
   {# U=, 0,  Gamma=, 0.1,  Delta=, 0,  t=, 0.}, 
   {#,  extrag=, 0,  extraomega=, 0.1}, 
   {# Gamma^(1/2)=, Sqrt[theta0]/Sqrt[Pi]}, {# Lambda=, 2.,  BAND=, flat}, 
   {# DISCRETIZATION=, Z,  z=, 1.}, 
   {# ops=, A_d n_d n_d^2 nph displ displ^2 acr}, {#!7}, 
   {# Number of channels, impurities, subspaces: }, {1, 1, 6}}
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.2"}
gammaPolCh={"0.252313252202016"}
checkdefinitions[] -> 6.509253008808065
calcgsenergy[]
diagvc[{-2, 1}]
Generating matrix and saving to ham.model..QS_-2.1
hamil={{(U - 2*coefzeta[1, 0])/2, -gep, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {-gep, omega + U/2 - coefzeta[1, 0], -(Sqrt[2]*gep), 0, 0, 0, 0, 0, 0, 0, 
    0}, {0, -(Sqrt[2]*gep), (4*omega + U - 2*coefzeta[1, 0])/2, 
    -(Sqrt[3]*gep), 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, -(Sqrt[3]*gep), (6*omega + U - 2*coefzeta[1, 0])/2, -2*gep, 0, 0, 
    0, 0, 0, 0}, {0, 0, 0, -2*gep, (8*omega + U - 2*coefzeta[1, 0])/2, 
    -(Sqrt[5]*gep), 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, -(Sqrt[5]*gep), 5*omega + U/2 - coefzeta[1, 0], 
    -(Sqrt[6]*gep), 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, -(Sqrt[6]*gep), 6*omega + U/2 - coefzeta[1, 0], 
    -(Sqrt[7]*gep), 0, 0, 0}, {0, 0, 0, 0, 0, 0, -(Sqrt[7]*gep), 
    7*omega + U/2 - coefzeta[1, 0], -2*Sqrt[2]*gep, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, -2*Sqrt[2]*gep, 8*omega + U/2 - coefzeta[1, 0], 
    -3*gep, 0}, {0, 0, 0, 0, 0, 0, 0, 0, -3*gep, 
    9*omega + U/2 - coefzeta[1, 0], -(Sqrt[10]*gep)}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[10]*gep), 
    (20*omega + U - 2*coefzeta[1, 0])/2}}
dim={11, 11}
det[vec]=1. 1-abs=0.
orthogonality check=11*Abs[-1 + 
      MyDot[{{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}, 
       {{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}]] + 
   110*Abs[MyDot[{{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}, 
      {{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}]]
diagvc[{-1, 2}]
Generating matrix and saving to ham.model..QS_-1.2
hamil={{U/2, gammaPolCh[1], -gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0}, {gammaPolCh[1], delta - coefzeta[1, 0], 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {-gep, 0, omega + U/2, gammaPolCh[1], -(Sqrt[2]*gep), 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, gammaPolCh[1], delta + omega - coefzeta[1, 0], 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, -(Sqrt[2]*gep), 0, (4*omega + U)/2, gammaPolCh[1], -(Sqrt[3]*gep), 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, gammaPolCh[1], delta + 2*omega - coefzeta[1, 0], 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, -(Sqrt[3]*gep), 0, (6*omega + U)/2, gammaPolCh[1], -2*gep, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, gammaPolCh[1], delta + 3*omega - coefzeta[1, 0], 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, -2*gep, 0, (8*omega + U)/2, gammaPolCh[1], 
    -(Sqrt[5]*gep), 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, gammaPolCh[1], delta + 4*omega - coefzeta[1, 0], 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[5]*gep), 0, 5*omega + U/2, gammaPolCh[1], 
    -(Sqrt[6]*gep), 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, gammaPolCh[1], 
    delta + 5*omega - coefzeta[1, 0], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[6]*gep), 0, 6*omega + U/2, 
    gammaPolCh[1], -(Sqrt[7]*gep), 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, gammaPolCh[1], 
    delta + 6*omega - coefzeta[1, 0], 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[7]*gep), 0, 7*omega + U/2, 
    gammaPolCh[1], -2*Sqrt[2]*gep, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, gammaPolCh[1], 
    delta + 7*omega - coefzeta[1, 0], 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -2*Sqrt[2]*gep, 0, 
    8*omega + U/2, gammaPolCh[1], -3*gep, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, gammaPolCh[1], 
    delta + 8*omega - coefzeta[1, 0], 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -3*gep, 0, 
    9*omega + U/2, gammaPolCh[1], -(Sqrt[10]*gep), 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, gammaPolCh[1], 
    delta + 9*omega - coefzeta[1, 0], 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[10]*gep), 
    0, (20*omega + U)/2, gammaPolCh[1]}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    gammaPolCh[1], delta + 10*omega - coefzeta[1, 0]}}
dim={22, 22}
det[vec]=-0.9999999999999997 1-abs=3.3306690738754696*^-16
orthogonality check=22*Abs[-1 + 
      MyDot[{{-0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., -0.7071067811865476, 0.7071067811865475, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., -0.7071067811865476, 0.7071067811865475, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0.}, {0.7071067811865475, 0.7071067811865476, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0.7071067811865475, 0.7071067811865476, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.7071067811865475, 0.7071067811865476, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0.7071067811865475, 0.7071067811865476, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., -0.7071067811865476, 0.7071067811865475, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., -0.7071067811865476, 0.7071067811865475}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0.7071067811865476, 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0.7071067811865476, 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0.7071067811865476, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865475, 0.7071067811865476, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865475, 0.7071067811865476}}, 
       {{-0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0.7071067811865475, 0., 0., 0., 0., 0., 0.7071067811865476, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 
         0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 
         0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 
         0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0.}\
, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0.7071067811865476}}]] + 
   462*Abs[MyDot[{{-0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., -0.7071067811865476, 0.7071067811865475, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., -0.7071067811865476, 0.7071067811865475, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0.7071067811865475, 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
        0., 0.}, {0., 0., 0.7071067811865475, 0.7071067811865476, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.7071067811865475, 0.7071067811865476, 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0.7071067811865475, 0.7071067811865476, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., -0.7071067811865476, 0.7071067811865475, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., -0.7071067811865476, 0.7071067811865475}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0.7071067811865476, 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0.7071067811865476, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865475, 0.7071067811865476, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.7071067811865475, 0.7071067811865476}}, 
      {{-0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0.7071067811865475, 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 
        0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 
        0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 
        0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 
        0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0.7071067811865476}}]]
diagvc[{0, 1}]
Generating matrix and saving to ham.model..QS_0.1
hamil={{U/2 + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], 0, -gep, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {Sqrt[2]*gammaPolCh[1], delta, Sqrt[2]*gammaPolCh[1], 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, Sqrt[2]*gammaPolCh[1], (4*delta + U - 2*coefzeta[1, 0])/2, 0, 0, gep, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0}, {-gep, 0, 0, omega + U/2 + coefzeta[1, 0], 
    Sqrt[2]*gammaPolCh[1], 0, -(Sqrt[2]*gep), 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, Sqrt[2]*gammaPolCh[1], delta + omega, Sqrt[2]*gammaPolCh[1], 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0}, {0, 0, gep, 0, Sqrt[2]*gammaPolCh[1], 
    2*delta + omega + U/2 - coefzeta[1, 0], 0, 0, Sqrt[2]*gep, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, -(Sqrt[2]*gep), 0, 0, 2*omega + U/2 + coefzeta[1, 0], 
    Sqrt[2]*gammaPolCh[1], 0, -(Sqrt[3]*gep), 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, Sqrt[2]*gammaPolCh[1], delta + 2*omega, 
    Sqrt[2]*gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, Sqrt[2]*gep, 0, 
    Sqrt[2]*gammaPolCh[1], (4*delta + 4*omega + U - 2*coefzeta[1, 0])/2, 0, 
    0, Sqrt[3]*gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0}, {0, 0, 0, 0, 0, 0, -(Sqrt[3]*gep), 0, 0, 
    3*omega + U/2 + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], 0, -2*gep, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[2]*gammaPolCh[1], delta + 3*omega, 
    Sqrt[2]*gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, Sqrt[3]*gep, 0, 
    Sqrt[2]*gammaPolCh[1], (4*delta + 6*omega + U - 2*coefzeta[1, 0])/2, 0, 
    0, 2*gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, -2*gep, 0, 0, 4*omega + U/2 + coefzeta[1, 0], 
    Sqrt[2]*gammaPolCh[1], 0, -(Sqrt[5]*gep), 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    Sqrt[2]*gammaPolCh[1], delta + 4*omega, Sqrt[2]*gammaPolCh[1], 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2*gep, 0, Sqrt[2]*gammaPolCh[1], 
    (4*delta + 8*omega + U - 2*coefzeta[1, 0])/2, 0, 0, Sqrt[5]*gep, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[5]*gep), 0, 0, 
    5*omega + U/2 + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], 0, 
    -(Sqrt[6]*gep), 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[2]*gammaPolCh[1], 
    delta + 5*omega, Sqrt[2]*gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[5]*gep, 0, 
    Sqrt[2]*gammaPolCh[1], 2*delta + 5*omega + U/2 - coefzeta[1, 0], 0, 0, 
    Sqrt[6]*gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[6]*gep), 0, 0, 
    6*omega + U/2 + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], 0, 
    -(Sqrt[7]*gep), 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    Sqrt[2]*gammaPolCh[1], delta + 6*omega, Sqrt[2]*gammaPolCh[1], 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[6]*gep, 0, 
    Sqrt[2]*gammaPolCh[1], 2*delta + 6*omega + U/2 - coefzeta[1, 0], 0, 0, 
    Sqrt[7]*gep, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -(Sqrt[7]*gep), 0, 
    0, 7*omega + U/2 + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], 0, 
    -2*Sqrt[2]*gep, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    Sqrt[2]*gammaPolCh[1], delta + 7*omega, Sqrt[2]*gammaPolCh[1], 0, 0, 0, 
    0, 0, 0, 0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, Sqrt[7]*gep, 0, Sqrt[2]*gammaPolCh[1], 
    2*delta + 7*omega + U/2 - coefzeta[1, 0], 0, 0, 2*Sqrt[2]*gep, 0, 0, 0, 
    0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, -2*Sqrt[2]*gep, 0, 0, 8*omega + U/2 + coefzeta[1, 0], 
    Sqrt[2]*gammaPolCh[1], 0, -3*gep, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    Sqrt[2]*gammaPolCh[1], delta + 8*omega, Sqrt[2]*gammaPolCh[1], 0, 0, 0, 
    0, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 2*Sqrt[2]*gep, 0, Sqrt[2]*gammaPolCh[1], 
    2*delta + 8*omega + U/2 - coefzeta[1, 0], 0, 0, 3*gep, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    -3*gep, 0, 0, 9*omega + U/2 + coefzeta[1, 0], Sqrt[2]*gammaPolCh[1], 0, 
    -(Sqrt[10]*gep), 0, 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[2]*gammaPolCh[1], delta + 9*omega, 
    Sqrt[2]*gammaPolCh[1], 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 3*gep, 0, Sqrt[2]*gammaPolCh[1], 
    2*delta + 9*omega + U/2 - coefzeta[1, 0], 0, 0, Sqrt[10]*gep}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, -(Sqrt[10]*gep), 0, 0, 10*omega + U/2 + coefzeta[1, 0], 
    Sqrt[2]*gammaPolCh[1], 0}, {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[2]*gammaPolCh[1], 
    delta + 10*omega, Sqrt[2]*gammaPolCh[1]}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, Sqrt[10]*gep, 0, Sqrt[2]*gammaPolCh[1], 
    (4*delta + 20*omega + U - 2*coefzeta[1, 0])/2}}
dim={33, 33}
det[vec]=0.9999999999999998 1-abs=2.220446049250313*^-16
orthogonality check=33*Abs[-1 + 
      MyDot[{{0.49999999999999994, -0.7071067811865475, 0.5, 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0.4999999999999999, -0.7071067811865475, 
         0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0.5, -0.7071067811865477, 0.5, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 
         -0.7071067811865477, 0.4999999999999999, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999998, -0.7071067811865479, 0.4999999999999999, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999999, -0.7071067811865475, 0.5, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0.7071067811865475, 2.6888213877640505*^-16, -0.7071067811865476, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.5000000000000001, -0.7071067811865477, 
         0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0.7071067811865476, 1.5959455978986625*^-16, 
         -0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.5000000000000004, -0.7071067811865475, 
         0.4999999999999997, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0.7071067811865477, 
         4.2500725161431764*^-16, -0.7071067811865476, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 
         -0.7071067811865475, 0.5000000000000002, 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865477, 
         3.4867941867133813*^-16, -0.7071067811865476, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 
         -0.7071067811865476, 0.49999999999999983, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, -9.194034422677074*^-17, -0.7071067811865474, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999999, -0.7071067811865475, 0.5000000000000001}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865474, -1.734723475976807*^-16, -0.7071067811865477, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0.49999999999999994, 0.7071067811865477, 0.4999999999999998, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865474, 5.20417042793042*^-18, 
         -0.7071067811865479, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0.4999999999999999, 0.7071067811865476, 
         0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.7071067811865475, 3.6949610038306006*^-16, 
         -0.7071067811865479, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0.4999999999999997, 0.7071067811865476, 
         0.5, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865479, 
         2.740863092043355*^-16, -0.7071067811865471, 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.4999999999999997, 
         0.7071067811865476, 0.5, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         3.4694469519536157*^-16, -0.7071067811865478, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.5000000000000001, 0.7071067811865474, 0.5000000000000001, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865478, 3.625572064791526*^-16, -0.7071067811865471}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999999, 0.7071067811865476, 0.49999999999999983, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.5000000000000001, 0.7071067811865475, 
         0.4999999999999998, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0.5, 0.7071067811865477, 
         0.4999999999999998, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0.49999999999999956, 
         0.7071067811865476, 0.5000000000000004, 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.49999999999999994, 
         -0.7071067811865475, -0.5000000000000002, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.49999999999999944, 0.7071067811865476, 0.5000000000000004}}, 
       {{0.49999999999999994, 0., 0., 0., 0., 0., 0.7071067811865475, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0.49999999999999994, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {-0.7071067811865475, 0., 0., 0., 0., 0., 2.6888213877640505*^-16, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865477, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0.5, 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.4999999999999998, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., -0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         1.5959455978986625*^-16, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 
         -0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0.5, 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865477, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999997, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., -0.7071067811865477, 0., 0., 0., 0., 0., 0., 0., 
         4.2500725161431764*^-16, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0.5, 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.}, {0., 0., 0., 0.5000000000000001, 0., 0., 0., 
         0., 0., 0., 0., 0., 0.7071067811865477, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.4999999999999997, 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0., -0.7071067811865477, 0., 0., 0., 0., 0., 0., 
         0., 0., 3.4867941867133813*^-16, 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5, 
         0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.4999999999999998, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.5000000000000001, 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., -0.7071067811865479, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., -9.194034422677074*^-17, 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865474, 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., -0.7071067811865474, 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.5000000000000001, 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865474, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., -0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., -1.734723475976807*^-16, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0.5, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865477, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.49999999999999983, 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0.7071067811865474, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.5000000000000001, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., -0.7071067811865477, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 5.20417042793042*^-18, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., -0.7071067811865479, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.4999999999999998, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000004, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 
         0., 0., 0., 0., 0., 0.5, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865475, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 3.6949610038306006*^-16, 0., 0., 
         0., 0., 0., 0., 0., 0., 0.7071067811865477, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.4999999999999997, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865479, 0., 0., 0., 
         0., 0., 0., 0., 0., 0.4999999999999998, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865479, 0., 
         0., 0., 0., 0., 0., 0., 0.49999999999999956, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865475, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 2.740863092043355*^-16, 
         0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000002, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865471, 0., 
         0., 0., 0., 0., 0., 0., 0.5000000000000004, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.5000000000000001, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., -0.49999999999999994, 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         3.4694469519536157*^-16, 0., 0., 0., 0., 0., 0., 
         -0.7071067811865475, 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.49999999999999983, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865478, 0., 0., 0., 0., 0., 0., -0.5000000000000002, 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865478, 0., 0., 0., 0., 0., 0.49999999999999944}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         3.625572064791526*^-16, 0., 0., 0., 0., 0., 0.7071067811865476}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.5000000000000001, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865471, 0., 0., 0., 0., 0., 0.5000000000000004}}]] + 
   1056*Abs[MyDot[{{0.49999999999999994, -0.7071067811865475, 0.5, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0.4999999999999999, -0.7071067811865475, 
        0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0.5, -0.7071067811865477, 0.5, 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.5000000000000001, -0.7071067811865477, 0.4999999999999999, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.4999999999999998, -0.7071067811865479, 0.4999999999999999, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.4999999999999999, -0.7071067811865475, 0.5, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0.7071067811865475, 2.6888213877640505*^-16, -0.7071067811865476, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.5000000000000001, -0.7071067811865477, 0.4999999999999999, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0.7071067811865476, 1.5959455978986625*^-16, 
        -0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0.5000000000000004, -0.7071067811865475, 
        0.4999999999999997, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0.7071067811865477, 
        4.2500725161431764*^-16, -0.7071067811865476, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 
        -0.7071067811865475, 0.5000000000000002, 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865477, 
        3.4867941867133813*^-16, -0.7071067811865476, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 
        -0.7071067811865476, 0.49999999999999983, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
        -9.194034422677074*^-17, -0.7071067811865474, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.4999999999999999, -0.7071067811865475, 0.5000000000000001}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865474, -1.734723475976807*^-16, -0.7071067811865477, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0.49999999999999994, 0.7071067811865477, 0.4999999999999998, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865474, 5.20417042793042*^-18, -0.7071067811865479, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0.4999999999999999, 0.7071067811865476, 
        0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0.7071067811865475, 3.6949610038306006*^-16, 
        -0.7071067811865479, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0.4999999999999997, 0.7071067811865476, 
        0.5, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0.7071067811865479, 
        2.740863092043355*^-16, -0.7071067811865471, 0., 0., 0., 0., 0., 0.}\
, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.4999999999999997, 
        0.7071067811865476, 0.5, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        3.4694469519536157*^-16, -0.7071067811865478, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 
        0.7071067811865474, 0.5000000000000001, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865478, 3.625572064791526*^-16, -0.7071067811865471}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.4999999999999999, 0.7071067811865476, 0.49999999999999983, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.5000000000000001, 0.7071067811865475, 0.4999999999999998, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0.5, 0.7071067811865477, 0.4999999999999998, 0., 
        0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0.49999999999999956, 
        0.7071067811865476, 0.5000000000000004, 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.49999999999999994, 
        -0.7071067811865475, -0.5000000000000002, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.49999999999999944, 0.7071067811865476, 0.5000000000000004}}, 
      {{0.49999999999999994, 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0.49999999999999994, 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {-0.7071067811865475, 0., 0., 0., 0., 0., 2.6888213877640505*^-16, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865477, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0.5, 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0.4999999999999998, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.4999999999999999, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., -0.7071067811865475, 0., 0., 0., 0., 0., 0., 
        1.5959455978986625*^-16, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.}, {0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 
        -0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.}, {0., 0., 0.5, 0., 0., 0., 0., 0., 0., 0., 0.7071067811865477, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.4999999999999997, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., -0.7071067811865477, 0., 0., 0., 0., 0., 0., 0., 
        4.2500725161431764*^-16, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0.5, 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.}, {0., 0., 0., 0.5000000000000001, 0., 0., 0., 0., 
        0., 0., 0., 0., 0.7071067811865477, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0.4999999999999997, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., -0.7071067811865477, 0., 0., 0., 0., 0., 0., 0., 0., 
        3.4867941867133813*^-16, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5, 
        0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.4999999999999998, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.5000000000000001, 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., -0.7071067811865479, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., -9.194034422677074*^-17, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865474, 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., -0.7071067811865474, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.5000000000000001, 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.7071067811865474, 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.4999999999999999, 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., -0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., -1.734723475976807*^-16, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0.7071067811865476, 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0.5, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865477, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.49999999999999983, 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0.7071067811865474, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0.5000000000000001, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., -0.7071067811865477, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 5.20417042793042*^-18, 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.7071067811865475, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0.4999999999999999, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., -0.7071067811865479, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0.4999999999999998, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000004, 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 
        0., 0., 0., 0.5, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865475, 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 3.6949610038306006*^-16, 0., 0., 0., 0., 
        0., 0., 0., 0., 0.7071067811865477, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.4999999999999997, 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., -0.7071067811865479, 0., 0., 0., 0., 0., 
        0., 0., 0., 0.4999999999999998, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000001, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865479, 0., 0., 0., 
        0., 0., 0., 0., 0.49999999999999956, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865475, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 2.740863092043355*^-16, 0., 0., 
        0., 0., 0., 0., 0., 0.7071067811865476, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.5000000000000002, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865471, 0., 0., 0., 
        0., 0., 0., 0., 0.5000000000000004, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.5000000000000001, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., -0.49999999999999994, 0.}\
, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 3.4694469519536157*^-16, 0., 
        0., 0., 0., 0., 0., -0.7071067811865475, 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.49999999999999983, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865478, 0., 0., 0., 0., 0., 0., -0.5000000000000002, 0.}\
, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.4999999999999999, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865478, 0., 0., 0., 0., 0., 0.49999999999999944}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        3.625572064791526*^-16, 0., 0., 0., 0., 0., 0.7071067811865476}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.5000000000000001, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865471, 0., 0., 0., 0., 0., 0.5000000000000004}}]]
diagvc[{0, 3}]
Generating matrix and saving to ham.model..QS_0.3
hamil={{delta, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, delta + omega, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, delta + 2*omega, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, delta + 3*omega, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, delta + 4*omega, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, delta + 5*omega, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, delta + 6*omega, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, delta + 7*omega, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, delta + 8*omega, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, delta + 9*omega, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, delta + 10*omega}}
dim={11, 11}
det[vec]=1. 1-abs=0.
orthogonality check=11*Abs[-1 + 
      MyDot[{{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}, 
       {{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}]] + 
   110*Abs[MyDot[{{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}, 
      {{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}]]
diagvc[{1, 2}]
Generating matrix and saving to ham.model..QS_1.2
hamil={{delta + coefzeta[1, 0], -gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {-gammaPolCh[1], (4*delta + U)/2, 0, gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0}, {0, 0, delta + omega + coefzeta[1, 0], 
    -gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, gep, -gammaPolCh[1], 2*delta + omega + U/2, 0, Sqrt[2]*gep, 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, delta + 2*omega + coefzeta[1, 0], -gammaPolCh[1], 0, 0, 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, Sqrt[2]*gep, -gammaPolCh[1], (4*delta + 4*omega + U)/2, 0, 
    Sqrt[3]*gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, delta + 3*omega + coefzeta[1, 0], -gammaPolCh[1], 0, 
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, Sqrt[3]*gep, -gammaPolCh[1], (4*delta + 6*omega + U)/2, 
    0, 2*gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, delta + 4*omega + coefzeta[1, 0], 
    -gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 2*gep, -gammaPolCh[1], (4*delta + 8*omega + U)/2, 
    0, Sqrt[5]*gep, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, delta + 5*omega + coefzeta[1, 0], 
    -gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[5]*gep, -gammaPolCh[1], 
    2*delta + 5*omega + U/2, 0, Sqrt[6]*gep, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, delta + 6*omega + coefzeta[1, 0], 
    -gammaPolCh[1], 0, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[6]*gep, -gammaPolCh[1], 
    2*delta + 6*omega + U/2, 0, Sqrt[7]*gep, 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    delta + 7*omega + coefzeta[1, 0], -gammaPolCh[1], 0, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[7]*gep, -gammaPolCh[1], 
    2*delta + 7*omega + U/2, 0, 2*Sqrt[2]*gep, 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    delta + 8*omega + coefzeta[1, 0], -gammaPolCh[1], 0, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2*Sqrt[2]*gep, 
    -gammaPolCh[1], 2*delta + 8*omega + U/2, 0, 3*gep, 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    delta + 9*omega + coefzeta[1, 0], -gammaPolCh[1], 0, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3*gep, 
    -gammaPolCh[1], 2*delta + 9*omega + U/2, 0, Sqrt[10]*gep}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
    delta + 10*omega + coefzeta[1, 0], -gammaPolCh[1]}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[10]*gep, 
    -gammaPolCh[1], (4*delta + 20*omega + U)/2}}
dim={22, 22}
det[vec]=-0.9999999999999997 1-abs=3.3306690738754696*^-16
orthogonality check=22*Abs[-1 + 
      MyDot[{{0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.7071067811865476, 0.7071067811865475, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0.7071067811865476, 0.7071067811865475, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0.7071067811865475, -0.7071067811865476, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.7071067811865475, -0.7071067811865476, 0., 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0.7071067811865475, -0.7071067811865476, 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865476, 0.7071067811865475, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.7071067811865476, 0.7071067811865475}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865475, -0.7071067811865476, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865475, -0.7071067811865476}}, 
       {{0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0.7071067811865475, 0., 0., 0., 0., 0., -0.7071067811865476, 0., 
         0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0., 0.}, {0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 
         -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
         0., 0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 
         0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 
         0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 
         0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 
         0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
         0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         0., 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0.}\
, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
         0., 0., 0., 0., 0., 0., -0.7071067811865476, 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
         0.7071067811865475, 0., 0., 0., 0., 0., -0.7071067811865476}}]] + 
   462*Abs[MyDot[{{0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.7071067811865476, 0.7071067811865475, 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0.7071067811865476, 0.7071067811865475, 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.7071067811865475, -0.7071067811865476, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0.7071067811865475, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0.7071067811865475, -0.7071067811865476, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865476, 0.7071067811865475, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}\
, {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0.7071067811865476, 0.7071067811865475}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, -0.7071067811865476, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865475, -0.7071067811865476, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.7071067811865475, -0.7071067811865476}}, 
      {{0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0.7071067811865475, 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.7071067811865475, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., -0.7071067811865476, 
        0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 
        -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 0., 
        0., -0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 0., 0., 
        0., 0., 0., -0.7071067811865476, 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 0., 0., 
        0., 0., 0., 0., -0.7071067811865476, 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 0., 0., 
        0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0., 0., 
        0., 0., 0., 0., 0., -0.7071067811865476, 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865476, 0., 
        0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0., 
        0., 0., 0., 0., 0., 0., -0.7071067811865476, 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0., 0., 0., 0., 0., 0., 0.7071067811865475, 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., 0., -0.7071067811865476, 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865476, 0., 0., 0., 0., 0., 0.7071067811865475}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 
        0.7071067811865475, 0., 0., 0., 0., 0., -0.7071067811865476}}]]
diagvc[{2, 1}]
Generating matrix and saving to ham.model..QS_2.1
hamil={{2*delta + U/2 + coefzeta[1, 0], gep, 0, 0, 0, 0, 0, 0, 0, 0, 0}, 
   {gep, 2*delta + omega + U/2 + coefzeta[1, 0], Sqrt[2]*gep, 0, 0, 0, 0, 0, 
    0, 0, 0}, {0, Sqrt[2]*gep, 2*delta + 2*omega + U/2 + coefzeta[1, 0], 
    Sqrt[3]*gep, 0, 0, 0, 0, 0, 0, 0}, 
   {0, 0, Sqrt[3]*gep, 2*delta + 3*omega + U/2 + coefzeta[1, 0], 2*gep, 0, 
    0, 0, 0, 0, 0}, {0, 0, 0, 2*gep, 
    2*delta + 4*omega + U/2 + coefzeta[1, 0], Sqrt[5]*gep, 0, 0, 0, 0, 0}, 
   {0, 0, 0, 0, Sqrt[5]*gep, 2*delta + 5*omega + U/2 + coefzeta[1, 0], 
    Sqrt[6]*gep, 0, 0, 0, 0}, {0, 0, 0, 0, 0, Sqrt[6]*gep, 
    2*delta + 6*omega + U/2 + coefzeta[1, 0], Sqrt[7]*gep, 0, 0, 0}, 
   {0, 0, 0, 0, 0, 0, Sqrt[7]*gep, 2*delta + 7*omega + U/2 + coefzeta[1, 0], 
    2*Sqrt[2]*gep, 0, 0}, {0, 0, 0, 0, 0, 0, 0, 2*Sqrt[2]*gep, 
    2*delta + 8*omega + U/2 + coefzeta[1, 0], 3*gep, 0}, 
   {0, 0, 0, 0, 0, 0, 0, 0, 3*gep, 2*delta + 9*omega + U/2 + coefzeta[1, 0], 
    Sqrt[10]*gep}, {0, 0, 0, 0, 0, 0, 0, 0, 0, Sqrt[10]*gep, 
    2*delta + 10*omega + U/2 + coefzeta[1, 0]}}
dim={11, 11}
det[vec]=1. 1-abs=0.
orthogonality check=11*Abs[-1 + 
      MyDot[{{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}, 
       {{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
        {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}]] + 
   110*Abs[MyDot[{{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}, 
      {{1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 1., 0., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 1., 0., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 1., 0., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 1., 0., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 1., 0., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0.}, 
       {0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1.}}]]
Lowest energies (absolute):{-0.5046265044040315, -0.4046265044040322, 
   -0.3046265044040317, -0.25231325220201606, -0.25231325220201606, 
   -0.20462650440403196, -0.15231325220201605, -0.15231325220201605, 
   -0.10462650440403215, -0.05231325220201605, -0.05231325220201605, 
   -0.004626504404032117, 0., 0., 0., 4.788107844239109*^-17, 
   0.047686747797983986, 0.047686747797983986, 0.09537349559596797, 0.1}
Lowest energies (GS shifted):{0., 0.09999999999999926, 0.19999999999999973, 
   0.2523132522020154, 0.2523132522020154, 0.2999999999999995, 
   0.35231325220201537, 0.35231325220201537, 0.3999999999999993, 
   0.4523132522020154, 0.4523132522020154, 0.49999999999999933, 
   0.5046265044040315, 0.5046265044040315, 0.5046265044040315, 
   0.5046265044040315, 0.5523132522020154, 0.5523132522020154, 
   0.5999999999999994, 0.6046265044040314}
Scale factor SCALE(Ninit):1.0201394465967895
Lowest energies (shifted and scaled):
  {0., 0.09802581434685399, 0.19605162869370918, 0.24733212017605893, 
   0.24733212017605893, 0.29407744304056366, 0.34535793452291363, 
   0.34535793452291363, 0.3921032573874182, 0.4433837488697684, 
   0.4433837488697684, 0.49012907173427295, 0.49466424035211853, 
   0.49466424035211853, 0.49466424035211853, 0.49466424035211853, 
   0.5414095632166231, 0.5414095632166231, 0.5881548860811278, 
   0.5926900546989733}
ireducTable: f[0]
Loading module operators.m
"operators.m $Id: operators.m 877 2009-12-16 16:34:11Z zitko $"
n_d op.model..QS.n_d nc[d[0, 0], d[1, 0]] + nc[d[0, 1], d[1, 1]]
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]]
ireducTable: d[#1, #2] & 
acr op.model..QS.acr nc[ket[1], bra[0]] + Sqrt[2]*nc[ket[2], bra[1]] + 
   Sqrt[3]*nc[ket[3], bra[2]] + 2*nc[ket[4], bra[3]] + 
   Sqrt[5]*nc[ket[5], bra[4]] + Sqrt[6]*nc[ket[6], bra[5]] + 
   Sqrt[7]*nc[ket[7], bra[6]] + 2*Sqrt[2]*nc[ket[8], bra[7]] + 
   3*nc[ket[9], bra[8]] + Sqrt[10]*nc[ket[10], bra[9]]
nph op.model..QS.nph nc[ket[1], bra[1]] + 2*nc[ket[2], bra[2]] + 
   3*nc[ket[3], bra[3]] + 4*nc[ket[4], bra[4]] + 5*nc[ket[5], bra[5]] + 
   6*nc[ket[6], bra[6]] + 7*nc[ket[7], bra[7]] + 8*nc[ket[8], bra[8]] + 
   9*nc[ket[9], bra[9]] + 10*nc[ket[10], bra[10]]
displ op.model..QS.displ nc[ket[0], bra[1]] + nc[ket[1], bra[0]] + 
   Sqrt[2]*nc[ket[1], bra[2]] + Sqrt[2]*nc[ket[2], bra[1]] + 
   Sqrt[3]*nc[ket[2], bra[3]] + Sqrt[3]*nc[ket[3], bra[2]] + 
   2*nc[ket[3], bra[4]] + 2*nc[ket[4], bra[3]] + Sqrt[5]*nc[ket[4], bra[5]] + 
   Sqrt[5]*nc[ket[5], bra[4]] + Sqrt[6]*nc[ket[5], bra[6]] + 
   Sqrt[6]*nc[ket[6], bra[5]] + Sqrt[7]*nc[ket[6], bra[7]] + 
   Sqrt[7]*nc[ket[7], bra[6]] + 2*Sqrt[2]*nc[ket[7], bra[8]] + 
   2*Sqrt[2]*nc[ket[8], bra[7]] + 3*nc[ket[8], bra[9]] + 
   3*nc[ket[9], bra[8]] + Sqrt[10]*nc[ket[9], bra[10]] + 
   Sqrt[10]*nc[ket[10], bra[9]]
displ^2 op.model..QS.displ^2 nc[ket[0], bra[0]] + 
   Sqrt[2]*nc[ket[0], bra[2]] + 3*nc[ket[1], bra[1]] + 
   Sqrt[6]*nc[ket[1], bra[3]] + Sqrt[2]*nc[ket[2], bra[0]] + 
   5*nc[ket[2], bra[2]] + 2*Sqrt[3]*nc[ket[2], bra[4]] + 
   Sqrt[6]*nc[ket[3], bra[1]] + 7*nc[ket[3], bra[3]] + 
   2*Sqrt[5]*nc[ket[3], bra[5]] + 2*Sqrt[3]*nc[ket[4], bra[2]] + 
   9*nc[ket[4], bra[4]] + Sqrt[30]*nc[ket[4], bra[6]] + 
   2*Sqrt[5]*nc[ket[5], bra[3]] + 11*nc[ket[5], bra[5]] + 
   Sqrt[42]*nc[ket[5], bra[7]] + Sqrt[30]*nc[ket[6], bra[4]] + 
   13*nc[ket[6], bra[6]] + 2*Sqrt[14]*nc[ket[6], bra[8]] + 
   Sqrt[42]*nc[ket[7], bra[5]] + 15*nc[ket[7], bra[7]] + 
   6*Sqrt[2]*nc[ket[7], bra[9]] + 2*Sqrt[14]*nc[ket[8], bra[6]] + 
   17*nc[ket[8], bra[8]] + 3*Sqrt[10]*nc[ket[8], bra[10]] + 
   6*Sqrt[2]*nc[ket[9], bra[7]] + 19*nc[ket[9], bra[9]] + 
   3*Sqrt[10]*nc[ket[10], bra[8]] + 10*nc[ket[10], bra[10]]
$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.077145`5.338852776872744}
{ham, 3.068527`6.160323692110611}
{maketable, 6.970705`7.294821697293583}
{xi, 1.000267`6.451660934645178}
{_, 0}
data
gammaPol=0.252313252202016

"Success!"


