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

QMC of Bose-Hubbard model with NNN interaction
by Tao Ying 24 Oct '15

24 Oct '15

Dear all, I want to use QMC in ALPS to simulate Bose-Hubbard model with NNN interaction, we know for NN interaction V, we can define V as following: <BONDTERM source="i" target="j"> <PARAMETER name="V" default="0"/> V*n(i)*n(j) </BONDTERM> Now I want to include NNN interaction, how can I define it? (source="i" target="*?*") Best regards, Tao Ying -- Tao Ying, Lecturer in department of physics, Harbin Institute of Technology

2 1

The problem of measuring stiffness in QMC with defining lattice using graph
by 郭怀明 20 Oct '15

20 Oct '15

Dear All, For the lattice in alps, we can use <LATTICES> <GRAPH>, in which can define an arbitrary lattice. However I find when I use this way, the stiffness measured in QMC is always zero. I wonder how to fix it. Thank you very much. Sincerely, Huaiming Guo

1 0

MPS
by Huaiming Guo 19 Oct '15

19 Oct '15

Dear All, I try the example script provided by alps. But the following information shows that there is an error. I attach the spin_one.py, which can also be found in the alps. What is the problem? Thank you very much. Sincerely, Huaiming Guo C:\alpsproject\mps>vispython spin_one.py C:\alpsproject\mps>REM Copyright Matthias Troyer 2010-2013. C:\alpsproject\mps>REM Distributed under the Boost Software License, Version 1. 0. C:\alpsproject\mps>REM (See accompanying file LICENSE_1_0.txt or copy at C:\alpsproject\mps>REM http://www.boost.org/LICENSE_1_0.txt) C:\alpsproject\mps>IF EXIST C:\"Program Files"\VisTrails\Python27\python.exe C:\ "Program Files"\VisTrails\Python27\python.exe spin_one.py C:\alpsproject\mps>IF EXIST C:\"Program Files"\VisTrails\Python27_64\python.exe C:\"Program Files"\VisTrails\Python27_64\python.exe spin_one.py mps_optim parm_spin_one.in.xml --write-xml Traceback (most recent call last): File "spin_one.py", line 57, in <module> for s in data[0]: IndexError: list index out of range C:\alpsproject\mps>IF EXIST C:\"Program Files (x86)"\VisTrails\Python27\python.e xe C:\"Program Files (x86)"\VisTrails\Python27\python.exe spin_one.py C:\alpsproject\mps>IF EXIST C:\"Program Files"\VisTrails\vistrails\Python27\pyth on.exe C:\"Program Files"\VisTrails\vistrails\Python27\python.exe spin_one.py C:\alpsproject\mps>IF EXIST C:\"Program Files"\VisTrails\vistrails\Python27_64\p ython.exe C:\"Program Files"\VisTrails\vistrails\Python27_64\python.exe spin_one .py C:\alpsproject\mps>IF EXIST C:\"Program Files (x86)"\VisTrails\vistrails\Python2 7\python.exe C:\"Program Files (x86)"\VisTrails\vistrails\Python27\python.exe sp in_one.py

2 1

Question on collect_result
by Leo Fang 19 Oct '15

19 Oct '15

Hello, Recently I am trying to understand how the hybridization qmc solver is implemented. In line 89 of the file hybmain.cpp, I noticed that the function "collect_result" is called and the returned stuff is assigned to a variable called"results". I know "results" is an "alps::results_type<hybridization>::type" object, aka an accumulator. What I don't understand is its data structure (how elements are stored and how to access them). For example, I find the usage of "results["Sign"]", "results[order_name.str()]", etc in hybevaluate.cpp confusing, because it doesn't seem to me that the accumulator stores information in the std::map type. Can someone explain it to me or point me to relevant documentation (if any), please? Thanks. Sincerely, Leo (Fang Yao-Lung) Duke Physics

1 0

XYZ interaction in alps-looper algorithm
by Saeed S. Jahromi 17 Oct '15

17 Oct '15

Dear All, I am trying to run the QMC "loop" code for the Kitaev honeycomb model which is actually an XYZ two-spin interaction on a honeycomb lattice with three different link types i.e. x,y and z. *The parameter file that I use is defined as:* parms = [] for t in [0.05, 0.1, 0.2, 0.3]: parms.append( { 'LATTICE' : "honeycomb kitaev lattice", 'T' : t, 'Jx' : .1 , 'Jz' : .1 , 'Jy' : .1 , 'ALGORITHM' : "loop", 'THERMALIZATION' : 65536/8, 'SWEEPS' : 65536, 'MODEL' : "K-H", 'L' : 4 } ) *The model that I use is defined as: * <HAMILTONIAN name="K-H"> <PARAMETER name="Jz" default="0"/> <PARAMETER name="Jx" default="0"/> <PARAMETER name="Jy" default="0"/> <PARAMETER name="hx" default="0"/> <PARAMETER name="hz" default="0"/> <BASIS ref="spin"/> <SITETERM site="i"> -hz*Sz(i)-hx*Sx(i) </SITETERM> <BONDTERM type="0" source="i" target="j"> Jz*Sz(i)*Sz(j) </BONDTERM> <BONDTERM type="1" source="i" target="j"> Jx*Sx(i)*Sx(j) </BONDTERM> <BONDTERM type="2" source="i" target="j"> Jy*Sy(i)*Sy(j) </BONDTERM> </HAMILTONIAN> *and the lattice is defined as**:* <UNITCELL name="honeycomb inhomogeneous" dimension="2"> <VERTEX><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX><COORDINATE>-0.288665771484375 0.5773468017578125</COORDINATE></VERTEX> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> <EDGE type="1"><SOURCE vertex="2"/><TARGET vertex="1" offset="0 1"/></EDGE> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2" offset="1 -1"/></EDGE> </UNITCELL> <LATTICEGRAPH name = "honeycomb kitaev lattice"> <FINITELATTICE> <LATTICE ref="triangular lattice"/> <PARAMETER name="W" default="L"/> <EXTENT dimension="1" size="L"/> <EXTENT dimension="2" size="W"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL ref="honeycomb inhomogeneous"/> </LATTICEGRAPH> *when I run the simulation with the "loop" code, the following message is produced and the simulation stops:* multi-cluster quantum Monte Carlo algorithms for spin systems available from http://wistaria.comp-phys.org/alps-looper/ copyright (c) 1997-2010 by Synge Todo <wistaria(a)comp-phys.org> using ALPS/parapack scheduler a Monte Carlo scheduler for multiple-level parallelization copyright (c) 1997-2013 by Synge Todo <wistaria(a)comp-phys.org> based on the ALPS libraries version 2.2.b3 available from http://alps.comp-phys.org/ copyright (c) 1994-2013 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. [2015-Oct-17 14:46:49]: starting scheduler on jahromi-Sharif master input file = /home/jahromi/Desktop/MyComputations/ALPS/Looper/parmlooper.in.xml master output file = /home/jahromi/Desktop/MyComputations/ALPS/Looper/parmlooper.out.xml termination file = [disabled] total number of thread(s) = 1 thread(s) per clone = 1 number of thread group(s) = 1 auto evaluation = yes time limit = unlimited interval between checkpointing = 3600 seconds interval between progress report = 600 seconds interval between vmusage report = infinity task range = all worker dump format = hdf5 worker dump policy = running workers only [2015-Oct-17 14:46:49]: task status: total number of tasks = 4 new = 4, running = 0, continuing = 0, suspended = 0, finished = 0, completed = 0, skipped = 0 [2015-Oct-17 14:46:49]: starting 1 threadgroup(s) [2015-Oct-17 14:46:49]: dispatching a new clone[1,1] on threadgroup[1] (Jx+Jy) and (Jx-Jy) may not have opposite signs FYI, I have enable the xyz bond support as instructed in this link and recompiled alps. https://lists.phys.ethz.ch/pipermail/comp-phys-alps-users/2013/002282.html I also tested the algorithm for the "XYZ spin" model available on the alps-looper source on wistaria GitHub and the results is the same and the problem exist there, either. I would be so thankful if you could help me on this matter Best -- Saeed S. Jahromi ========================================================== Postdoctoral Researcher Physics Department, Sharif University of Technology P.O. Box 11155-9161, Tehran , Iran. Mailing Adrress: Azadi Ave., 14588-89694, Tehran, Iran Homepage: http://www.saeedjahromi.com ==========================================================

2 1

inhomogeneous Bose-Hubbard model with MPS
by Dominique Delande 09 Oct '15

09 Oct '15

Dear all, Since at least r7593 (end of July 2015), it is possible to study with MPS inhomogeneous systems, for example an 1D inhomogeneous Bose-Hubbard model with an harmonic trap by using 'inhomogeneous chain lattice' for the lattice and the following definition for the Hamiltonian <HAMILTONIAN name="harm boson Hubbard"> <PARAMETER name="mu" default="0"/> <PARAMETER name="t" default="1"/> <PARAMETER name="V" default="0"/> <PARAMETER name="t'" default="0"/> <PARAMETER name="V'" default="0"/> <PARAMETER name="U" default="0"/> <PARAMETER name="t0" default="t"/> <PARAMETER name="t1" default="t'"/> <PARAMETER name="V0" default="V"/> <PARAMETER name="V1" default="V'"/> <PARAMETER name="K" default="0"/> <BASIS ref="boson"/> <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> <PARAMETER name="K#" default="K"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2+K#*n(i)*(x-0.5*(L-1))^2 </SITETERM> <BONDTERM source="i" target="j"> <PARAMETER name="t#" default="0"/> <PARAMETER name="V#" default="0"/> -t#*(bdag(i)*b(j)+bdag(j)*b(i)) + V#*n(i)*n(j) </BONDTERM> </HAMILTONIAN> I would like to use an inhomogeneous BONDTERM. I tried things like -(t#+tharm#*(i+j-L)^2)*(bdag(i)*b(j)+bdag(j)*b(i)) or -(t#+tharm#*(x-0.5*(L-1))^2)*(bdag(i)*b(j)+bdag(j)*b(i)) in the definition of the Hamiltonian, but they all produce error messages indicating that mps_optim cannot evaluate such a Hamiltonian. I thus have two questions: * Is inhomogeneous BONDTERM possible with MPS, and how to proceed? * MORE IMPORTANTLY: How can I figure out whether it is possible or not? I find nothing in the documentation, nothing in the Wiki. Is the only solution to look in the guts of the code, or to ask on the mailing list? Thanks a lot in advance for your help. Dominique -- Dominique Delande (Dominique.Delande(a)lkb.upmc.fr) Laboratoire Kastler-Brossel - Case 74 - Universite P. et M. Curie 4, place Jussieu, F-75252 Paris Cedex 05, FRANCE Phone : +33 1 44 27 27 97 - Fax : +33 1 44 27 38 45 Acces : Tour 13, Couloir 12-13, 3eme etage - Bureau 316

2 2

A question about the DMRG of alps
by 郭怀明 06 Oct '15

06 Oct '15

Dear all, Recently I use dmrg of alps. I run three dmrg programs simultaneously and sometimes I get the following information, which stop the program. I raise a question: since the dmrg write files while it is running, is there a problem when several programs run? These program may wirte the same files to the same folder, which will cause a problem. (I use the alps in windows 7). Thank you very much. Sincerely, Huaiming Guo The group '/spectrum' does not exist. Traceback (most recent call last): File "C:\Program Files\VisTrails\Python27_64\lib\site-packages\pyalps\load.py" , line 269, in ReadDiagDataFromFile if 'results' in self.h5f.list_children(respath): PathNotFound: The group '/spectrum' does not exist.

3 6

Defining lattice using graph
by Huaiming Guo 02 Oct '15

02 Oct '15

Dear All, I want to ask how to use graph to define the lattice. I write a graph file named as "a.xml". In the parameter file, I use: GRAPH="C:\alpsproject\lattice\a.xml"; . The a.xml file is as the following. However the output is "could not find". I want to ask how to use it. Thank you very much. Sincerely, Huaiming Guo <GRAPH dimension="2" vertices="8" edges="12"> <VERTEX id="1" type="0"><COORDINATE>0 0</COORDINATE></VERTEX> <VERTEX id="2" type="0"><COORDINATE>0 1</COORDINATE></VERTEX> <VERTEX id="3" type="0"><COORDINATE>1 0</COORDINATE></VERTEX> <VERTEX id="4" type="0"><COORDINATE>1 1</COORDINATE></VERTEX> <VERTEX id="5" type="0"><COORDINATE>2 0</COORDINATE></VERTEX> <VERTEX id="6" type="0"><COORDINATE>2 1</COORDINATE></VERTEX> <VERTEX id="7" type="0"><COORDINATE>3 0</COORDINATE></VERTEX> <VERTEX id="8" type="0"><COORDINATE>3 1</COORDINATE></VERTEX> <EDGE source="1" target="3" id="1" type="0" vector="1 0"/> <EDGE source="1" target="2" id="2" type="1" vector="0 1"/> <EDGE source="2" target="4" id="3" type="0" vector="1 0"/> <EDGE source="3" target="5" id="4" type="0" vector="1 0"/> <EDGE source="3" target="4" id="5" type="1" vector="0 1"/> <EDGE source="4" target="6" id="6" type="0" vector="1 0"/> <EDGE source="5" target="7" id="7" type="0" vector="1 0"/> <EDGE source="5" target="6" id="8" type="1" vector="0 1"/> <EDGE source="6" target="8" id="9" type="0" vector="1 0"/> <EDGE source="7" target="1" id="10" type="0" vector="1 0"/> <EDGE source="7" target="8" id="11" type="1" vector="0 1"/> <EDGE source="8" target="2" id="12" type="0" vector="1 0"/> </GRAPH>

2 2

Re: [ALPS-users] Complex terms in DMRG
by Michele Dolfi 02 Oct '15

02 Oct '15

Hi, The error that you see is not an ALPS problem, but an error from the Python interpreter saying that the variable “I” is not defined. You just have to define the parameter as a string: ‘t0' :-1.0*I, becomes: ‘t0’ : '-1.0*I', Michele -- ETH Zurich Michele Dolfi Institute for Theoretical Physics HIT G 32.4 Wolfgang-Pauli-Str. 27 8093 Zurich Switzerland dolfim(a)phys.ethz.ch www.itp.phys.ethz.ch +41 44 633 78 56 phone +41 44 633 11 15 fax > On Oct 2, 2015, at 10:42 AM, Huaiming Guo <hmguo(a)buaa.edu.cn> wrote: > > Dear Michele Dolfi, > I use DMRG of alps with complex value. I use Complex=1 in the command line and that works. However when I use the expression in the python way, it return that “I is not defined”. I wonder how to write it in the python way. Thank you very much. > > Sincerely, > Huaiming Guo > <parm.py>

1 0

Problem running ALPS executables in parallel over Mac OS X and Ubuntu
by Andre Vieira 30 Sep '15

30 Sep '15

Hello, I would be glad if ALPS users could help me with a problem I've been experiencing. I have ALPS 2.2.b3 installed in an iMac (i5 mid 2011) and an i7 desktop running Ubuntu 15.04 64 bits, with openmpi (1.6.5) also present. I have no problem running ALPS executables in either computer, nor running basic programs using openmpi (for instance `mpirun -host <host> hostname`) from either machine to the other. I can also successfully run ALPS executables in either machine using openmpi locally. However, when I try to run an ALPS executable (e.g. loop) via openmpi using both machines simultaneously, the jobs are initiated, keep running indefinitely, but no output is produced, not even to the screen. Using mpirun with the -v flag gives no additional information. All I get is the following: dirac-2:qmc apvieira$ mpirun -v -np 7 loop --mpi heisloop20.in.xml --Tmin 5 ALPS/looper version 3.2b12-20100128 (2010/01/28) multi-cluster quantum Monte Carlo algorithms for spin systems available from http://wistaria.comp-phys.org/alps-looper/ copyright (c) 1997-2010 by Synge Todo <wistaria(a)comp-phys.org> using ALPS/parapack scheduler a Monte Carlo scheduler for multiple-level parallelization copyright (c) 1997-2013 by Synge Todo <wistaria(a)comp-phys.org> based on the ALPS libraries version 2.2.b3 available from http://alps.comp-phys.org/ copyright (c) 1994-2013 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. Did anyone have this kind of problem before? Thanks, Andre Vieira

2 5