Thank you very much Synge. Yes, that was the problem. I deleted all the files and then only run the script and it worked now.
Many thanks,
Santu
On 7 August 2018 at 21:08, Synge Todo wistaria@phys.s.u-tokyo.ac.jp wrote:
Dear Santu,
It seems that there are two data points at each temperature. Please make sure that there remain no old files from your previous run for a different parameter set.
Please find attached the plot produced from your script. A clear Curie law is observed.
Best, Synge
2018/08/07 17:03、S Baidya santubaidya2009@gmail.comのメール:
Thank you Rongyang for your reply.
My script is :
import pyalps import matplotlib.pyplot as plt import pyalps.plot
#prepare the input parameters parms = [] for t in [0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.25, 1.5, 1.75, 2.0]: parms.append( { 'LATTICE' : "Kagome lattice", 'MODEL' : "spin", 'local_S' : 1.0, 'T' : t, 'J' : -5 , 'THERMALIZATION' : 5000, 'SWEEPS' : 50000, 'L' : 24, 'ALGORITHM' : "loop" } )
#write the input file and run the simulation input_file = pyalps.writeInputFiles('parm2d',parms) pyalps.runApplication('loop',input_file)
#load the susceptibility and collect it as function of temperature T data = pyalps.loadMeasurements(pyalps.getResultFiles(prefix=' parm2d'),'Susceptibility') susceptibility = pyalps.collectXY(data,x='T',y='Susceptibility')
#make plot plt.figure() pyalps.plot.plot(susceptibility) plt.xlabel('Temperature $T/J$') plt.ylabel('Susceptibility $\chi J$') plt.ylim(0,1000.00) plt.title('Kagome Heisenberg lattice') plt.show()
The code run without problem but there can be some problem with J directions related to Kagome.
Any help would be appreciated.
Thank you,
Santu
On 7 August 2018 at 16:58, Rongyang Sun sun-rongyang@outlook.com wrote: Dear Santu,
I am not familiar with the MC algorithm, but can you show your python script which can get this result?
BTW, you can define your own lattice library in an XML file and use it in your python script without recompile your ALPS source code.
Best, Rongyang
From: Comp-phys-alps-users <comp-phys-alps-users-bounces@ lists.phys.ethz.ch> on behalf of S Baidya santubaidya2009@gmail.com Sent: Tuesday, August 7, 2018 15:43 To: comp-phys-alps-users@lists.phys.ethz.ch Subject: Re: [ALPS-users] ALPS code error with Kagome spin lattice to calculate susceptibility
Dear Alps users,
I just resolved now to get the susceptibility vs temp plot for Kagome FM lattice. But, The plot appears to be unacceptable.
I would ask if there is any problem with the definition of Kagome lattice already in the lattice library. Should I define by myself and recompile the code?
<Kagome-FM-chi-vs-temp.png>
Please give me one solution for this problem.
Thanks,
Santu
On 4 August 2018 at 14:08, S Baidya santubaidya2009@gmail.com wrote: Thank you for your reply. I recheked again. I did not have problem in getting the original tutorial mc-02 plots although there was some matplotlib problem.
I already reproduced also tutorial mc-03 and there was no problem in getting the plots.
So, is there any other solution possible.
Thanks,
Santu
On 3 August 2018 at 18:33, Synge Todo <wistaria@phys.s.u-tokyo.ac.jp
wrote:
Dear Santu,
It seems to me that the QMC calculation it self finishes normally, and there is a problem in installation of matplotlib on your system. Have you succeeded in plotting the susceptibility for the original tutorial2d on the spin ladder system?
Best, Synge
2018/08/03 16:57、S Baidya santubaidya2009@gmail.comのメール:
Dear Synge and alps user,
This time I did a refresh calculation for FM Kagome lattice with
S=3/2 and J=-4.2 using tutorial2d.py but still got no data as output.
/usr/local/lib/python2.7/dist-packages/matplotlib/style/core.py:203:
UserWarning: In /usr/share/matplotlib/stylelib/classic.mplstyle: text.dvipnghack is obsolete. Please remove it from your matplotlibrc and/or style files.
warnings.warn(message) loop parm2d.in.xml ALPS/looper version 3.2b12-20100128 (2010/01/28) multi-cluster quantum Monte Carlo algorithms for spin systems available from http://wistaria.comp-phys.org/alps-looper/ copyright (c) 1997-2010 by Synge Todo wistaria@comp-phys.org
using ALPS/parapack scheduler a Monte Carlo scheduler for multiple-level parallelization copyright (c) 1997-2016 by Synge Todo wistaria@comp-phys.org
based on the ALPS libraries version 2.3.0 available from http://alps.comp-phys.org/ copyright (c) 1994-2016 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.
[2018-Aug-03 16:05:55]: starting scheduler on cces-System-Product-Name master input file = /home/cces/kagome-alps/mc-02/parm2d.in.xml master output file = /home/cces/kagome-alps/mc-02/parm2d.out.xml termination file = [disabled] total number of thread(s) = 1 thread(s) per clone = 1 number of thread group(s) = 1 auto evaluation = yes time limit = unlimited interval between checkpointing = 3600 seconds interval between progress report = 600 seconds interval between vmusage report = infinity task range = all worker dump format = hdf5 worker dump policy = running workers only [2018-Aug-03 16:05:55]: task status: total number of tasks = 15 new = 15, running = 0, continuing = 0, suspended = 0, finished = 0,
completed = 0, skipped = 0
[2018-Aug-03 16:05:55]: starting 1 threadgroup(s) [2018-Aug-03 16:05:55]: dispatching a new clone[1,1] on threadgroup[1] [2018-Aug-03 16:11:55]: checkpointing task files save task [2018-Aug-03 16:11:55]: task status: total number of tasks = 15 new = 14, running = 1, continuing = 0, suspended = 0, finished = 0,
completed = 0, skipped = 0
[2018-Aug-03 16:15:55]: progress report: clone[1,1] is running (50.9%
done)
[2018-Aug-03 16:25:54]: clone[1,1] finished on threadgroup[1] [2018-Aug-03 16:25:54]: dispatching a new clone[2,1] on threadgroup[1] [2018-Aug-03 16:34:18]: clone[2,1] finished on threadgroup[1] [2018-Aug-03 16:34:18]: dispatching a new clone[3,1] on threadgroup[1] [2018-Aug-03 16:38:29]: clone[3,1] finished on threadgroup[1] [2018-Aug-03 16:38:29]: dispatching a new clone[4,1] on threadgroup[1] [2018-Aug-03 16:41:19]: clone[4,1] finished on threadgroup[1] [2018-Aug-03 16:41:19]: dispatching a new clone[5,1] on threadgroup[1] [2018-Aug-03 16:43:27]: clone[5,1] finished on threadgroup[1] [2018-Aug-03 16:43:27]: dispatching a new clone[6,1] on threadgroup[1] [2018-Aug-03 16:45:09]: clone[6,1] finished on threadgroup[1] [2018-Aug-03 16:45:09]: dispatching a new clone[7,1] on threadgroup[1] [2018-Aug-03 16:46:33]: clone[7,1] finished on threadgroup[1] [2018-Aug-03 16:46:33]: dispatching a new clone[8,1] on threadgroup[1] [2018-Aug-03 16:47:46]: clone[8,1] finished on threadgroup[1] [2018-Aug-03 16:47:46]: dispatching a new clone[9,1] on threadgroup[1] [2018-Aug-03 16:48:50]: clone[9,1] finished on threadgroup[1] [2018-Aug-03 16:48:50]: dispatching a new clone[10,1] on threadgroup[1] [2018-Aug-03 16:49:47]: clone[10,1] finished on threadgroup[1] [2018-Aug-03 16:49:47]: dispatching a new clone[11,1] on threadgroup[1] [2018-Aug-03 16:50:40]: clone[11,1] finished on threadgroup[1] [2018-Aug-03 16:50:40]: dispatching a new clone[12,1] on threadgroup[1] [2018-Aug-03 16:51:22]: clone[12,1] finished on threadgroup[1] [2018-Aug-03 16:51:22]: dispatching a new clone[13,1] on threadgroup[1] [2018-Aug-03 16:51:58]: clone[13,1] finished on threadgroup[1] [2018-Aug-03 16:51:58]: dispatching a new clone[14,1] on threadgroup[1] [2018-Aug-03 16:52:30]: clone[14,1] finished on threadgroup[1] [2018-Aug-03 16:52:30]: dispatching a new clone[15,1] on threadgroup[1] [2018-Aug-03 16:52:58]: clone[15,1] finished on threadgroup[1] [2018-Aug-03 16:52:58]: all tasks have been finished [2018-Aug-03 16:52:58]: task status: total number of tasks = 15 new = 0, running = 0, continuing = 0, suspended = 0, finished = 0,
completed = 15, skipped = 0
[2018-Aug-03 16:52:58]: all threads halted [2018-Aug-03 16:52:58]: starting evaluation on cces-System-Product-Name evaluating parm2d.task1.out.xml loading clones: 1 evaluating parm2d.task2.out.xml loading clones: 1 evaluating parm2d.task3.out.xml loading clones: 1 evaluating parm2d.task4.out.xml loading clones: 1 evaluating parm2d.task5.out.xml loading clones: 1 evaluating parm2d.task6.out.xml loading clones: 1 evaluating parm2d.task7.out.xml loading clones: 1 evaluating parm2d.task8.out.xml loading clones: 1 evaluating parm2d.task9.out.xml loading clones: 1 evaluating parm2d.task10.out.xml loading clones: 1 evaluating parm2d.task11.out.xml loading clones: 1 evaluating parm2d.task12.out.xml loading clones: 1 evaluating parm2d.task13.out.xml loading clones: 1 evaluating parm2d.task14.out.xml loading clones: 1 evaluating parm2d.task15.out.xml loading clones: 1 [2018-Aug-03 16:52:58]: all tasks evaluated Cannot open file '/usr/share/matplotlib/images/matplotlib.svg',
because: No such file or directory
Cannot open file '/usr/share/matplotlib/images/matplotlib.svg',
because: No such file or directory
Can anyone please tell me if there is any modification is needed for two
dimensional Kagome spin lattice.
Thnaks,
Santu
On 3 August 2018 at 14:56, S Baidya santubaidya2009@gmail.com wrote: Dear Synge,
I just removed param2d and used tutorial2d.py as
import pyalps import matplotlib.pyplot as plt import pyalps.plot
#prepare the input parameters parms = [] for t in [0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.25,
1.5, 1.75, 2.0]:
parms.append( { 'LATTICE' : "Kagome lattice", 'MODEL' : "spin", 'local_S' : 1.5, 'T' : t, 'J' : -4.2 , 'THERMALIZATION' : 5000, 'SWEEPS' : 50000, 'L' : 60, 'ALGORITHM' : "loop" } )#write the input file and run the simulation input_file = pyalps.writeInputFiles('parm2d',parms) pyalps.runApplication('loop',input_file)
#load the susceptibility and collect it as function of temperature T data = pyalps.loadMeasurements(pyalps.getResultFiles(prefix='
parm2d'),'Susceptibility')
susceptibility = pyalps.collectXY(data,x='T',y='Susceptibility')
#make plot plt.figure() pyalps.plot.plot(susceptibility) plt.xlabel('Temperature $T/J$') plt.ylabel('Susceptibility $\chi J$') plt.ylim(0,0.22) plt.title('Quantum Heisenberg ladder') plt.show()
Then typed alpspython tutorial2d.py and it shows the same situation
again.....
/usr/local/lib/python2.7/dist-packages/matplotlib/style/core.py:203:
UserWarning: In /usr/share/matplotlib/stylelib/classic.mplstyle: text.dvipnghack is obsolete. Please remove it from your matplotlibrc and/or style files.
warnings.warn(message) loop parm2d.in.xml ALPS/looper version 3.2b12-20100128 (2010/01/28) multi-cluster quantum Monte Carlo algorithms for spin systems available from http://wistaria.comp-phys.org/alps-looper/ copyright (c) 1997-2010 by Synge Todo wistaria@comp-phys.org
using ALPS/parapack scheduler a Monte Carlo scheduler for multiple-level parallelization copyright (c) 1997-2016 by Synge Todo wistaria@comp-phys.org
based on the ALPS libraries version 2.3.0 available from http://alps.comp-phys.org/ copyright (c) 1994-2016 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.
[2018-Aug-03 14:54:16]: starting scheduler on cces-System-Product-Name master input file = /home/cces/kagome-alps/mc-02/parm2d.in.xml master output file = /home/cces/kagome-alps/mc-02/parm2d.out.xml termination file = [disabled] total number of thread(s) = 1 thread(s) per clone = 1 number of thread group(s) = 1 auto evaluation = yes time limit = unlimited interval between checkpointing = 3600 seconds interval between progress report = 600 seconds interval between vmusage report = infinity task range = all worker dump format = hdf5 worker dump policy = running workers only [2018-Aug-03 14:54:16]: task status: total number of tasks = 15 new = 15, running = 0, continuing = 0, suspended = 0, finished = 0,
completed = 0, skipped = 0
[2018-Aug-03 14:54:16]: starting 1 threadgroup(s) [2018-Aug-03 14:54:16]: dispatching a new clone[1,1] on threadgroup[1]
Thanks,
Santu
On 3 August 2018 at 14:51, Synge Todo wistaria@phys.s.u-tokyo.ac.jp
wrote:
Dear Santu,
Strange. Could you remove parm2d.* in your working directory and run
tutorial2d.py again?
BTW, 2d lattice with L=60 is a bit large as a starting point. I would
recommend you to begin with
smaller value, say, L=8.
Best, Synge
2018/08/03 14:39、S Baidya santubaidya2009@gmail.comのメール:
Dear Synge,
Thank you for your reply. I then made a small mistake. The "J" value
should be "-4.2" as it is FM. If it is now FM it should work like ising model in any 2d lattice. But I am still getting the same error.
Is there any other solution for dealing this situation.
Thanking you,
Santu
On 3 August 2018 at 14:30, Synge Todo wistaria@phys.s.u-tokyo.ac.jp
wrote:
Dear Santu Baidya,
The antiferromagnetic (positive J) Heisenberg model on the kagome
lattice has frustration.
This causes the so-called negative sign problem in the quantum Monte
Carlo methods.
Although the loop code can simulate such a system and produce some
numbers, the errorbar
becomes very large especially at low temperatures.
You can use the exact diagonalization for the frustrated spin models,
but the lattice size is
very limited in this case due to the spin size S=3/2.
Best, Synge
2018/08/03 14:10、S Baidya santubaidya2009@gmail.comのメール:
Dear Alps users,
I am totally new to ALPS code using quantum monte carlo. I have
started to use this code for few small purpose like magnetic suceptibility vs temperature for various spin lattice models.
Currently I have used ALPS code for calculating susceptibility for
ising type S=3/2 Kagome lattice. I found in the lattice library that Kagome lattice is defined. So in the tutorial mc-02-susceptibilities/ I changed tutorial2d.py by editing
#prepare the input parameters parms = [] for t in [0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0,
1.25, 1.5, 1.75, 2.0]:
parms.append( { 'LATTICE' : "Kagome lattice", 'MODEL' : "spin", 'local_S' : 1.5, 'T' : t, 'J' : 4.2 , 'THERMALIZATION' : 5000, 'SWEEPS' : 50000, 'L' : 60, 'ALGORITHM' : "loop" } )Then I changed the same way param2d file and ran the alpspython
tutorial2d.py and I got strange error
/usr/local/lib/python2.7/dist-packages/matplotlib/style/core.py:203:
UserWarning: In /usr/share/matplotlib/stylelib/classic.mplstyle: text.dvipnghack is obsolete. Please remove it from your matplotlibrc and/or style files.
warnings.warn(message) loop parm2d.in.xml ALPS/looper version 3.2b12-20100128 (2010/01/28) multi-cluster quantum Monte Carlo algorithms for spin systems available from http://wistaria.comp-phys.org/alps-looper/ copyright (c) 1997-2010 by Synge Todo wistaria@comp-phys.org
using ALPS/parapack scheduler a Monte Carlo scheduler for multiple-level parallelization copyright (c) 1997-2016 by Synge Todo wistaria@comp-phys.org
based on the ALPS libraries version 2.3.0 available from http://alps.comp-phys.org/ copyright (c) 1994-2016 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.
[2018-Aug-03 14:02:51]: starting scheduler on
cces-System-Product-Name
master input file = /home/cces/kagome-alps/mc-02/parm2d.in.xml master output file = /home/cces/kagome-alps/mc-02/parm2d.out.xml termination file = [disabled] total number of thread(s) = 1 thread(s) per clone = 1 number of thread group(s) = 1 auto evaluation = yes time limit = unlimited interval between checkpointing = 3600 seconds interval between progress report = 600 seconds interval between vmusage report = infinity task range = all worker dump format = hdf5 worker dump policy = running workers only [2018-Aug-03 14:02:51]: task status: total number of tasks = 15 new = 15, running = 0, continuing = 0, suspended = 0, finished =
0, completed = 0, skipped = 0
[2018-Aug-03 14:02:51]: starting 1 threadgroup(s) [2018-Aug-03 14:02:51]: dispatching a new clone[1,1] on
threadgroup[1]
WARNING: model is classically frustrated WARNING: model has negative signs
I donot understand the actual reason and how to solve this
situation. I used this tutorial before and it worked for ladder spin lattice.
Can anyone please tell me what should I do to solve this situation.
Thanking you,
Santu Baidya
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@
lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@
lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@
lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@
lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@
lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@ lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@ lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@ lists.phys.ethz.ch.
Comp-phys-alps-users Mailing List for the ALPS Project http://alps.comp-phys.org/
List info: https://lists.phys.ethz.ch//listinfo/comp-phys-alps-users Archive: https://lists.phys.ethz.ch//pipermail/comp-phys-alps-users
Unsubscribe by writing a mail to comp-phys-alps-users-leave@ lists.phys.ethz.ch.
