Permalink
Browse files

Dont now what the hell is happening with majorana case

  • Loading branch information...
cifu9502 committed Nov 13, 2017
1 parent 497bb56 commit c6b4944615e5abaf2ef19708be7f3609c0df0d76
@@ -877,8 +877,6 @@ complex<double> OneChNupPdn_H0DQD_MatEl(vector<double> Params,
//state came from
int stcfj=pAbasis->StCameFrom[jst];
// Get c1_up matrix element from MatArray[0]
// Get c1_dn matrix element from MatArray[1]
// Get c2_up matrix element from MatArray[2]
@@ -902,10 +900,12 @@ complex<double> OneChNupPdn_H0DQD_MatEl(vector<double> Params,
// if zero, try h.c.
int help = 0.0;
if ((dEqual(fabs(OldEl[icounter].real()),0.0))&&(dEqual(fabs(OldEl[icounter].imag()),0.0))){
OldEl[icounter]=MatArray[icounter].cGetMatEl(stcfj,stcfi);
typep=typej;
type=typei;
help=1.0;
}
// I NEED TO CHANGE THIS RULE NOW Get SingleSite QNumbers
@@ -922,29 +922,45 @@ complex<double> OneChNupPdn_H0DQD_MatEl(vector<double> Params,
)
FermiSign=-1.0;
// I HAVE TO CHANGE THIS TABLE?
double FullMatEl=OneCh_fd_table(sigma,typep,type)*FermiSign;
double FullMatEl= OneCh_fd_table(sigma,typep,type)*FermiSign;
cMatEl+=chi_N[idot-1]*OldEl[icounter]*FullMatEl;
//MatEl+=chi_N[idot-1]*FullMatEl;
//if (( (ist==8)||(ist==10) )&&( (jst==8)||(jst==10) ))
if (dEqual(Nupi,0.0)&&dEqual(Pdni,-1.0))
cout << " Oiiiiiiiii"
<< " ist= " << ist
if (dEqual(Nupi,0.0)&&dEqual(Pdni,-1.0)){
cout << " ist= " << ist
<< " jst = " << jst
<< " idot = " << idot
<< " FermiSign = " << FermiSign
<< " Chi = " << chi_N[idot-1]
<< " sigma = " << sigma
//<< " Szold = " << Szold
<< " OldEl = " << OldEl[icounter]
<< " FullMatEl = " << FullMatEl
<< " MatEl = " << cMatEl
<< "chi_N[idot-1]*OldEl[icounter]" << chi_N[idot-1]*OldEl[icounter]
<< " Prefactor" << chi_N[idot-1]*OldEl[icounter]
<< " Help" << help
<< endl
<< endl;
}
if (dEqual(Nupi,0.0)&&dEqual(Pdni,1.0)){
cout << " ist= " << ist
<< " jst = " << jst
<< " idot = " << idot
<< " Chi = " << chi_N[idot-1]
<< " sigma = " << sigma
//<< " Szold = " << Szold
<< " OldEl = " << OldEl[icounter]
<< " FullMatEl = " << FullMatEl
<< " MatEl = " << cMatEl
<< " Prefactor" << chi_N[idot-1]*OldEl[icounter]
<< " Help" << help
<< endl
<< endl;
}
icounter++;
}
// end loop in sigma
View
BIN +1.87 KB (100%) src/Main/Main_OpMatRules.o
Binary file not shown.
View
@@ -2681,9 +2681,10 @@ void OneChNupPdn_SetH0_AndersonMajorana(vector<double> Params,
//NEW CODE, I START BY ADDING SOME MISSING PARAMETERS, THE PARAMETERS OF THE SECOND QD WILL BE INITIALIZED EQUAL TO THE FIRST ONE
//double tdots = t1+t2;
double t11 = t1;
double t12 = 0;
//double t11 = t1;
//double t12 = 0;
double t11 = 0;
double t12 = t1;
// double U2 = Utilde;
//double ed2 = edtilde;
@@ -3307,7 +3308,7 @@ void OneChNupPdn_SetH0_AndersonMajorana(vector<double> Params,
CommonQNs.push_back(1); // pos of Parity QN
//exit(0);
// totSpos.push_back(1); // SU(2) symmetry in position 1
// totSpos.push_back(1); // SU(2) symmetry in position
// AbasisHm1 becomes "AcutHm1"
AbasisHm1.FalseCut(&AeigHm1);
View
BIN +0 Bytes (100%) src/Main/Main_SetH0.o
Binary file not shown.
View
BIN +4 KB (100%) src/Main/NRG_main
Binary file not shown.
@@ -0,0 +1,58 @@
2
1000
5.E-1
2.82691E-2
-0.25
2.5
0
0
0
0
9999999
0.0
0.02
0.0
0.0
0.0
0.0
0.0
0.0
0.5
2.82691E-2
-2.5E-1
Nsitesmax+1 (watch out! this is Nsitesmax+1 no Nsitesmax...)
Ncutoff (Cutoff states)
U_d
Gamma1 (J1)
e_d
Lambda
Dband
tela (Display messages on screen? y - 1; n - 0)
setGSzero/UpdateBefCut (in old/new code y - 1; n - 0)
calc_density (calculate the spectral density? y - 1; n - 0; 2- Calc Suscep; 3- calc chain suscep)
9999999-------- cont for DQD+majorana
h
t1
t2
U2
Gamma2 (J2) lambda_ph
e2 alpha
tdots
Sample parameters for 1ch calc:
U/D=0.5 ; U/(Pi Gamma) = 5.63
5.E-1
2.82691E-2
-2.5E-1
2.5
1.0
@@ -1,8 +1,8 @@
2
1000
0
0
0
5.E-1
2.82691E-2
-0.25
2.5
0
0
@@ -12,9 +12,9 @@
0.0
0.0
0.02
5.E-1
2.82691E-2
-0.25
0.0
0.0
0.0
0.0
0.0
@@ -20,12 +20,11 @@
0.5
2.82691E-2
-2.5E-1
Nsitesmax+1 (watch out! this is Nsitesmax+1 no Nsitesmax...)
Ncutoff (Cutoff states)
U_d
Oops, something went wrong.

0 comments on commit c6b4944

Please sign in to comment.