Dear ALPS team,
I am trying to obtain magnetization curve of anisotropic XY model on frustrated diamond spin chain in transverse magnetic field (m^x) using full diagonalization. For this purpose I have prepared such input file for delta=1.5;
LATTICE="lattice_diamond" MODEL="spin" LATTICE_LIBRARY="lattice-diamond.xml" local_S=1/2
Jxy1 = 0.85 Jxy2 = 3.0 Jxy3 = 1.15
Jz1 = 1.275 Jz2 = 4.5 Jz3 = 1.725
L = 2 T = 0.001
{Gamma =4;} {Gamma =4.05;} {Gamma =4.1;} {Gamma =4.15;} {Gamma =4.2;} {Gamma =4.25;} {Gamma =4.3;} {Gamma =4.35;} {Gamma =4.4;} {Gamma =4.45;} {Gamma =4.5;} {Gamma =4.55;} {Gamma =4.6;} {Gamma =4.65;} {Gamma =4.7;} {Gamma =4.75;} {Gamma =4.8;} {Gamma =4.85;} {Gamma =4.9;} {Gamma =4.95;} {Gamma =5;} {Gamma =5.05;} {Gamma =5.1;} {Gamma =5.15;} {Gamma =5.2;} {Gamma =5.25;} {Gamma =5.3;} {Gamma =5.35;} {Gamma =5.4;} {Gamma =5.45;} {Gamma =5.5;} {Gamma =5.55;} {Gamma =5.6;} {Gamma =5.65;} {Gamma =5.7;} {Gamma =5.75;} {Gamma =5.8;} {Gamma =5.85;} {Gamma =5.9;} {Gamma =5.95;} {Gamma =6;}
lattice-diamond.xml is:
<LATTICES> <LATTICEGRAPH name="lattice_diamond"> <FINITELATTICE> <LATTICE dimension="1"> <BASIS> <VECTOR>1</VECTOR> </BASIS> </LATTICE> <PARAMETER name="L"/> <EXTENT size="L"/> <BOUNDARY type="periodic"/> </FINITELATTICE> <UNITCELL dimension="1" vertices="3"> <VERTEX id="1" type="0"> </VERTEX> <VERTEX id="2" type="0"> </VERTEX> <VERTEX id="3" type="0"> </VERTEX> <EDGE type="2"> <SOURCE vertex="1"/> <TARGET vertex="2"/> </EDGE> <EDGE type="3"> <SOURCE vertex="1"/> <TARGET vertex="3"/> </EDGE> <EDGE type="1"> <SOURCE vertex="2"/> <TARGET vertex="3"/> </EDGE> <EDGE type="3"> <SOURCE vertex="3"/> <TARGET vertex="2" offset="1"/> </EDGE> <EDGE type="1"> <SOURCE vertex="3"/> <TARGET vertex="1" offset="1"/> </EDGE> </UNITCELL> </LATTICEGRAPH> </LATTICES>
It works perfect for L=2 (6 sites), but I cant obtain correct magnetization curve for L=3 (9 sites) or L= 4 (12 sites) etc. Would you be so kind to explain me why it works well for small chain and worse for slightly larger chain?
Best regards,
Olesia Krupnitska
comp-phys-alps-users@lists.phys.ethz.ch