FermionSiteU1xU1.h

Go to the documentation of this file.

#ifndef FERMIONSITEU1xU1_H_
#define FERMIONSITEU1xU1_H_
 
#include "DmrgTypedefs.h"
 
#include "symmetry/S1xS2.h"
#include "symmetry/SU2.h"
#include "symmetry/U1.h"
 
#include "sites/U1Helper.h"
 
#include "tensors/SiteOperatorQ.h"
 
template <typename Symmetry> class FermionSite;
 
// template <>
// class FermionSite<Sym::S1xS2<Sym::U1<Sym::SpinU1>,Sym::U1<Sym::ChargeU1> > >
// {
//  typedef double Scalar;
//  typedef Sym::S1xS2<Sym::U1<Sym::SpinU1>,Sym::U1<Sym::ChargeU1> > Symmetry;
//  typedef SiteOperatorQ<Symmetry,Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic> > OperatorType;
// public:
//  FermionSite() {};
//  FermionSite(bool U_IS_INFINITE, bool UPH_IS_INFINITE);
    
//  OperatorType Id_1s() const {return Id_1s_;}
//  OperatorType F_1s() const {return F_1s_;}
    
//  OperatorType c_1s(SPIN_INDEX sigma) const { if (sigma == UP) {return cup_1s_;} return cdn_1s_;}
//  OperatorType cdag_1s(SPIN_INDEX sigma) const { if (sigma == UP) {return cup_1s_.adjoint();} return cdn_1s_.adjoint();}
 
//  OperatorType n_1s() const {return n_1s(UP) + n_1s(DN);}
//  OperatorType n_1s(SPIN_INDEX sigma) const { if (sigma == UP) {return nup_1s_;} return ndn_1s_;}
//  OperatorType ns_1s() const {return n_1s() - 2.*d_1s();}
//  OperatorType nh_1s() const {return 2.*d_1s() - n_1s() + Id_1s();}
//  OperatorType d_1s() const {return d_1s_;}
 
//  OperatorType Sz_1s() const {return Sz_1s_;}
//  OperatorType Sp_1s() const {return Sp_1s_;}
//  OperatorType Sm_1s() const {return Sm_1s_;}
 
//  OperatorType Tz_1s() const {return 0.5*(n_1s() - Id_1s());}
//  OperatorType cc_1s() const {return cc_1s_;}
//  OperatorType cdagcdag_1s() const {return cdagcdag_1s_;}
 
//  Qbasis<Symmetry> basis_1s() const {return basis_1s_;}
 
//  template<typename Symmetry> friend void fill_SiteOps (FermionSite<Symmetry> &Site, bool U_IS_INFINITE, bool UPH_IS_INFINITE);
// protected:
    
//  Qbasis<Symmetry> basis_1s_;
 
//  OperatorType Id_1s_; //identity
//  OperatorType F_1s_; //Fermionic sign
//  OperatorType cup_1s_; //annihilation
//  OperatorType cdn_1s_; //annihilation
    
//  OperatorType n_1s_; //particle number
//  OperatorType nup_1s_; //particle number
//  OperatorType ndn_1s_; //particle number
//  OperatorType d_1s_; //double occupancy
 
//  OperatorType Sz_1s_; //orbital spin
//  OperatorType Sp_1s_; //orbital spin
//  OperatorType Sm_1s_; //orbital spin
    
//  OperatorType Tz_1s_; //orbital pseude spin
//  OperatorType cc_1s_; //pairing
//  OperatorType cdagcdag_1s_; //pairing adjoint
// };
 
// FermionSite<Sym::S1xS2<Sym::U1<Sym::SpinU1>,Sym::U1<Sym::ChargeU1> > >::
// FermionSite(bool U_IS_INFINITE, bool UPH_IS_INFINITE)
// {
//  //create basis for one Fermionic Site
//  typename Symmetry::qType Q={0,0}; //empty occupied state
//  Eigen::Index inner_dim = 1;
//  std::vector<std::string> ident;
    
//  if (!UPH_IS_INFINITE)
//  {
//      ident.push_back("empty");
//      this->basis_1s_.push_back(Q,inner_dim,ident);
//      ident.clear();
//  }
    
//  Q={+1,1}; //up spin state
//  inner_dim = 1;
//  ident.push_back("up");
//  this->basis_1s_.push_back(Q,inner_dim,ident);
//  ident.clear();
 
//  Q={-1,1}; //down spin state
//  inner_dim = 1;
//  ident.push_back("dn");
//  this->basis_1s_.push_back(Q,inner_dim,ident);
//  ident.clear();
    
//  if (!U_IS_INFINITE and !UPH_IS_INFINITE)
//  {
//      Q={0,2}; //doubly occupied state
//      inner_dim = 1;
//      ident.push_back("double");
//      this->basis_1s_.push_back(Q,inner_dim,ident);
//      ident.clear();
//  }
 
//  cout << "single site basis" << endl << this->basis_1s_ << endl;
 
//  ::fill_SiteOps(*this, U_IS_INFINITE, UPH_IS_INFINITE);
// }
 
#endif //FERMIONSITESU2xU1_H_