February 2011 - Comp-phys-alps-users

the test is killed without any prompt message
by Chen Ahai 01 Mar '11

01 Mar '11

dear all, I run a DMRG code to calculate the correlation function with about 50 sweeps and 500 maxstates. Some tests is killed without any prompt message, while the others are successful. I can also get the results with lower sweeps and maxstates. I wonder if there are some ways to solve this problem or just a better computer. . thank you very much -- -------------------------------------------------------------------------- Chen Ahai Department of Physics, Zhejiang Normal University Jinhua, Zhejiang Province, 321004, China e-mail: chenahaiphysics(a)gmail.com; chenahai(a)yahoo.cn

2 2

Re: [ALPS-users] continue or extend a run using directloop_sse
by Kubila 27 Feb '11

27 Feb '11

Dear Prof. Troyer Thank you very much for your response. When I resubmit the run, I always got a message saying " segmentation fault ". Followings are some relevant files. Sincerely, Kubila The submit.file (qsub submit.file) #!/bin/sh #$ -V #$ -N square_2.5196_L24 #$ -S /bin/sh #$ -cwd #$ -j y #echo $SGE_TASK_ID /home/kubila/ALPS/bin/parameter2xml -f square_2.5196_L24 /home/kubila/ALPS/bin/dirloop_sse --write-xml square_2.5196_L24.in.xml The resubmit.file (qsub resubmit.file) #!/bin/sh #$ -V #$ -N square_2.5196_L24 #$ -S /bin/sh #$ -cwd #$ -j y #echo $SGE_TASK_ID /home/kubila/ALPS/bin/dirloop_sse --write-xml square_2.5196_L24.out.xml The square_2.5196_L24.out.xml file <?xml version="1.0" encoding="UTF-8"?> <?xml-stylesheet type="text/xsl" href="ALPS.xsl"?> <JOB xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/job.xsd"> <TASK status="running"> <INPUT file="square_2.5196_L24.task1.out.xml"/> </TASK> </JOB> The parameter file LATTICE="staggered anisotropic square lattice" MODEL = "anisotropic_spin" LATTICE_LIBRARY="./../lattices1.xml" MODEL_LIBRARY="./../models.xml" local_S=1/2 Q=1.0 Q0=2.5196 beta = 24.0 T = 1/beta EPSILON = 0.0001 THERMALIZATION=10000 SWEEPS=450000 {L=24;W=24;} --- On Fri, 2/18/11, Matthias Troyer <troyer(a)phys.ethz.ch> wrote: From: Matthias Troyer <troyer(a)phys.ethz.ch> Subject: Re: [ALPS-users] continue or extend a run using directloop_sse To: comp-phys-alps-users(a)lists.phys.ethz.ch Date: Friday, February 18, 2011, 12:20 AM Hi Kubila, can you send the input file and tell me excatly what you did? Is it the chaining of jobs that fails? And if so, what is the problem? Matthias On Feb 17, 2011, at 4:42 PM, Kubila wrote: Hi ALPS 2.0 users I was wondering did anyone succeed in extending or continuing a run using directloop_sse in ALPS 2.0? Using job.out.xml seems to work with loop (./loop job.out.xml works), but not with dirctloop_sse (./directloop_sse job.out.xml fails).... Any help will appreciate very much! Sincerely, Kubila

2 3

TEBD application: How to access the Schmidt eigenvalues ?
by Grégoire Misguich 24 Feb '11

24 Feb '11

Dear alps users, I would like to know if one can simply have access the to the Schmidt eigenvalues which describe the wave-function in (or at the end of) the TEBD application. It there a "measurement" for that ? Thank you ! Gregoire Misguich

1 0

continue or extend a run using directloop_sse
by Kubila 18 Feb '11

18 Feb '11

Hi ALPS 2.0 users I was wondering did anyone succeed in extending or continuing a run using directloop_sse in ALPS 2.0? Using job.out.xml seems to work with loop (./loop job.out.xml works), but not with dirctloop_sse (./directloop_sse job.out.xml fails).... Any help will appreciate very much! Sincerely, Kubila

2 1

quantum XY model with an uniform magnetic field in the x-direction
by Kubila 11 Feb '11

11 Feb '11

Hi ALPS users I was wondering can one use looper to simulate quantum XY model with an uniform magnetic field in the x-direction. Thanks.

3 6

