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

mps_evolve to calculate the time evolution of Fermi-Hubbard
by xiaofanlaohu 14 Jun '16

14 Jun '16

To whom it may concern, When I use the mps_evolve to calculate the time evolution of Fermi-Hubbard, there are some problems that I can not solve. The python file and the error information are shown as follow. best, Xiaofan Zhou ------------------------------------------------------------------------------------------ The python file is: import pyalps import matplotlib.pyplot as plt import pyalps.plot import numpy as np from collections import OrderedDict from copy import deepcopy basename = 'fermi' ## common model parameters model = OrderedDict() model['LATTICE_LIBRARY' ] = 'lattices.xml' model['LATTICE'] = 'open chain lattice' model['L' ] = 20 model['MODEL_LIBRARY' ] = 'models.xml' model['MODEL' ] = 'fermion Hubbard' model['CONSERVED_QUANTUMNUMBERS'] = 'Nup,Ndown' model['Nup_total' ] = 5 model['Ndown_total' ] = 5 model['t' ] = 1 model['U' ] = 1 model['MAXSTATES'] = 100 model['MEASURE_LOCAL[Local Nup]'] = 'n_up' model['MEASURE_LOCAL[Local Ndown]'] = 'n_down' #prepare the input parameters parms = [] for tau in [20.0]: ns = 50 dt = tau / ns p = deepcopy(model) p['init_state' ] = 'local_quantumnumbers' p['initial_local_Nup'] = ','.join((['0.5']+ ['0.0'])* 10) p['initial_local_Ndown'] = ','.join((['0.5']+ ['0.0'])* 10) p['te_order' ] = 'second' p['DT' ] = dt p['TIMESTEPS'] = ns p['tau' ] = tau # not used in the simulation, but useful in the evaluation below p['ALWAYS_MEASURE'] = 'Local Nup' p['chkp_each' ] = ns p['measure_each' ] = 5 p['COMPLEX' ] = 1 parms.append(p) ## write input files and run application input_file = pyalps.writeInputFiles(basename, parms) res = pyalps.runApplication('mps_evolve', input_file) -------------------------------------------------------------------------------------------------------- The error is, ALPS/MPS Time Evolution version MAQUIS repo r4143 (2013-2014) Density Matrix Renormalization Group algorithm available from http://alps.comp-phys.org/ copyright (c) 2013 Institute for Theoretical Physics, ETH Zurich copyright (c) 2010-2011 by Bela Bauer copyright (c) 2011-2013 by Michele Dolfi for details see the publication: M. Dolfi et al., Computer Physics Communications 185, 3430 (2014). doi: 10.1016/j.cpc.2014.08.019 based on the ALPS libraries version 2.2.b5 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. Running task 1. This binary contains symmetries: nu1 MAQUIS repo r4143 Temporary storage is disabled site_types: 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 phys[0]: |( <2,2,0,0,0,0>: 1 )( <2,0,0,0,0,0>: 1 )( <0,2,0,0,0,0>: 1 )( <0,0,0,0,0,0>: 1 )| Abort by signal 11 In /home/sxdxcg/alps-2.2.b5-src-with-boost/alps/src/alps/ngs/lib/signal.cpp on 104 in segfault /home/sxdxcg/alps_b5/lib/libalps.so.2(_ZN4alps3ngs6signal8segfaultEi+0x3b) [0x7f8858184d6b] /lib/x86_64-linux-gnu/libpthread.so.0(+0xfcb0) [0x7f88567cfcb0] mps_evolve(_ZNK9ALPSModelIN4alps7numeric6matrixISt7complexIdESt6vectorIS4_SaIS4_EEEE12NU1_templateILi6EiEE11initializerERK7LatticeR14BaseParameters+0x7fb) [0x5e5c5b] mps_evolve(_ZN3simIN4alps7numeric6matrixISt7complexIdESt6vectorIS4_SaIS4_EEEE12NU1_templateILi6EiEEC1ERK14DmrgParameters+0xf63) [0x4fe1f3] mps_evolve(_ZN9tevol_simIN4alps7numeric6matrixISt7complexIdESt6vectorIS4_SaIS4_EEEE12NU1_templateILi6EiE25nearest_neighbors_evolverIS8_SA_EEC1ERK14DmrgParametersb+0x26) [0x4fe9a6] mps_evolve(_ZN10simulationI12NU1_templateILi6EiEE3runER14DmrgParametersb+0x32c) [0x4ff7cc] mps_evolve(_Z7run_simRKN5boost10filesystem4pathES3_bd+0x1b20) [0x49bf20] mps_evolve(_ZN9Scheduler3runEv+0x1d8) [0x55cef8] mps_evolve(main+0xbf) [0x49667f] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f8855c1176d] mps_evolve() [0x4973d5] Aborted (core dumped) 2016-06-14 xiaofanlaohu

2 1

String order parameters measurements
by 梁霄 14 Jun '16

14 Jun '16

ALPS users and developers, I am wondering in ALPS how to construct the measurement of String Order Parameter in paper: 2012 New J. Phys. 14 065012. I have tried to measure the operators independently but the result was wrong. Xiao Liang

2 4

compile error, svn version 7708
by Amir M.Aghaei 08 Jun '16

08 Jun '16

