Detailed Description

Heisenberg Model.

MPO representation of

$H = J \sum_{<ij>} \left(\mathbf{S_i}\mathbf{S_j}\right) +J' \sum_{<<ij>>} \left(\mathbf{S_i}\mathbf{S_j}\right)$

Note: Makes use of the spin-SU(2) symmetry, which implies no magnetic fields. For B-fields see VMPS::HeisenbergU1.; The default variable settings can be seen in HeisenbergSU2::defaults.; is antiferromagnetic

Definition at line 29 of file HeisenbergSU2.h.

#include <HeisenbergSU2.h>

Inheritance diagram for VMPS::HeisenbergSU2:

[legend]

Public Types
typedef Sym::SU2< Sym::SpinSU2 >	Symmetry

typedef Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic >	MatrixType

typedef SiteOperatorQ< Symmetry, MatrixType >	OperatorType

Public Types inherited from Mpo< Sym::SU2< Sym::SpinSU2 >, double >
typedef SparseMatrixXd	SparseMatrixType

typedef SiteOperator< Sym::SU2< Sym::SpinSU2 >, double >	OperatorType

typedef Symmetry::qType	qType

typedef Matrix< double, Dynamic, Dynamic >	MatrixType

typedef double	Scalar_

typedef Mps< Sym::SU2< Sym::SpinSU2 >, double >	StateXd

typedef Umps< Sym::SU2< Sym::SpinSU2 >, double >	StateUd

typedef Mps< Sym::SU2< Sym::SpinSU2 >, std::complex< double > >	StateXcd

typedef Umps< Sym::SU2< Sym::SpinSU2 >, std::complex< double > >	StateUcd

typedef MpsCompressor< Sym::SU2< Sym::SpinSU2 >, double, double >	CompressorXd

typedef MpsCompressor< Sym::SU2< Sym::SpinSU2 >, std::complex< double >, double >	CompressorXcd

typedef Mpo< Sym::SU2< Sym::SpinSU2 > >	Operator

Public Types inherited from MpoTerms< Symmetry, Scalar >
typedef SiteOperator< Symmetry, Scalar >	OperatorType

typedef Symmetry::qType	qType

typedef Eigen::SparseMatrix< Scalar, Eigen::ColMajor, EIGEN_DEFAULT_SPARSE_INDEX_TYPE >	MatrixType

Public Member Functions
bool	validate (qarray< 1 > qnum) const


	HeisenbergSU2 ()

	HeisenbergSU2 (const size_t &L, const vector< Param > &params={}, const BC &boundary=BC::OPEN, const DMRG::VERBOSITY::OPTION &VERB=DMRG::VERBOSITY::ON_EXIT)

	HeisenbergSU2 (Mpo< Symmetry > &Mpo_input, const vector< Param > &params)

Public Member Functions inherited from Mpo< Sym::SU2< Sym::SpinSU2 >, double >
	Mpo ()

	Mpo (MpoTerms< Sym::SU2< Sym::SpinSU2 >, double > &Terms)

	Mpo (size_t L_input)

	Mpo (std::size_t L_input, qType Qtot_input, std::string label_input="Mpo", bool HERMITIAN_input=false, bool UNITARY_input=false, BC BC_input=BC::OPEN, DMRG::VERBOSITY::OPTION VERB_input=DMRG::VERBOSITY::OPTION::SILENT)

void	construct_from_pushlist (const PushType< OperatorType, CouplScalar > &pushlist, const std::vector< std::vector< std::string > > &labellist, size_t Lcell)

void	calc_reversedData_from_pushlist (const PushType< OperatorType, CouplScalar > &pushlist, double tolerance=::mynumeric_limits< double >::epsilon())

void	push_qpath (const std::size_t loc, const std::vector< OperatorType > &opList, const std::vector< qType > &qList, const double lambda=1.0)

void	setLocal (std::size_t loc, const OperatorType &op)

void	setLocal (std::size_t loc, const OperatorType &op, const OperatorType &signOp)

void	setLocal (std::size_t loc, const OperatorType &op, const std::vector< OperatorType > &signOp)

void	setLocal (const std::vector< std::size_t > &locs, const std::vector< OperatorType > &ops)

void	setLocal (const std::vector< std::size_t > &locs, const std::vector< OperatorType > &ops, const OperatorType &signOp)