TEBD tutorial "2a" (kink & XX chain)
by Grégoire Misguich 09 Feb '11

09 Feb '11

Dear ALPS community, As my first attempt to use the TEBD code provided in alps-2.0, I ran a tutorial concerning the simulation of the XX spin chain (tutorial_2a.py), which simulates the dynamics starting from a domain wall spin configuration: http://alps.comp-phys.org/mediawiki/index.php/ALPS_2_Tutorials:TEBD-02_kink I got the following figures (attached .png files). The Wiki says "We see that the magnetization agrees very well to visual accuracy", but in the plot I obtained, it is not quite the case (very "noisy"). I played a bit with some parameters (increasing the number of time steps and the max Schmidt dimension), but it still get some ugly curves. So, I suspect that something might we wrong... Any help or suggestion would be very much appreciated ! Gregoire Misguich.

1 1

RESTRICT_MEASUREMENTS[Sz] dirloop_sse alps 1.35
by Marcelo Guimarães 08 Feb '11

08 Feb '11

Hello, I would like to know if the parameter RESTRICT_MEASUREMENTS[Sz] works for dirloop_sse as indicated at http://alps.comp-phys.org/mediawiki/index.php/Documentation:dirloop_sse. I'm trying to obtain the gap for (3d), H=J( S+(i)S(i+1)- + S-(i)S+(i+1) ) +D(Sz)^2 at D > Dc by calculating energy difference from E(Mz=1) , E(Mz=0) at low temperature. Since <Sz(i)Sz(i+1)> ~ 0 , and dirloop does not measure off-diagonal correlations , I was wandering if I could use this procedure in dirloop to estimate the gap. Thanks, Marcelo Guimaraes.

2 1

Re: [ALPS-users] Extracting dates
by Xiao Quan Yu 07 Feb '11

07 Feb '11

