Dear Synge,
thank you for your hints. It works nice.
I think, it would be useful to update the tutorial on that subject, to make
it more clear,
that parameter file is not the lattice description one.
I knew it was the two-dimensional model. I put it on purpose, just as an
example.
Can someone experienced in QMC-loop simulations comment on my bigger
problem, please.
>
> Message: 10
> Date: Mon, 29 Sep 2014 11:03:36 +0900
> From: Synge Todo <wistaria(a)comp-phys.org>
> To: comp-phys-alps-users(a)lists.phys.ethz.ch
> Subject: Re: [ALPS-users] complexity estimation for the 3D
> ferromagnetic Heisenberg model
> Message-ID: <9E0990A4-3582-45B9-B8C6-F704BE452A27(a)comp-phys.org>
> Content-Type: text/plain; charset=windows-1252
>
> Dear Oleh,
>
> It seems that your lattice ?p_lat? is not defined in THREE dimensions, but
> in TWO dimensions.
> Are you really simulating three-dimensional model?
>
> Best,
> Synge
>
>
> On 2014/09/29, at 5:09, Menchyshyn Oleh <oleh.menchyshyn(a)gmail.com> wrote:
>
> > Dear ALPS community,
> >
> > I am trying to obtain the critical temperature for the 3D ferromagnetic
> Heisenberg model S=1/2 on the cubic lattice with a unit cell which contains
> 16 atoms. I have made the QMC "loop" simulations for the lattice with
> dimensions L = 8, i.e. 8*8*8 unit cells, also for L=10, and L=12.
> > But there was just a little bend on the magnetisation curve as a sign of
> the phase transition which should occur for the ferromagnetic model for
> certain. I know I should use the Binder cumulant and the finite size
> scaling to locate the phase transition point correctly.
> > I have enlarged my lattice to L=24, but simulation goes very-very
> slowly. As I have limited resources to just tens of cores and I have a
> feeling I would need to take L=32(48?) at least I want to ask:
> >
> > Based on your experience what order should be the THERMALIZATION and
> SWEEPS parameters?
> > How can I estimate the computational complexity of my problem and the
> time it would take?
> > Or maybe, I do anything wrong?
> >
> > A one more technical issue. The tutorial on the ALPS says that
> correctness of a lattice definition can be checked with "printgraph" tool.
> I used it with a couple definition files but all resulted in:
> >
> > Caught exception: parameter parse error at "<LATTICES> <LATTICEGRAPH
> name="p"
> >
> > Just for a reference I put a definition of a lattice I know worked well
> with simulation tools, but failed with "printgraph"
> >
> > <LATTICES>
> > <LATTICEGRAPH name="p_lat">
> > <FINITELATTICE>
> > <LATTICE dimension="2">
> > <BASIS> <VECTOR>1 0</VECTOR><VECTOR>0 1</VECTOR></BASIS>
> > </LATTICE>
> > <PARAMETER name="L"/>
> > <PARAMETER name="M"/>
> > <EXTENT dimension="1" size="L"/>
> > <EXTENT dimension="2" size="M"/>
> > <BOUNDARY type="periodic"/>
> > </FINITELATTICE>
> > <UNITCELL dimension="2" vertices="6">
> > <VERTEX id="1"><COORDINATE> 0.6 0.2 </COORDINATE></VERTEX>
> > <VERTEX id="2"><COORDINATE> 0.6 0.6 </COORDINATE></VERTEX>
> > <VERTEX id="3"><COORDINATE> 0.2 0.6 </COORDINATE></VERTEX>
> > <VERTEX id="4"><COORDINATE> 0.2 0.2 </COORDINATE></VERTEX>
> > <VERTEX id="5"><COORDINATE> 0.8 0.4 </COORDINATE></VERTEX>
> > <VERTEX id="6"><COORDINATE> 0.4 0.8 </COORDINATE></VERTEX>
> > <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE>
> > <EDGE type="2"><SOURCE vertex="2"/><TARGET vertex="3"/></EDGE>
> > <EDGE type="2"><SOURCE vertex="3"/><TARGET vertex="4"/></EDGE>
> > <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="4"/></EDGE>
> > <EDGE type="3"><SOURCE vertex="1"/><TARGET vertex="5"/></EDGE>
> > <EDGE type="3"><SOURCE vertex="2"/><TARGET vertex="6"/></EDGE>
> > <EDGE type="1"><SOURCE vertex="2"/><TARGET vertex="5"/></EDGE>
> > <EDGE type="1"><SOURCE vertex="3"/><TARGET vertex="6"/></EDGE>
> > <EDGE type="3"><SOURCE vertex="3"/><TARGET vertex="5" offset="-1
> 0"/></EDGE>
> > <EDGE type="1"><SOURCE vertex="4"/><TARGET vertex="5" offset="-1
> 0"/></EDGE>
> > <EDGE type="3"><SOURCE vertex="4"/><TARGET vertex="6" offset="0
> -1"/></EDGE>
> > <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="6" offset="0
> -1"/></EDGE>
> > </UNITCELL>
> > </LATTICEGRAPH>
> > </LATTICES>
> >
> >
> > Best regards,
> > Oleh Menchyshyn
>
>
>
> ------------------------------
>
> Message: 11
> Date: Mon, 29 Sep 2014 11:44:07 +0900
> From: Synge Todo <wistaria(a)comp-phys.org>
> To: comp-phys-alps-users(a)lists.phys.ethz.ch
> Subject: Re: [ALPS-users] complexity estimation for the 3D
> ferromagnetic Heisenberg model
> Message-ID: <1DFD019C-D86C-4BB5-8152-974E17F7B366(a)comp-phys.org>
> Content-Type: text/plain; charset=windows-1252
>
> Dear Oleh,
>
> As for ?printgraph? tool, you have to specify the name of parameter file
> (not lattice file).
> Please prepare a parameter file (say, ?parms?) that includes, e.g.
>
> LATTICE_LIBRARY = "lattice.xml"
> LATTICE = "p_lat"
> L = 4
> M = 4
>
> and execute printgraph tool with the following command line option:
>
> printgraph parms
>
> Best,
> Synge
>
>
> On 2014/09/29, at 11:03, Synge Todo <wistaria(a)comp-phys.org> wrote:
>
> > Dear Oleh,
> >
> > It seems that your lattice ?p_lat? is not defined in THREE dimensions,
> but in TWO dimensions.
> > Are you really simulating three-dimensional model?
> >
> > Best,
> > Synge
> >
> >
> > On 2014/09/29, at 5:09, Menchyshyn Oleh <oleh.menchyshyn(a)gmail.com>
> wrote:
> >
> >> Dear ALPS community,
> >>
> >> I am trying to obtain the critical temperature for the 3D ferromagnetic
> Heisenberg model S=1/2 on the cubic lattice with a unit cell which contains
> 16 atoms. I have made the QMC "loop" simulations for the lattice with
> dimensions L = 8, i.e. 8*8*8 unit cells, also for L=10, and L=12.
> >> But there was just a little bend on the magnetisation curve as a sign
> of the phase transition which should occur for the ferromagnetic model for
> certain. I know I should use the Binder cumulant and the finite size
> scaling to locate the phase transition point correctly.
> >> I have enlarged my lattice to L=24, but simulation goes very-very
> slowly. As I have limited resources to just tens of cores and I have a
> feeling I would need to take L=32(48?) at least I want to ask:
> >>
> >> Based on your experience what order should be the THERMALIZATION and
> SWEEPS parameters?
> >> How can I estimate the computational complexity of my problem and the
> time it would take?
> >> Or maybe, I do anything wrong?
> >>
> >> A one more technical issue. The tutorial on the ALPS says that
> correctness of a lattice definition can be checked with "printgraph" tool.
> I used it with a couple definition files but all resulted in:
> >>
> >> Caught exception: parameter parse error at "<LATTICES> <LATTICEGRAPH
> name="p"
> >>
> >> Just for a reference I put a definition of a lattice I know worked well
> with simulation tools, but failed with "printgraph"
> >>
> >> <LATTICES>
> >> <LATTICEGRAPH name="p_lat">
> >> <FINITELATTICE>
> >> <LATTICE dimension="2">
> >> <BASIS> <VECTOR>1 0</VECTOR><VECTOR>0 1</VECTOR></BASIS>
> >> </LATTICE>
> >> <PARAMETER name="L"/>
> >> <PARAMETER name="M"/>
> >> <EXTENT dimension="1" size="L"/>
> >> <EXTENT dimension="2" size="M"/>
> >> <BOUNDARY type="periodic"/>
> >> </FINITELATTICE>
> >> <UNITCELL dimension="2" vertices="6">
> >> <VERTEX id="1"><COORDINATE> 0.6 0.2 </COORDINATE></VERTEX>
> >> <VERTEX id="2"><COORDINATE> 0.6 0.6 </COORDINATE></VERTEX>
> >> <VERTEX id="3"><COORDINATE> 0.2 0.6 </COORDINATE></VERTEX>
> >> <VERTEX id="4"><COORDINATE> 0.2 0.2 </COORDINATE></VERTEX>
> >> <VERTEX id="5"><COORDINATE> 0.8 0.4 </COORDINATE></VERTEX>
> >> <VERTEX id="6"><COORDINATE> 0.4 0.8 </COORDINATE></VERTEX>
> >> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE>
> >> <EDGE type="2"><SOURCE vertex="2"/><TARGET vertex="3"/></EDGE>
> >> <EDGE type="2"><SOURCE vertex="3"/><TARGET vertex="4"/></EDGE>
> >> <EDGE type="2"><SOURCE vertex="1"/><TARGET vertex="4"/></EDGE>
> >> <EDGE type="3"><SOURCE vertex="1"/><TARGET vertex="5"/></EDGE>
> >> <EDGE type="3"><SOURCE vertex="2"/><TARGET vertex="6"/></EDGE>
> >> <EDGE type="1"><SOURCE vertex="2"/><TARGET vertex="5"/></EDGE>
> >> <EDGE type="1"><SOURCE vertex="3"/><TARGET vertex="6"/></EDGE>
> >> <EDGE type="3"><SOURCE vertex="3"/><TARGET vertex="5" offset="-1
> 0"/></EDGE>
> >> <EDGE type="1"><SOURCE vertex="4"/><TARGET vertex="5" offset="-1
> 0"/></EDGE>
> >> <EDGE type="3"><SOURCE vertex="4"/><TARGET vertex="6" offset="0
> -1"/></EDGE>
> >> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="6" offset="0
> -1"/></EDGE>
> >> </UNITCELL>
> >> </LATTICEGRAPH>
> >> </LATTICES>
> >>
> >>
> >> Best regards,
> >> Oleh Menchyshyn
>