- Comp-phys-alps-users - lists.phys.ethz.ch

DMRG Su(n) model
by Moreno William 08 Jan '14

08 Jan '14

Dear developers, I have written a model (who can be used as a model library) who generalize to classical Heisenberg model for fermion with spin up and down to fermion with ’n’ degrees of freedom and two particles per sites. Here is my code <MODELS> <SITEBASIS name ="SU5_m2"> <QUANTUMNUMBER name="N_1" min ="0" max ="1" type="fermionic"/> <QUANTUMNUMBER name="N_2" min ="0" max ="1" type="fermionic"/> <QUANTUMNUMBER name="N_3" min ="0" max ="(1 - N_1*N_2)" type ="fermionic"/> <QUANTUMNUMBER name="N_4" min ="(1 - N_1)*(1 - N_2)*(1 - N_3)" max ="(1 - N_1*N_2)*(1 - N_1*N_3)*(1 - N_2*N_3)" type ="fermionic"/> <QUANTUMNUMBER name="N_5" min ="2 - N_1 - N_2 - N_3 - N_4" max ="2 - N_1 - N_2 - N_3 - N_4" type ="fermionic"/> <OPERATOR name="S_1_1" matrixelement="N_1"/> <OPERATOR name="S_2_2" matrixelement="N_2"/> <OPERATOR name="S_3_3" matrixelement="N_3"/> <OPERATOR name="S_4_4" matrixelement="N_4"/> <OPERATOR name="S_5_5" matrixelement="N_5"/> <OPERATOR name="S_1_2" matrixelement="1"> <CHANGE quantumnumber="N_2" change="-1"/> <CHANGE quantumnumber="N_1" change="+1"/> </OPERATOR> <OPERATOR name="S_1_3" matrixelement="1"> <CHANGE quantumnumber="N_3" change="-1"/> <CHANGE quantumnumber="N_1" change="+1"/> </OPERATOR> <OPERATOR name="S_1_4" matrixelement="1"> <CHANGE quantumnumber="N_4" change="-1"/> <CHANGE quantumnumber="N_1" change="+1"/> </OPERATOR> <OPERATOR name="S_1_5" matrixelement="1"> <CHANGE quantumnumber="N_5" change="-1"/> <CHANGE quantumnumber="N_1" change="+1"/> </OPERATOR> <OPERATOR name="S_2_1" matrixelement="1"> <CHANGE quantumnumber="N_1" change="-1"/> <CHANGE quantumnumber="N_2" change="+1"/> </OPERATOR> <OPERATOR name="S_2_3" matrixelement="1"> <CHANGE quantumnumber="N_3" change="-1"/> <CHANGE quantumnumber="N_2" change="+1"/> </OPERATOR> <OPERATOR name="S_2_4" matrixelement="1"> <CHANGE quantumnumber="N_4" change="-1"/> <CHANGE quantumnumber="N_2" change="+1"/> </OPERATOR> <OPERATOR name="S_2_5" matrixelement="1"> <CHANGE quantumnumber="N_5" change="-1"/> <CHANGE quantumnumber="N_2" change="+1"/> </OPERATOR> <OPERATOR name="S_3_1" matrixelement="1"> <CHANGE quantumnumber="N_1" change="-1"/> <CHANGE quantumnumber="N_3" change="+1"/> </OPERATOR> <OPERATOR name="S_3_2" matrixelement="1"> <CHANGE quantumnumber="N_2" change="-1"/> <CHANGE quantumnumber="N_3" change="+1"/> </OPERATOR> <OPERATOR name="S_3_4" matrixelement="1"> <CHANGE quantumnumber="N_4" change="-1"/> <CHANGE quantumnumber="N_3" change="+1"/> </OPERATOR> <OPERATOR name="S_3_5" matrixelement="1"> <CHANGE quantumnumber="N_5" change="-1"/> <CHANGE quantumnumber="N_3" change="+1"/> </OPERATOR> <OPERATOR name="S_4_1" matrixelement="1"> <CHANGE quantumnumber="N_1" change="-1"/> <CHANGE quantumnumber="N_4" change="+1"/> </OPERATOR> <OPERATOR name="S_4_2" matrixelement="1"> <CHANGE quantumnumber="N_2" change="-1"/> <CHANGE quantumnumber="N_4" change="+1"/> </OPERATOR> <OPERATOR name="S_4_3" matrixelement="1"> <CHANGE quantumnumber="N_3" change="-1"/> <CHANGE quantumnumber="N_4" change="+1"/> </OPERATOR> <OPERATOR name="S_4_5" matrixelement="1"> <CHANGE quantumnumber="N_5" change="-1"/> <CHANGE quantumnumber="N_4" change="+1"/> </OPERATOR> <OPERATOR name="S_5_1" matrixelement="1"> <CHANGE quantumnumber="N_1" change="-1"/> <CHANGE quantumnumber="N_5" change="+1"/> </OPERATOR> <OPERATOR name="S_5_2" matrixelement="1"> <CHANGE quantumnumber="N_2" change="-1"/> <CHANGE quantumnumber="N_5" change="+1"/> </OPERATOR> <OPERATOR name="S_5_3" matrixelement="1"> <CHANGE quantumnumber="N_3" change="-1"/> <CHANGE quantumnumber="N_5" change="+1"/> </OPERATOR> <OPERATOR name="S_5_4" matrixelement="1"> <CHANGE quantumnumber="N_4" change="-1"/> <CHANGE quantumnumber="N_5" change="+1"/> </OPERATOR> </SITEBASIS> <BASIS name = "SU5_m2" > <SITEBASIS ref="SU5_m2"/> <CONSTRAINT quantumnumber="N_1" value="N_total"/> <CONSTRAINT quantumnumber="N_2" value="N_total"/> <CONSTRAINT quantumnumber="N_3" value="N_total"/> <CONSTRAINT quantumnumber="N_4" value="N_total"/> <CONSTRAINT quantumnumber="N_5" value="N_total"/> </BASIS> <HAMILTONIAN name ="SU5_m2" > <BASIS ref="SU5_m2"/> <BONDTERM source="i" target="j"> S_1_1(i)*S_1_1(j) + S_1_2(i)*S_2_1(j) + S_1_3(i)*S_3_1(j) + S_1_4(i)*S_4_1(j) + S_1_5(i)*S_5_1(j) + S_2_1(i)*S_1_2(j) + S_2_2(i)*S_2_2(j) + S_2_3(i)*S_3_2(j) + S_2_4(i)*S_4_2(j) + S_2_5(i)*S_5_2(j) + S_3_1(i)*S_1_3(j) + S_3_2(i)*S_2_3(j) + S_3_3(i)*S_3_3(j) + S_3_4(i)*S_4_3(j) + S_3_5(i)*S_5_3(j) + S_4_1(i)*S_1_4(j) + S_4_2(i)*S_2_4(j) + S_4_3(i)*S_3_4(j) + S_4_4(i)*S_4_4(j) + S_4_5(i)*S_5_4(j) + S_5_1(i)*S_1_5(j) + S_5_2(i)*S_2_5(j) + S_5_3(i)*S_3_5(j) + S_5_4(i)*S_4_5(j) + S_5_5(i)*S_5_5(j) </BONDTERM> </HAMILTONIAN> </MODELS> The problem is that the ground state energy is -622 and the first excited energy is -670 in a open chain lattice with 5 sites. Obviously somethings is wrong because the ground state energy is greater than the first excited energy. I don’t know if the min and max condition in the quantum number are correctly implemented or maybe this is not the correct way to implement the bond term? Can you help me ? Best regards, William Moreno

