Comp-phys-alps-users

comp-phys-alps-users@lists.phys.ethz.ch

1244 discussions

ALPS 2.0.1 released
by Matthias Troyer 19 Sep '12

19 Sep '12

Dear ALPS users, ALPS 2.0.1 has been released. This is the version included with the ALPS paper: B. Bauer et al., J. Stat. Mech. (2011) P05001, which is available at : http://www.iop.org/EJ/abstract/1742-5468/2011/05/P05001 .Besides references to this paper, this release contains many minor bug fixes. The MacOS X 10.5 binary release will be uploaded shortly. The others are ready for download. Best regards Matthias Troyer

4 33

questions about worm code
by Lev 18 Sep '12

18 Sep '12

Dear Dr. Troyer and users of ALPS, I am trying to measure local magnetization with MEASURE_LOCAL[Local magnetization]=Sx for the following model: <HAMILTONIAN name="my spin model 1"> <BASIS ref="spin"/> <SITETERM site="i"> -2*xi*Sz(i)-2*h*Sx(i) </SITETERM> <BONDTERM source="i" target="j"> -g*(Splus(i)*Sminus(j)+Sminus(i)*Splus(j)) </BONDTERM> </HAMILTONIAN> Value of parameter xi depends on the site in an actual problem, but this is not important now. Since dirloop_sse says that "Offdiagonal site terms are not implemented in this SSE code", I use the worm code. Worm code says the following for this: "Will not measure "Local magnetization" since it is off-diagonal or not a site operator". As a trick to overcome this problem, I tried to rotate the coordinate system in order to measure the required values with MEASURE_LOCAL[Local magnetization]=Sz In this case the rotated model is: <SITEOPERATOR name="Sy" site="x"> 1/2/I*(Splus(x)-Sminus(x)) </SITEOPERATOR> <HAMILTONIAN name="my spin model 2"> <BASIS ref="spin"/> <SITETERM site="i"> -2*xi*Sy(i)-2*h*Sz(i) </SITETERM> <BONDTERM source="i" target="j"> -2*g*Sz(i)*Sz(j)-2*g*Sx(i)*Sx(j) </BONDTERM> </HAMILTONIAN> Now the worm code says the following: "can not convert complex number into real one". Therefore, the worm code does not seem to understand Sy(i) as a site term in Hamiltonian (?). Finally, I tried another version of rotation, which also would allow to measure with MEASURE_LOCAL[Local magnetization]=Sz The rotated model is: <SITEOPERATOR name="Sy" site="x"> 1/2/I*(Splus(x)-Sminus(x)) </SITEOPERATOR> <HAMILTONIAN name="my spin model 3"> <BASIS ref="spin"/> <SITETERM site="i"> 2*xi*Sx(i)-2*h*Sz(i) </SITETERM> <BONDTERM source="i" target="j"> -2*g*Sz(i)*Sz(j)-2*g*Sy(i)*Sy(j) </BONDTERM> </HAMILTONIAN> Now the worm code says the following: "Cannot simulate this bond term with the worm code". Therefore, the worm code does not seem to understand Sy(i)*Sy(j) as a bond term in Hamiltonian (?). Could you let me know whether there are ways to perform the measurement ? Sincerely yours, Lev Barash

2 1

Dynamical quantities with ALPS ED
by Mykola Maksymenko 13 Sep '12

13 Sep '12

Dear everybody, I am curious I it is possible to write a function in ALPS ED code to calculate the dynamical responses in a spectral representation ? It seems like in a simplest case one can compute it directly by having access to the matrix elements of the operator and corresponding energy eigenvalues of the Hamiltonian. Isn't it ? Will be thankful for the answer! With best regards, Mykola

1 0

Binary installation of ALPS at an alternative location in windows
by Chao Shen 09 Sep '12

09 Sep '12

Hi all, I try to install ALPS binary in a position other than 'C:\program files' because there I would need administrator privilege to create files, etc. I noticed on ALPS website it says I just need to set environment variable ALPS_XML_DIR to point to the position of the actual XML directory, e.g. 'E:\ALPS\lib\xml'. However this doesn't really have any effect. When I run 'tutorial1a.py' which calls 'spinmc.exe' in the tutorial folder, 'spinmc.exe' complains about missing 'lattice.xml' at location 'C:\program files\ALPS\lib\xml'. Apparently the environment variable ALPS_XML_DIR didn't take effect at all. I am wondering if this is a bug or there are actually other environment variables to set? Thanks in advance. Yours, Chao

1 0

Re: [ALPS-users] Comp-phys-alps-users Digest, Vol 77, Issue 5
by 陈巧妮 30 Aug '12

30 Aug '12

