October 2016 - Comp-phys-alps-users

reduced density matrices eigenvalues
by Robertson Esperanza 03 Nov '16

03 Nov '16

Hello, Thank you for replying to my previous inquiries. For now, I'm looking at the von neumann entropy of some variation of the Hubbard Hamiltonian, and I'm using the mps_optim application. Since the von Neumann entropy can be calculated, will I be able to access the eigenvalues of the reduced density matrices? If yes, how can I do that? I would like to also do the same thing using the exact diagonalization, but there is no function, as far as I know, to calculate the non Neumann entropy. How can I access the eigenvalued of the reduced density matrices using the ALPS application for exact diagonalization? Regards, Robbie

3 7

DMRG code
by Julien 02 Nov '16

02 Nov '16

Hello ALPS community, First, sorry for my English and for the trivial question that I'm going to ask you. I already made the tutorial on the S=(1/2, 1) antiferromagnetic Heisenberg chain where I found the correlation function and the gap in the thermodynamic limit (for S = 1). Now, I'm trying to do the same thing for the S=1/2 ferromagnetic Heisenberg chain but I guess I made a mistake when we have to specify the different parameters. Here is what I wrote in my computer code : parms.append({ 'LATTICE' : 'open chain lattice', 'MODEL' : 'spin', 'CONSERVED_QUANTUMNUMBERS' : 'N,Sz', 'Sz_total' : 0, 'J' : -1, 'SWEEPS' : 4, 'NUMBER_EIGENVALUES' : 1, 'L' : 32, 'MAXSTATES' : 20, 'MEASURE_CORRELATIONS[Diagonal spin correlations]' : 'Sz', }) but the result is completely wrong. I only modify the sign of the coupling constant compared to the antiferromagnetic case. Thank you in advance for your help, Best regards, Julien

4 11

MPI problem
by Tadeusz Wasiutyński 31 Oct '16

31 Oct '16

Hello, I met the problem on my CENTOS 7 installed on station with 24 threads. After: ~/opt/bin/cmake -D Boost_ROOT_DIR:PATH=~/alps-2.2.b4-src-with-boost/boost/ -D MPI_C_LIBRARIES=/usr/lib64/openmpi/libmpi.so -D MPI_C_INCLUDE_PATH=/usr/lib64/openmpi/include/ -D MPI_CXX_INCLUDE_PATH=/usr/lib64/openmpi/include -D LPSolve_LIBRARY=/usr/lib64/liblpsolve55.so -D LPSolve_INCLUDE_DIR=/usr/include/lpsolve/ -D HDF5_LIBRARIES=~/opt/lib/libhdf5.so -D HDF5_INCLUDE_DIR=~/opt/include/ -D SZIP_LIBRARIES=/usr/lib64/lib/libsz.so -D SZIP_INCLUDE_DIRS=/usr/lib64/include/ ~/alps-2.2.b4-src-with-boost/alps/ I receive: ---- Compiler version: c++ (GCC) 4.8.5 20150623 (Red Hat 4.8.5-4) -- Build type: Release -- Python interpreter /usr/bin/python -- Python interpreter ok : version 2. PYTHON_INCLUDE_DIRS = /usr/include/python2.7 -- PYTHON_NUMPY_INCLUDE_DIR = /usr/lib64/python2.7/site-packages/numpy/core/include -- PYTHON_SITE_PKG = /usr/lib/python2.7/site-packages -- PYTHON_LIBRARY = /usr/lib64/python2.7/config/libpython2.7.so -- PYTHON_EXTRA_LIBS =-lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -- ALPS version: 2.2.b4 -- Looking for Boost Source -- Found Boost Source: /home/twasiutynsk/alps-2.2.b4-src-with-boost/boost -- Boost Version: 1_58_0 -- Adding Boost dir: /home/twasiutynsk/alps-2.2.b4-src-with-boost/boost -- MPI compiler was /usr/lib64/openmpi/bin/mpicxx -- Falling back to CMake provied LAPACK/BLAS detection. -- A library with BLAS API found. -- A library with BLAS API found. -- A library with LAPACK API found. -- SQLite Library: not found -- Could NOT find SZIP (missing: SZIP_LIBRARIES SZIP_INCLUDE_DIRS) -- HDF5 without THREADSAFE mode. ALPS will ensure thread safety by HDF5 running sequentially. -- Python interpreter /usr/bin/python -- Python interpreter ok : version 2.7.5 -- PYTHON_INCLUDE_DIRS = /usr/include/python2.7 -- PYTHON_NUMPY_INCLUDE_DIR = /usr/lib64/python2.7/site-packages/numpy/core/include -- PYTHON_SITE_PKG = /usr/lib/python2.7/site-packages -- PYTHON_LIBRARY = /usr/lib64/python2.7/config/libpython2.7.so -- PYTHON_EXTRA_LIBS =-lpthread -ldl -lutil -- PYTHON_LINK_FOR_SHARED = -Xlinker -export-dynamic -- Numpy include in /usr/lib64/python2.7/site-packages/numpy/core/include -- ALPS XML dir is /opt/alps/lib/xml -- HDF5 without THREADSAFE mode. ALPS will ensure thread safety by HDF5 running sequentially. -- MPS: enabling NU1 symmetry. -- HDF5 without THREADSAFE mode. ALPS will ensure thread safety by HDF5 running sequentially. -- tebd will not be built -- HDF5 without THREADSAFE mode. ALPS will ensure thread safety by HDF5 running sequentially. -- Configuring done -- Generating done -- Build files have been written to: /home/twasiutynsk/build in ccmake I see MPI ON. make and make tests (all passed) go smoothly. In tutorial runs I see however mess with MPI=24 while everything goes OK with MPI=1. Is something wrong in my HDF5? In docker hub I found dolfim/alps but could not run probably because of some path problems. Anyone did it successfully? Regards -- Tadeusz Wasiutyński -- Tadeusz Wasiutyński