1 0

Hybridization expansion Solver runs longer than specified
by backes＠th.physik.uni-frankfurt.de 06 Jan '14

06 Jan '14

Dear developers, I recently checked out the newest version of the Alps code rev. 7291 and compiled the Hybridization expansion Impurity solver with the Intel compiler composer_xe_2013.3.163 on Ubuntu 12.04 LTS. My problem is that the code takes longer to exit than specified in MAX_TIME in the parameter file (SWEEPS much higher than possible in this time). I.e. MAX_TIME = 30 --> runs for 60sec MAX_TIME = 60 --> runs for 180sec MAX_TIME = 120 --> runs for 180sec MAX_TIME = 180 --> runs for 420sec MAX_TIME = 240 --> runs for 420sec Also for longer time and mpi-parallelization the behaviour is similar. With an older version checked out about 2 months ago this problem did not occur. Could you check if this problem is caused by some change in the code? Thanks, Steffen

3 4

Re: [ALPS-users] Constraints for a unit cell (Elton Nascimento)
by Elton Nascimento 02 Jan '14

02 Jan '14

Dear Matthias Troyer You are right, but the lattice we are using shows hop energy between the two sites of the unit cell (edge type = 0) and an exchange energy between sites in adjacent cells (edge type = 1), as bellow <UNITCELL name="Nascimento" dimension="1"> <VERTEX type="0"/> <VERTEX type="0"/> <EDGE type="1"><SOURCE vertex="1"/><TARGET vertex="1" offset="1"/></EDGE> <EDGE type="1"><SOURCE vertex="2"/><TARGET vertex="2" offset="1"/></EDGE> <EDGE type="0"><SOURCE vertex="1"/><TARGET vertex="2"/></EDGE> </UNITCELL> Because of the hop energy, the particles have transition probability between the sites of the unit cell, but not between adjacent cells. This is why we need ensure the total particles of the unit cell be constant. Elton M. Nascimento