Dear Prof. Troyer, Thank you for your replying. One of my out put file is the following: <?xml version="1.0" encoding="UTF-8"?> <?xml-stylesheet type="text/xsl" href="file:///u/he/warter/ALPS/lib/xml/ALPS.xsl"?> <SIMULATION xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2002/10/ALPS.xsd"> <PARAMETERS> <PARAMETER name="LATTICE_LIBRARY">lattices.xml</PARAMETER> <PARAMETER name="LATTICE">inhomogeneous simple cubic lattice</PARAMETER> <PARAMETER name="L">5</PARAMETER> <PARAMETER name="MODEL_LIBRARY">models.xml</PARAMETER> <PARAMETER name="MODEL">hardcore boson</PARAMETER> <PARAMETER name="t">1.0</PARAMETER> <PARAMETER name="T">0.1</PARAMETER> <PARAMETER name="delta">10.0</PARAMETER> <PARAMETER name="SWEEPS">500000</PARAMETER> <PARAMETER name="THERMALIZATION">10000</PARAMETER> <PARAMETER name="EPSILON">0.25</PARAMETER> <PARAMETER name="DISORDERSEED">65897</PARAMETER> <PARAMETER name="mu">delta*2*(random()-0.5)</PARAMETER> <PARAMETER name="SEED">668368703</PARAMETER> </PARAMETERS> <AVERAGES> <SCALAR_AVERAGE name="Density"> <COUNT>497664</COUNT> <MEAN method="jackknife">0.52646836</MEAN> <ERROR converged="yes" method="jackknife">2.09e-05</ERROR> <VARIANCE method="simple">5.45e-05</VARIANCE> <AUTOCORR method="jackknife">1.5</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Density^2"> <COUNT>497664</COUNT> <MEAN method="jackknife">34.652926</MEAN> <ERROR converged="yes" method="jackknife">0.00276</ERROR> <VARIANCE method="simple">0.944</VARIANCE> <AUTOCORR method="jackknife">1.5</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Energy"> <COUNT>497664</COUNT> <MEAN method="jackknife">-251.9492</MEAN> <ERROR converged="yes" method="jackknife">0.0221</ERROR> <VARIANCE method="simple">64.1</VARIANCE> <AUTOCORR method="jackknife">1.4</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Energy Density"> <COUNT>497664</COUNT> <MEAN method="jackknife">-2.0155936</MEAN> <ERROR converged="yes" method="jackknife">0.000177</ERROR> <VARIANCE method="simple">0.0041</VARIANCE> <AUTOCORR method="jackknife">1.4</AUTOCORR> </SCALAR_AVERAGE> <SCALAR_AVERAGE name="Stiffness"> <COUNT>497664</COUNT> <MEAN method="jackknife">2.17262</MEAN> <ERROR converged="yes" method="jackknife">0.00305</ERROR> <VARIANCE method="simple">3.16</VARIANCE> 1,1 Top ..................... In the mean time, I will try to use ALPS 2.0. Thank you. Best, Xiaoquan > Today's Topics: > > 1. Re: Extracting dates (Matthias Troyer) > 2. Re: Dzyaloshinskii-Moriya in spinmc (Matthias Troyer) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 6 Feb 2011 23:03:46 +0200 > From: Matthias Troyer <troyer(a)phys.ethz.ch> > Subject: Re: [ALPS-users] Extracting dates > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: <ACE5D63D-7EEA-43D9-BDAC-862222FFEBD5(a)phys.ethz.ch> > Content-Type: text/plain; charset=us-ascii > > Can you send your output files? The problem is probably that mu is not > just a number. You might also look into using ALPS 2 instead which has > much better plotting and evaluation tools. > > Matthias > > On Feb 4, 2011, at 12:33 PM, Xiao Quan Yu wrote: > >> Dear ALPS community, >> >> I met a trouble for extracting the calculated results from the >> XML output files. My setting is the following: >> >> plot2.xml file: >> >> <?xml version="1.0" encoding="UTF-8"?> >> >> <plot name="Stiffness versus disorder for three-dimesnional >> HardcoreBoson >> models"> >> >> <legend show="true"/> >> <xaxis label="Disorder" type="PARAMETER" name="mu"/> >> <yaxis label="Stiffness" type="SCALAR_AVERAGE"/> >> >> <set label="Inhomogeneous simple cubic lattice"/> >> >> </plot> >> >> and my "input" file is >> >> LATTICE_LIBRARY="lattices.xml"; >> LATTICE="inhomogeneous simple cubic lattice"; >> L=5; >> MODEL_LIBRARY="models.xml"; >> MODEL="hardcore boson"; >> t = 1.0; >> T = 0.1; >> delta=10.0; >> SWEEPS=500000; >> THERMALIZATION=10000; >> EPSILON=0.25; >> { DISORDERSEED = 34275; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 49802; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 82529; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 13564; mu=delta*2*(random()-0.5); } >> { DISORDERSEED = 65897; mu=delta*2*(random()-0.5); }. >> >> When I typed the command "extracttext plot2.xml input.task*.out.xml", >> the >> following thing happend: >> >> >> plot2.xml:4: parser error : XML declaration allowed only at the start of >> the document >> <?xml version="1.0" encoding="UTF-8"?> >> ^ >> unable to parse plot2.xml >> /tmp/archive2plot3652.xsl:1: parser error : Document is empty >> >> ^ >> /tmp/archive2plot3652.xsl:1: parser error : Start tag expected, '<' not >> found >> >> ^ >> cannot parse /tmp/archive2plot3652.xsl >> /tmp/plot3652:1: parser error : Document is empty >> >> ^ >> /tmp/plot3652:1: parser error : Start tag expected, '<' not found >> >> ^ >> unable to parse /tmp/plot3652 >> >> >> I don't know what is going on here, could you help me figure it out? >> Thank >> you in advance! >> >> >> Best, >> Xiaoquan >> >> >> >> ---------------------------------------------------------------- >> SISSA Webmail https://webmail.sissa.it/ >> Powered by SquirrelMail http://www.squirrelmail.org/ >> > > > > ------------------------------ > > Message: 2 > Date: Sun, 6 Feb 2011 23:11:48 +0200 > From: Matthias Troyer <troyer(a)phys.ethz.ch> > Subject: Re: [ALPS-users] Dzyaloshinskii-Moriya in spinmc > To: comp-phys-alps-users(a)lists.phys.ethz.ch > Message-ID: <574253B4-92C0-4A0F-97AB-0ADF062B8A7A(a)phys.ethz.ch> > Content-Type: text/plain; charset=iso-8859-1 > > Dear Francesco, > > The source code needs to be modified, but it is very easy to write such a > code. Anyone in my group at ETH can teach you how to do that. I'm on > vacation and without good internet this week but we can meet next week, or > you can ask others on my group whether they can help you. Writing a local > update code for your model should not take more than half a day at the > most. > > Matthias > > > > > On Feb 2, 2011, at 6:32 PM, Francesco Casola wrote: > >> Dear ALPS community, >> >> could you please let me know about the possibility of including >> Dzyaloshinskii-Moriya in spinmc, possibly through the presence of an >> off-diagonal exchange coupling constant? >> >> Is it necessary to modify the source code for such a case, or can I >> manage it in some other way not explicitly mentioned in the ALPS >> tutorial? >> >> Best wishes >> Francesco >> >> -- >> Francesco Casola, Ph. D. Student >> Laboratory for Solid State Physics >> ETH Z?rich H?nggerberg >> HPF E15.1 Schafmattstrasse 16 >> CH-8093 Switzerland >> >> Tel. Office: +41 (0)44 633 22 96 / 03 >> Tel. Lab: +41 (0)44 633 23 01 >> E-mail: fcasola(a)phys.ethz.ch >> >> Group website: https://protos.ethz.ch/dscm >> > > > > End of Comp-phys-alps-users Digest, Vol 59, Issue 4 > *************************************************** > ---------------------------------------------------------------- SISSA Webmail https://webmail.sissa.it/ Powered by SquirrelMail http://www.squirrelmail.org/