void	setLocal (const std::vector< std::size_t > &locs, const std::vector< OperatorType > &ops, const std::vector< OperatorType > &signOps)

void	setLocalStag (std::size_t loc, const OperatorType &op, const std::vector< OperatorType > &stagSignOps)

void	setLocalSum (const OperatorType &op, double(*f)(int)=localSumTrivial)

void	setLocalSum (const std::vector< OperatorType > &op, std::vector< double > coeffs)

void	setProductSum (const OperatorType &op1, const OperatorType &op2)

void	precalc_TwoSiteData (bool FORCE=false)

std::string	info (bool REDUCED=false) const

int	get_dAux_max (int power=1) const

std::size_t	length () const

std::size_t	volume () const

int	locality () const

void	set_locality (std::size_t LocalSite_input)

OperatorType	localOperator () const

void	set_localOperator (OperatorType LocalOp_input)

Mpo< Sym::SU2< Sym::SpinSU2 >, double >	Identity () const

Mpo< Sym::SU2< Sym::SpinSU2 >, double >	Zero () const

bool	IS_UNITARY () const

bool	IS_HERMITIAN () const

bool	IS_HAMILTONIAN () const

bool	HAS_TWO_SITE_DATA () const

bool	HAS_W () const

boost::multi_array< double, 4 >	H2site (std::size_t loc, bool HALF_THE_LOCAL_TERM=false) const

void	generate_label (std::size_t Lcell)

void	push_width (const std::size_t n, const std::size_t loc, const double lambda, const OperatorType &outOp, const std::vector< OperatorType > &trans, const OperatorType &inOp)

void	push_local (const std::size_t loc, const double lambda, const OperatorType &op)

void	push_tight (const std::size_t loc, const double lambda, const OperatorType &op1, const OperatorType &op2)

void	push_nextn (const std::size_t loc, const double lambda, const OperatorType &op1, const OperatorType &trans, const OperatorType &op2)

Public Member Functions inherited from MpoTerms< Symmetry, Scalar >
void	reconstruct (const std::vector< std::map< std::array< qType, 2 >, std::vector< std::vector< std::map< qType, OperatorType > > > > > &O_in, const std::vector< Qbasis< Symmetry > > &qAux_in, const std::vector< std::vector< qType > > &qPhys_in, const bool FINALIZED_IN, const BC boundary_condition_in, const qType &qTot_in=Symmetry::qvacuum())

	MpoTerms (const std::size_t L=0, const BC boundary_condition_in=BC::OPEN, const qType &qTot_in=Symmetry::qvacuum(), const DMRG::VERBOSITY::OPTION &VERB_in=DMRG::VERBOSITY::OPTION::SILENT)

void	initialize (const std::size_t L, const BC boundary_condition_in, const qType &qTot_in)

virtual void	push (const std::size_t loc, const std::vector< OperatorType > &opList, const std::vector< qType > &qList, const Scalar lambda=1.0)

void	push (const std::size_t loc, const std::vector< OperatorType > &opList, const Scalar lambda=1.0)

void	show () const

void	save_label (const std::size_t loc, const std::string &info_label)

void	set_name (const std::string &label_in)

const std::string	get_name () const

std::vector< std::string >	get_info () const

void	scale (const Scalar factor, const Scalar offset=0., const std::size_t power=0ul, const double tolerance=1.e-14)

template<typename OtherScalar >
MpoTerms< Symmetry, OtherScalar >	cast ()

void	set_qPhys (const std::size_t loc, const std::vector< qType > &qPhys_in)

void	finalize (const bool COMPRESS=true, const std::size_t power=1, const double tolerance=::mynumeric_limits< double >::epsilon())

void	calc (const std::size_t power)

bool	is_finalized () const

const std::vector< std::map< std::array< qType, 2 >, std::vector< std::vector< std::map< qType, OperatorType > > > > > &	get_O () const

const std::vector< std::vector< std::vector< std::vector< Biped< Symmetry, MatrixType > > > > > &	get_W () const

const std::vector< Qbasis< Symmetry > > &	get_qAux () const

const std::vector< std::vector< qType > > &	get_qOp () const

const std::vector< std::vector< qType > > &	get_qPhys () const

const std::vector< std::vector< std::vector< std::vector< Biped< Symmetry, MatrixType > > > > > &	get_W_power (std::size_t power) const

const std::vector< Qbasis< Symmetry > > &	get_qAux_power (std::size_t power) const

