- Comp-phys-alps-users - lists.phys.ethz.ch

Fulldiag
by zhian asadzadeh 14 Apr '09

14 Apr '09

Hi Dear, I am using fulldiag application for my problem but when I use this command: *fulldiag_evaluate parm1.task1.out.xml* I got this error: *bash: fulldiag_evaluate: command not found* And also I have not, for example, Magnetization and susceptebility and ... in the parm1.task.out.xml. What is wrong? Thanks, Zhian.

1 0

Specific heat and Susceptibility with DMRG
by zhian asadzadeh 12 Apr '09

12 Apr '09

Hello All, Let me know could we calculate Specific and Susceptibility with the DMRG's application of current's ALPS version? Thanks, Zhian

2 1

Energy Gap
by zhian asadzadeh 10 Apr '09

10 Apr '09

Hello Dear, yes it worked. If we don't have a good symmetry in the Hamiltonian for example anisotropic spin chain in transversed magnetic field, we don't know easily, in what sector ground state and first exited stated lie ? in this case is it true that I trust to the results which have given with setting just NUM_EIGENVALUES=2 in parm file? Thank you very much, Zhian. On Tue, Mar 24, 2009 at 2:30 PM, <comp-phys-alps-users-request(a)phys.ethz.ch>wrote: > Send Comp-phys-alps-users mailing list submissions to > comp-phys-alps-users(a)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(a)phys.ethz.ch > > You can reach the person managing the list at > comp-phys-alps-users-owner(a)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. dynamical correlation functions (jessica.alfonsi(a)unipd.it) > 2. Energy Gap (zhian asadzadeh) > 3. Re: Energy Gap (Adrian E. Feiguin) > 4. uniform magnetic susceptibility (Guilin Zhuang) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 23 Mar 2009 12:48:14 +0100 (CET) > From: jessica.alfonsi(a)unipd.it > Subject: [ALPS-users] dynamical correlation functions > To: comp-phys-alps-users(a)phys.ethz.ch > Message-ID: > <25a2823199cb7c2f2d21169a68d1a5b0.squirrel(a)webmail.unipd.it> > Content-Type: text/plain;charset=iso-8859-1 > > Hello everybody, > I'm interested in starting using ALPS project, in particular for computing > optical spectral functions starting from exact diagonalization and > dynamical correlation function for optical conductivity. However, I read > in papers about ALPS 1.3 release that these will be implemented in release > 1.4 . When do you plan to release at least the first beta of ALPS 1.4 ? Is > it possible to use a patched SVN version ? > > Thank you and best wishes > > Jessica Alfonsi > > > > > ------------------------------ > > Message: 2 > Date: Tue, 24 Mar 2009 09:57:55 +0330 > From: zhian asadzadeh <zhian.asadzadeh(a)gmail.com> > Subject: [ALPS-users] Energy Gap > To: comp-phys-alps-users(a)phys.ethz.ch > Message-ID: > <89e5d6760903232327p5fe2f942x8fde9c344cceb0d2(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Hi Dear, > > I tried to calculate energy gap for ferrimagnetic ladder. in this paper ( > Phy Rev B, 63, 054432) it has been done which you can find in > table I. I used the periodic boundary condition because in the case of open > boundary condition the energy converges very slowly. > I have achieved the following result. > E0= -0.757036 > E1= -0.757032 > > I have used the following input file: I have tried for inter chain coupling > J' = 0.2 with the notation of this paper. > > LATTICE="ferrimagnetic ladder" > SWEEPS=4 > MAXSTATES=300 > MODEL="spin", L=30 > local_S0=1/2,local_S1=1,Jxy0=1,h=0.0,Jxy1=0.2,Jz1=0.2 > BC1=periodic,BC2=periodic > NUMBER_EIGENVALUES=2 > Jz0=1 > {Gamma=0.0;} > > > which L is the number of unit cell and h, Gamma are longitudinal and > transverse magnetic field respectively. Jz0 , Jxy0 are chain couplings and > Jz1, and Jxy1 are inter chain couplings. > we see that the result is not concur with this paper for energy gap. > I see in this paper lanczos method is used but I expect that the result be > same in 4 digit. > please let me know what is wrong? > > Thanks, > Mohammad. > > > > > > > On Wed, Mar 18, 2009 at 2:30 PM, > <comp-phys-alps-users-request(a)phys.ethz.ch>wrote: > > > Send Comp-phys-alps-users mailing list submissions to > > comp-phys-alps-users(a)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(a)phys.ethz.ch > > > > You can reach the person managing the list at > > comp-phys-alps-users-owner(a)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. Energy Gap (zhian asadzadeh) > > 2. Re: Energy Gap (Matthias Troyer) > > > > > > ---------------------------------------------------------------------- > > > > Message: 1 > > Date: Tue, 17 Mar 2009 14:43:46 +0330 > > From: zhian asadzadeh <zhian.asadzadeh(a)gmail.com> > > Subject: [ALPS-users] Energy Gap > > To: comp-phys-alps-users(a)phys.ethz.ch > > Message-ID: > > <89e5d6760903170413w7e72ef85i873b8d47564f4320(a)mail.gmail.com> > > Content-Type: text/plain; charset="iso-8859-1" > > > > Hi Dear, > > > > Thank you for your consideration. > > My model is anisotropic Ferrimagnetic chain in transverse magnetic field. > > when I say the result is wrong I mean when I have tried to obtain energy > > gap > > for some problems and models which have been solved earlier I saw the > > result doesn't agree with those earlier results. for example spin one and > > half-one chain and ferrimagnetic chain and anisotropic ferrimagnetic > > ladder. > > Something is not very obvius for me, how DMRG calculate ground state and > > first exited state from one wave function ?! > > > > Thanks, > > Zhian. > > > > > > On Tue, Mar 17, 2009 at 2:30 PM, > > <comp-phys-alps-users-request(a)phys.ethz.ch>wrote: > > > > > Send Comp-phys-alps-users mailing list submissions to > > > comp-phys-alps-users(a)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(a)phys.ethz.ch > > > > > > You can reach the person managing the list at > > > comp-phys-alps-users-owner(a)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: Energy Gap (Matthias Troyer) > > > > > > > > > ---------------------------------------------------------------------- > > > > > > Message: 1 > > > Date: Mon, 16 Mar 2009 12:33:11 -0700 > > > From: Matthias Troyer <troyer(a)phys.ethz.ch> > > > Subject: Re: [ALPS-users] Energy Gap > > > To: comp-phys-alps-users(a)phys.ethz.ch > > > Message-ID: <43E825E5-7BE3-4EA0-9C2A-20A1FC8C21AF(a)phys.ethz.ch> > > > Content-Type: text/plain; charset=US-ASCII; format=flowed; delsp=yes > > > > > > > > > On Mar 16, 2009, at 2:01 AM, zhian asadzadeh wrote: > > > > > > > Hi All, > > > > > > > > I want to calculate energy gap for spin chain with DMRG. I set > > > > NUM_EIGENVALUES=2 in my parm file but the result which achieved is > > > > wrong. > > > > Is it necessary to determine something else for example the sector > > > > which DMRG work in it? in my problem I dont have any good symmetry > > > > and I can not easily determine the sector which the first exited > > > > state is in it. I would be very appreciate if someone could please > > > > help me. > > > > > > > > Thanks, > > > > Zhian. > > > > > > Can you please send more details about your model, or better yet the > > > input file? An what do you mean exactly by saying that the result is > > > wrong? > > > > > > Matthias > > > > > > > > > > > > End of Comp-phys-alps-users Digest, Vol 36, Issue 10 > > > **************************************************** > > > > >

