SpinSiteXxSU2.h

Go to the documentation of this file.

#ifndef SPINSITEXxSU2_H_
#define SPINSITEXxSU2_H_
 
#include "symmetry/kind_dummies.h"
 
#include "DmrgTypedefs.h"
 
#include "symmetry/S1xS2.h"
#include "symmetry/SU2.h"
#include "symmetry/U1.h"
 
#include "tensors/SiteOperatorQ.h"
 
template <typename Symmetry, size_t order> class SpinSite;
 
template <typename OtherSym>
class SpinSite<Sym::S1xS2<OtherSym, Sym::SU2<Sym::SpinSU2> >, 1ul>
{
    typedef double Scalar;
    typedef typename Sym::S1xS2<OtherSym, Sym::SU2<Sym::SpinSU2> > Symmetry;
    typedef SiteOperatorQ<Symmetry,Eigen::Matrix<Scalar,Eigen::Dynamic,Eigen::Dynamic> > OperatorType;
public:
    SpinSite() {};
    SpinSite(std::size_t D_input);
    
    OperatorType Id_1s() const {return Id_1s_;}
    
    OperatorType S_1s() const {return S_1s_;}
 
    OperatorType Q_1s() const {return Q_1s_;}
 
    Qbasis<Symmetry> basis_1s() const {return basis_1s_;}
protected:
    std::size_t D;
    
    Qbasis<Symmetry> basis_1s_;
 
    OperatorType Id_1s_; //identity
    OperatorType S_1s_; //spin
    OperatorType Sdag_1s_; //spin
    OperatorType Q_1s_; //quadrupole moment
};
 
template<typename OtherSym>
SpinSite<Sym::S1xS2<OtherSym, Sym::SU2<Sym::SpinSU2> >, 1ul>::
SpinSite(std::size_t D_input)
:D(D_input)
{
    //create basis for one Spin Site
    typename Symmetry::qType Q = join(OtherSym::qvacuum(),qarray<1>{static_cast<int>(D)});
    Eigen::Index inner_dim = 1;
    std::vector<std::string> ident;
    ident.push_back("spin");
    
    basis_1s_.push_back(Q,inner_dim,ident);
    
    Id_1s_ = OperatorType(Symmetry::qvacuum(),basis_1s_,"id");
    S_1s_ = OperatorType(join(OtherSym::qvacuum(),qarray<1>{3}),basis_1s_,"S");
 
    Scalar locS = 0.5*static_cast<double>(D-1);
    S_1s_( "spin", "spin" ) = std::sqrt(locS*(locS+1.));
    Sdag_1s_ = S_1s_.adjoint();
    Q_1s_ = OperatorType::prod(S_1s_,S_1s_,join(OtherSym::qvacuum(),qarray<1>{5}));
    Id_1s_.setIdentity();
}
 
#endif //SPINSITESU2xX_H_