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
Dear Dominique and Michele,
Could the problem be that DMRG codes need a truncation of the local basis to a finite number of states?
Matthias
On 24 Jul 2015, at 14:08, Dominique Delande Dominique.Delande@lkb.upmc.fr wrote:
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
-- Dominique Delande (Dominique.Delande@lkb.upmc.fr) Laboratoire Kastler-Brossel - Case 74 - Universite P. et M. Curie 4, place Jussieu, F-75252 Paris Cedex 05, FRANCE Phone : +33 1 44 27 27 97 - Fax : +33 1 44 27 38 45 Acces : Tour 13, Couloir 12-13, 3eme etage - Bureau 316
Dear Dominique,
The problem in the MPS codes has been fixed in revision 7593. Just upgrade to the latest /trunk version to get the fix.
Thanks a lot for reporting it. During the implementation we forgot about inhomogeneous lattices.
Note that your parameters seem to work with the sparsediag application, once one actually specifies an inhomogeneous lattice. In your XML parameters you specified the conventional “open chain lattice”.
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
+41 44 633 78 56 phone +41 44 633 11 15 fax
On 24 Jul 2015, at 22:08, Dominique Delande Dominique.Delande@lkb.upmc.fr wrote:
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
-- Dominique Delande (Dominique.Delande@lkb.upmc.fr) Laboratoire Kastler-Brossel - Case 74 - Universite P. et M. Curie 4, place Jussieu, F-75252 Paris Cedex 05, FRANCE Phone : +33 1 44 27 27 97 - Fax : +33 1 44 27 38 45 Acces : Tour 13, Couloir 12-13, 3eme etage - Bureau 316
comp-phys-alps-users@lists.phys.ethz.ch