2 1

Install the library without g77
by Guilin Zhuang 10 Apr '09

10 Apr '09

Dear all users, This question have confused me for a long time: Last week, I reinstall the new version Linux OS (Ubuntu 8.10),in which gcc package does not include the g77 compiler. I also tried to install but failed. Do I compile the library without the support of fortran 77,but with the support the gfortran or intel fortran? If yes, how should I do? Best wishes! Thank you for your attention to my question! Zhuang

2 1

basis states along with corresponding energies
by khalid hassan 09 Apr '09

09 Apr '09

Dear alps team, I asked for the basis staties and hamiltonian matrix for the system described by parameters in case of exact diagonalization. Dr Matthias Troyer told me of a program in alps/model directory. This program is print_numeric and it displays the basis states and hamiltonian matrix. Now I am interested with the basis states of lowest energy, as in case case of sparsediag which only has lowest energy states or fulldiag with firse states. I want that as sparse diag or full diag give the output energies it also gives the corresponding states along with those energies. this is what I want, then I can use the lowest energy basis states for my research as the professor will told me. I have look in the alps code it hard for me as i dont know c++ of the level it needs and learning c++ is taking much of the time. If the author of the program write a small routine just as print_numeric, it will serve my purpose, it will be great help of me. Otherwise please tell me to look in what files or modify what programs that i get the basis state of the model corresponding to its energy. Best Regards Khalid Loane

