Thank you very much for your reply.
I have solved the problem with my own DMRG code.
However, I will keep tracking ALPS project. It is indeed very great.
Thank you for your work.
Send Comp-phys-alps-users mailing list submissions to
comp-phys-alps-users@lists.phys.ethz.ch
To subscribe or unsubscribe via the World Wide Web, visit
https://lists.phys.ethz.ch/listinfo/comp-phys-alps-users
or, via email, send a message with subject or body 'help' to
comp-phys-alps-users-request@lists.phys.ethz.ch
You can reach the person managing the list at
comp-phys-alps-users-owner@lists.phys.ethz.ch
When replying, please edit your Subject line so it is more specific
than "Re: Contents of Comp-phys-alps-users digest..."
Today's Topics:
- Re: How to define an operator A in "models.xml" (Matthias Troyer)
- multiple eigenvalues (degenerate eigenstates) in sparsediag
(Andriy Zhugayevych)
Message: 1
Date: Thu, 24 Mar 2011 12:28:28 -0500
From: Matthias Troyer troyer@phys.ethz.ch
Subject: Re: [ALPS-users] How to define an operator A in "models.xml"
To: comp-phys-alps-users@lists.phys.ethz.ch
Message-ID: DAC79DA4-5BC5-49DA-BAAE-1B885B6736BB@phys.ethz.ch
Content-Type: text/plain; charset="utf-8"
Dear Sunzhaoyu,
This is not implemented yet. For now you will have to decompose it into
more elementary operators.
Matthias
On Mar 18, 2011, at 12:38 AM, ??? wrote:
Dear ALPS users and developers
I want to re-produce the result of a 12x2 spin ladder with four-site
interaction.
The four-site interaction is between two nearest-neighboring rungs.
Please see:
[ Phys. Rev. B 74, 155119 (2006),
http://prb.aps.org/abstract/PRB/v74/i15/e155119 ]
However, it seems that ALPS 2.0 do not support 4-site terms.
Thus I define every rung as a single site,
then the ladder becomes a one-dimensional chain (local degree of freedom
is 4),
and the four-site interaction becomes a two-site interaction.
In original ladder, the SITEBASIS is just {sz} = {-1/2, 1/2}.
In this basis operators such as Sminus have very compact expressions
thus they can be simply define as:
<OPERATOR name="Sminus" matrixelement="sqrt(S*(S+1)-Sz*(Sz-1))">
However, in the new model, a possible choice of SITEBASIS is {J,Jz}={00,
1-1, 10, 11}.
In this new basis, operators do not have very compact expressions,
however, the elements of the operator can be easily obtained through
unitary transformation.
The problem is, I can figure out all the elements of an operator in
basis{J, Jz}, but I do not know how to define the operator in "models.xml".
In other words,
suppose an operator X in basis {J, Jz} is explicitly found to be
[ 1, 0, 0, 0
0, 2, 1, 0
0, 1, 3, 0
0, 0, 0, 4]
then in ALPS 2.0, is it possible to define X by simply writing down every
elements Xij in "models.xml" ?
Any suggestion would be appreciated.
sunzhaoyu