I'm calculating the one particle reduced density matrix, so it's the cdag_up:c_up, and cdag_down:c_down. 2012/8/30 <comp-phys-alps-users-request(a)lists.phys.ethz.ch> > Send Comp-phys-alps-users mailing list submissions to > comp-phys-alps-users(a)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(a)lists.phys.ethz.ch > > You can reach the person managing the list at > comp-phys-alps-users-owner(a)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: > > 1. using sparsediag (???) > 2. Re: using sparsediag (Matthias Troyer) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 29 Aug 2012 21:35:13 -0400 > From: ??? <qiaoni2233(a)gmail.com> > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Subject: [ALPS-users] using sparsediag > Message-ID: > <CAEM5Jm67=pHMzhVxKwZuPKFJjiNpw8A3+LSSa3hG0MJ9zEAC= > Q(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Hi, everyone > I'm using sparsediag to get the ground state of Hubbard model, e.g. L=12, > with the particle number constraint Nup_total=6, Ndown_total=6. Using the > command line. If I calculate the correlation function , it takes 900 sec, > if not it takes 100 sec, does anyone know how to speed up this calculation? > > Btw if I use python to run the result, I can't see the result from the > task.out.xml file, how to look up the result > > -- > Qiaoni Chen > Department of Chemistry, > Princeton University > Princeton, NJ, 08544 > Tel +1-609-865-0817 >

1 1

using sparsediag
by 陈巧妮 30 Aug '12

30 Aug '12

Hi, everyone I'm using sparsediag to get the ground state of Hubbard model, e.g. L=12, with the particle number constraint Nup_total=6, Ndown_total=6. Using the command line. If I calculate the correlation function , it takes 900 sec, if not it takes 100 sec, does anyone know how to speed up this calculation? Btw if I use python to run the result, I can't see the result from the task.out.xml file, how to look up the result -- Qiaoni Chen Department of Chemistry, Princeton University Princeton, NJ, 08544 Tel +1-609-865-0817

2 1

ngs mcresult operators
by august＠fzu.cz 21 Aug '12

21 Aug '12

Hi all, I found a strange behavior of ngs mcresult operators "+" and "*". I made a fresh instalation of ALPS-2.1.1-r6193, edited /alps/applications/dmft/qmc/hybridization2/main.C and added the following code between lines 52 and 53: alps::mcresult img, temp; img = results["gw_im"]; temp = img + img; std::vector<double> temp_mean, img_mean; temp_mean = temp.mean<std::vector<double> >(); img_mean = img.mean<std::vector<double> >(); I would expect temp_mean to be just twice the img_mean. However, it seems to be 2*img_mean*img.bin_size()/img_mean.size(). Similarly, the "*" oerator yields wrongly scalled result as well (I do not know what the scaling factor is). On the contrary, "img-img" and "img/img" gives correct answers. Is it a bug or just my misunderstanding of what these operators should do? Besides that, I would like to use the mcresult::covariance method, but it is not implemented yet. Could I ask anyone to change that? Thank you for any help or comments! Best regards, Pavel

1 0

pyalps.loadMeasurements & MPI -> Floating point exception
by Kevin Slagle 17 Aug '12

17 Aug '12

Hello When I run the attached python file, the simulation runs way too fast (ie it isn't running correctly) and then crashes at the pyalps.loadMeasurements function since the data output isn't valid. The cause seems to be the MPI=4 parameter in the pyalps.runApplication function. I realize that MPI=4 is silly for running over only one set of parameters (eg different temperatures). However, I have similar issues when I run over multiple parameters. Any ideas what I need to do to fix this? Or is this a reproducible bug (seems unlikely..)? my system: cpu = 64bit Ivy Bridge i7 distro = Arch Linux (packages up to date) $ ./Heisenburg.py mpirun -np 4 dirloop_sse --mpi Heisenburg.in.xml --Tmin 1 Quantum Monte Carlo simulations using the SSE algorithm v. 4.1.1 available from http://alps.comp-phys.org/ copyright (c) 2003-2010 by Sergei Isakov <isakov(a)itp.phys.ethz.ch> using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer <troyer(a)comp-phys.org>. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998). based on the ALPS libraries version 2.1.1 available from http://alps.comp-phys.org/ copyright (c) 1994-2012 by the ALPS collaboration. Consult the web page for license details. For details see the publication: B. Bauer et al., J. Stat. Mech. (2011) P05001. parsing task files ... Creating a new simulation: 1 on 4 processes Created run 1 locally Created run 2 remote on Host ID: 1 Created run 3 remote on Host ID: 2 Created run 4 remote on Host ID: 3 All processes have been assigned Checking if Simulation 1 is finished: not yet, next check in 1 seconds ( 0% done). Checking if Simulation 1 is finished: not yet, next check in 1 seconds ( 58% done). NUMBER_OF_WORMS_PER_SWEEP: 9 NUMBER_OF_WORMS_PER_SWEEP: 9 NUMBER_OF_WORMS_PER_SWEEP: 9 NUMBER_OF_WORMS_PER_SWEEP: 9 Checking if Simulation 1 is finished: Finished Halted Simulation 1 Checkpointing Simulation 1 Finished with all tasks. loading measurements: Floating point exception thanks! Kevin

2 1

non-standard constraint: "NL + nR" value="N_total"
by Mateusz Łącki 04 Aug '12

04 Aug '12

Dear All, I want to define a model for boson-boson mixture to be run by a worm code. I have two types of bosons: of type L and type R. I would like to be able to consider a conservation of total number of bosons that is : nL + nR = const. What should I put into "contraint": <CONSTRAINT quantumnumber="N" value="N_total"/> ? I would like to have something like <CONSTRAINT quantumnumber="NL + nR" value="N_total"/> An option is probably to redefine quantum numbers from: nR and nL into nA and nB by nA = (nR + nL) and nB = (nL-nR) (now nA is >=0 and nB is arbitrary) of course one would have to change creation and annihilation operators, all observables... this i why I would like to avoid it. Regards, Mateusz Łącki

2 2

Question
by Majid Abedi 30 Jul '12

30 Jul '12

Hi I run the script '$ alpspython build_lattice.py 32 > my_lattice.xml'. The coordinate for site "1" is "0". Other sites are different. For example site "2" have coordination "2". Am I making a mistake or is there a reason behind that? Sincerely Majid Abedi Student of 7-year PhD Course Institute for Advanced Studies in Basic Sciences (IASBS) Gavazang road Zanjan, Iran P.O.Box: 45195 - 1159 Tel: +98-241-415-2023 (IASBS) Fax: +98-241-415-2104 W: www.iasbs.ac.ir~abedi

2 1

Jump to page:

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Comp-phys-alps-users