1 0

DMRG question
by Justin David Peel 07 Apr '09

07 Apr '09

Hi all, I hope this is okay to ask here, but I didn't know where else to ask it. I've been trying to do a different periodic DMRG algorithm that splits the first site so that it is in both the left and right blocks. This is in trying to do a DMRG algorithm kind of like the matrix product state sort that works so much better. I've been having some problems with doing this and was wondering if anyone had any ideas. I'm only considering a 1D Heisenberg S=1/2 model for now. First of all, I keep track of the Sz operator of the first site of the left block and the last site of the right block (which are really the same site). I make the superblock hamiltonian and diagonalize within the states that have Sz_total=0 for lowest energy AND Sz of first site of left block equal to the last site of the right block (they should have the same Sz spin since they are essentially the same site). When making the superblock hamiltonian, the Sz terms don't require any special consideration. However, the S_plus/S_minus ones do because when the 1st site of the left block is raised, the last site of the right block must also be raised, but the energy should be that of only raising one site. This is where there are problems. I have tried several ways but been unsuccessful. What must be made is an operator that raises both sites as if they are one. Any ideas out there? Is this possible to do at all? Is it impossible because of the minus sign problem with the electrons? Let me know if I wasn't clear enough. Also, when I rotate the a block after finding the highest eigenvalued eigenvectors of the reduced density matrix, I have to find the eigenvectors that preserve both Sz_total and Sz of the first site (in the case of the left block). I'm thinking that this will make the method less effective, but do you think that this will make the method ineffective? Thanks for your time and thoughts, Justin Peel

1 0

configure flags for ALPS applications
by Andrea Taroni 07 Apr '09

07 Apr '09

Dear All, I am having trouble installing ALPS applications on a Ubuntu 8.04 system. Installing the libraries is not a problem. Installing the applications wouldn't be, if I were not interested in including the lapack libraries in the installation. I've followed the instructions for the setup and installation of the applications, where it states that adding the option --with-lapack=-lLIB when running configure "specifies the library files for the LAPACK libraries". When I do this, I get an output to the effect of "configure: WARNING: Unrecognized options: --with-lapack". The same occurs if I include the option "--with-lapack-dir". Essentially, my question is what am I doing wrong here? Surely the option should be recognised by the configure file? Any help on this issue is greatly appreciated. Kind Regards, Andrea ps I suspect my problem has something to do with this: https://webmail.phys.ethz.ch/pipermail/comp-phys-alps-users/2008-May/000253… but I can't find much other information out there.

2 2

Classical Heisenberg model
by Gianni Profeta 06 Apr '09

06 Apr '09

Dear users, I'm trying to perform classical MC simulation for the 3D Heisenberg model. The model is formed by random magnetic impurities in periodic lattice interacting with different Jij for various neighboring distances r between the magnetic sistes (i and j). Any suggestions on how to define different J's between different sites? Thanks a lot, Gianni

4 7

Re: [ALPS-users] working on alps code
by Jakub Zakrzewski 28 Mar '09

28 Mar '09