2 1

ALPS 2.2 now available in Archlinux AUR repository
by Mohammadreza Abdollahzadeh 02 Jan '14

02 Jan '14

Hi everyone, for those who use Archlinux like us! you can now install alps 2.2 from AUR repositoty: https://aur.archlinux.org/packages/alps/

1 0

cmake error in configuring alps 2.1.1
by Mohammadreza Abdollahzadeh 02 Jan '14

02 Jan '14

2 2

Re: [ALPS-users] Comp-phys-alps-users Digest, Vol 93, Issue 11
by Elton Nascimento 31 Dec '13

31 Dec '13

Dear Matthias Troyer Thanks for you attention. You are right on your first option, that is exactly I have: a complex lattice with two sites per unit cell and I must ensure exactly two particles per unit cell. 2013/12/30 <comp-phys-alps-users-request(a)lists.phys.ethz.ch> > Send Comp-phys-alps-users mailing list submissions to > comp-phys-alps-users(a)lists.phys.ethz.ch > > To subscribe or unsubscribe via the World Wide Web, visit > https://lists.phys.ethz.ch/listinfo/comp-phys-alps-users > or, via email, send a message with subject or body 'help' to > comp-phys-alps-users-request(a)lists.phys.ethz.ch > > You can reach the person managing the list at > comp-phys-alps-users-owner(a)lists.phys.ethz.ch > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of Comp-phys-alps-users digest..." > > > Today's Topics: > > 1. Re: Constraints for a unit cell (Matthias Troyer) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 29 Dec 2013 15:37:02 +0100 > From: Matthias Troyer <troyer(a)phys.ethz.ch> > To: ALPS user list <comp-phys-alps-users(a)lists.phys.ethz.ch> > Subject: Re: [ALPS-users] Constraints for a unit cell > Message-ID: <1FC5704D-82EB-4A08-A16F-08376885EA58(a)phys.ethz.ch> > Content-Type: text/plain; charset=iso-8859-1 > > Hi, > > what do you mean by 2-site unit cell? Do you have a complex lattice with > two sites per unit cell and you want to have exactly two particles per unit > cell, or do you want to have two particles in total on the whole lattice? > > > On 28 Dec 2013, at 16:59, Elton Nascimento <emn(a)fis.ufal.br> wrote: > > > Dear All > > > > How can I constrain the particle numbers of a 2-sites unit cell for a > constant N=2? I am trying to obtain local magnetization using DMRG > application and t-J model. > > > > Thanks everybody. > > > > Elton M. Nascimento > > Instituto de F?sica - UFAL-Brazil > > > > End of Comp-phys-alps-users Digest, Vol 93, Issue 11 > **************************************************** >

2 1

Fw: cmake error in configuring alps 2.1.1
by Mohammadreza Abdollahzadeh 31 Dec '13

31 Dec '13

hi every one I try to install alps, but I get following error from cmake: CMake Error: The following variables are used in this project, but they are set to NOTFOUND. Please set them or make sure they are set and tested correctly in the CMake files: HDF5_dl_LIBRARY_DEBUG (ADVANCED) HDF5_hdf5_LIBRARY_DEBUG (ADVANCED) HDF5_hdf5_hl_LIBRARY_DEBUG (ADVANCED) HDF5_m_LIBRARY_DEBUG (ADVANCED) HDF5_pthread_LIBRARY_DEBUG (ADVANCED) HDF5_z_LIBRARY_DEBUG (ADVANCED) help me please! ps: I use archlinux distribution with following updated packages installed: hdf5 1.8.11-1 boost 1.55.0-4 cmake 2.8.12.1-3 python2 2.7.6-1 python2-h5py 2.2.0-1 python2-matplotlib 1.3.1-1 python2-numpy 1.8.0-1 python2-scipy 0.13.2-1 blas 3.4.2-2 lapack 3.4.2-2

