Hi All,

I have asked this question earlier, how I can change the code for measuring advanced correlation functions? I mean I must refer to what file to modify ?

Thanks,
zhian.

On Sun, May 31, 2009 at 2:30 PM, <comp-phys-alps-users-request@phys.ethz.ch> wrote:
Send Comp-phys-alps-users mailing list submissions to
       comp-phys-alps-users@phys.ethz.ch

To subscribe or unsubscribe via the World Wide Web, visit
       https://webmail.phys.ethz.ch/mailman/listinfo/comp-phys-alps-users
or, via email, send a message with subject or body 'help' to
       comp-phys-alps-users-request@phys.ethz.ch

You can reach the person managing the list at
       comp-phys-alps-users-owner@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:

  1. Re: current correlation function for fermionic hubbard    model
     (Matthias Troyer)
  2. Re: Setting up polaritonic model in ALPS (Matthias Troyer)
  3. Re: Advanced correlation functions (Matthias Troyer)
  4. Re: Advanced correlation functions (Matthias Troyer)
  5. Re: Problems with worms?? (Matthias Troyer)
  6. Re: any news about next  alps release ? (Matthias Troyer)


----------------------------------------------------------------------

Message: 1
Date: Sat, 30 May 2009 14:22:22 -0600
From: Matthias Troyer <troyer@phys.ethz.ch>
Subject: Re: [ALPS-users] current correlation function for fermionic
       hubbard model
To: comp-phys-alps-users@phys.ethz.ch, khalid hassan
       <mkhloane@yahoo.com>
Message-ID: <6EAA816E-F4FC-498A-ABCF-C2B8A1BB080A@phys.ethz.ch>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes


On 26 Feb 2009, at 14:25, khalid hassan wrote:

>
> Can alps can calculate the current correlation function for the
> fermionic hubbard model model for say two site. I have been using
> alps to calculate cahrge coorelation function
>
> the charge correlation function is
>
> < n_i n_j > = < c_i^+ c_i  c_i^+ c_i >

This can be done

>
>
> where as the current correlation function is
>
> < j_i j_j > = < c_i^+ c_j  c_k^+ c_l >

We have not implemented bond-bond correlation functions yet, but it
can be done if you are interested.

Matthias



------------------------------

Message: 2
Date: Sat, 30 May 2009 16:49:02 -0600
From: Matthias Troyer <troyer@phys.ethz.ch>
Subject: Re: [ALPS-users] Setting up polaritonic model in ALPS
To: comp-phys-alps-users@phys.ethz.ch
Message-ID: <E831C875-FB09-40A6-9FC9-68CAE96C4CAA@phys.ethz.ch>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes

On 29 May 2009, at 19:09, Su Chun-Hsu wrote:
>
> Dear all,
> I'm new to using ALPS and not entirely confident with setting up my
> model,
> which didn't appear to work at the moment. I'm interested in
> polaritonic
> system -- evanescently coupled-cavity system, each cavity containing a
> two-level atom. I have set up my model in the attached xml, Sminus
> and Splus
> are atomic projection operators and b are photon annihilation
> operator.
> The total number of excitations (atomic and photonic) is conserved.
> When I try running it, there is an bad_alloc error. Will anyone
> be kindly give me some pointers? Much appreciated and thank you very
> much
> in advance.

Not knowing exactly what you have in mind I can only say that the
model looks reasonable. You have to be very careful though as you have
a large number of local states and already a 2x2 lattice without
constraints has a Hilbert space of 202^4, which is far too large to be
solved by exact diagonalization.

With your constraints the Hilbert space is smaller but there is still
the problem that the ALPS program wants to build the full Hamiltonian
matrix for two sites. The 2-site Hilbert space in your model is 202^2
= 40804 states. The matrix for two sites thus 40804 x 40804 = 1 664
966 416 elements large, which needs 12.5 GB memory. The allocation of
this matrix fails and thus there is a bad_alloc exception. Such a
model is to big to be studied by exact diagonalization. If you reduce
the maximum occupancy t0 10 or 20 it will work

Matthias



------------------------------