const std::vector< std::vector< qType > > &	get_qOp_power (std::size_t power) const

std::size_t	get_hilbert_dimension (const std::size_t loc) const

BC	get_boundary_condition () const

const qType &	get_qVac () const

const qType &	get_qTot () const

const std::size_t	get_pos_qTot () const

bool	check_qPhys () const

bool	check_power (std::size_t power) const

std::size_t	maxPower () const

std::size_t	size () const

void	transform_base (const qType &qShift, const bool PRINT=false, const int factor=-1, const std::size_t powre=0ul)

std::vector< std::vector< TwoSiteData< Symmetry, Scalar > > >	calc_TwoSiteData () const

std::vector< std::pair< qType, std::size_t > >	base_order_IBC (const std::size_t power=1) const

double	memory (MEMUNIT memunit=kB) const

double	sparsity (const std::size_t power=1, const bool PER_MATRIX=false) const

void	set_Identity (const typename Symmetry::qType &Q=Symmetry::qvacuum())

void	set_Zero ()

void	set_verbosity (const DMRG::VERBOSITY::OPTION VERB_in)

DMRG::VERBOSITY::OPTION	get_verbosity () const

void	save (std::string filename)

void	load (std::string filename)

const std::vector< std::vector< std::vector< Biped< Symmetry, MatrixType > > > > &	W_at (const std::size_t loc) const

const std::vector< std::vector< std::vector< std::vector< Biped< Symmetry, MatrixType > > > > > &	W_full () const

const std::vector< std::vector< qType > > &	locBasis () const

const std::vector< qType > &	locBasis (const std::size_t loc) const

const std::vector< Qbasis< Symmetry > > &	auxBasis () const

const Qbasis< Symmetry > &	auxBasis (const std::size_t loc) const

const Qbasis< Symmetry > &	inBasis (const std::size_t loc) const

const Qbasis< Symmetry > &	outBasis (const std::size_t loc) const

const std::vector< std::vector< qType > > &	opBasis () const

const std::vector< qType > &	opBasis (const std::size_t loc) const

const qType &	Qtarget () const

void	setLocBasis (const std::vector< std::vector< qType > > &q)

void	setLocBasis (const std::vector< qType > &q, std::size_t loc)

const std::vector< std::vector< std::vector< Biped< Symmetry, MatrixType > > > > &	Wsq_at (const std::size_t loc) const

const std::vector< qType > &	opBasisSq (const std::size_t loc) const

const std::vector< std::vector< qType > > &	opBasisSq () const

const bool	check_SQUARE () const

double	sparsity (bool USE_SQUARE, bool PER_MATRIX) const

void	setQtarget (const qType &q)

Public Member Functions inherited from HeisenbergObservables< Sym::SU2< Sym::SpinSU2 > >
std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, complex< double > > >::type	S_ky (vector< complex< double > > phases) const

std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, complex< double > > >::type	Sdag_ky (vector< complex< double > > phases, double factor=sqrt(3.)) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, complex< double > > >::type	Sa_ky (SPINOP_LABEL Sa, vector< complex< double > > phases) const

std::vector< Mpo< Sym::SU2< Sym::SpinSU2 >, double > >	make_spinPermutation (const Permutation &permutations) const

MpoTerms< Sym::SU2< Sym::SpinSU2 >, double >	spin_swap_operator_D2 (const std::size_t locx1, const std::size_t locx2, const std::size_t locy1=0ul, const std::size_t locy2=0ul) const

MpoTerms< Sym::SU2< Sym::SpinSU2 >, double >	Psinglet (const std::size_t locx1, const std::size_t locx2, const std::size_t locy1=0ul, const std::size_t locy2=0ul) const

MpoTerms< Sym::SU2< Sym::SpinSU2 >, double >	spin_swap_operator_D3 (const std::size_t locx1, const std::size_t locx2, const std::size_t locy1=0ul, const std::size_t locy2=0ul) const

	HeisenbergObservables ()

	HeisenbergObservables (const size_t &L)

	HeisenbergObservables (const size_t &L, const vector< Param > &params, const std::map< string, std::any > &defaults)

std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	S (size_t locx, size_t locy=0, double factor=1.) const