1 0

Constraints for a unit cell
by Elton Nascimento 29 Dec '13

29 Dec '13

Dear All How can I constrain the particle numbers of a 2-sites unit cell for a constant N=2? I am trying to obtain local magnetization using DMRG application and t-J model. Thanks everybody. Elton M. Nascimento Instituto de Física - UFAL-Brazil

2 1

How can DMRG in ALPS calculate this unparallel edge magnetic fields model
by cuishuai_757 27 Dec '13

27 Dec '13

Dear all ALPS users, I want to calculate ground state energy of XXX Heisenberg chain of open boundary conditon with unparallel edge magnetic fields,i.e.,H=\sum_{j=1}^{N-1}\sigma_j \sigma_{j+1}+\frac{1}{p}\sigma_N^z+\frac{1}{q}(\sigma_1^z +\xi\sigma_1^x)it have sigma^z term in site N, sigma^z and sigma^x terms in site 1. but the dmrg in ALPS, can not give correct result. I think the reason is sigma^x term, which lead to the quantity "s_z" doesn't conserve, as I set the coeffecient of sigma^x to 0, the energy agrees to the data using exact diagonalization very well. Is DMRG in ALPS able to calculate this unparallel edge magnetic fields model? I believe DMRG can solve this model, should I modify some parameter or code? My files are below: LATTICE_LIBRARY="my_lattice.xml"LATTICE="open chain lattice with special edges 16 "MODEL="spin"J=4SWEEPS=4NUMBER_EIGENVALUES=1local_S=0.5h0=-0.5Gamma0=-0.0625h2=0.25{MAXSTATES=40} <LATTICES> <GRAPH name="open chain lattice with special edges 16" dimension="1" vertices=" 16 " edges="15"> <VERTEX id="1" type="0"><COORDINATE>0</COORDINATE></VERTEX> <VERTEX id="2" type="1"><COORDINATE>2</COORDINATE></VERTEX> <VERTEX id="3" type="1"><COORDINATE>3</COORDINATE></VERTEX> <VERTEX id="4" type="1"><COORDINATE>4</COORDINATE></VERTEX> <VERTEX id="5" type="1"><COORDINATE>5</COORDINATE></VERTEX> <VERTEX id="6" type="1"><COORDINATE>6</COORDINATE></VERTEX> <VERTEX id="7" type="1"><COORDINATE>7</COORDINATE></VERTEX> <VERTEX id="8" type="1"><COORDINATE>8</COORDINATE></VERTEX> <VERTEX id="9" type="1"><COORDINATE>9</COORDINATE></VERTEX> <VERTEX id="10" type="1"><COORDINATE>10</COORDINATE></VERTEX> <VERTEX id="11" type="1"><COORDINATE>11</COORDINATE></VERTEX> <VERTEX id="12" type="1"><COORDINATE>12</COORDINATE></VERTEX> <VERTEX id="13" type="1"><COORDINATE>13</COORDINATE></VERTEX> <VERTEX id="14" type="1"><COORDINATE>14</COORDINATE></VERTEX> <VERTEX id="15" type="1"><COORDINATE>15</COORDINATE></VERTEX> <VERTEX id="16" type="2"><COORDINATE>16</COORDINATE></VERTEX> <EDGE source="1" target="2" id="1" type="0" vector="1"/> <EDGE source="2" target="3" id="2" type="0" vector="1"/> <EDGE source="3" target="4" id="3" type="0" vector="1"/> <EDGE source="4" target="5" id="4" type="0" vector="1"/> <EDGE source="5" target="6" id="5" type="0" vector="1"/> <EDGE source="6" target="7" id="6" type="0" vector="1"/> <EDGE source="7" target="8" id="7" type="0" vector="1"/> <EDGE source="8" target="9" id="8" type="0" vector="1"/> <EDGE source="9" target="10" id="9" type="0" vector="1"/> <EDGE source="10" target="11" id="10" type="0" vector="1"/> <EDGE source="11" target="12" id="11" type="0" vector="1"/> <EDGE source="12" target="13" id="12" type="0" vector="1"/> <EDGE source="13" target="14" id="13" type="0" vector="1"/> <EDGE source="14" target="15" id="14" type="0" vector="1"/> <EDGE source="15" target="16" id="15" type="0" vector="1"/> </GRAPH></LATTICES> the energy of ALPS is : -6.4655612295266689numerical dialoge data: -27.684359423813586the difference is large! why? Best wishes! Shuai Cui