Dear Zhenya, Many thanks!! It is not ugly it is beautiful and it works !! All the best jakub > Message: 3 > Date: Fri, 27 Mar 2009 21:48:00 +0100 > From: Evgeny Burovskiy <burovski(a)alumni.umass.edu> > Subject: Re: [ALPS-users] working on alps code > To: comp-phys-alps-users(a)phys.ethz.ch > Message-ID: <49CD3B80.9080803(a)physik.uni-greifswald.de> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > > Dear Jakub, > > > So I guess from examples .. but I cannot guess everything. To be specific. > > I run QMC codes for Bose-Hubbard. I know how to get occupation at different > > sides (for inhomogenous model, e.g. with harmonic trap and/or disorder), > > call them <n_i>. But I am interested in second moment <n_i^2> > > > > > > There is a solution which looks ugly, but works : > > patch the definition of the site basis of the model you're working > with. > In your example, Bose-Hubbard Hamiltonian uses the basis called > "boson". > Go to models.xml file and then add a line to the SITEBASIS definition: > > <SITEBASIS name="boson"> > <PARAMETER name="Nmax" default="infinity"/> > <QUANTUMNUMBER name="N" min="0" max="Nmax"/> > <OPERATOR name="bdag" matrixelement="sqrt(N+1)"> > <CHANGE quantumnumber="N" change="1"/> > </OPERATOR> > <OPERATOR name="b" matrixelement="sqrt(N)"> > <CHANGE quantumnumber="N" change="-1"/> > </OPERATOR> > <OPERATOR name="n" matrixelement="N"/> > <OPERATOR name="n2" matrixelement="N*N"/>  > </SITEBASIS> > > Like I said, it's ugly, but it works. > And, like you, I'd be very much interested in learning how to add > measurements properly, without dirty tricks :). > > > Best, > > > Zhenya > > > End of Comp-phys-alps-users Digest, Vol 36, Issue 16 > ****************************************************

1 0

Re: [ALPS-users] working on alps code
by Jakub Zakrzewski 27 Mar '09

27 Mar '09

Dear Matthias, 1. I congratulate you and your colleagues for creating such a fantastic tool as ALPS is. 2. See point 1. 3. Still there is a lot of truth in what Khalid Hassan wrote. To be specific: - I cannot go to Zurich for personal consultation. Both from Paris where I am on sabbatical and from Krakow (POland) it is a bit too far just for asking a couple of simple questions. Instead I go to Tutorial and Documentation link on https://alps.ethz.ch/mediawiki/index.php/Main_Page and want to read about: Definition of custom measurements or Common hardcoded measurements to know what I can get. But these sections of the Documentation are EMPTY. So I guess from examples .. but I cannot guess everything. To be specific. I run QMC codes for Bose-Hubbard. I know how to get occupation at different sides (for inhomogenous model, e.g. with harmonic trap and/or disorder), call them <n_i>. But I am interested in second moment <n_i^2> What I can do? a/ I wrote a long email to you b/ 2 times I posed questions on the very same subject on this list without getting any answer. From more experienced colleague, I learn that I may modify qmc.h so I get <n_i^2> INSTEAD of <n_i> - but this requires recompilation of applications and running the whole simulation twice. There MUST be a more intelligent way of getting the info I need. Or correlations <a_i a_j^\dagger> or <n_i n_j> .... May somebody give an advice on that? 4. Simon Trebst was happy that you reached 100 citations to ALPS (I guess this is both to old and new reference). I bet you, were there more user friendly documentation available and more active list of users, the number of citations would increase tremendously... I am sure that at least by 1 since the things we work on now (with ALPS) could be finished in a more reasonable ammount of time... All the best Jakub > > > > ------------------------------ > > Message: 2 > Date: Thu, 26 Mar 2009 20:19:22 -0700 > From: Matthias Troyer <troyer(a)phys.ethz.ch> > > Dear Khalid Loane, > > The ALPS libraries are meant to be used in user code, and are > documented. The applications had been designed as black box codes. > However, if you would like to learn more about them and work on one of > > them, we will be delighted to help. Best would be a personal visit > though to explain the code structure - and clean it up while doing that. > > Matthias > > > > > On 26 Mar 2009, at 14:12, khalid hassan wrote: > > > > > alps is a program that works well with a comprehensive help. but if > > > one has use its code for one self research, the code provides little > > > help. it has very few comments and difficult to understand. Can > > there can be any help for the code of full diagonalization for > > example logical structure, flow chart, algorithm available. If this > > > is in personal files of the developers can they send by email or if > > > someone visits them they can give him. > > > > Best Regards > > Khalid Loane > > ITCP, Graz Austria > > > > > > > > > > End of Comp-phys-alps-users Digest, Vol 36, Issue 15 > ****************************************************

2 1