2 1

ALPS install
by Safir 14 Oct '16

14 Oct '16

*Sir, I tried to install Alps from source code. But while running make install it giving errors. What are correct steps to install alps in ubuntu 16.04, whether VisTrails is necessary for alps?. If So how we can install Vistrails. How we can check Alps is installed or not. I am attaching screen shot of terminal* * I look forward to hearing from you.* -- *Many Thanks* *----------------------------------------------------------------------------------------------------------------------------------------------------*Safir TK Research Scholar Theoretical Physics Lab(TPL) Department Of physics National Institute of Technology Karnataka (NITK) Surathkal, P.O: Srinivasnagar, Mangalore, India - 575 025 Tel: +91 824 2473289 (Office) Mob: +91 - 7411797881, +91 -9946635799 Email : safirmon(a)gmail.com

3 9

EBH in DWA
by Mateusz Łącki 13 Oct '16

13 Oct '16

Dear All, could I ask if the DWA code actually supports the extended BH model (as defined in “boson Hubbard” hamiltonian in a standard issue of models.xml)? the results I obtain seem not to depend at all on the value of the V parameter. With sparsediag I observe no such behaviour. thebondterm contains “hopping” (proportional to t#) and interaction (proportional to V#). If I wanted to split (as the hopping is treated separately in worm-type algorithms I guess…) it into two bond-terms it seems that only the last bondterm is taken into account: this results in t=0 result <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2 </SITETERM> <BONDTERM source="i" target="j"> <PARAMETER name="t#" default="0"/> -t#*(bdag(i)*b(j)+bdag(j)*b(i)) </BONDTERM> <BONDTERM source="i" target="j"> <PARAMETER name="V#" default="0"/> -V#*(n(j)*n(i)) </BONDTERM> this in V=0 result: <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2 </SITETERM> <BONDTERM source="i" target="j"> <PARAMETER name="V#" default="0"/> -V#*(n(j)*n(i)) </BONDTERM> <BONDTERM source="i" target="j"> <PARAMETER name="t#" default="0"/> -t#*(bdag(i)*b(j)+bdag(j)*b(i)) </BONDTERM> Best, Mateusz Łącki

3 6