Message: 3
Date: Sat, 30 May 2009 16:50:51 -0600
From: Matthias Troyer <troyer@phys.ethz.ch>
Subject: Re: [ALPS-users] Advanced correlation functions
To: comp-phys-alps-users@phys.ethz.ch
Message-ID: <CB7794C7-1A03-4A39-BD73-24CF100175D2@phys.ethz.ch>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes

Just define Sx, which is real:

<SITEOPERATOR name="Sx" site="i">
1/2*(Splus(i)+Sminus(i))
</SITEOPERATOR>

On 28 Apr 2009, at 12:44, Justin David Peel wrote:

> You're correct. It isn't quite right. It should be divided by i:
>
> <SITEOPERATOR name="Sy" site="x">
> 1/2*(Splus(x)-Sminus(x))/sqrt(-1)
> </SITEOPERATOR>
>
>
> ________________________________________
> From: comp-phys-alps-users-bounces@phys.ethz.ch [comp-phys-alps-users-bounces@phys.ethz.ch
> ] On Behalf Of zhian asadzadeh [zhian.asadzadeh@gmail.com]
> Sent: Tuesday, April 28, 2009 5:56 AM
> To: comp-phys-alps-users@phys.ethz.ch
> Subject: [ALPS-users] Advanced correlation functions
>
> Hi Dear,
>
> Thank you for your guidance, but the mentioned definition for Sy
> doesn't work.
> And for the first question, how I can change the code for measuring
> those  two kind of  correlations? I mean the correlation of first
> site with others and middle site with others too.
>
> Thanks,
> Zhian.
>
>
>
>
> On Tue, Apr 28, 2009 at 2:30 PM, <comp-phys-alps-users-request@phys.ethz.ch
> <mailto:comp-phys-alps-users-request@phys.ethz.ch>> wrote:
> Send Comp-phys-alps-users mailing list submissions to
>       comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >
>
> To subscribe or unsubscribe via the World Wide Web, visit
>       https://webmail.phys.ethz.ch/mailman/listinfo/comp-phys-alps-users
> or, via email, send a message with subject or body 'help' to
>       comp-phys-alps-users-request@phys.ethz.ch<mailto:comp-phys-alps-users-request@phys.ethz.ch
> >
>
> You can reach the person managing the list at
>       comp-phys-alps-users-owner@phys.ethz.ch<mailto:comp-phys-alps-users-owner@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:
>
>  1. Re: lapack support in ALPS (Synge Todo)
>  2. Re: Advanced correlation functions (Justin David Peel)
>  3. Re: Comp-phys-alps-users Digest, Vol 37, Issue 18
>     (Justin David Peel)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 27 Apr 2009 19:51:14 +0900
> From: Synge Todo <wistaria@comp-phys.org<mailto:wistaria@comp-
> phys.org>>
> Subject: Re: [ALPS-users] lapack support in ALPS
> To: comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >
> Message-ID: <04B17733-4F07-4867-8B22-4343337160BB@comp-phys.org<mailto:04B17733-4F07-4867-8B22-4343337160BB@comp-phys.org
> >>
> Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes
>
> Dear Fu-Jiun Jiang,
>
> On 2009/04/23, at 23:39, Fu-Jiun Jiang ITP wrote:
>> When install the ALPS with lapack being supported, what command
>> should be used? In the instruction, it is said
>> --with-lapack-dir=/path/to/my/lapack/ --with-lapack="-l my_lapack".
>> Is it true that first one should point to where the lapack library
>> is.
>> How about the second one? Thanks.
>
> In many cases, lapack installation will be found automatically without
> any options.  If this is not the case, please your
> system administrator about the linker options to link your program
> against lapack.  For example, if you are linking
> liblapack.a and libblas.a in /path/to/my/lapack, then the configure
> options should be
>
>  --with-lapack-dir=/path/to/my/lapack --with-lapack="-llapack -lblas"
>
> Synge
>
>
>
> ------------------------------
>
> Message: 2
> Date: Mon, 27 Apr 2009 11:35:24 -0600
> From: Justin David Peel <justin.peel@utah.edu<mailto:justin.peel@utah.edu
> >>
> Subject: Re: [ALPS-users] Advanced correlation functions
> To: "comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >"
>       <comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >>
> Message-ID:
>       <3E5F591958142A40B693A5FC6672D4942624680106@C8V1.xds.umail.utah.edu
> <mailto:3E5F591958142A40B693A5FC6672D4942624680106@C8V1.xds.umail.utah.edu
> >>
> Content-Type: text/plain; charset="us-ascii"
>
> I'm not sure about only measuring those two correlations (you might
> have to change the code for that), but, as to your last question,
> you should be able to use Sx:Sx because the Sx operator is defined
> in the models.xml file in terms of the raising and lowering
> operators. If you want to measure Sy, you would have to add the
> following lines to the models.xml file:
>
> <SITEOPERATOR name="Sy" site="x">
> 1/2*(Splus(x)-Sminus(x))
> </SITEOPERATOR>
>
> because it hasn't been defined as far as I can tell.
>
> Good luck,
> Justin
> ________________________________________
> From: comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> > [comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> >] On Behalf Of zhian asadzadeh [zhian.asadzadeh@gmail.com<mailto:zhian.asadzadeh@gmail.com
> >]
> Sent: Sunday, April 26, 2009 1:54 AM
> To: Comp-phys-alps-users@phys.ethz.ch<mailto:Comp-phys-alps-users@phys.ethz.ch
> >
> Subject: [ALPS-users] Advanced correlation functions
>
> Hi All,
>
> I am using DMRG for my problem and I am going to measure diagonal
> and offdiagonal correlation functions.
> I am interest to obtain the correlation of first site with other
> sites and  the middle site with other sites too.
> The output of DMRG consist of   all the combinations which are
> possible but I just want the two correlations mentioned above
> because the output will be very large ( I am studying the
> ferrimagnetic chain with 100 sites).
> can I do that? how?
> And also can I use Sx:Sx for measuring off diagonal correlation
> functions?
>
> Thanks,
> Zhian.
>
>
>
> ------------------------------
>
> Message: 3
> Date: Mon, 27 Apr 2009 11:42:56 -0600
> From: Justin David Peel <justin.peel@utah.edu<mailto:justin.peel@utah.edu
> >>
> Subject: Re: [ALPS-users] Comp-phys-alps-users Digest, Vol 37, Issue
>       18
> To: "comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >"
>       <comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >>
> Message-ID:
>       <3E5F591958142A40B693A5FC6672D4942624680107@C8V1.xds.umail.utah.edu
> <mailto:3E5F591958142A40B693A5FC6672D4942624680107@C8V1.xds.umail.utah.edu
> >>
> Content-Type: text/plain; charset="us-ascii"
>
> Someone correct me if I'm wrong, but I believe this isn't very
> directly feasible. The ground state vector from the exact ground
> state will be in a different basis and have a different length than
> the ground state vector from DMRG (it's a truncated basis after all).
> ________________________________________
> From: comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> > [comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> >] On Behalf Of Gang Chen [chggst@physics.ucsb.edu<mailto:chggst@physics.ucsb.edu
> >]
> Sent: Saturday, April 25, 2009 9:55 AM
> To: comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >
> Subject: Re: [ALPS-users] Comp-phys-alps-users Digest, Vol 37, Issue
> 18
>
> Thanks, I want to compute the overlap between the exact ground state
> and the state I generated from other method.
>
> best,
> Gang
>
> On Sat, Apr 25, 2009 at 3:00 AM, <comp-phys-alps-users-request@phys.ethz.ch
> <mailto:comp-phys-alps-users-request@phys.ethz.ch><mailto:comp-phys-alps-users-request@phys.ethz.ch
> <mailto:comp-phys-alps-users-request@phys.ethz.ch>>> wrote:
> Send Comp-phys-alps-users mailing list submissions to
>      comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> ><mailto:comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >>
>
> To subscribe or unsubscribe via the World Wide Web, visit
>      https://webmail.phys.ethz.ch/mailman/listinfo/comp-phys-alps-
> users
> or, via email, send a message with subject or body 'help' to
>      comp-phys-alps-users-request@phys.ethz.ch<mailto:comp-phys-alps-users-request@phys.ethz.ch
> ><mailto:comp-phys-alps-users-request@phys.ethz.ch<mailto:comp-phys-alps-users-request@phys.ethz.ch
> >>
>
> You can reach the person managing the list at
>      comp-phys-alps-users-owner@phys.ethz.ch<mailto:comp-phys-alps-users-owner@phys.ethz.ch
> ><mailto:comp-phys-alps-users-owner@phys.ethz.ch<mailto:comp-phys-alps-users-owner@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:
>
> 1. extract the ground state from DMRG or Exact       diagonalization
>    (Gang Chen)
> 2. Re: extract the ground state from DMRG or Exact
>    diagonalization (Justin David Peel)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 24 Apr 2009 14:57:02 -0700
> From: Gang Chen <chggst@physics.ucsb.edu<mailto:chggst@physics.ucsb.edu
> ><mailto:chggst@physics.ucsb.edu<mailto:chggst@physics.ucsb.edu>>>
> Subject: [ALPS-users] extract the ground state from DMRG or Exact
>      diagonalization
> To: comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> ><mailto:comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >>
> Message-ID:
>      <3924ab970904241457k17ab51b1tbd002d09ae7287e@mail.gmail.com<mailto:3924ab970904241457k17ab51b1tbd002d09ae7287e@mail.gmail.com
> ><mailto:3924ab970904241457k17ab51b1tbd002d09ae7287e@mail.gmail.com<mailto:3924ab970904241457k17ab51b1tbd002d09ae7287e@mail.gmail.com
> >>>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi All,
>
> If I use DMRG or Exact diagonalization to find the ground state of my
> Hamiltonian. How can
> I extract the ground state vector (the numerical vector) from my
> results?
>
> best,
>
> --
> GANG CHEN
> Department of Physics
> University of California
> Santa Barbara, CA 93106
> Office: Broida 6216
> Phone: 805-893-5260
> Fax: 805-893-3378
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: https://webmail.phys.ethz.ch/pipermail/comp-phys-alps-users/attachments/20090424/1340d4ea/attachment.html
>
> ------------------------------
>
> Message: 2
> Date: Fri, 24 Apr 2009 16:17:36 -0600
> From: Justin David Peel <justin.peel@utah.edu<mailto:justin.peel@utah.edu
> ><mailto:justin.peel@utah.edu<mailto:justin.peel@utah.edu>>>
> Subject: Re: [ALPS-users] extract the ground state from DMRG or Exact
>      diagonalization
> To: "comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> ><mailto:comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >>"
>      <comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> ><mailto:comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >>>
> Message-ID:
>      <3E5F591958142A40B693A5FC6672D4942624680104@C8V1.xds.umail.utah.edu
> <mailto:3E5F591958142A40B693A5FC6672D4942624680104@C8V1.xds.umail.utah.edu
> ><mailto:3E5F591958142A40B693A5FC6672D4942624680104@C8V1.xds.umail.utah.edu
> <mailto:3E5F591958142A40B693A5FC6672D4942624680104@C8V1.xds.umail.utah.edu
> >>>
> Content-Type: text/plain; charset="us-ascii"
>
> Hi,
>
> Perhaps you can tell us why you want the vector - what you plan on
> doing with it because there isn't much point of knowing the ground
> state vector from DMRG unless you have the operator matrices to go
> along with it.
>
> Justin
> ________________________________________
> From: comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> ><mailto:comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> >> [comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> ><mailto:comp-phys-alps-users-bounces@phys.ethz.ch<mailto:comp-phys-alps-users-bounces@phys.ethz.ch
> >>] On Behalf Of Gang Chen [chggst@physics.ucsb.edu<mailto:chggst@physics.ucsb.edu
> ><mailto:chggst@physics.ucsb.edu<mailto:chggst@physics.ucsb.edu>>]
> Sent: Friday, April 24, 2009 3:57 PM
> To: comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> ><mailto:comp-phys-alps-users@phys.ethz.ch<mailto:comp-phys-alps-users@phys.ethz.ch
> >>
> Subject: [ALPS-users] extract the ground state from DMRG or
> Exact       diagonalization
>
> Hi All,
>
> If I use DMRG or Exact diagonalization to find the ground state of
> my Hamiltonian. How can
> I extract the ground state vector (the numerical vector) from my
> results?
>
> best,
>
> --
> GANG CHEN
> Department of Physics
> University of California
> Santa Barbara, CA 93106
> Office: Broida 6216
> Phone: 805-893-5260
> Fax: 805-893-3378
>
>
> End of Comp-phys-alps-users Digest, Vol 37, Issue 18
> ****************************************************
>
>
>
>
> --
> GANG CHEN
> Department of Physics
> University of California
> Santa Barbara, CA 93106
> Office: Broida 6216
> Phone: 805-893-5260
> Fax: 805-893-3378
>
>
> End of Comp-phys-alps-users Digest, Vol 37, Issue 20
> ****************************************************



------------------------------

Message: 4
Date: Sat, 30 May 2009 16:51:33 -0600
From: Matthias Troyer <troyer@phys.ethz.ch>
Subject: Re: [ALPS-users] Advanced correlation functions
To: comp-phys-alps-users@phys.ethz.ch
Message-ID: <EF110C15-B270-463F-9DC8-3A2D3D42FCD3@phys.ethz.ch>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes


On 30 Apr 2009, at 01:04, Ryo IGARASHI wrote:

> Hi Zhian,
>
> On Thu, Apr 30, 2009 at 3:59 PM, zhian asadzadeh
> <zhian.asadzadeh@gmail.com> wrote:
>> Please let me know can I change the code for measuring just two
>> kinds of
>> correlations?
>> (Correlation of first site with other sites and the middle site
>> with other
>> sites).
>
>> From the licencing point of view,  there is no restriction for
> modifying the ALPS code.
> Just go ahead. Of course, patches are always welcome.

To modify the DMRG code you can ask Adrian Feiguin for help

Matthias



------------------------------

Message: 5
Date: Sat, 30 May 2009 16:54:33 -0600
From: Matthias Troyer <troyer@phys.ethz.ch>
Subject: Re: [ALPS-users] Problems with worms??
To: comp-phys-alps-users@phys.ethz.ch
Message-ID: <84A3C6A6-16BD-40D7-B4D3-B8B5DB710764@phys.ethz.ch>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes

We are aware of this problem and will fix it soon.

Matthias

On 12 May 2009, at 08:39, Jakub Zakrzewski wrote:

> Hello,
>
> Sometimes QMC worm code (used on Bose-Hubbard) terminates with the
> following error:
>
> here come first earlier succesful simulations then
> ...
> Created run 85 remote on Host ID: 84
> Created run 86 remote on Host ID: 85
> Created run 87 remote on Host ID: 86
> Created run 88 remote on Host ID: 87
> All processes have been assigned
> Checking if Simulation 1 is finished: not yet, next check in 120
> seconds ( 0% done).
> q = -0 state1 = 0 state2 = 0 bond_type = 0
> zero matrix element in remove_jump
> application called MPI_Abort(MPI_COMM_WORLD, -2) - process 67
>
> Could somebody advice what I can do?
>
> Best regards  Kuba



------------------------------

Message: 6
Date: Sat, 30 May 2009 16:55:32 -0600
From: Matthias Troyer <troyer@phys.ethz.ch>
Subject: Re: [ALPS-users] any news about next  alps release ?
To: jessica.alfonsi@unipd.it, comp-phys-alps-users@phys.ethz.ch
Message-ID: <067AA75A-302C-447C-BAF9-393177A1ABDF@phys.ethz.ch>
Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes


On 27 May 2009, at 02:45, jessica.alfonsi@unipd.it wrote:

>
> Hello everybody,
> I asked some time ago, but I got no reply. When do you plan to
> release the
> first beta of ALPS 1.4 ? Is there any hope to see it by this
> summer ? I
> ask because I would be interested into the new exact diagonalization
> features...

We have not implemented these features yet because none of us needed
it. What exactly do you require?

Matthias



End of Comp-phys-alps-users Digest, Vol 38, Issue 11
****************************************************