That looks ok if you use a site type of 1 in the middle and 0 at the boundary.
Matthias
On 31 Jan 2012, at 12:09, Javad Vahedi wrote:
Dear M. Troyer
Thanks for your immediate comments, I have prepared the following input file, but I am not very sure that it is correct or not. I would be grateful If you check and find my mistakes.
######################################### LATTICE_LIBRARY="my_lattices.xml" LATTICE="open chain lattice with especial edges 32" MODEL="spin" local_S=1 SWEEPS=6 MAXSTATES=40 Jz=-1 Jxy=0 %h0=boundary field h0=0.5 Gamma=1 D=1 NUMBER_EIGENVALUES=2 MEASURE_AVERAGE[Magnetization]=Sz MEASURE_LOCAL[Local magnetization]=Sz %h1=bulk field { h1=0 } { h1=0.1 } { h1=0.2 } { h1=0 .3} { h1=0 4} { h1=0.5} ... %and so on ######################################
-- Best Regards Javad Vahedi, Dept. of Physics, Sari Branch of Islamic Azad University, P.O.Box 48164-194, Sari, Iran. Tel: (+98) 911 1554504 Fax: (+98)151 2256679 E-mail: javahedi@gmail.com javahedi@iausari.ac.ir