VMPS++
Loading...
Searching...
No Matches
VMPS::Heisenberg Class Reference

Detailed Description

Heisenberg Model.

MPO representation of

\[
H =  J \sum_{<ij>} \left(\mathbf{S_i} \cdot \mathbf{S_j}\right) 
     J' \sum_{<<ij>>} \left(\mathbf{S_i} \cdot \mathbf{S_j}\right)
    -B_z \sum_i S^z_i
    -B_x \sum_i S^x_i
    +K_z \sum_i \left(S^z_i\right)^2
    +K_x \sum_i \left(S^x_i\right)^2
    -D_y \sum_{<ij>} \left(\mathbf{S_i} \times \mathbf{S_j}\right)_y
    -D_y' \sum_{<<ij>>} \left(\mathbf{S_i} \times \mathbf{S_j}\right)_y
\]

Parameters
D: $D=2S+1$ where $S$ is the spin
Note
Uses no symmetries. Any parameter constellations are allowed. For variants with symmetries, see VMPS::HeisenbergU1 or VMPS::HeisenbergSU2.
The default variable settings can be seen in Heisenberg::defaults.
$J>0$ is antiferromagnetic
This is the real version of the Heisenbergmodel without symmetries, so $J_x = J_y$ is mandatory. For general couplings use VMPS::HeisenbergXYZ.

Definition at line 33 of file Heisenberg.h.

#include <Heisenberg.h>

Inheritance diagram for VMPS::Heisenberg:

Public Types

typedef Sym::U0 Symmetry
 
- Public Types inherited from Mpo< Sym::U0, double >
typedef SparseMatrixXd SparseMatrixType
 
typedef SiteOperator< Sym::U0, double > OperatorType
 
typedef Symmetry::qType qType
 
typedef Matrix< double, Dynamic, Dynamic > MatrixType
 
typedef double Scalar_
 
typedef Mps< Sym::U0, double > StateXd
 
typedef Umps< Sym::U0, double > StateUd
 
typedef Mps< Sym::U0, std::complex< double > > StateXcd
 
typedef Umps< Sym::U0, std::complex< double > > StateUcd
 
typedef MpsCompressor< Sym::U0, double, double > CompressorXd
 
typedef MpsCompressor< Sym::U0, std::complex< double >, double > CompressorXcd
 
typedef Mpo< Sym::U0Operator
 
- 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_TYPEMatrixType
 

Public Member Functions

 Heisenberg ()
 
 Heisenberg (const size_t &L, const vector< Param > &params, const BC &boundary=BC::OPEN, const DMRG::VERBOSITY::OPTION &VERB=DMRG::VERBOSITY::OPTION::ON_EXIT)
 
 Heisenberg (Mpo< Symmetry > &Mpo_input, const vector< Param > &params)
 