std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Sdag (size_t locx, size_t locy=0, double factor=std::sqrt(3.)) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Scomp (SPINOP_LABEL Sa, size_t locx, size_t locy=0, double factor=1.) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Qcomp (SPINOP_LABEL Sa, size_t locx, size_t locy=0, double factor=1.) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Sz (size_t locx, size_t locy=0) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Sx (size_t locx, size_t locy=0) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	iSy (size_t locx, size_t locy=0) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Sp (size_t locx, size_t locy=0) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Sm (size_t locx, size_t locy=0) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, complex< double > > >::type	exp_ipiSz (size_t locx, size_t locy=0, double factor=1.) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	ScompScomp (SPINOP_LABEL Sa1, SPINOP_LABEL Sa2, size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0, double fac=1.) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	QcompQcomp (SPINOP_LABEL Sa1, SPINOP_LABEL Sa2, size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0, double fac=1.) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	SpSm (size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0, double fac=1.) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	SmSp (size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0, double fac=1.) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	SzSz (size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0) const

std::enable_if< Dummy::NO_SPIN_SYM(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	SxSx (size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0) const

std::conditional< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double >, vector< Mpo< Sym::SU2< Sym::SpinSU2 >, double > > >::type	SdagS (size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0, double factor=1.) const

std::conditional< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double >, vector< Mpo< Sym::SU2< Sym::SpinSU2 >, double > > >::type	QdagQ (size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0) const

std::conditional< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double >, vector< Mpo< Sym::SU2< Sym::SpinSU2 >, double > > >::type	SdagSxS (size_t locx1, size_t locx2, size_t locx3, size_t locy1=0, size_t locy2=0, size_t locy3=0) const

std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Stot (size_t locy1=0, double factor=1., int dLphys=1) const

std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Sdagtot (size_t locy1=0, double factor=std::sqrt(3.), int dLphys=1) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Scomptot (SPINOP_LABEL Sa, size_t locy1=0, double factor=1., int dLphys=1) const

std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Q (size_t locx, size_t locy=0, double factor=1.) const

std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	Qdag (size_t locx, size_t locy=0, double factor=std::sqrt(5.)) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, complex< double > > >::type	Rcomp (SPINOP_LABEL Sa, size_t locx, size_t locy=0) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	String (STRING STR, size_t locx, size_t locy=0) const

std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::SU2< Sym::SpinSU2 >, double > >::type	StringCorr (STRING STR, size_t locx1, size_t locx2, size_t locy1=0, size_t locy2=0) const

Public Member Functions inherited from ParamReturner
	ParamReturner ()

	ParamReturner (const std::map< string, std::any > &defaults_input)

DMRG::CONTROL::DYN	get_DmrgDynParam (const vector< Param > &params={}) const

DMRG::CONTROL::GLOB	get_DmrgGlobParam (const vector< Param > &params={}) const

VUMPS::CONTROL::DYN	get_VumpsDynParam (const vector< Param > &params={}) const

VUMPS::CONTROL::GLOB	get_VumpsGlobParam (const vector< Param > &params={}) const

Static Public Member Functions
static qarray< 1 >	singlet (int N=0)

static void	set_operators (const std::vector< SpinBase< Symmetry > > &B, const ParamHandler &P, PushType< SiteOperator< Symmetry, double >, double > &pushlist, std::vector< std::vector< std::string > > &labellist, const BC boundary=BC::OPEN)

Static Public Member Functions inherited from Mpo< Sym::SU2< Sym::SpinSU2 >, double >
static std::vector< T >	get_N_site_interaction (T const &Op0, Operator const &... Ops)

static Mpo< Sym::SU2< Sym::SpinSU2 >, double >	Identity (const std::vector< std::vector< qType > > &qPhys, const qType &Q=Symmetry::qvacuum())

static Mpo< Sym::SU2< Sym::SpinSU2 >, double >	Zero (const std::vector< std::vector< qType > > &qPhys)

static Mpo< Sym::SU2< Sym::SpinSU2 >, double >	cast_Terms_to_Mpo (const MpoTerms< Sym::SU2< Sym::SpinSU2 >, double > &input)

Static Public Member Functions inherited from MpoTerms< Symmetry, Scalar >
static MpoTerms< Symmetry, Scalar >	prod (const MpoTerms< Symmetry, Scalar > &top, const MpoTerms< Symmetry, Scalar > &bottom, const qType &qTot, const double tolerance=::mynumeric_limits< double >::epsilon())

static MpoTerms< Symmetry, Scalar >	sum (const MpoTerms< Symmetry, Scalar > &top, const MpoTerms< Symmetry, Scalar > &bottom, const double tolerance=::mynumeric_limits< double >::epsilon())

Static Public Attributes
static constexpr MODEL_FAMILY	FAMILY = HEISENBERG

static const std::map< string, std::any >	defaults

static const std::map< string, std::any >	sweep_defaults

Static Public Attributes inherited from Mpo< Sym::SU2< Sym::SpinSU2 >, double >
static constexpr size_t	Nq

Private Types
typedef Eigen::Index	Index

typedef Symmetry::qType	qType

typedef Eigen::SparseMatrix< double >	SparseMatrixType

Additional Inherited Members
Public Attributes inherited from Mpo< Sym::SU2< Sym::SpinSU2 >, double >
Symmetry::qType	qVac

bool	GOT_TWO_SITE_DATA

std::vector< std::vector< TwoSiteData< Sym::SU2< Sym::SpinSU2 >, double > > >	TSD

bool	UNITARY

bool	HERMITIAN

bool	HAMILTONIAN

bool	GOT_SEMIOPEN_LEFT

bool	GOT_SEMIOPEN_RIGHT

OperatorType	LocalOp

int	LocalSite

Public Attributes inherited from MpoTerms< Symmetry, Scalar >
std::vector< std::vector< qType > >	qOp

std::vector< Qbasis< Symmetry > >	qAux

std::vector< std::vector< qType > >	qPhys

std::vector< std::vector< std::vector< std::vector< Biped< Symmetry, MatrixType > > > > >	W

std::vector< std::vector< std::vector< qType > > >	qOp_powers

std::vector< std::vector< Qbasis< Symmetry > > >	qAux_powers

std::vector< std::vector< std::vector< std::vector< std::vector< Biped< Symmetry, MatrixType > > > > > >	W_powers

std::size_t	N_sites

std::size_t	N_phys = 0ul

reversedData	reversed

std::string	label = "MPO"

Protected Member Functions inherited from MpoTerms< Symmetry, Scalar >
void	increment_auxdim (const std::size_t loc, const qType &q)

void	increment_first_auxdim_OBC (const qType &qIn)

void	increment_last_auxdim_OBC (const qType &qOut)

void	decrement_auxdim (const std::size_t loc, const qType &q)

void	decrement_first_auxdim_OBC (const qType &qIn)

void	decrement_last_auxdim_OBC (const qType &qOut)

void	add (const std::size_t loc, const OperatorType &op, const qType &qIn, const qType &qOut, const std::size_t IndexIn, const std::size_t IndexOut)

std::size_t	get_auxdim (const std::size_t loc, const qType &q) const

void	assert_hilbert (const std::size_t loc, const std::size_t dim)

bool	eliminate_linearlyDependent_rows (const std::size_t loc, const qType &qIn, const double tolerance)

bool	eliminate_linearlyDependent_cols (const std::size_t loc, const qType &qOut, const double tolerance)

std::map< qType, std::vector< std::map< qType, OperatorType > > >	delete_row (const std::size_t loc, const qType &qIn, const std::size_t row_to_delete, const bool SAMESITE)

std::map< qType, std::vector< std::map< qType, OperatorType > > >	delete_col (const std::size_t loc, const qType &qOut, const std::size_t col_to_delete, const bool SAMESITE)

void	add_to_row (const std::size_t loc, const qType &qIn, const std::size_t row, const std::map< qType, std::vector< std::map< qType, OperatorType > > > &ops, const Scalar factor)

void	add_to_col (const std::size_t loc, const qType &qOut, const std::size_t col, const std::map< qType, std::vector< std::map< qType, OperatorType > > > &ops, const Scalar factor)

void	calc_qOp ()

void	calc_qAux ()

void	calc_W ()

void	got_update ()

void	initialize ()

void	compress (const double tolerance)

void	renormalize ()

void	clear_opLabels ()

std::tuple< std::size_t, std::size_t, std::size_t, std::size_t >	auxdim_infos () const

std::vector< qType >	calc_qList (const std::vector< OperatorType > &opList)

void	fill_O_from_W ()

Protected Member Functions inherited from HeisenbergObservables< Sym::SU2< Sym::SpinSU2 > >
Mpo< Sym::SU2< Sym::SpinSU2 >, double >	make_local (size_t locx, size_t locy, const OperatorType &Op, double factor=1., bool HERMITIAN=false, STRING STR=NOSTRING) const

Mpo< Sym::SU2< Sym::SpinSU2 >, double >	make_localSum (const vector< OperatorType > &Op, vector< double > factor, bool HERMITIAN) const

Mpo< Sym::SU2< Sym::SpinSU2 >, double >	make_corr (size_t locx1, size_t locx2, size_t locy1, size_t locy2, const OperatorType &Op1, const OperatorType &Op2, qarray< Symmetry::Nq > Qtot, double factor, bool HERMITIAN, STRING STR=NOSTRING) const

Mpo< Sym::SU2< Sym::SpinSU2 >, complex< double > >	make_FourierYSum (string name, const vector< OperatorType > &Ops, double factor, bool HERMITIAN, const vector< complex< double > > &phases) const

Static Protected Member Functions inherited from MpoTerms< Symmetry, Scalar >
static void	prod_swap_IBC (std::vector< std::map< std::array< qType, 2 >, std::vector< std::vector< std::map< qType, OperatorType > > > > > &O_out, std::vector< std::size_t > &row_qVac, std::vector< std::size_t > &col_qVac, std::vector< std::size_t > &row_qTot, std::vector< std::size_t > &col_qTot)

static void	prod_delZeroCols_OBC (std::map< std::array< qType, 2 >, std::vector< std::vector< std::map< qType, OperatorType > > > > &O_last, Qbasis< Symmetry > &qAux_last, Qbasis< Symmetry > &qAux_prev, const qType &qTot, const std::size_t col_qTot)

static std::string	power_to_string (std::size_t power)

static std::pair< std::string, std::size_t >	detect_and_remove_power (const std::string &name_w_power)

Protected Attributes inherited from MpoTerms< Symmetry, Scalar >
std::vector< std::map< std::array< qType, 2 >, std::vector< std::vector< std::map< qType, OperatorType > > > > >	O

std::vector< std::map< qType, std::size_t > >	auxdim

qType	qTot

qType	qVac = Symmetry::qvacuum()

std::size_t	pos_qTot

std::size_t	pos_qVac

BC	boundary_condition

std::vector< std::vector< std::string > >	info

std::vector< int >	hilbert_dimension

MPO_STATUS::OPTION	status = MPO_STATUS::ALTERABLE

DMRG::VERBOSITY::OPTION	VERB = DMRG::VERBOSITY::OPTION::SILENT

std::string	before_verb_set

std::size_t	current_power = 1

bool	GOT_QOP = false

bool	GOT_QAUX = false

bool	GOT_W = false

std::vector< bool >	GOT_QPHYS

Protected Attributes inherited from HeisenbergObservables< Sym::SU2< Sym::SpinSU2 > >
vector< SpinBase< Sym::SU2< Sym::SpinSU2 > > >	B

Member Typedef Documentation

◆ Index

typedef Eigen::Index VMPS::HeisenbergSU2::Index

private

Definition at line 42 of file HeisenbergSU2.h.

◆ MatrixType

typedef Eigen::Matrix<double,Eigen::Dynamic,Eigen::Dynamic> VMPS::HeisenbergSU2::MatrixType

Definition at line 34 of file HeisenbergSU2.h.

◆ OperatorType

typedef SiteOperatorQ<Symmetry,MatrixType> VMPS::HeisenbergSU2::OperatorType

Definition at line 35 of file HeisenbergSU2.h.

◆ qType

typedef Symmetry::qType VMPS::HeisenbergSU2::qType

private

Definition at line 43 of file HeisenbergSU2.h.

◆ SparseMatrixType

typedef Eigen::SparseMatrix<double> VMPS::HeisenbergSU2::SparseMatrixType

private

Definition at line 44 of file HeisenbergSU2.h.

◆ Symmetry

typedef Sym::SU2<Sym::SpinSU2> VMPS::HeisenbergSU2::Symmetry

Definition at line 32 of file HeisenbergSU2.h.

Constructor & Destructor Documentation

◆ HeisenbergSU2() [1/3]

VMPS::HeisenbergSU2::HeisenbergSU2 ( )

inline

Do nothing.

Definition at line 52 of file HeisenbergSU2.h.

◆ HeisenbergSU2() [2/3]

VMPS::HeisenbergSU2::HeisenbergSU2	(	const size_t &	L,
		const vector< Param > &	params = `{}`,
		const BC &	boundary = `BC::OPEN`,
		const DMRG::VERBOSITY::OPTION &	VERB = `DMRG::VERBOSITY::ON_EXIT`
	)

Parameters

L	: chain length
params	: Parmaterpack....
boundary	: Boundary condition. Either BC::OPEN, BC::PERIODIC or BC::INFINITE.

Definition at line 115 of file HeisenbergSU2.h.

◆ HeisenbergSU2() [3/3]

VMPS::HeisenbergSU2::HeisenbergSU2	(	Mpo< Symmetry > &	Mpo_input,
		const vector< Param > &	params
	)

inline

Do nothing.

Definition at line 61 of file HeisenbergSU2.h.

Member Function Documentation

◆ set_operators()

void VMPS::HeisenbergSU2::set_operators	(	const std::vector< SpinBase< Symmetry > > &	B,
		const ParamHandler &	P,
		PushType< SiteOperator< Symmetry, double >, double > &	pushlist,
		std::vector< std::vector< std::string > > &	labellist,
		const BC	boundary = `BC::OPEN`
	)

static

This function merges the different local operators together to a PushType instance (See DmrgTypedefs.h for the type definition.). This pushlist is used for the construction of the Mpo.

Parameters

B	: Base class from which the local operators are received
P	: The parameters
pushlist	: All the local operators for the Mpo will be pushed into `pushlist`.
labellist	: All the labels for the Mpo will be put into `labellist`. Mpo::generate_label will produce a nice label from the data in labellist.
boundary	: Boundary condition. Either BC::OPEN, BC::PERIODIC or BC::INFINITE.

Definition at line 153 of file HeisenbergSU2.h.

◆ singlet()

static qarray< 1 > VMPS::HeisenbergSU2::singlet ( int N = 0 )

inlinestatic

Definition at line 37 of file HeisenbergSU2.h.

◆ validate()

bool VMPS::HeisenbergSU2::validate ( qarray< 1 > qnum ) const

Validates whether a given total quantum number qnum is a possible target quantum number for an Mps.

Returns: true if valid, false if not

Definition at line 143 of file HeisenbergSU2.h.

Member Data Documentation

◆ defaults

const std::map< string, std::any > VMPS::HeisenbergSU2::defaults

static

Initial value:

= 
{
    {"J",1.}, {"Jprime",0.}, {"Jprimeprime",0.}, {"Jrung",1.},
    {"R",0.}, {"Offset",0.},
    {"D",2ul}, {"maxPower",2ul}, {"CYLINDER",false}, {"Ly",1ul}, {"mfactor",1}
}

Definition at line 93 of file HeisenbergSU2.h.

◆ FAMILY

constexpr MODEL_FAMILY VMPS::HeisenbergSU2::FAMILY = HEISENBERG

staticconstexpr

Definition at line 38 of file HeisenbergSU2.h.

◆ sweep_defaults

const std::map< string, std::any > VMPS::HeisenbergSU2::sweep_defaults

static

Initial value:

= 
{
    {"max_alpha",100.}, {"min_alpha",1.e-11}, {"lim_alpha",10ul}, {"eps_svd",1.e-7},
    {"Mincr_abs", 50ul}, {"Mincr_per", 2ul}, {"Mincr_rel", 1.1},
    {"min_Nsv",0ul}, {"max_Nrich",-1},
    {"max_halfsweeps",50ul}, {"min_halfsweeps",20ul},
    {"Minit",1ul}, {"Qinit",1ul}, {"Mlimit",1000ul},
    {"tol_eigval",1e-7}, {"tol_state",1e-6},
    {"savePeriod",0ul}, {"CALC_S_ON_EXIT", true}, {"CONVTEST", DMRG::CONVTEST::VAR_HSQ}
}

Definition at line 94 of file HeisenbergSU2.h.

The documentation for this class was generated from the following file:

HeisenbergSU2.h

Detailed Description

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Private Types

Additional Inherited Members

Member Typedef Documentation

◆ Index

◆ MatrixType

◆ OperatorType

◆ qType

◆ SparseMatrixType

◆ Symmetry

Constructor & Destructor Documentation

◆ HeisenbergSU2() [1/3]

◆ HeisenbergSU2() [2/3]

◆ HeisenbergSU2() [3/3]

Member Function Documentation

◆ set_operators()

◆ singlet()

◆ validate()

Member Data Documentation

◆ defaults

◆ FAMILY

◆ sweep_defaults