Hi,

DWA will only work with bosonic Hubbard models. For spin chains use the loop code when you have no field, and the directed loop code when you have a field

Matthias

On 22 Oct 2014, at 07:55, Rene John Kerkdyk rene.kerkdyk@theorie.physik.uni-goettingen.de wrote:

Dear ALPS developers and users,

while using the dwa application I stumbled upon a segmentation fault. This happens when I try to simulate a S=1/2 Heisenberg chain. Below is the input file "dwa.parm" used and the output generated by ALPS. Is it not possible to calculate this spin model with the dwa application?

Best, René

Reference: I am using the alps-2.2.b3-r7462-src to build ALPS and work on the following operating system: Linux version 2.6.32-431.20.5.el6.x86_64 (mockbuild@c6b8.bsys.dev.centos.org) (gcc version 4.4.7 20120313 (Red Hat 4.4.7-4) (GCC) ) #1 SMP Fri Jul 25 08:34:44 UTC 2014 CentOS release 6.5 (Final) Kernel \r on an \m

Input & Output: -bash-4.1$ cat dwa.parm LATTICE="chain lattice"; MODEL="spin"; LATTICE_LIBRARY="./lattices.xml" MODEL_LIBRARY="./models.xml" local_S = 1/2 J = 1 h=0 L=4; T = 0.001; SWEEPS=5000000; THERMALIZATION=100000; SKIP=500; MEASURE[Winding Number]=1 { t=0.01; }

-bash-4.1$ parameter2xml dwa.parm Converting parameter file dwa.parm to dwa.parm.in.xml -bash-4.1$ dwa --write-xml dwa.parm.in.xml /***************************************************************************** /SNIP/ *****************************************************************************/

using the ALPS parallelizing scheduler copyright (c) 1994-2006 by Matthias Troyer troyer@comp-phys.org. see Lecture Notes in Computer Science, Vol. 1505, p. 191 (1998).

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.

parsing task files ...

Initialization stage 1 (Site states) ... starting ... finished.

Initialization stage 2 (Hamiltonian) ... starting i. State minimum/maximum ... starting ... done. i. Onsite matrix ... starting a. Step 1: partial evaluate everything up to x,y,z terms ... starting ... done. b. Step 2: Initialize onsite matrix ... starting: ... done. ... done. iii. Bond strength matrix ... starting ... done. ... finished.

Initialization stage 3 (Lookups) ... starting ... finished.

Initialization stage 4 (Measurements) ... starting ... finished. Parameters ==========

LATTICE = "chain lattice"; MODEL = spin; LATTICE_LIBRARY = ./lattices.xml; MODEL_LIBRARY = ./models.xml; local_S = 1/2; J = 1; h = 0; L = 4; T = 0.001; SWEEPS = 5000000; THERMALIZATION = 100000; SKIP = 500; MEASURE[Winding Number] = 1; t = 0.01; SEED = 856848597;

Simulation

    Sweeps                              : 5000000
    Thermalization sweeps               : 100000
    Skip (Sweeps per measurement)       : 500
    Sweep counter                       : 0
    Sweep failure counter               : 0
    Propagation counter                 : 0
    Propagation failure counter         : 0

    Interested in time-of-flight images : 0
    Waist correction                    : 0

Lattice

Name : chain lattice Dimension : 1 Total Sites/Bonds : 4 / 4 Periodic : 1

Model

Name : spin Charge model : 0 Spin model : 1 Site Basis : spin states : -0.5 0.5

Measurements

MEASURE : 1 MEASURE[Simulation Speed] : 1 MEASURE[Total Particle Number^2] : 0 MEASURE[Energy^2] : 0 MEASURE[Density^2] : 0 MEASURE[Energy Density^2] : 0 MEASURE[Winding Number] : 1 MEASURE[Local Kink:Number] : 0 MEASURE[Local Density] : 0 MEASURE[Local Density^2] : 0 MEASURE[Green Function] : 0

Created run 1 locally Starting task 1. Segmentation fault

---

René John Kerkdyk Institute for Theoretical Physics University of Goettingen Friedrich-Hund-Platz 1 D-37077 Goettingen Germany Office: A4.126 Phone: +49 551 39 7687