- Public Member Functions inherited from Mpo< Sym::U0, double >
 Mpo ()
 
 Mpo (MpoTerms< Sym::U0, 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::U0, double > Identity () const
 
Mpo< Sym::U0, 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 qTypeget_qVac () const
 
const qTypeget_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 qTypeQtarget () 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::U0 >
std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, complex< double > > >::type S_ky (vector< complex< double > > phases) const
 
std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, complex< double > > >::type Sdag_ky (vector< complex< double > > phases, double factor=sqrt(3.)) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, complex< double > > >::type Sa_ky (SPINOP_LABEL Sa, vector< complex< double > > phases) const
 
std::vector< Mpo< Sym::U0, double > > make_spinPermutation (const Permutation &permutations) const
 
MpoTerms< Sym::U0, 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::U0, 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::U0, 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::U0, double > >::type S (size_t locx, size_t locy=0, double factor=1.) const
 
std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, 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::U0, 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::U0, 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::U0, double > >::type Sz (size_t locx, size_t locy=0) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, double > >::type Sx (size_t locx, size_t locy=0) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, double > >::type iSy (size_t locx, size_t locy=0) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, double > >::type Sp (size_t locx, size_t locy=0) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, double > >::type Sm (size_t locx, size_t locy=0) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, 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::U0, 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::U0, 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::U0, 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::U0, 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::U0, 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::U0, 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::U0, double >, vector< Mpo< Sym::U0, 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::U0, double >, vector< Mpo< Sym::U0, 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::U0, double >, vector< Mpo< Sym::U0, 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::U0, double > >::type Stot (size_t locy1=0, double factor=1., int dLphys=1) const
 
std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, 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::U0, 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::U0, double > >::type Q (size_t locx, size_t locy=0, double factor=1.) const
 
std::enable_if< Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, 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::U0, complex< double > > >::type Rcomp (SPINOP_LABEL Sa, size_t locx, size_t locy=0) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, double > >::type String (STRING STR, size_t locx, size_t locy=0) const
 
std::enable_if<!Dummy::IS_SPIN_SU2(), Mpo< Sym::U0, 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< 0 > singlet (int N=0)
 
static void add_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 refEnergy ref (const vector< Param > &params, double L=numeric_limits< double >::infinity())
 
- Static Public Member Functions inherited from Mpo< Sym::U0, double >
static std::vector< T > get_N_site_interaction (T const &Op0, Operator const &... Ops)
 
static Mpo< Sym::U0, double > Identity (const std::vector< std::vector< qType > > &qPhys, const qType &Q=Symmetry::qvacuum())
 
static Mpo< Sym::U0, double > Zero (const std::vector< std::vector< qType > > &qPhys)
 
static Mpo< Sym::U0, double > cast_Terms_to_Mpo (const MpoTerms< Sym::U0, 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::U0, double >
static constexpr size_t Nq
 

Private Types

typedef Symmetry::qType qType
 

Additional Inherited Members

- Public Attributes inherited from Mpo< Sym::U0, double >
Symmetry::qType qVac
 
bool GOT_TWO_SITE_DATA
 
std::vector< std::vector< TwoSiteData< Sym::U0, 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< qTypecalc_qList (const std::vector< OperatorType > &opList)
 
void fill_O_from_W ()
 
- Protected Member Functions inherited from HeisenbergObservables< Sym::U0 >
Mpo< Sym::U0, double > make_local (size_t locx, size_t locy, const OperatorType &Op, double factor=1., bool HERMITIAN=false, STRING STR=NOSTRING) const
 
Mpo< Sym::U0, double > make_localSum (const vector< OperatorType > &Op, vector< double > factor, bool HERMITIAN) const
 
Mpo< Sym::U0, 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::U0, 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::U0 >
vector< SpinBase< Sym::U0 > > B
 

Member Typedef Documentation

◆ qType

Definition at line 43 of file Heisenberg.h.

◆ Symmetry

Definition at line 36 of file Heisenberg.h.

Constructor & Destructor Documentation

◆ Heisenberg() [1/3]

VMPS::Heisenberg::Heisenberg ( )
inline

Definition at line 48 of file Heisenberg.h.

◆ Heisenberg() [2/3]

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

Definition at line 109 of file Heisenberg.h.

◆ Heisenberg() [3/3]

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

Definition at line 52 of file Heisenberg.h.

Member Function Documentation

◆ add_operators()

void VMPS::Heisenberg::add_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

Definition at line 148 of file Heisenberg.h.

◆ ref()

refEnergy VMPS::Heisenberg::ref ( const vector< Param > &  params,
double  L = numeric_limits<double>::infinity() 
)
static

Definition at line 298 of file Heisenberg.h.

◆ singlet()

static qarray< 0 > VMPS::Heisenberg::singlet ( int  N = 0)
inlinestatic

Definition at line 39 of file Heisenberg.h.

Member Data Documentation

◆ defaults

const std::map< string, std::any > VMPS::Heisenberg::defaults
static
Initial value:
=
{
{"J",0.}, {"Jprime",0.}, {"Jrung",0.},
{"Jxy",0.}, {"Jxyprime",0.}, {"Jxyrung",0.},
{"Jz",0.}, {"Jzprime",0.}, {"Jzrung",0.},
{"Jx",0.}, {"Jxrung",0.},
{"Jy",0.}, {"Jyrung",0.},
{"Jz",0.}, {"Jzrung",0.},
{"Jw",0.}, {"Jwrung",0.},
{"betaKT",0.},
{"R",0.},
{"Bz",0.}, {"Bx",0.},
{"Kz",0.}, {"Kx",0.},
{"Dy",0.}, {"Dyprime",0.}, {"Dyrung",0.},
{"t",0.}, {"mu",0.}, {"Delta",0.},
{"mu",0.}, {"nu",0.},
{"D",2ul}, {"maxPower",2ul}, {"CYLINDER",false}, {"Ly",1ul}, {"mfactor",1}
}

Definition at line 71 of file Heisenberg.h.

◆ FAMILY

constexpr MODEL_FAMILY VMPS::Heisenberg::FAMILY = HEISENBERG
staticconstexpr

Definition at line 40 of file Heisenberg.h.

◆ sweep_defaults

const std::map< string, std::any > VMPS::Heisenberg::sweep_defaults
static
Initial value:
=
{
{"max_alfa",100.}, {"min_alfa",1.e-11}, {"lim_alfa",10ul}, {"eps_svd",1.e-7},
{"Mincr_abs", 10ul}, {"Mincr_per", 2ul}, {"Mincr_rel", 1.1},
{"min_Nsv",0ul}, {"max_Nrich",-1},
{"max_halfsweeps",40ul}, {"min_halfsweeps",1ul},
{"Minit",10ul}, {"Qinit",1ul}, {"Mlimit",500ul},
{"tol_eigval",1.e-5}, {"tol_state",1.e-5},
{"savePeriod",0ul}, {"CALC_S_ON_EXIT", true}, {"CONVTEST", DMRG::CONVTEST::VAR_2SITE}
}

Definition at line 72 of file Heisenberg.h.


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