bugfix

Prog/Hamiltonian_Examples.f90

@@ -377,6 +377,7 @@
              Do nf = 1,N_FL
                 nc = 0
                 X = 1.d0
+                if (nf == 2) X = -1.d0 ! This induces the asymmetry to the other spin sector
                 
                 Do i = 1,Ndim
                    nc = nc + 1
@@ -390,10 +391,9 @@
                    Call Op_set( Op_V(nc,nf,1) )
                    
                    ! Here comes the S_z rep part (non-default)
-                   if (nf == 2) X = -1.d0 ! This induces the asymmetry to the other spin sector
                    Op_V(nc,nf,2)%P(1) = I
                    Op_V(nc,nf,2)%O(1,1) = cmplx(1.d0, 0.d0, kind(0.D0))
-                   Op_V(nc,nf,2)%g      = X*SQRT(CMPLX(DTAU*ham_U/2.d0, 0.D0, kind(0.D0))) ! Note the sign
+                   Op_V(nc,nf,2)%g      = X*SQRT(CMPLX(-DTAU*ham_U/2.d0, 0.D0, kind(0.D0))) ! changed the sign...
                    Op_V(nc,nf,2)%alpha  = cmplx(0.d0, 0.d0, kind(0.D0)) ! Note the alpha
                    Op_V(nc,nf,2)%type   = 2
                    Call Op_set( Op_V(nc,nf,2) )