1 0

DMRG in ALPS unable to calculate the unparallel edge magnetic fields model
by cuishuai_757 25 Dec '13

25 Dec '13

Dear all ALPS users, I want to calculate ground state energy of XXX Heisenberg chain of open boundary conditon with unparallel edge magnetic fields,i.e.,H=\sum_{j=1}^{N-1}\sigma_j \sigma_{j+1}+\frac{1}{p}\sigma_N^z+\frac{1}{q}(\sigma_1^z +\xi\sigma_1^x)it have sigma^z term in site N, sigma^z and sigma^x terms in site 1. but the dmrg in ALPS, can not give the correct result. I think the reason is sigma^x term, which lead to the quantity "s_z" doesn't conserve, as I set the coeffecient of sigma^x to 0, the energy agrees to the data using exact diagonalization very well. Is DMRG in ALPS able to calculate this unparallel edge magnetic fields model? My files are below: LATTICE_LIBRARY="my_lattice.xml"LATTICE="open chain lattice with special edges 16 "MODEL="spin"J=4SWEEPS=4NUMBER_EIGENVALUES=1local_S=0.5h0=-0.5Gamma0=-0.0625h2=0.25{MAXSTATES=40} <LATTICES> <GRAPH name="open chain lattice with special edges 16" dimension="1" vertices=" 16 " edges="15"> <VERTEX id="1" type="0"><COORDINATE>0</COORDINATE></VERTEX> <VERTEX id="2" type="1"><COORDINATE>2</COORDINATE></VERTEX> <VERTEX id="3" type="1"><COORDINATE>3</COORDINATE></VERTEX> <VERTEX id="4" type="1"><COORDINATE>4</COORDINATE></VERTEX> <VERTEX id="5" type="1"><COORDINATE>5</COORDINATE></VERTEX> <VERTEX id="6" type="1"><COORDINATE>6</COORDINATE></VERTEX> <VERTEX id="7" type="1"><COORDINATE>7</COORDINATE></VERTEX> <VERTEX id="8" type="1"><COORDINATE>8</COORDINATE></VERTEX> <VERTEX id="9" type="1"><COORDINATE>9</COORDINATE></VERTEX> <VERTEX id="10" type="1"><COORDINATE>10</COORDINATE></VERTEX> <VERTEX id="11" type="1"><COORDINATE>11</COORDINATE></VERTEX> <VERTEX id="12" type="1"><COORDINATE>12</COORDINATE></VERTEX> <VERTEX id="13" type="1"><COORDINATE>13</COORDINATE></VERTEX> <VERTEX id="14" type="1"><COORDINATE>14</COORDINATE></VERTEX> <VERTEX id="15" type="1"><COORDINATE>15</COORDINATE></VERTEX> <VERTEX id="16" type="2"><COORDINATE>16</COORDINATE></VERTEX> <EDGE source="1" target="2" id="1" type="0" vector="1"/> <EDGE source="2" target="3" id="2" type="0" vector="1"/> <EDGE source="3" target="4" id="3" type="0" vector="1"/> <EDGE source="4" target="5" id="4" type="0" vector="1"/> <EDGE source="5" target="6" id="5" type="0" vector="1"/> <EDGE source="6" target="7" id="6" type="0" vector="1"/> <EDGE source="7" target="8" id="7" type="0" vector="1"/> <EDGE source="8" target="9" id="8" type="0" vector="1"/> <EDGE source="9" target="10" id="9" type="0" vector="1"/> <EDGE source="10" target="11" id="10" type="0" vector="1"/> <EDGE source="11" target="12" id="11" type="0" vector="1"/> <EDGE source="12" target="13" id="12" type="0" vector="1"/> <EDGE source="13" target="14" id="13" type="0" vector="1"/> <EDGE source="14" target="15" id="14" type="0" vector="1"/> <EDGE source="15" target="16" id="15" type="0" vector="1"/> </GRAPH></LATTICES> the energy of ALPS is : -6.4655612295266689numerical dialoge data: -27.684359423813586the difference is large! why? Best wishes1 Shuai Cui

1 0