Dear Michele,
The model is the following one: <HAMILTONIAN name="harm boson Hubbard"> <PARAMETER name="mu" default="0"/> <PARAMETER name="t" default="1"/> <PARAMETER name="V" default="0"/> <PARAMETER name="t'" default="0"/> <PARAMETER name="V'" default="0"/> <PARAMETER name="U" default="0"/> <PARAMETER name="t0" default="t"/> <PARAMETER name="t1" default="t'"/> <PARAMETER name="V0" default="V"/> <PARAMETER name="V1" default="V'"/> <PARAMETER name="K" default="0"/> <BASIS ref="boson"/> <SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> <PARAMETER name="K#" default="K"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2+K#*n(i)*(x-0.5*(L-1))^2 </SITETERM> <BONDTERM source="i" target="j"> <PARAMETER name="t#" default="0"/> <PARAMETER name="V#" default="0"/> -t#*(bdag(i)*b(j)+bdag(j)*b(i)) + V#*n(i)*n(j) </BONDTERM> </HAMILTONIAN>
and the parameters: <?xml version="1.0" encoding="UTF-8"?> <?xml-stylesheet type="text/xsl" href="ALPS.xsl"?> <SIMULATION xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://xml.comp-phys.org/2003/8/QMCXML.xsd"> <PARAMETERS> <PARAMETER name="Nmax">3</PARAMETER> <PARAMETER name="MEASURE_LOCAL[N]">n</PARAMETER> <PARAMETER name="init_state">thin</PARAMETER> <PARAMETER name="MAXSTATES">60</PARAMETER> <PARAMETER name="K">0.1</PARAMETER> <PARAMETER name="SWEEPS">16</PARAMETER> <PARAMETER name="NONLOCAL">1</PARAMETER> <PARAMETER name="L">10</PARAMETER> <PARAMETER name="LATTICE">open chain lattice</PARAMETER> <PARAMETER name="SEED">745388772</PARAMETER> <PARAMETER name="U">1.0</PARAMETER> <PARAMETER name="t">0.1</PARAMETER> <PARAMETER name="MEASURE[EnergyVariance]">1</PARAMETER> <PARAMETER name="MODEL">harm boson Hubbard</PARAMETER> <PARAMETER name="N_total">5</PARAMETER> <PARAMETER name="MODEL_LIBRARY">modelkub.xml</PARAMETER> <PARAMETER name="CONSERVED_QUANTUMNUMBERS">N</PARAMETER> </PARAMETERS> </SIMULATION>
It works perfectly with dmrg. Fulldiag and sparsediag complain: Cannot evaluate expression 0.1 * (-4.5 + x)^2 which seems to mean that the parabolic trap term K#*n(i)*(x-0.5*(L-1))^2 in the model is not understood.
Thanks for your help
Dear Dominique,
If you send the parameters and model files we could help you more. I?m currently away, I won?t be able to try it out before next week.
Did you try already if the model works with the fulldiag, sparsediag and dmrg applications?
Best regards, Michele
-- ETH Zurich Michele Dolfi Institute for Theoretical Physics HIT G 32.4 Wolfgang-Pauli-Str. 27 8093 Zurich Switzerland
dolfim@phys.ethz.ch www.itp.phys.ethz.ch