2 1

Dzyaloshinskii-Moriya in spinmc
by Francesco Casola 06 Feb '11

06 Feb '11

Dear ALPS community, could you please let me know about the possibility of including Dzyaloshinskii-Moriya in spinmc, possibly through the presence of an off-diagonal exchange coupling constant? Is it necessary to modify the source code for such a case, or can I manage it in some other way not explicitly mentioned in the ALPS tutorial? Best wishes Francesco -- Francesco Casola, Ph. D. Student Laboratory for Solid State Physics ETH Zürich Hönggerberg HPF E15.1 Schafmattstrasse 16 CH-8093 Switzerland Tel. Office: +41 (0)44 633 22 96 / 03 Tel. Lab: +41 (0)44 633 23 01 E-mail: fcasola(a)phys.ethz.ch Group website: https://protos.ethz.ch/dscm

2 1

Extracting dates
by Xiao Quan Yu 06 Feb '11

06 Feb '11

Dear ALPS community, I met a trouble for extracting the calculated results from the XML output files. My setting is the following: plot2.xml file: <?xml version="1.0" encoding="UTF-8"?> <plot name="Stiffness versus disorder for three-dimesnional HardcoreBoson models"> <legend show="true"/> <xaxis label="Disorder" type="PARAMETER" name="mu"/> <yaxis label="Stiffness" type="SCALAR_AVERAGE"/> <set label="Inhomogeneous simple cubic lattice"/> </plot> and my "input" file is LATTICE_LIBRARY="lattices.xml"; LATTICE="inhomogeneous simple cubic lattice"; L=5; MODEL_LIBRARY="models.xml"; MODEL="hardcore boson"; t = 1.0; T = 0.1; delta=10.0; SWEEPS=500000; THERMALIZATION=10000; EPSILON=0.25; { DISORDERSEED = 34275; mu=delta*2*(random()-0.5); } { DISORDERSEED = 49802; mu=delta*2*(random()-0.5); } { DISORDERSEED = 82529; mu=delta*2*(random()-0.5); } { DISORDERSEED = 13564; mu=delta*2*(random()-0.5); } { DISORDERSEED = 65897; mu=delta*2*(random()-0.5); }. When I typed the command "extracttext plot2.xml input.task*.out.xml", the following thing happend: plot2.xml:4: parser error : XML declaration allowed only at the start of the document <?xml version="1.0" encoding="UTF-8"?> ^ unable to parse plot2.xml /tmp/archive2plot3652.xsl:1: parser error : Document is empty ^ /tmp/archive2plot3652.xsl:1: parser error : Start tag expected, '<' not found ^ cannot parse /tmp/archive2plot3652.xsl /tmp/plot3652:1: parser error : Document is empty ^ /tmp/plot3652:1: parser error : Start tag expected, '<' not found ^ unable to parse /tmp/plot3652 I don't know what is going on here, could you help me figure it out? Thank you in advance! Best, Xiaoquan ---------------------------------------------------------------- SISSA Webmail https://webmail.sissa.it/ Powered by SquirrelMail http://www.squirrelmail.org/

2 1