Hi, I tried compiling the latest svn version of ALPS (7708) last night and I got the following compile error: */alps/applications/dmrg/dmrg/dmtk/system.h:594:48: error: lvalue required as unary ‘&’ operand s.write((const char *)&_target_weight[i], sizeof(double));* I looked it up and sounds like this line of code was written back in 2006. So that can be a compiler issue but I didn't find anything that I could relate to the problem online. I am using gcc 6.1.1. Bests, Amir

1 1

a question of external potential in MPS
by 北尘心 28 May '16

28 May '16

Previously, I used DMRG to calculate the model containing the external potential, but now I would like to try to use MPS to calculate it.When i use MPS, I meet a question of potential. My potential setting like this <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> <PARAMETER name="Vout" default="0"/> -mu#*n(i)+U#*n_up(i)*n_down(i)+Vout*((x-(L-1)/2.0)^2)*n(i) </SITETERM> I use the model ,lattice and command like I used in DMRG, but the computer prompt like : Exception thrown: Cannot evaluate x By using DMRG, it can be calculated correctly, but MPS will prompt the error, i wonder what should I do to correct it? It is highly grateful for your any suggestion. Yours sincerely, Xingbo Wei Zhejiang Normal University Jinhua,230014, China

2 1

infinite lattice size and initial states definitions.
by 梁霄 28 May '16

28 May '16

ALPS users and developers, I have two questions about the definitions of lattice sizes and the initial states for mps_optim. 1, how to create a infinite size of lattice? And how to define filling of the lattice for DMRG or MPS? 2, if the lattice size is finite and is very large, how to quickly set the quantumnumbers of the initial state for all sites? Thanks. Xiao Liang

2 1

strange chemical potential of extended boson hubbard model
by 梁霄 27 May '16

27 May '16

ALPS users and developers, I am trying to calculate the phase diagram of extended boson hubbard using dmrg. The lattice is the finite one-dimension lattice with open boundaries. The chemical potentials are: u-=E(N)-E(N-1), u+=E(N+1)-E(N). E(N) is the ground energy for a total particle number N. However I found that there is something unusual when calculating the chemical potential: E(N+1)-E(N). When V is not zero, the u+ calculated by DMRG, sparsediag and mps_optim is less than the expected value. For example when L=100, N=50, t=0 and V=0.4, the correct value of u+ is 0.8, however using dmrg(mps and sparsediag) the value is 0.4 . I can obtain the correct value using u+=E(N+2)-E(N+1) when choosing L=41 and N=21 instead. However within L=101 and N=101, even choose u+=E(N+2)-E(N+1) the result is still 0.4. It seems that when V does not equal to zero the calculated chemical potential is related to specific total particle numbers. Xiao Liang

2 2

Help with install
by Thornton, Curtislee 26 May '16

26 May '16

Hello my name is Curtis, (My setup: Windows 10 pro 64bit, Python(x,y) which is python 2.7 32bit) I am trying to install ALPS and I am running into the problem with the PATH being too long. I read everything on the website and searched google and youtube for anything I can find but still nothing helps. I tried version 2.1 and 2.2 (32 and 64 bit versions) binary installer. I tried installing it to PATH from the beginning but get the error, I then tried installing it and adding to the PATH (through the windows setting) later (I tried adding every folder to the PATH seperatly as well) and still does not 'see' it when I try to "import pyalps" and tells me there is no module with that name. I also tried setting the path inside the program 'Spyder' (IDE) which I run python with. I also tried to copy and paste the files from the ALPS install into the C:\python27 folders that hold the same file types (Dlls, Lib, libs, and include). When I copy them over I then tried to "import pyalps" again and it finds the module but then throws an error saying that some other module was not found. Any help with this will be great. Thanks for your time.

2 1

Measuring a sum (or sum squared) of site operators
by Larry Engelhardt 14 May '16

14 May '16

Is there a way to directly measure sums of site operators and squares of sums of site operators? (I am using the directed loop SSE code, but I think this question also applied to other algorithms.) I could accomplish this indirectly by measuring the local site operators and also measuring correlations between the site operators and then adding them together (to manually expand the square), but if there is a way to do this directly, that would seem to be a more elegant/pretty solution. The motivation for this questions is that I need to compute the “Magnetization" and “Magnetization Squared” — which already exists as pre-created measurements — but my system has site dependent g values, so instead of measuring the sum of Sz_i (the square of the sum), I need to measure the sum of g_i***Sz_i (and the square of the sum), where the sum extends over all sites, i. Is this possible? If so, can you provide an example? Thanks! Larry

1 0

Constraints from two types of sitebasis
by 梁霄 08 May '16

08 May '16

Dear all, Recently I am using ALPS simulating JCH model. To construct the basis for atom-light interaction, two types of basis are defined on each site. My question is how to set a constraint of the summation from two types of basis? For example, <BASIS> <SITEBASIS type="0" ref="M"/> <SITEBASIS type="1" ref="N"/> <CONSTRAINT quanumnumber="A+B" value="C"/> </BASIS> A and B are from sitebasis M and N. Does the method above work? Thanks.

2 1

Define a J-C model
by 梁霄 06 May '16

06 May '16

Dear all, In alps models all the basis are defined on site positions. When there are interactions between atoms and photons, an atom on each site can absorb or emit photons, how to define such basis for photons? Thanks.

3 5